title: A Census of the 32 Ori Association with Gaia

abstract: I have used high-precision photometry and astrometry from the third data release of Gaia (DR3) to identify candidate members of the 32 Ori association. Spectral types and radial velocities have been measured for subsets of the candidates using new and archival spectra. For the candidates that have radial velocity measurements, I have used UVW velocities to further constrain their membership, arriving at a final catalog of 169 candidates. I estimate that the completeness of the survey is ~90% for spectral types of <=M7 (>=0.06 Msun). The histogram of spectral types for the 32 Ori candidates exhibits a maximum at M5 (~0.15 Msun), resembling the distributions measured for other young clusters and associations in the solar neighborhood. The available UVW velocities indicate that the association is expanding, but they do not produce a well-defined kinematic age. Based on their sequences of low-mass stars in color-magnitude diagrams, the 32 Ori association and Upper Centaurus-Lupus/Lower Centaurus-Crux (UCL/LCC) are coeval to within +/-1.2 Myr, and they are younger than the Beta Pic moving group by ~3 Myr, which agrees with results from previous analysis based on the second data release of Gaia. Finally, I have used mid-infrared (IR) photometry from the Wide-field Infrared Survey Explorer to check for excess emission from circumstellar disks among the 32 Ori candidates. Disks are detected for 18 candidates, half of which are reported for the first time in this work. The fraction of candidates at <=M6 that have full, transitional, or evolved disks is 10/149=0.07+0.03/-0.02, which is consistent with the value for UCL/LCC.


title: Universal and dynamic ridge filter for pencil beam scanning particle therapy: novel concept for ultra-fast treatment delivery

abstract: Purpose In PBS particle therapy, short treatment delivery time is paramount for the efficient treatment of moving targets with motion mitigation techniques (such as breath-hold, rescanning, and gating). Energy and spot position change time are limiting factors in reducing treatment time. In this study, we designed a universal and dynamic energy modulator (ridge filter, RF) to broaden the Bragg peak, to reduce the number of energies and spots required to cover the target volume, thus lowering the treatment time. Methods Our RF unit comprises two identical RFs placed just before the isocenter. Both RFs move relative to each other, changing the Bragg peaks characteristics dynamically. We simulated different Bragg peak shapes with the RF in TOPAS and validated them experimentally. We then delivered single-field plans with 1Gy/fraction to different geometrical targets in water, to measure the dose delivery time using the RF and compare it with the clinical settings. Results Aligning the RFs in different positions produces different broadening in the Bragg peak; we achieved a maximum broadening of 2 cm. With RF we reduced the number of energies in a field by more than 60%, and the dose delivery time by 50%, for all geometrical targets investigated, without compromising the dose distribution transverse and distal fall-off. Conclusions Our novel universal and dynamic RF allows for the adaptation of the Bragg peak broadening for a spot and/or energy layer based on the requirement of dose shaping in the target volume. It significantly reduces the number of energy layers and spots to cover the target volume, and thus the treatment time. This RF design is ideal for ultra-fast treatment delivery within a single breath-hold (5-10 sec), efficient delivery of motion mitigation techniques, and small animal irradiation with ultra-high dose rates (FLASH).


title: An Empirical Comparison of Explainable Artificial Intelligence Methods for Clinical Data: A Case Study on Traumatic Brain Injury

abstract: A longstanding challenge surrounding deep learning algorithms is unpacking and understanding how they make their decisions. Explainable Artificial Intelligence (XAI) offers methods to provide explanations of internal functions of algorithms and reasons behind their decisions in ways that are interpretable and understandable to human users. . Numerous XAI approaches have been developed thus far, and a comparative analysis of these strategies seems necessary to discern their relevance to clinical prediction models. To this end, we first implemented two prediction models for short- and long-term outcomes of traumatic brain injury (TBI) utilizing structured tabular as well as time-series physiologic data, respectively. Six different interpretation techniques were used to describe both prediction models at the local and global levels. We then performed a critical analysis of merits and drawbacks of each strategy, highlighting the implications for researchers who are interested in applying these methodologies. The implemented methods were compared to one another in terms of several XAI characteristics such as understandability, fidelity, and stability. Our findings show that SHAP is the most stable with the highest fidelity but falls short of understandability. Anchors, on the other hand, is the most understandable approach, but it is only applicable to tabular data and not time series data.


title: Metasurface-Stabilized Optical Microcavities

abstract: We demonstrate stable optical microcavities by counteracting the phase evolution of the cavity modes using an amorphous silicon metasurface as one of the two cavity end mirrors. Careful design allows us to limit the metasurface scattering losses at telecom wavelengths to less than 2% and using a distributed Bragg reflector as metasurface substrate ensures high reflectivity. Our first demonstration experimentally achieves telecom-wavelength microcavities with quality factors of up to 4600, spectral resonance linewidths below 0.4 nm, and mode volumes down to below 2.7$\lambda ^3$. We then show that the method introduces unprecedented freedom to stabilize modes with arbitrary transverse intensity profiles and design cavity-enhanced hologram modes. Our approach introduces the nanoscopic light control capabilities of dielectric metasurfaces to cavity electrodynamics and is directly industrially scalable using widespread semiconductor manufacturing processes.


title: StraKLIP: A novel pipeline for detection and characterization of close-in faint companions through Karhunen-Lo\^eve Image Processing algorithm

abstract: We present a new pipeline developed to detect and characterize faint astronomical companions at small angular separation from the host star using sets of wide-field imaging observations not specifically designed for High Contrast Imaging analysis. The core of the pipeline relies on Karhunen-Lo\^eve truncated transformation of the reference PSF library to perform PSF subtraction and identify candidates. Tests of reliability of detections and characterization of companions are made through simulation of binaries and generation of Receiver Operating Characteristic curves for false positive/true positive analysis. The algorithm has been successfully tested on large HST/ACS and WFC3 datasets acquired for two HST Treasury Programs on the Orion Nebula Cluster. Based on these extensive numerical experiments we find that, despite being based on methods designed for observations of single star at a time, our pipeline performs very well on mosaic space based data. In fact, we are able to detect brown dwarf-mass companions almost down to the planetary mass limit. The pipeline is able to reliably detect signals at separations as close as $\gtrsim 0.1 "$ with a completeness of $\gtrsim 10\%$, or $\sim 0.2"$ with a completeness of $\sim 30\%$. This approach can potentially be applied to a wide variety of space based imaging surveys, starting with data in the existing HST archive, near-future JWST mosaics, and future wide-field Roman images.


title: A Near-Optimal Algorithm for Univariate Zeroth-Order Budget Convex Optimization

abstract: This paper studies a natural generalization of the problem of minimizing a univariate convex function $f$ by querying its values sequentially. At each time-step $t$, the optimizer can invest a budget $b_t$ in a query point $X_t$ of their choice to obtain a fuzzy evaluation of $f$ at $X_t$ whose accuracy depends on the amount of budget invested in $X_t$ across times. This setting is motivated by the minimization of objectives whose values can only be determined approximately through lengthy or expensive computations. We design an any-time parameter-free algorithm called Dyadic Search, for which we prove near-optimal optimization error guarantees. As a byproduct of our analysis, we show that the classical dependence on the global Lipschitz constant in the error bounds is an artifact of the granularity of the budget. Finally, we illustrate our theoretical findings with numerical simulations.


title: Lyapunov Design for Robust and Efficient Robotic Reinforcement Learning

abstract: Recent advances in the reinforcement learning (RL) literature have enabled roboticists to automatically train complex policies in simulated environments. However, due to the poor sample complexity of these methods, solving RL problems using real-world data remains a challenging problem. This paper introduces a novel cost-shaping method which aims to reduce the number of samples needed to learn a stabilizing controller. The method adds a term involving a Control Lyapunov Function (CLF) -- an `energy-like' function from the model-based control literature -- to typical cost formulations. Theoretical results demonstrate the new costs lead to stabilizing controllers when smaller discount factors are used, which is well-known to reduce sample complexity. Moreover, the addition of the CLF term `robustifies' the search for a stabilizing controller by ensuring that even highly sub-optimal polices will stabilize the system. We demonstrate our approach with two hardware examples where we learn stabilizing controllers for a cartpole and an A1 quadruped with only seconds and a few minutes of fine-tuning data, respectively. Furthermore, simulation benchmark studies show that obtaining stabilizing policies by optimizing our proposed costs requires orders of magnitude less data compared to standard cost designs.


title: A hands-on gaze on HTTP/3 security through the lens of HTTP/2 and a public dataset

abstract: Following QUIC protocol ratification on May 2021, the third major version of the Hypertext Transfer Protocol, namely HTTP/3, was published around one year later in RFC 9114. In light of these consequential advancements, the current work aspires to provide a full-blown coverage of the following issues, which to our knowledge have received feeble or no attention in the literature so far. First, we provide a complete review of attacks against HTTP/2, and elaborate on if and in which way they can be migrated to HTTP/3. Second, through the creation of a testbed comprising the at present six most popular HTTP/3-enabled servers, we examine the effectiveness of a quartet of attacks, either stemming directly from the HTTP/2 relevant literature or being entirely new. This scrutiny led to the assignment of at least one CVE ID with a critical base score by MITRE. No less important, by capitalizing on a realistic, abundant in devices testbed, we compiled a voluminous, labeled corpus containing traces of ten diverse attacks against HTTP and QUIC services. An initial evaluation of the dataset mainly by means of machine learning techniques is included as well. Given that the 30 GB dataset is made available in both pcap and CSV formats, forthcoming research can easily take advantage of any subset of features, contingent upon the specific network topology and configuration.


title: Jacobi Set Driven Search for Flexible Fiber Surface Extraction

abstract: Isosurfaces are an important tool for analysis and visualization of univariate scalar fields. Earlier works have demonstrated the presence of interesting isosurfaces at isovalues close to critical values. This motivated the development of efficient methods for computing individual components of isosurfaces restricted to a region of interest. Generalization of isosurfaces to fiber surfaces and critical points to Jacobi sets has resulted in new approaches for analyzing bivariate scalar fields. Unlike isosurfaces, there exists no output sensitive method for computing fiber surfaces. Existing methods traverse through all the tetrahedra in the domain. In this paper, we propose the use of the Jacobi set to identify fiber surface components of interest and present an output sensitive approach for its computation. The Jacobi edges are used to initiate the search towards seed tetrahedra that contain the fiber surface, thereby reducing the search space. This approach also leads to effective analysis of the bivariate field by supporting the identification of relevant fiber surfaces near Jacobi edges.


title: CollComm: Enabling Efficient Collective Quantum Communication Based on EPR buffering

abstract: The noisy and lengthy nature of quantum communication hinders the development of distributed quantum computing. The inefficient design of existing compilers for distributed quantum computing worsens the situation. Previous compilation frameworks couple communication hardware with the implementation of expensive remote gates. However, we discover that the efficiency of quantum communication, especially collective communication, can be significantly boosted by decoupling communication resources from remote operations, that is, the communication hardware would be used only for preparing remote entanglement, and the computational hardware, the component used to store program information, would be used for conducting remote gates. Based on the observation, we develop a compiler framework to optimize the collective communication happening in distributed quantum programs. In this framework, we decouple the communication preparation process in communication hardware from the remote gates conducted in computational hardware by buffering EPR pairs generated by communication hardware in qubits of the computational hardware. Experimental results show that the proposed framework can halve the communication cost of various distributed quantum programs, compared to state-of-the-art compilers for distributed quantum computing.


title: Propagation of singularities for the wave equation

abstract: This expository note gives a digest version of Hormander's propagation of singularities theorem for the wave equation.


title: Predicting the cascading dynamics in complex networks via the bimodal failure size distribution

abstract: Cascading failure as a systematic risk occurs in a wide range of real-world networks. Cascade size distribution is a basic and crucial characteristic of systemic cascade behaviors. Recent research works have revealed that the distribution of cascade sizes is a bimodal form indicating the existence of either very small cascades or large ones. In this paper, we aim to understand the properties and formation of such bimodal distribution of cascade sizes in complex networks, and further predict the final cascade size. We first find that the bimodal distribution of cascade sizes is ubiquitous in both synthetic and real networks. Moreover, the large cascade sizes distributed in the right peak of bimodal distribution are resulted from either the failure of nodes with high load at the first step of the cascade or multiple rounds of cascades triggered by the initial failure. Accordingly, we propose a hybrid load metric (HLM), which combines the load of the initial broken node and the load of failed nodes triggered by the initial failure, to predict the final size of cascading failures. Finally, we validate the effectiveness of HLM by computing the accuracy of identifying the cascades belonging to the right and left peaks of the bimodal distribution. The results show that HLM is a better predictor than commonly used network centrality metrics in both synthetic and real-world networks.


title: Reliable emulation of complex functionals by active learning with error control

abstract: A statistical emulator can be used as a surrogate of complex physics-based calculations to drastically reduce the computational cost. Its successful implementation hinges on an accurate representation of the nonlinear response surface with a high-dimensional input space. Conventional "space-filling" designs, including random sampling and Latin hypercube sampling, become inefficient as the dimensionality of the input variables increases, and the predictive accuracy of the emulator can degrade substantially for a test input distant from the training input set. To address this fundamental challenge, we develop a reliable emulator for predicting complex functionals by active learning with error control (ALEC). The algorithm is applicable to infinite-dimensional mapping with high-fidelity predictions and a controlled predictive error. The computational efficiency has been demonstrated by emulating the classical density functional theory (cDFT) calculations, a statistical-mechanical method widely used in modeling the equilibrium properties of complex molecular systems. We show that ALEC is much more accurate than conventional emulators based on the Gaussian processes with "space-filling" designs and alternative active learning methods. Besides, it is computationally more efficient than direct cDFT calculations. ALEC can be a reliable building block for emulating expensive functionals owing to its minimal computational cost, controllable predictive error, and fully automatic features.


title: On making physics relevant to society in general and to scientists in particular: Closing the epistemic gap

abstract: Physics has a bad press: it is seen by a majority of people as a boring discipline ever since their High School days. There is no glamour to it, just toil and pain, and for many who engaged in it, the end sight is often unemployment. Could it be that physicists don't know how to communicate what their discipline entails to? I will be tackling the problematic of making physics relevant to society and to the scientists in general. I will also be dealing with the methodological and educational aspects of teaching and practicing physics, and the need to close the epistemic gap between physics teaching and the physicist's understanding ... By the way, do physicists understand physics?


title: Optimal Recovery for Causal Inference

abstract: Problems in causal inference can be fruitfully addressed using signal processing techniques. As an example, it is crucial to successfully quantify the causal effects of an intervention to determine whether the intervention achieved desired outcomes. We present a new geometric signal processing approach to classical synthetic control called ellipsoidal optimal recovery (EOpR), for estimating the unobservable outcome of a treatment unit. EOpR provides policy evaluators with both worst-case and typical outcomes to help in decision making. It is an approximation-theoretic technique that relates to the theory of principal components, which recovers unknown observations given a learned signal class and a set of known observations. We show EOpR can improve pre-treatment fit and mitigate bias of the post-treatment estimate relative to other methods in causal inference. Beyond recovery of the unit of interest, an advantage of EOpR is that it produces worst-case limits over the estimates produced. We assess our approach on artificially-generated data, on datasets commonly used in the econometrics literature, and in the context of the COVID-19 pandemic, showing better performance than baseline techniques


title: A simple micro-swimmer model inspired by the general equation for nonequilibrium reversible-irreversible coupling

abstract: A simple mean-field micro-swimmer model is presented. The model is inspired by the nonequilibrium thermodynamics of multi-component fluids that undergo chemical reactions. These thermodynamics can be rigorously described in the context of the GENERIC (general equation for the nonequilibrium reversible-irreversible coupling) framework. More specifically this approach was recently applied to non-ideal polymer solutions (Indei and Schieber, J. Chem. Phys, 146, 184902, 2017). One of the species of the solution is an unreactive polymer chain represented by the bead-spring model. Using this detailed description as inspiration we then make several simplifying assumptions to obtain a mean-field model for a Janus micro-swimmer. The swimmer model considered here consists of a polymer dumbbell in a sea of reactants. One of the beads of the dumbbell is allowed to act as a catalyst for a chemical reaction between the reactants. We show that the mean-squared displacement, MSD, of the center of mass of this Janus dumbbell exhibits ballistic behavior at time scales at which the concentration of reactant is large. The time scales at which the ballistic behavior is observed in the MSD coincides with the time scales at which the cross-correlation between the swimmer's orientation and the direction of its displacement exhibit a maximum. Since the swimmer model was inspired by the GENERIC framework it is possible to ensure that the entropy generation is always positive and therefore the second law of thermodynamics is obeyed.


title: All identically conserved gravitational tensors are metric variations of an action

abstract: I prove an old unsolved conjecture, the hard part of its obvious sufficiency, namely that all identically conserved symmetric two-tensors are necessarily metric variations of invariant actions, thus sparing the world more alternative gravity theories. The proof is reasonably simple, if perhaps a physicist's.


title: An integrated photonic engine for programmable atomic control

abstract: Solutions for scalable, high-performance optical control are important for the development of scaled atom-based quantum technologies. Modulation of many individual optical beams is central to the application of arbitrary gate and control sequences on arrays of atoms or atom-like systems. At telecom wavelengths, miniaturization of optical components via photonic integration has pushed the scale and performance of classical and quantum optics far beyond the limitations of bulk devices. However, these material platforms for high-speed telecom integrated photonics are not transparent at the short wavelengths required by leading atomic systems. Here, we propose and implement a scalable and reconfigurable photonic architecture for multi-channel quantum control using integrated, visible-light modulators based on thin-film lithium niobate. Our approach combines techniques in free-space optics, holography, and control theory together with a sixteen-channel integrated photonic device to stabilize temporal and cross-channel power deviations and enable precise and uniform control. Applying this device to a homogeneous constellation of silicon-vacancy artificial atoms in diamond, we present techniques to spatially and spectrally address a dynamically-selectable set of these stochastically-positioned point emitters. We anticipate that this scalable and reconfigurable optical architecture will lead to systems that could enable parallel individual programmability of large many-body atomic systems, which is a critical step towards universal quantum computation on such hardware.


title: Automated Conversion of Axiomatic to Operational Models: Theory and Practice

abstract: A system may be modelled as an operational model (which has explicit notions of state and transitions between states) or an axiomatic model (which is specified entirely as a set of invariants). Most formal methods techniques (e.g., IC3, invariant synthesis, etc) are designed for operational models and are largely inaccessible to axiomatic models. Furthermore, no prior method exists to automatically convert axiomatic models to operational ones, so operational equivalents to axiomatic models had to be manually created and proven equivalent. In this paper, we advance the state-of-the-art in axiomatic to operational model conversion. We show that general axioms in the $\mu$spec axiomatic modelling framework cannot be translated to equivalent finite-state operational models. We also derive restrictions on the space of $\mu$spec axioms that enable the feasible generation of equivalent finite-state operational models for them. As for practical results, we develop a methodology for automatically translating $\mu$spec axioms to equivalent finite-state automata-based operational models. We demonstrate the efficacy of our method by using the models generated by our procedure to prove the correctness of ordering properties on three RTL designs.


title: An Answer Verbalization Dataset for Conversational Question Answerings over Knowledge Graphs

abstract: We introduce a new dataset for conversational question answering over Knowledge Graphs (KGs) with verbalized answers. Question answering over KGs is currently focused on answer generation for single-turn questions (KGQA) or multiple-tun conversational question answering (ConvQA). However, in a real-world scenario (e.g., voice assistants such as Siri, Alexa, and Google Assistant), users prefer verbalized answers. This paper contributes to the state-of-the-art by extending an existing ConvQA dataset with multiple paraphrased verbalized answers. We perform experiments with five sequence-to-sequence models on generating answer responses while maintaining grammatical correctness. We additionally perform an error analysis that details the rates of models' mispredictions in specified categories. Our proposed dataset extended with answer verbalization is publicly available with detailed documentation on its usage for wider utility.


title: Positivity of the exterior power of the tangent bundles

abstract: Let $X$ be a complex smooth projective variety such that the exterior power of the tangent bundle $\bigwedge^{r} T_X$ is nef for some $1\leq r<\dim X$. We prove that, up to an \'etale cover, $X$ is a Fano fiber space over an Abelian variety. This gives generalizations of the structure theorem of varieties with nef tangent bundle by Demailly, Peternell and Schneider and that of varieties with nef $\bigwedge^{2} T_X$ by the author. Our result also gives an answer to a question raised by Li, Ou and Yang for varieties with strictly nef $\bigwedge^{r} T_X$ when $r < \dim X$.


title: Linear Stability of liquid Lane-Emden stars

abstract: We establish various qualitative properties of liquid Lane-Emden stars in $\mathbb{R}^d$, including bounds for its density profile $\rho$ and radius $R$. Using them we prove that against radial perturbations, the liquid Lane-Emden stars are linearly stable when $\gamma\geq 2(d-1)/d$; linearly stable when $\gamma<2(d-1)/d$ for stars with small relative central density $\rho(0)-\rho(R)$; and linearly unstable when $\gamma<2(d-1)/d$ for stars with large central density. Such dependence on central density is not seen in the gaseous Lane-Emden stars.


title: The Effects of the Interplay Between Motor and Brownian Forces on the Rheology of Active Gels

abstract: Active gels perform key mechanical roles inside the cell, such as cell division, motion and force sensing. The unique mechanical properties required to perform such functions arise from the interactions between molecular motors and semi-flexible polymeric filaments. Molecular motors can convert the energy released in the hydrolysis of ATP into forces of up to pico-Newton magnitudes. Moreover, the polymeric filaments that form active gels are flexible enough to respond to Brownian forces, but also stiff enough to support the large tensions induced by the motor-generated forces. Brownian forces are expected to have a significant effect especially at motor activities at which stable non-contractile in vitro active gels are prepared for rheological measurements. Here, a microscopic mean-field theory of active gels originally formulated in the limit of motor-dominated dynamics is extended to include Brownian forces. In the model presented here Brownian forces are included accurately, at real room temperature, even in systems with high motor activity. It is shown that a subtle interplay, or competition, between motor-generated forces and Brownian forces has an important impact in the mass transport and rheological properties of active gels. The model predictions show that at low frequencies the dynamic modulus of active gels is determined mostly by motor protein dynamics. However, Brownian forces significantly increase the breadth of the relaxation spectrum and can affect the shape of the dynamic modulus over a wide frequency range even for ratios of motor to Brownian forces of more than a hundred. Since the ratio between motor and Brownian forces is sensitive to ATP concentration, the results presented here shed some light on how the transient mechanical response of active gels changes with varying ATP concentration.


title: Structure induced by a multiple membership transformation on the Conditional Autoregressive model

abstract: The objective of disease mapping is to model data aggregated at the areal level. In some contexts, however, (e.g. residential histories, general practitioner catchment areas) when data is arising from a variety of sources, not necessarily at the same spatial scale, it is possible to specify spatial random effects, or covariate effects, at the areal level, by using a multiple membership principle (MM) (Petrof et al. 2020, Gramatica et al. 2021). A weighted average of conditional autoregressive (CAR) spatial random effects embeds spatial information for a spatially-misaligned outcome and estimate relative risk for both frameworks (areas and memberships). In this paper we investigate the theoretical underpinnings of these application of the multiple membership principle to the CAR prior, in particular with regard to parameterisation, properness and identifiability. We carry out simulations involving different numbers of memberships as compared to number of areas and assess impact of this on estimating parameters of interest. Both analytical and simulation study results show under which conditions parameters of interest are identifiable, so that we can offer actionable recommendations to practitioners. Finally, we present the results of an application of the multiple membership model to diabetes prevalence data in South London, together with strategic implications for public health considerations


title: Machine Learning Based Radiomics for Glial Tumor Classification and Comparison with Volumetric Analysis

abstract: Purpose; The purpose of this study is to classify glial tumors into grade II, III and IV categories noninvasively by application of machine learning to multi-modal MRI features in comparison with volumetric analysis. Methods; We retrospectively studied 57 glioma patients with pre and postcontrast T1 weighted, T2 weighted, FLAIR images, and ADC maps acquired on a 3T MRI. The tumors were segmented into enhancing and nonenhancing portions, tumor necrosis, cyst and edema using semiautomated segmentation of ITK-SNAP open source tool. We measured total tumor volume, enhancing-nonenhancing tumor, edema, necrosis volume and the ratios to the total tumor volume. Training of a support vector machine (SVM) classifier and artificial neural network (ANN) was performed with labeled data designed to answer the question of interest. Specificity, sensitivity, and AUC of the predictions were computed by means of ROC analysis. Differences in continuous measures between groups were assessed by using Kruskall Wallis, with post hoc Dunn correction for multiple comparisons. Results; When we compared the volume ratios between groups, there was statistically significant difference between grade IV and grade II-III glial tumors. Edema and tumor necrosis volume ratios for grade IV glial tumors were higher than that of grade II and III. Volumetric ratio analysis could not distinguish grade II and III tumors successfully. However, SVM and ANN correctly classified each group with accuracies up to 98% and 96%. Conclusion; Application of machine learning methods to MRI features can be used to classify brain tumors noninvasively and more readily in clinical settings.


title: Lipschitz decompositions of domains with bilaterally flat boundaries

abstract: We study classes of domains in $\mathbb{R}^{d+1},\ d \geq 2$ with sufficiently flat boundaries which admit a decomposition or covering of bounded overlap by Lipschitz graph domains with controlled total surface area. This study is motivated by the following result proved by Peter Jones as a piece of his proof of the Analyst's Traveling Salesman Theorem in the complex plane: Any simply connected domain $\Omega\subseteq\mathbb{C}$ with finite boundary length $\mathcal{H}^1(\partial\Omega)$ can be decomposed into Lipschitz graph domains with total boundary length bounded above by $M\mathcal{H}^1(\partial\Omega)$ for some $M$ independent of $\Omega$. In this paper, we prove an analogous Lipschitz decomposition result in higher dimensions for domains with Reifenberg flat boundaries satisfying a uniform beta-squared sum bound. We use similar techniques to show that domains with general Reifenberg flat or uniformly rectifiable boundaries admit similar Lipschitz decompositions while allowing the constituent domains to have bounded overlaps rather than be disjoint.


title: Indecomposable solutions of the Yang-Baxter equation with permutation group of sizes $pq$ and $p^2q$

abstract: In this paper we study the problem of classification of indecomposable solutions of the Yang-Baxter equation. Using a scheme proposed by Bachiller, Ced\'o, and Jespers, and recent advances in the classification of braces we classify all indecomposable solutions with some particular permutation groups. We do this for all groups of size $pq$, all abelian groups of size $p^2q$ and all dihedral groups of size $p^2q$.


title: Feasibility Layer Aided Machine Learning Approach for Day-Ahead Operations

abstract: Day-ahead operations involves a complex and computationally intensive optimization process to determine the generator commitment schedule and dispatch. The optimization process is a mixed-integer linear program (MILP) also known as security-constrained unit commitment (SCUC). Independent system operators (ISOs) run SCUC daily and require state-of-the-art algorithms to speed up the process. Existing patterns in historical information can be leveraged for model reduction of SCUC, which can provide significant time savings. In this paper, machine learning (ML) based classification approaches, namely logistic regression, neural networks, random forest and K-nearest neighbor, were studied for model reduction of SCUC. The ML was then aided with a feasibility layer (FL) and post-process technique to ensure high-quality solutions. The proposed approach is validated on several test systems namely, IEEE 24-Bus system, IEEE-73 Bus system, IEEE 118-Bus system, 500-Bus system, and Polish 2383-Bus system. Moreover, model reduction of a stochastic SCUC (SSCUC) was demonstrated utilizing a modified IEEE 24-Bus system with renewable generation. Simulation results demonstrate a high training accuracy to identify commitment schedule while FL and post-process ensure ML predictions do not lead to infeasible solutions with minimal loss in solution quality.


title: Enhancing Graph Contrastive Learning with Node Similarity

abstract: Graph Neural Networks (GNNs) have achieved great success in learning graph representations and thus facilitating various graph-related tasks. However, most GNN methods adopt a supervised learning setting, which is not always feasible in real-world applications due to the difficulty to obtain labeled data. Hence, graph self-supervised learning has been attracting increasing attention. Graph contrastive learning (GCL) is a representative framework for self-supervised learning. In general, GCL learns node representations by contrasting semantically similar nodes (positive samples) and dissimilar nodes (negative samples) with anchor nodes. Without access to labels, positive samples are typically generated by data augmentation, and negative samples are uniformly sampled from the entire graph, which leads to a sub-optimal objective. Specifically, data augmentation naturally limits the number of positive samples that involve in the process (typically only one positive sample is adopted). On the other hand, the random sampling process would inevitably select false-negative samples (samples sharing the same semantics with the anchor). These issues limit the learning capability of GCL. In this work, we propose an enhanced objective that addresses the aforementioned issues. We first introduce an unachievable ideal objective that contains all positive samples and no false-negative samples. This ideal objective is then transformed into a probabilistic form based on the distributions for sampling positive and negative samples. We then model these distributions with node similarity and derive the enhanced objective. Comprehensive experiments on various datasets demonstrate the effectiveness of the proposed enhanced objective under different settings.


title: Self-avoiding walks and polygons crossing a domain on the square and hexagonal lattices

abstract: We have analysed the recently extended series for the number of self-avoiding walks (SAWs) $C_L(1)$ that cross an $L \times L$ square between diagonally opposed corners. The number of such walks is known to grow as $\lambda_S^{L^2}.$ We have made more precise the estimate of $\lambda_S,$ based on additional series coefficients provided by several authors, and refined analysis techniques. We estimate that $\lambda_S = 1.7445498 \pm 0.0000012.$ We have also studied the subdominant behaviour, and conjecture that $$ C_L(1) \sim \lambda_S^{L^2+bL+c}\cdot L^g,$$ where $b=-0.04354 \pm 0.0001,$ $c=0.5624 \pm 0.0005,$ and $g=0.000 \pm 0.005.$ We implemented a very efficient algorithm for enumerating paths on the square and hexagonal lattices making use of a minimal perfect hash function and in-place memory updating of the arrays for the counts of the number of paths. Using this algorithm we extended and then analysed series for SAWs spanning the square lattice and self-avoiding polygons (SAPs) crossing the square lattice. These are known to also grow as $\lambda_S^{L^2}.$ The sub-dominant term $\lambda^b$ is found to be the same as for SAWs crossing the square, while the exponent $g = 1.75\pm 0.01$ for spanning SAWs and $g = -0.500 \pm 0.005$ for SAPs. We have also studied the analogous problems on the hexagonal lattice, and generated series for a number of geometries. In particular, we study SAWs and SAPs crossing rhomboidal, triangular and square domains on the hexagonal lattice, as well as SAWs spanning a rhombus. We estimate that the analogous growth constant $\lambda_H=1.38724951 \pm 0.00000005,$ so an even more precise estimate than found for the square lattice. We also give estimates of the sub-dominant terms.


title: 2D-XY ferromagnetism with high transition temperature in Janus monolayer V$_{2}$XN (X = P, As)

abstract: Two-dimensional (2D) XY magnets with easy magnetization planes support the nontrivial topological spin textures whose dissipationless transport is highly desirable for 2D spintronic devices. Here, we predicted that Janus monolayer V$_{2}$XN (X = P, As) with a square lattice are 2D-XY ferromagnets by first-principles calculations. Both the magnetocrystalline anisotropy and magnetic shape anisotropy favor an in-plane magnetization, leading to an easy magnetization $xy$-plane in Janus monolayer V$_{2}$XN. Resting on the Monte Carlo simulations, we observed the Berezinskii-Kosterlitz-Thouless (BKT) phase transition in monolayer V$_{2}$XN with transition temperature $T_{\rm BKT}$ being above the room temperature. Especially, monolayer V$_{2}$AsN has a magnetic anisotropy energy (MAE) of 292.0 $\mu$eV per V atom and a $T_{\rm BKT}$ of 434 K, which is larger than that of monolayer V$_{2}$PN. Moreover, a tensile strain of 5\% can further improve the $T_{\rm BKT}$ of monolayer V$_{2}$XN to be above 500 K. Our results indicated that Janus monolayer V$_{2}$XN (X = P, As) were candidate materials to realize high-temperature 2D-XY ferromagnetism for spintronics applications.


title: Contrastive Counterfactual Learning for Causality-aware Interpretable Recommender Systems

abstract: The field of generating recommendations within the framework of causal inference has seen a recent surge, with recommendations being likened to treatments. This approach enhances insights into the influence of recommendations on user behavior and helps in identifying the underlying factors. Existing research has often leveraged propensity scores to mitigate bias, albeit at the risk of introducing additional variance. Others have explored the use of unbiased data from randomized controlled trials, although this comes with assumptions that may prove challenging in practice. In this paper, we first present the causality-aware interpretation of recommendations and reveal how the underlying exposure mechanism can bias the maximum likelihood estimation (MLE) of observational feedback. Recognizing that confounders may be elusive, we propose a contrastive self-supervised learning to minimize exposure bias, employing inverse propensity scores and expanding the positive sample set. Building on this foundation, we present a novel contrastive counterfactual learning method (CCL) that incorporates three unique positive sampling strategies grounded in estimated exposure probability or random counterfactual samples. Through extensive experiments on two real-world datasets, we demonstrate that our CCL outperforms the state-of-the-art methods.


title: Volume effects on the QCD critical end point from thermal fluctuations within the super statistics framework

abstract: We investigate the impact of the finite volume and the thermal fluctuations on the Critical End Point of the QCD phase diagram. To do so, we implement the super statistics framework with Gamma, $F$, and log-normal distributions and their relation with the Tsallis non-extensive thermodynamics. We compute an effective thermodynamic potential as a function of the inverse temperature fluctuations and explicit dependence on the system volume. To find an analytic expression for the effective potential, we expand the modified Boltzmann factor by using the equilibrium thermodynamic potential computed in the Linear Sigma Model coupled to quarks. We find that the pseudo-critical temperature of transition at vanishing baryon chemical potential is modified by the size of the system being about $7\%$ lower for small volumes. Additionally, the critical endpoint moves to higher densities and lower temperatures (about $12\%$ in both cases). Interestingly, the results are quantitatively the same when the parameter that models the out-of-equilibrium situation is modified, indicating that the chiral symmetry restoration is robust against the thermal fluctuations in this approximation.


title: Learning to Infer Counterfactuals: Meta-Learning for Estimating Multiple Imbalanced Treatment Effects

abstract: We regularly consider answering counterfactual questions in practice, such as "Would people with diabetes take a turn for the better had they choose another medication?". Observational studies are growing in significance in answering such questions due to their widespread accumulation and comparatively easier acquisition than Randomized Control Trials (RCTs). Recently, some works have introduced representation learning and domain adaptation into counterfactual inference. However, most current works focus on the setting of binary treatments. None of them considers that different treatments' sample sizes are imbalanced, especially data examples in some treatment groups are relatively limited due to inherent user preference. In this paper, we design a new algorithmic framework for counterfactual inference, which brings an idea from Meta-learning for Estimating Individual Treatment Effects (MetaITE) to fill the above research gaps, especially considering multiple imbalanced treatments. Specifically, we regard data episodes among treatment groups in counterfactual inference as meta-learning tasks. We train a meta-learner from a set of source treatment groups with sufficient samples and update the model by gradient descent with limited samples in target treatment. Moreover, we introduce two complementary losses. One is the supervised loss on multiple source treatments. The other loss which aligns latent distributions among various treatment groups is proposed to reduce the discrepancy. We perform experiments on two real-world datasets to evaluate inference accuracy and generalization ability. Experimental results demonstrate that the model MetaITE matches/outperforms state-of-the-art methods.


title: Tensor Algebra on an Optoelectronic Microchip

abstract: Tensor algebra lies at the core of computational science and machine learning. Due to its high usage, entire libraries exist dedicated to improving its performance. Conventional tensor algebra performance boosts focus on algorithmic optimizations, which in turn lead to incremental improvements. In this paper, we describe a method to accelerate tensor algebra a different way: by outsourcing operations to an optical microchip. We outline a numerical programming language developed to perform tensor algebra computations that is designed to leverage our optical hardware's full potential. We introduce the language's current grammar and go over the compiler design. We then show a new way to store sparse rank-n tensors in RAM that outperforms conventional array storage (used by C++, Java, etc.). This method is more memory-efficient than Compressed Sparse Fiber (CSF) format and is specifically tuned for our optical hardware. Finally, we show how the scalar-tensor product, rank-$n$ Kronecker product, tensor dot product, Khatri-Rao product, face-splitting product, and vector cross product can be compiled into operations native to our optical microchip through various tensor decompositions.


title: Medium induced jet broadening in a quantum computer

abstract: QCD jets provide one of the best avenues to extract information about the quark-gluon plasma produced in the aftermath of ultra relativistic heavy ions collisions. The structure of jets is determined by multiparticle quantum interference hard to tackle using perturbative methods. When jets evolve in a QCD medium this interference pattern is modified, adding another layer of complexity. By taking advantage of the recent developments in quantum technologies, such effects might be better understood via direct quantum simulation of jet evolution. In this work, we introduce a precursor to such simulations. Based on the light-front Hamiltonian formalism, we construct a digital quantum circuit that tracks the evolution of a single hard probe in the presence of a stochastic color background. In terms of the jet quenching parameter $\hat q$, the results obtained using classical simulators of ideal quantum computers agree with known analytical results. With this study, we hope to provide a baseline for future in-medium jet physics studies using quantum computers.


title: A note on another approach on power sums

abstract: In this note, we first review the novel approach to power sums put forward recently by Muschielok in arXiv:2207.01935v1, which can be summarized by the formula $S_m^{(a)}(n) = \sum_{k} c_{mk} \psi_k^{(a)}(n)$, where the $c_{mk}$'s are the expansion coefficients and where the basis functions $\psi_m^{(a)}(n)$ fulfil the recursive property $\psi_m^{(a+1)}(n)= \sum_{i=1}^n \psi_m^{(a)}(i)$. Then, we point out a number of supplementary facts concerning the said approach not contemplated explicitly in Muschielok's paper. In particular, we show that, for any given $m$, the values of the $c_{mk}$'s can be obtained by inverting a matrix involving only binomial coefficients. This may be compared with the original approach of Muschielok, where the values of the $c_{mk}$'s can be obtained by inverting a lower triangular matrix involving the Stirling numbers of the first kind. Also, we make a conjecture about the functional form of the coefficients $c_{m\, m-k}$.


title: DAOS as HPC Storage, a view from Numerical Weather Prediction

abstract: Object storage solutions potentially address long-standing performance issues with POSIX file systems for certain I/O workloads, and new storage technologies offer promising performance characteristics for data-intensive use cases. In this work, we present a preliminary assessment of Intel's Distributed Asynchronous Object Store (DAOS), an emerging high-performance object store, in conjunction with non-volatile storage and evaluate its potential use for HPC storage. We demonstrate DAOS can provide the required performance, with bandwidth scaling linearly with additional DAOS server nodes in most cases, although choices in configuration and application design can impact achievable bandwidth. We describe a new I/O benchmark and associated metrics that address object storage performance from application-derived workloads.


title: Sharp Frequency Bounds for Sample-Based Queries

abstract: A data sketch algorithm scans a big data set, collecting a small amount of data -- the sketch, which can be used to statistically infer properties of the big data set. Some data sketch algorithms take a fixed-size random sample of a big data set, and use that sample to infer frequencies of items that meet various criteria in the big data set. This paper shows how to statistically infer probably approximately correct (PAC) bounds for those frequencies, efficiently, and precisely enough that the frequency bounds are either sharp or off by only one, which is the best possible result without exact computation.


title: Impact of curvature based geometric constraints on $F(R)$ theory

abstract: Theories of gravity are fundamentally a relation between matter and the geometric structure of the underlying spacetime. So once we put some additional restrictions on the spacetime geometry, the theory of gravity is bound to get the impact, irrespective of whether it is general relativity or the modified theories of gravity. In the present article, we consider two curvature-based constraints, namely the almost pseudo-Ricci symmetric and weakly Ricci symmetric condition. As a novel result, such spacetimes with non-null associated vectors are entirely classified, and then applying the obtained results, we investigate these spacetimes as solutions of the $F(R)$-gravity theory. The modified Friedmann equations are derived and analysed in a model-independent way first. Finally, two $F(R)$ gravity models are examined for recent observational constrained values of the deceleration, jerk, and Hubble parameters. We further discuss the behavior of energy conditions.


title: Left-invariant almost para-K\"{a}hler structures on six-dimensional nilpotent Lie groups

abstract: In this paper, we consider left-invariant para-complex structures on six-dimensional nilpotent Lie groups. A complete list of six-dimensional nilpotent Lie groups that admit para-K\"{a}hler structures is obtained, explicit expressions for para-complex structures are found, and curvature properties of associated para-K\"{a}hler metrics are investigated. It is shown that paracomplex structures are nilpotent and the corresponding para-K\"{a}hler metrics are Ricci-flat.


title: Predicting skull fractures via CNN with classification algorithms

abstract: Computer Tomography (CT) images have become quite important to diagnose diseases. CT scan slice contains a vast amount of data that may not be properly examined with the requisite precision and speed using normal visual inspection. A computer-assisted skull fracture classification expert system is needed to assist physicians. Convolutional Neural Networks (CNNs) are the most extensively used deep learning models for image categorization since most often time they outperform other models in terms of accuracy and results. The CNN models were then developed and tested, and several convolutional neural network (CNN) architectures were compared. ResNet50, which was used for feature extraction combined with a gradient boosted decision tree machine learning algorithm to act as a classifier for the categorization of skull fractures from brain CT scans into three fracture categories, had the best overall F1-score of 96%, Hamming Score of 95%, Balanced accuracy Score of 94% & ROC AUC curve of 96% for the classification of skull fractures.


title: A Preliminary Study on the Potential Usefulness of Open Domain Model for Missing Software Requirements Recommendation

abstract: Completeness is one of the most important attributes of software requirement specifications. Unfortunately, incompleteness is meanwhile one of the most difficult problems to detect. Some approaches have been proposed to detect missing requirements based on the requirement-oriented domain model. However, this kind of models are lacking for lots of domains. Fortunately, the domain models constructed for different purposes can usually be found online. This raises a question: whether or not these domain models are helpful in finding the missing functional information in requirement specification? To explore this question, we design and conduct a preliminary study by computing the overlapping rate between the entities in domain models and the concepts of natural language software requirements and then digging into four regularities of the occurrence of these entities(concepts) based on two example domains. The usefulness of these regularities, especially the one based on our proposed metric AHME (with F2 gains of 146% and 223% on the two domains than without any regularity), has been shown in experiments.


title: The ultraviolet CII lines as a diagnostic of kappa-distributed electrons in planetary nebulae

abstract: Non-Maxwellian $\kappa$ electron energy distributions (EEDs) have been proposed in recent years to resolve the so-called ``electron temperature and abundance discrepancy problem'' in the study of planetary nebulae (PNe). Thus the need to develop diagnostic tools to determine from observations the EED of PNe is raised. Arising from high energy levels, the ultraviolet (UV) emission lines from PNe present intensities that depend sensitively on the high-energy tail of the EED. In this work, we investigate the feasibility of using the \ion{C}{2}]$\lambda$2326/\ion{C}{2}$\lambda$1335 intensity ratios as a diagnostic of the deviation of the EED from the Maxwellian distribution (as represented by the $\kappa$ index). We use a Maxwellian decomposition approach to derive the theoretical $\kappa$-EED-based collisionally excited coefficients of \ion{C}{2}, and then compute the \ion{C}{2} UV intensity ratio as a function of the $\kappa$ index. We analyze the archival spectra acquired by the {\it International Ultraviolet Explorer} and measure the intensities of \ion{C}{2} UV lines from 12 PNe. By comparing the observed line ratios and the theoretical predictions, we can infer their $\kappa$ values. With the Maxwellian-EED hypothesis, the observed \ion{C}{2}]$\lambda$2326/\ion{C}{2}$\lambda$1335 ratios are found to be generally lower than those predicted from the observed optical spectra. This discrepancy can be explained in terms of the $\kappa$ EED. Our results show that the $\kappa$ values inferred range from 15 to infinity, suggesting a mild or modest deviation from the Maxwellian distribution. However, the $\kappa$-distributed electrons are unlikely to exist throughout the whole nebulae. A toy model shows that if just about 1--5 percent of the free electrons in a PN had a $\kappa$-EED as small as $\kappa=3$, it would be sufficient to account for the observations.


title: On variational principles of metric mean dimension on subset in Feldman-Katok metric

abstract: In this paper, we studied the metric mean dimension in Feldman-Katok(FK for short) metric. We introduced the notions of FK-Bowen metric mean dimension and FK-Packing metric mean dimension on subset. And we established two variational principles.


title: Electroweak precision tests for triplet scalars

abstract: Electroweak precision observables are fundamentally important for testing the standard model (SM) or its extensions. The influences to observables from new physics within the electroweak sector can be expressed in terms of oblique parameters S, T, U. The recently reported W mass excess anomaly by CDF modifies these parameters in a significant way. By performing the global fit with the CDF new W mass measurement data, we obtain $S=0.03 \pm 0.03$, $T=0.06 \pm 0.02$ and $U=0.16 \pm 0.03$ (or $S=0.14 \pm 0.03$, $T=0.24 \pm 0.02$ with $U=0$) which is significantly away from zero as SM would predict. The CDF excess strongly indicates the need of new physics beyond SM. We carry out a global fit to study the influence of two different cases of simple extensions by a hyper-charge $Y=1$ triplet scalar $\Delta$ (corresponding to the type-II seesaw model) and a $Y=0$ real triplet scalar $\Sigma$ on electroweak precision tests to determine parameter space in these models to solve the W mass anomaly and discuss the implications. We find that these triplets can affect the oblique parameters significantly at the tree and loop levels.


title: MAFNet: A Multi-Attention Fusion Network for RGB-T Crowd Counting

abstract: RGB-Thermal (RGB-T) crowd counting is a challenging task, which uses thermal images as complementary information to RGB images to deal with the decreased performance of unimodal RGB-based methods in scenes with low-illumination or similar backgrounds. Most existing methods propose well-designed structures for cross-modal fusion in RGB-T crowd counting. However, these methods have difficulty in encoding cross-modal contextual semantic information in RGB-T image pairs. Considering the aforementioned problem, we propose a two-stream RGB-T crowd counting network called Multi-Attention Fusion Network (MAFNet), which aims to fully capture long-range contextual information from the RGB and thermal modalities based on the attention mechanism. Specifically, in the encoder part, a Multi-Attention Fusion (MAF) module is embedded into different stages of the two modality-specific branches for cross-modal fusion at the global level. In addition, a Multi-modal Multi-scale Aggregation (MMA) regression head is introduced to make full use of the multi-scale and contextual information across modalities to generate high-quality crowd density maps. Extensive experiments on two popular datasets show that the proposed MAFNet is effective for RGB-T crowd counting and achieves the state-of-the-art performance.


title: Determination of the asymptotic limits of adaptive photon counting measurements for coherent-state optical phase estimation

abstract: Physical realizations of the canonical phase measurement for the optical phase are unknown. Single-shot phase estimation, which aims to determine the phase of an optical field in a single shot, is critical in quantum information processing and metrology. Here we present a family of strategies for single-shot phase estimation of coherent states based on adaptive non-Gaussian, photon counting, measurements with coherent displacements that maximize information gain as the measurement progresses, which have higher sensitivities over the best known adaptive Gaussian strategies. To gain understanding about their fundamental characteristics and demonstrate their superior performance, we develop a comprehensive statistical analysis based on the Bayesian optimal design of experiments, which provides a natural description of these non-Gaussian strategies. This mathematical framework, together with numerical analysis and Monte Carlo methods, allows us to determine the asymptotic limits in sensitivity of strategies based on photon counting designed to maximize information gain, which up to now had been a challenging problem. Moreover, we show that these non-Gaussian phase estimation strategies have the same functional form as the canonical phase measurement in the asymptotic limit differing only by a scaling factor, thus providing the highest sensitivity among physically-realizable measurements for single-shot phase estimation of coherent states known to date. This work shines light into the potential of optimized non-Gaussian measurements based on photon counting for optical quantum metrology and phase estimation.


title: Efficient quantum algorithms for solving quantum linear system problems

abstract: We transform the problem of solving linear system of equations $A\mathbf{x}=\mathbf{b}$ to a problem of finding the right singular vector with singular value zero of an augmented matrix $C$, and present two quantum algorithms for solving this problem. The first algorithm solves the problem directly by applying the quantum eigenstate filtering algorithm with query complexity of $O\left( s\kappa \log \left( 1/\epsilon \right) \right) $ for a $s$-sparse matrix $C$, where $\kappa $ is the condition number of the matrix $A$, and $\epsilon $ is the desired precision. The second algorithm uses the quantum resonant transition approach, the query complexity scales as $O\left[s\kappa + \log\left( 1/\epsilon \right)/\log \log \left( 1/\epsilon \right) \right] $. Both algorithms meet the optimal query complexity in $\kappa $, and are simpler than previous algorithms.


title: Triply-Resonant Sum Frequency Conversion with Gallium Phosphide Ring Resonators

abstract: We demonstrate quasi-phase matched, triply-resonant sum frequency conversion in 10.6-um-diameter integrated gallium phosphide ring resonators. A small-signal, waveguide-to-waveguide power conversion efficiency of 8%/mW is measured for conversion from telecom (1536 nm) and near infrared (1117 nm) to visible (647 nm) wavelengths with an absolute power conversion efficiency of 6.3% measured at saturation pump power. For the complementary difference frequency generation process, a single photon conversion efficiency of 7.2%/mW from visible to telecom is projected for resonators with optimized coupling. Efficient conversion from visible to telecom will facilitate long-distance transmission of spin-entangled photons from solid-state emitters such as the diamond NV center, allowing long-distance entanglement for quantum networks.


title: Extended-localized transition in diffusive quasicrystals

abstract: Compared to periodic systems, quasicrystals without translational invariance exhibit unexpected localization properties. The extended-localized transition in quasicrystals has been observed in both quantum and classical wave systems. However, its manifestation in diffusion systems, which serve as novel platforms for exploring phases of matter in condensed matter physics, remains unexplored. Here, we present the implementation of the extended-localized transition in a diffusive quasicrystal based on the coupled ring chain structure. By modulating the thermal conductivities of rings, we obtain the diffusive one-dimensional Aubry-Andr\'e-Harper (AAH) model, which exhibits an extended-localized transition. Thanks to the ring-shaped chain, we clearly demonstrate the extended-localized transition under the uniform excitation through temperature field simulations. For the localized state, the temperature field clearly demonstrates a multiple localization centers phenomenon, which has no counterpart in wave systems. We also quantitatively investigate the temperature evolution and size effect of this transition. Furthermore, the local excitation has been adopted to demonstrate the temperature field for both the extended and localized states. Besides, we implement the non-Hermitian diffusive AAH model by rotating rings, whose temperature field shows a moving multiple localization centers phenomenon in the localized phase. Finally, we give the experimental suggestions for the diffusive AAH model and propose a potential application named as double-trace distributed generator. Our results can facilitate the design of flexible thermal devices and efficient heat management.


title: A reconstruction method for binary limited-data tomography using a dictionary-based sparse shape recovery

abstract: Binary tomography is concerned with reconstructing a binary image from a very small number or other limited CT projection data. This problem itself not only possesses several medical imaging applications but also can be considered a model of general inverse problems to recover the object shape from limited measured data. Several approaches such as the Mumford-Shah method and various level-set methods have been investigated, but most of them lead to a non-convex optimization due to the difficulty to handle the binary constraint. We propose a new method based on a convex optimization inspired by dictionary-based shape recovery. In the proposed method, the object boundary of the binary image is represented by a level set of linear combinations of basis vectors in the dictionary. Using the dictionary, the object boundary is reconstructed by finding weights of the linear combination that best match the measured data. We create the dictionary by using the Gaussian radial basis function (GRBF). More concretely, we use Gaussian functions as a basis function placed at sparse grid points to represent the parametric level-set function and provide more flexibility in the binary representation of the reconstructed image. The simulation results of CT image reconstruction from only four projection data demonstrate that the proposed method can recover the object boundary more accurately compared with other competitive methods. The significance of our approach is the formulation with a tractable convex program while keeping moderate mathematical rigorousness.


title: IPvSeeYou: Exploiting Leaked Identifiers in IPv6 for Street-Level Geolocation

abstract: We present IPvSeeYou, a privacy attack that permits a remote and unprivileged adversary to physically geolocate many residential IPv6 hosts and networks with street-level precision. The crux of our method involves: 1) remotely discovering wide area (WAN) hardware MAC addresses from home routers; 2) correlating these MAC addresses with their WiFi BSSID counterparts of known location; and 3) extending coverage by associating devices connected to a common penultimate provider router. We first obtain a large corpus of MACs embedded in IPv6 addresses via high-speed network probing. These MAC addresses are effectively leaked up the protocol stack and largely represent WAN interfaces of residential routers, many of which are all-in-one devices that also provide WiFi. We develop a technique to statistically infer the mapping between a router's WAN and WiFi MAC addresses across manufacturers and devices, and mount a large-scale data fusion attack that correlates WAN MACs with WiFi BSSIDs available in wardriving (geolocation) databases. Using these correlations, we geolocate the IPv6 prefixes of $>$12M routers in the wild across 146 countries and territories. Selected validation confirms a median geolocation error of 39 meters. We then exploit technology and deployment constraints to extend the attack to a larger set of IPv6 residential routers by clustering and associating devices with a common penultimate provider router. While we responsibly disclosed our results to several manufacturers and providers, the ossified ecosystem of deployed residential cable and DSL routers suggests that our attack will remain a privacy threat into the foreseeable future.


title: Flow-Guided Transformer for Video Inpainting

abstract: We propose a flow-guided transformer, which innovatively leverage the motion discrepancy exposed by optical flows to instruct the attention retrieval in transformer for high fidelity video inpainting. More specially, we design a novel flow completion network to complete the corrupted flows by exploiting the relevant flow features in a local temporal window. With the completed flows, we propagate the content across video frames, and adopt the flow-guided transformer to synthesize the rest corrupted regions. We decouple transformers along temporal and spatial dimension, so that we can easily integrate the locally relevant completed flows to instruct spatial attention only. Furthermore, we design a flow-reweight module to precisely control the impact of completed flows on each spatial transformer. For the sake of efficiency, we introduce window partition strategy to both spatial and temporal transformers. Especially in spatial transformer, we design a dual perspective spatial MHSA, which integrates the global tokens to the window-based attention. Extensive experiments demonstrate the effectiveness of the proposed method qualitatively and quantitatively. Codes are available at https://github.com/hitachinsk/FGT.


title: Partial Differential Equations of Mixed Type: Analysis and Applications

abstract: Partial differential equations (PDEs) are at the heart of many mathematical and scientific advances. While great progress has been made on the theory of PDEs of standard types during the last eight decades, the analysis of nonlinear PDEs of mixed type is still in its infancy. The aim of this expository paper is to show, through several longstanding fundamental problems in fluid mechanics, differential geometry, and other areas, that many nonlinear PDEs arising in these areas are no longer of standard types, but lie at the boundaries of the classification of PDEs or, indeed, go beyond the classification to be of mixed type. Some interrelated connections, historical perspectives, recent developments, and current trends in the analysis of nonlinear PDEs of mixed type are also presented.

categories: math.AP math-ph math.DG math.FA math.MP physics.flu-dyn


title: Joint User Association and Resource Pricing for Metaverse: Distributed and Centralized Approaches

abstract: Metaverse as the next-generation Internet provides users with physical-virtual world interactions. To improve the quality of immersive experience, users access to Metaverse service providers (MSPs) and purchase bandwidth resource to reduce the communication latency of the Metaverse services. The MSPs decide selling price of the bandwidth resource to maximize the revenue. This leads to a joint user association and resource pricing problem between all users and MSPs. To tackle the problem, we formulate a Stackelberg game where the MSPs are game leaders and users are game followers. We resolve the Stackelberg equilibrium via the distributed and centralized approaches, according to different privacy requirements. In the distributed approach, the MSPs compete against each other to maximize the individual revenue, and a user selects an MSP in a probabilistic manner. The Stackelberg equilibrium is achieved in a privacy-friendly way. In the centralized approach, all MSPs and users accept the unified management and their strategies are instructed. The centralized approach acquires the superior decision-making performance but sacrifices the privacy of the game players. Finally, we provide numerical results to demonstrate the effectiveness and efficiency of our schemes.


title: A 100% Renewable Energy System: Enabling Zero CO2 Emission Offshore Platforms

abstract: The total electricity consumption from offshore oil/gas platforms is around 16 TWh worldwide in 2019. The majority offshore platforms are powered by the diesel generators while the rest mainly uses gas turbines, which emits large amounts of CO2 per year. The fast development of offshore wind turbines (WT) can potentially replace traditional fossil fuel based resources to power offshore loads. Thus, a novel offshore hybrid renewable energy sources (OHRES) system is proposed to enable a zero CO2 emission offshore platform mitigating climate change. Battery energy storage system (BESS) and hydrogen energy storage system (HESS) are considered to mitigate the fluctuation of wind power in the proposed OHRES system. Resilience models are designed to enhance the resilience of the proposed OHRES system with extra energy stored in BESS and/or HESS. Case studies demonstrate the feasibility of the proposed OHRES system to power offshore platforms. The economic analysis reports the planning cost for the proposed OHRES system under different resilience levels, which may benefit the decision to balance the carbon emission and investment cost.


title: Sum uncertainty relations based on $(\alpha,\beta,\gamma)$ weighted Wigner-Yanase-Dyson skew information

abstract: We introduce ($\alpha,\beta,\gamma$) weighted Wigner-Yanase-Dyson (($\alpha,\beta,\gamma$) WWYD) skew information and ($\alpha,\beta,\gamma$) modified weighted Wigner-Yanase-Dyson (($\alpha,\beta,\gamma$) MWWYD) skew information. We explore the sum uncertainty relations for arbitrary $N$ mutually noncommutative observables based on ($\alpha,\beta,\gamma$) WWYD skew information. A series of uncertainty inequalities are derived. We show by detailed example that our results cover and improve the previous ones based on the original Wigner-Yanase (WY) skew information. Finally, we establish new sum uncertainty relations in terms of the ($\alpha,\beta,\gamma$) MWWYD skew information for arbitrary $N$ quantum channels.


title: TL;DW? Summarizing Instructional Videos with Task Relevance & Cross-Modal Saliency

abstract: YouTube users looking for instructions for a specific task may spend a long time browsing content trying to find the right video that matches their needs. Creating a visual summary (abridged version of a video) provides viewers with a quick overview and massively reduces search time. In this work, we focus on summarizing instructional videos, an under-explored area of video summarization. In comparison to generic videos, instructional videos can be parsed into semantically meaningful segments that correspond to important steps of the demonstrated task. Existing video summarization datasets rely on manual frame-level annotations, making them subjective and limited in size. To overcome this, we first automatically generate pseudo summaries for a corpus of instructional videos by exploiting two key assumptions: (i) relevant steps are likely to appear in multiple videos of the same task (Task Relevance), and (ii) they are more likely to be described by the demonstrator verbally (Cross-Modal Saliency). We propose an instructional video summarization network that combines a context-aware temporal video encoder and a segment scoring transformer. Using pseudo summaries as weak supervision, our network constructs a visual summary for an instructional video given only video and transcribed speech. To evaluate our model, we collect a high-quality test set, WikiHow Summaries, by scraping WikiHow articles that contain video demonstrations and visual depictions of steps allowing us to obtain the ground-truth summaries. We outperform several baselines and a state-of-the-art video summarization model on this new benchmark.


title: Cryptanalyzing an Image Encryption Algorithm Underpinned by 2D Lag-Complex Logistic Map

abstract: This paper analyzes security performance of an image encryption algorithm using 2D lag-complex Logistic map (LCLM), which adopts it as a pseudo-random number generator, and uses the sum of all pixel values of the plain-image as its initial value to control the random combination of the basic encryption operations. However, multiple factors make the final pseudo-random sequences controlling the encryption process may be the same for different plain-images. Based on this point, we proposed a chosen-plaintext attack by attacking the six encryption steps with a strategy of divide and conquer. Using the pitfalls of 2D-LCLM, the number of required chosen plain-images is further reduced to $5\cdot\log_2(MN)+95$, where $\mathit{MN}$ is the number of pixels of the plain-image.


title: The high-resolution in vivo measurement of replication fork velocity and pausing by lag-time analysis

abstract: An important step towards understanding the mechanistic basis of the central dogma is the quantitative characterization of the dynamics of nucleic-acid-bound molecular motors in the context of the living cell, where a crowded cytoplasm as well as competing and potentially antagonistic processes may significantly affect their rapidity and reliability. To capture these dynamics, we develop a novel method, lag-time analysis, for measuring in vivo dynamics. The approach uses exponential growth as the stopwatch to resolve dynamics in an asynchronous culture and therefore circumvents the difficulties and potential artifacts associated with synchronization or fluorescent labeling. Although lag-time analysis has the potential to be widely applicable to the quantitative analysis of in vivo dynamics, we focus on an important application: characterizing replication dynamics. To benchmark the approach, we analyze replication dynamics in three different species and a collection of mutants. We provide the first quantitative locus-specific measurements of fork velocity, in units of kb per second, as well as replisome-pause durations, some with the precision of seconds. The measured fork velocity is observed to be both locus and time dependent, even in wild-type cells. In addition to quantitatively characterizing known phenomena, we detect brief, locus-specific pauses at rDNA in wild-type cells for the first time. We also observe temporal fork velocity oscillations in three highly-divergent bacterial species. Lag-time analysis not only has great potential to offer new insights into replication, as demonstrated in the paper, but also has potential to provide quantitative insights into other important processes.


title: Link-Backdoor: Backdoor Attack on Link Prediction via Node Injection

abstract: Link prediction, inferring the undiscovered or potential links of the graph, is widely applied in the real-world. By facilitating labeled links of the graph as the training data, numerous deep learning based link prediction methods have been studied, which have dominant prediction accuracy compared with non-deep methods. However,the threats of maliciously crafted training graph will leave a specific backdoor in the deep model, thus when some specific examples are fed into the model, it will make wrong prediction, defined as backdoor attack. It is an important aspect that has been overlooked in the current literature. In this paper, we prompt the concept of backdoor attack on link prediction, and propose Link-Backdoor to reveal the training vulnerability of the existing link prediction methods. Specifically, the Link-Backdoor combines the fake nodes with the nodes of the target link to form a trigger. Moreover, it optimizes the trigger by the gradient information from the target model. Consequently, the link prediction model trained on the backdoored dataset will predict the link with trigger to the target state. Extensive experiments on five benchmark datasets and five well-performing link prediction models demonstrate that the Link-Backdoor achieves the state-of-the-art attack success rate under both white-box (i.e., available of the target model parameter)and black-box (i.e., unavailable of the target model parameter) scenarios. Additionally, we testify the attack under defensive circumstance, and the results indicate that the Link-Backdoor still can construct successful attack on the well-performing link prediction methods. The code and data are available at https://github.com/Seaocn/Link-Backdoor.


title: An extension of the Fukaya-Kato method

abstract: In a groundbreaking paper, T. Fukaya and K. Kato proved a slight weakening of a conjecture of the author's relating modular symbols and cup products of cyclotomic units under an assumption that a Kubota-Leopoldt p-adic L-function has no multiple zeros. This article describes a refinement of their method that sheds light on the role of the p-adic L-function.


title: $\beta$-Divergence-Based Latent Factorization of Tensors model for QoS prediction

abstract: A nonnegative latent factorization of tensors (NLFT) model can well model the temporal pattern hidden in nonnegative quality-of-service (QoS) data for predicting the unobserved ones with high accuracy. However, existing NLFT models' objective function is based on Euclidean distance, which is only a special case of $\beta$-divergence. Hence, can we build a generalized NLFT model via adopting $\beta$-divergence to achieve prediction accuracy gain? To tackle this issue, this paper proposes a $\beta$-divergence-based NLFT model ($\beta$-NLFT). Its ideas are two-fold 1) building a learning objective with $\beta$-divergence to achieve higher prediction accuracy, and 2) implementing self-adaptation of hyper-parameters to improve practicability. Empirical studies on two dynamic QoS datasets demonstrate that compared with state-of-the-art models, the proposed $\beta$-NLFT model achieves the higher prediction accuracy for unobserved QoS data.


title: An outlook on the estimate of the solar quadrupole moment from relativistic gravitation contributions

abstract: Of all the solar fundamental parameters (mass, diameter, gravity at the surface,...), the gravitational moments have been quite often ignored in the past, mainly due to the great difficulty to get a reliable estimate. Even though the order of magnitude of the solar quadrupole moment $J_2$ is now known to be $10^{-7}$, its accurate value is still discussed. Indeed, the expansion in multipoles $J_{(l,~ l = 2, ...)}$ of the gravitational potential of a rotating body affects the orbital motion of planets at a relativistic level. We will recall here the recent progresses made in testing General Relativity through the contribution of the first solar quadrupole moment. Using the Eddington-Robertson parameters, we recall the constraints both on a theoretical and experimental point of view. Together with $\gamma$, which encodes the amount of curvature of space-time per unit rest-mass, the Post--Newtonian Parameter $\beta$ contributes to the relativistic precession of planets. The latter parameter encodes the amount of non-linearity in the superposition law of gravitation. Even though in principle, it would be possible to extract $J_2$ from planetary ephemerides, we observe that it is significantly correlated with other solution parameters (semi-major axis of planets, mass of asteroids...). Focusing on the $J_2$ correlations, we show that in general, when ~$\beta$ and ~$\gamma$ are freed, the correlations ~[$\beta, J_2$] and ~[$\gamma, J_2$] are $\approx$ 45\% and $\approx$ 55\% respectively. Moreover, all the planetary dynamics-based values are biased by the Lense--Thiring effect, which has never been modeled and solved for so far, but can be estimated to $\approx$ 7\%. It is thus possible to get a good estimate of the solar quadrupole moment:$1.66\times10^{-7}$$\leq$$J_2$$\leq$$2.32\times10^{-7}$.


title: Uncertainty of quantum channels via modified generalized variance and modified generalized Wigner-Yanase-Dyson skew information

abstract: Uncertainty relation is a fundamental issue in quantum mechanics and quantum information theory. By using modified generalized variance (MGV), and modified generalized Wigner-Yanase-Dyson skew information (MGWYD), we identify the total and quantum uncertainty of quantum channels. The elegant properties of the total uncertainty of quantum channels are explored in detail. In addition, we present a trade-off relation between the total uncertainty of quantum channels and the entanglement fidelity and establish the relationships between the total uncertainty and entropy exchange/coherent information. Detailed examples are given to the explicit formulas of the total uncertainty and the quantum uncertainty of quantum channels. Moreover, utilizing a realizable experimental measurement scheme by using the Mach-Zehnder interferometer proposed in Nirala et al. (Phys Rev A 99:022111, 2019), we discuss how to measure the total/quantum uncertainty of quantum channels for pure states.


title: Joint Channel Estimation and Data Detection for Hybrid RIS aided Millimeter Wave OTFS Systems

abstract: For high mobility communication scenario, the recently emerged orthogonal time frequency space (OTFS) modulation introduces a new delay-Doppler domain signal space, and can provide better communication performance than traditional orthogonal frequency division multiplexing system. This article focuses on the joint channel estimation and data detection (JCEDD) for hybrid reconfigurable intelligent surface (HRIS) aided millimeter wave (mmWave) OTFS systems. Firstly, a new transmission structure is designed. Within the pilot durations of the designed structure, partial HRIS elements are alternatively activated. The time domain channel model is then exhibited. Secondly, the received signal model for both the HRIS over time domain and the base station over delay-Doppler domain are studied. Thirdly, by utilizing channel parameters acquired at the HRIS, an HRIS beamforming design strategy is proposed. For the OTFS transmission, we propose a JCEDD scheme over delay-Doppler domain. In this scheme, message passing (MP) algorithm is designed to simultaneously obtain the equivalent channel gain and the data symbols. On the other hand, the channel parameters, i.e., the Doppler shift, the channel sparsity, and the channel variance, are updated through expectation-maximization (EM) algorithm. By iteratively executing the MP and EM algorithm, both the channel and the unknown data symbols can be accurately acquired. Finally, simulation results are provided to validate the effectiveness of our proposed JCEDD scheme.


title: Performance Analysis of Charging Infrastructure Sharing in UAV and EV-involved Networks

abstract: Electric vehicles (EVs) and unmanned aerial vehicles (UAVs) show great potential in modern transportation and communication networks, respectively. However, with growing demands for such technologies, the limited energy infrastructure becomes the bottleneck for their future growth. It might be of high cost and low energy efficiency for all the operators to each have their own dedicated energy infrastructure, such as charging stations. In this paper, we analyze a wireless charging infrastructure sharing strategy in UAV and EV-involved networks. We consider a scenario where UAVs can charge in EV charging stations and pay for the sharing fee. On the EVs' side, sharing infrastructure can earn extra profit but their service quality, such as waiting time, might slightly reduce. On the UAVs' side, if renting EV charging stations can achieve an acceptable system performance, say high coverage probability, while considering the cost, they may not need to build their dedicated charging stations. In this case, we use tools from stochastic geometry to model the locations and propose an optimization problem that captures the aforementioned trade-offs between cost or profit and quality of service. Our numerical results show that sharing infrastructure slightly increases the waiting time of EVs, say within $5$ min, but dramatically decreases the waiting time of drones, say more than $50$ min, and deploying more charging stations do achieve better performances, but all these better performances are expected to cost more.


title: Stabilizing unstable periodic orbit of unknown fractional-order systems via adaptive delayed feedback control

abstract: This article presents an adaptive nonlinear delayed feedback control scheme for stabilizing the unstable periodic orbit of unknown fractional-order chaotic systems. The proposed control framework uses the Lyapunov approach and sliding mode control technique to guarantee that the closed-loop system is asymptotically stable on a periodic trajectory sufficiently close to the unstable periodic orbit of the system. The proposed method has two significant advantages. First, it employs a direct adaptive control method, making it easy to implement this method on systems with unknown parameters. Second, the framework requires only the period of the unstable periodic orbit. The robustness of the closed-loop system against system uncertainties and external disturbances with unknown bounds is guaranteed. Simulations on fractional-order duffing and gyro systems are used to illustrate the effectiveness of the theoretical results. The simulation results demonstrate that our approach outperforms the previously developed linear feedback control method for stabilizing unstable periodic orbits in fractional-order chaotic systems, particularly in reducing steady-state error and achieving faster convergence of tracking error.


title: On the Gasca-Maeztu conjecture for $n=6$

abstract: A two-dimensional $n$-correct set is a set of nodes admitting unique bivariate interpolation with polynomials of total degree at most ~$n$. We are interested in correct sets with the property that all fundamental polynomials are products of linear factors. In 1982, M.~Gasca and J.~I.~Maeztu conjectured that any such set necessarily contains $n+1$ collinear nodes. So far, this had only been confirmed for $n\leq 5.$ In this paper, we make a step for proving the case $n=6.$


title: A probabilistic view on predictive constructions for Bayesian learning

abstract: Given a sequence $X=(X_1,X_2,\ldots)$ of random observations, a Bayesian forecaster aims to predict $X_{n+1}$ based on $(X_1,\ldots,X_n)$ for each $n\ge 0$. To this end, in principle, she only needs to select a collection $\sigma=(\sigma_0,\sigma_1,\ldots)$, called ``strategy" in what follows, where $\sigma_0(\cdot)=P(X_1\in\cdot)$ is the marginal distribution of $X_1$ and $\sigma_n(\cdot)=P(X_{n+1}\in\cdot\mid X_1,\ldots,X_n)$ the $n$-th predictive distribution. Because of the Ionescu-Tulcea theorem, $\sigma$ can be assigned directly, without passing through the usual prior/posterior scheme. One main advantage is that no prior probability is to be selected. In a nutshell, this is the predictive approach to Bayesian learning. A concise review of the latter is provided in this paper. We try to put such an approach in the right framework, to make clear a few misunderstandings, and to provide a unifying view. Some recent results are discussed as well. In addition, some new strategies are introduced and the corresponding distribution of the data sequence $X$ is determined. The strategies concern generalized P\'olya urns, random change points, covariates and stationary sequences.


title: Numerical simulations of a continuously injected relativistic electron beam relaxation into a plasma with large-scale density gradients

abstract: In this paper the influence of large-scale decreasing and increasing gradients of the density of magnetized plasma on the relaxation process of a continuously injected relativistic electron beam with an energy of 611 keV ($v_b=0.9c$) and a pitch-angle distribution is studied using particle-in-cell numerical simulations. It is found that for the selected parameters in the case of a smoothly decreasing gradient and in a homogeneous plasma the formation of spatially limited plasma oscillations of large amplitude occurs. In such cases, modulation instability develops and a long-wave longitudinal modulation of the ion density is formed. In addition, the large amplitude of plasma waves accelerates plasma electrons to energies on the order of the beam energy. In the case of increasing and sharply decreasing gradients, a significant decrease in the amplitude of plasma oscillations and the formation of a turbulent ion density spectrum are observed. The possibility of acceleration of beam electrons to energies more than 2 times higher than the initial energy of the beam particles is also demonstrated. This process takes place not only during beam propagation in growing plasma density, but also in homogeneous plasma due to interaction of beam particles with plasma oscillations of large amplitude.


title: HDR-Plenoxels: Self-Calibrating High Dynamic Range Radiance Fields

abstract: We propose high dynamic range (HDR) radiance fields, HDR-Plenoxels, that learn a plenoptic function of 3D HDR radiance fields, geometry information, and varying camera settings inherent in 2D low dynamic range (LDR) images. Our voxel-based volume rendering pipeline reconstructs HDR radiance fields with only multi-view LDR images taken from varying camera settings in an end-to-end manner and has a fast convergence speed. To deal with various cameras in real-world scenarios, we introduce a tone mapping module that models the digital in-camera imaging pipeline (ISP) and disentangles radiometric settings. Our tone mapping module allows us to render by controlling the radiometric settings of each novel view. Finally, we build a multi-view dataset with varying camera conditions, which fits our problem setting. Our experiments show that HDR-Plenoxels can express detail and high-quality HDR novel views from only LDR images with various cameras.


title: Blowups for a class of second order nonlinear hyperbolic equations: A reduced model of nonlinear Jeans instability

abstract: Due to the importance of understanding the formation of the nonlinear structures in the universe and stellar systems, the nonlinear Jeans instability which plays the key role in these formation problems is a long-standing open problem presented by astrophysics more than one hundred years. In this article, we focus on a reduced model of the nonlinear Jeans instability for the expanding Newtonian universe which is characterized by a class of second order nonlinear hyperbolic equations \begin{equation*} \Box \varrho(x^\mu) +\frac{\mathcal{a} }{t} \partial_{t}\varrho(x^\mu) - \frac{\mathcal{b}}{t^2} \varrho(x^\mu) (1+ \varrho(x^\mu) ) -\frac{\mathcal{c}-\mathcal{k}}{1+\varrho(x^\mu)} (\partial_{t}\varrho(x^\mu))^2= \mathcal{k} F(t). \end{equation*} We establish a family of nonlinear self-increase blowup solution (the solution itself becomes infinite, stable ODE-type blowup) to this equation and give estimates on the growth rate of $\varrho$ which may explain the much faster growth rate of the nonlinear structures of universe in astrophysical observations than the rate given by the classical Jeans instability.


title: Entanglement, CPT and neutral kaons

abstract: Direct tests of T, CP, CPT symmetries in transitions using entangled neutral kaons produced at a $\phi$-factory are briefly reviewed.


title: Hybrid explicit-implicit learning for multiscale problems with time dependent source

abstract: The splitting method is a powerful method for solving partial differential equations. Various splitting methods have been designed to separate different physics, nonlinearities, and so on. Recently, a new splitting approach has been proposed where some degrees of freedom are handled implicitly while other degrees of freedom are handled explicitly. As a result, the scheme contains two equations, one implicit and the other explicit. The stability of this approach has been studied. It was shown that the time step scales as the coarse spatial mesh size, which can provide a significant computational advantage. However, the implicit solution part can still be expensive, especially for nonlinear problems. In this paper, we introduce modified partial machine learning algorithms to replace the implicit solution part of the algorithm. These algorithms are first introduced in arXiv:2109.02147, where a homogeneous source term is considered along with the Transformer, which is a neural network that can predict future dynamics. In this paper, we consider time-dependent source terms which is a generalization of the previous work. Moreover, we use the whole history of the solution to train the network. As the implicit part of the equations is more complicated to solve, we design a neural network to predict it based on training. Furthermore, we compute the explicit part of the solution using our splitting strategy. In addition, we use Proper Orthogonal Decomposition based model reduction in machine learning. The machine learning algorithms provide computational saving without sacrificing accuracy. We present three numerical examples which show that our machine learning scheme is stable and accurate.


title: On the recent-$k$-record of discrete random variables

abstract: Let $X_1,~X_2,\cdots$ be a sequence of i.i.d random variables which are supposed to be observed in sequence. The $n$th value in the sequence is a $k-record~value$ if exactly $k$ of the first $n$ values (including $X_n$) are at least as large as it. Let ${\bf R}_k$ denote the ordered set of $k$-record values. The famous Ignatov's Theorem states that the random sets ${\bf R}_k(k=1,2,\cdots)$ are independent with common distribution. We introduce one new record named $recent-k-record$ (RkR in short) in this paper: $X_n$ is a $j$-RkR if there are exactly $j$ values at least as large as $X_n$ in $X_{n-k},~X_{n-k+1},\cdots,~X_{n-1}$. It turns out that RkR brings many interesting problems and some novel properties such as prediction rule and Poisson approximation which are proved in this paper. One application named "No Good Record" via the Lov{\'a}sz Local Lemma is also provided. We conclude this paper with some possible connection with scan statistics.


title: Improving Phishing Detection Via Psychological Trait Scoring

abstract: Phishing emails exhibit some unique psychological traits which are not present in legitimate emails. From empirical analysis and previous research, we find three psychological traits most dominant in Phishing emails - A Sense of Urgency, Inducing Fear by Threatening, and Enticement with Desire. We manually label 10% of all phishing emails in our training dataset for these three traits. We leverage that knowledge by training BERT, Sentence-BERT (SBERT), and Character-level-CNN models and capturing the nuances via the last layers that form the Phishing Psychological Trait (PPT) scores. For the phishing email detection task, we use the pretrained BERT and SBERT model, and concatenate the PPT scores to feed into a fully-connected neural network model. Our results show that the addition of PPT scores improves the model performance significantly, thus indicating the effectiveness of PPT scores in capturing the psychological nuances. Furthermore, to mitigate the effect of the imbalanced training dataset, we use the GPT-2 model to generate phishing emails (Radford et al., 2019). Our best model outperforms the current State-of-the-Art (SOTA) model's F1 score by 4.54%. Additionally, our analysis of individual PPTs suggests that Fear provides the strongest cue in detecting phishing emails.