title: Deep Learning of Phase Transitions for Quantum Spin Chains from Correlation Aspects

abstract: Using machine learning (ML) to recognize different phases of matter and to infer the entire phase diagram has proven to be an effective tool given a large dataset. In our previous proposals, we have successfully explored phase transitions for topological phases of matter at low dimensions either in a supervised or an unsupervised learning protocol with the assistance of quantum information related quantities. In this work, we adopt our previous ML procedures to study quantum phase transitions of magnetism systems such as the XY and XXZ spin chains by using spin-spin correlation functions as the input data. We find that our proposed approach not only maps out the phase diagrams with accurate phase boundaries, but also indicates some new features that have not observed before. In particular, we define so-called relevant correlation functions to some corresponding phases that can always distinguish between those and their neighbors. Based on the unsupervised learning protocol we proposed [Phys. Rev. B 104, 165108 (2021)], the reduced latent representations of the inputs combined with the clustering algorithm show the connectedness or disconnectedness between neighboring clusters (phases), just corresponding to the continuous or disrupt quantum phase transition, respectively.


title: Investigation of the laser-induced lineshape change in attosecond transient absorption spectra by employing a time-dependent generalized Floquet approach

abstract: We introduce a time-dependent generalized Floquet (TDGF) approach to calculate attosecond transient absorption spectra of helium atoms subjected to the combination of an attosecond extreme ultraviolet (XUV) pulse and a delayed few-cycle infrared (IR) laser pulse. This TDGF approach provides a Floquet understanding of the laser-induced change of resonant absorption lineshape. It is analytically demonstrated that, the phase shift of the time-dependent dipole moment that results in the lineshape changes consists of the \emph{adiabatic} laser-induced phase (LIP) due to the IR-induced stark shifts of adiabatic Floquet states and the \emph{non-adiabatic} phase correction due to the non-adiabatic IR-induced coupling between adiabatic Floquet states. Comparisons of the spectral lineshape calculated based on the TDGF approach with the results obtained with the LIP model [S. Chen \emph{et al.}, Phys. Rev. A \textbf{88}, 033409(2013)] and the rotating-wave approximation (RWA) are made in several typical cases. It is suggested in the picture of adiabatic Floquet states that, the LIP model works as long as the generalized adiabatic theorem [A. Dodin \emph{et al.}, Phys. Rev. X Quantum \textbf{2}, 030302(2021)] fulfils, and the RWA works when the higher-order IR-coupling effect in the formation of adiabatic Floquet states is neglectable.


title: Boosting Light Absorption of a Therapeutic Microcapsule by Means of Auxiliary Solid Nanoparticles

abstract: Multilayer microparticles with a liquid core and a polycomposite light-absorbing shell are important components of modern bio- and medical technologies, in particular, as transport microcontainers in the system of targeted delivery of therapeutic nanodoses to the desired region of the biological tissues. For reliably opening the microcapsule shell by an optical radiation and releasing the payload, it is necessary to dramatically increase the light absorption of such a microcontainer. To this end, we propose surrounding the microcapsule with specially added auxiliary nanoparticles, which can accumulate optical energy near its surface and direct a concentrated photonic flux to the target microcapsule thus leading to its booster heating. Using numerical finite-difference time-domain (FDTD) calculations, we simulate and examine the absorption dynamics of the near-infrared optical radiation in a spherical microcapsule surrounded by solid nanoparticles of different optical properties (metal, biocompatible dielectric). We show that due to light scattering on nanoparticles, the optical field superlocalization in the "hot regions" on the microcapsule surface occurs. The three-fold light absorption enhancement can be achieved due to the addition of nanoparticles.


title: Custom 8-bit floating point value format for reducing shared memory bank conflict in approximate nearest neighbor search

abstract: The k-nearest neighbor search is used in various applications such as machine learning, computer vision, database search, and information retrieval. While the computational cost of the exact nearest neighbor search is enormous, an approximate nearest neighbor search (ANNS) has been attracting much attention. IVFPQ is one of the ANNS methods. Although we can leverage the high bandwidth and low latency of shared memory to compute the search phase of the IVFPQ on NVIDIA GPUs, the throughput can degrade due to shared memory bank conflict. To reduce the bank conflict and improve the search throughput, we propose a custom 8-bit floating point value format. This format doesn't have a sign bit and can be converted from/to FP32 with a few instructions. We use this format for IVFPQ on GPUs and achieved better performance without significant recall loss compared to FP32 and FP16.


title: Multi Kernel Positional Embedding ConvNeXt for Polyp Segmentation

abstract: Medical image segmentation is the technique that helps doctor view and has a precise diagnosis, particularly in Colorectal Cancer. Specifically, with the increase in cases, the diagnosis and identification need to be faster and more accurate for many patients; in endoscopic images, the segmentation task has been vital to helping the doctor identify the position of the polyps or the ache in the system correctly. As a result, many efforts have been made to apply deep learning to automate polyp segmentation, mostly to ameliorate the U-shape structure. However, the simple skip connection scheme in UNet leads to deficient context information and the semantic gap between feature maps from the encoder and decoder. To deal with this problem, we propose a novel framework composed of ConvNeXt backbone and Multi Kernel Positional Embedding block. Thanks to the suggested module, our method can attain better accuracy and generalization in the polyps segmentation task. Extensive experiments show that our model achieves the Dice coefficient of 0.8818 and the IOU score of 0.8163 on the Kvasir-SEG dataset. Furthermore, on various datasets, we make competitive achievement results with other previous state-of-the-art methods.


title: Multi-fidelity surrogate modeling for temperature field prediction using deep convolution neural network

abstract: Temperature field prediction is of great importance in the thermal design of systems engineering, and building the surrogate model is an effective way for the task. Generally, large amounts of labeled data are required to guarantee a good prediction performance of the surrogate model, especially the deep learning model, which have more parameters and better representational ability. However, labeled data, especially high-fidelity labeled data, are usually expensive to obtain and sometimes even impossible. To solve this problem, this paper proposes a pithy deep multi-fidelity model (DMFM) for temperature field prediction, which takes advantage of low-fidelity data to boost the performance with less high-fidelity data. First, a pre-train and fine-tune paradigm are developed in DMFM to train the low-fidelity and high-fidelity data, which significantly reduces the complexity of the deep surrogate model. Then, a self-supervised learning method for training the physics-driven deep multi-fidelity model (PD-DMFM) is proposed, which fully utilizes the physics characteristics of the engineering systems and reduces the dependence on large amounts of labeled low-fidelity data in the training process. Two diverse temperature field prediction problems are constructed to validate the effectiveness of DMFM and PD-DMFM, and the result shows that the proposed method can greatly reduce the dependence of the model on high-fidelity data.


title: Artificial intelligence as a gateway to scientific discovery: Uncovering features in retinal fundus images

abstract: Purpose: Convolutional neural networks can be trained to detect various conditions or patient traits based on retinal fundus photographs, some of which, such as the patient sex, are invisible to the expert human eye. Here we propose a methodology for explainable classification of fundus images to uncover the mechanism(s) by which CNNs successfully predict the labels. We used patient sex as a case study to validate our proposed methodology. Approach: First, we used a set of 4746 fundus images, including training, validation and test partitions, to fine-tune a pre-trained CNN on the sex classification task. Next, we utilized deep learning explainability tools to hypothesize possible ways sex differences in the retina manifest. We measured numerous retinal properties relevant to our hypotheses through image segmentation to identify those significantly different between males and females. To tackle the multiple comparisons problem, we shortlisted the parameters by testing them on a set of 100 fundus images distinct from the images used for fine-tuning. Finally, we used an additional 400 images, not included in any previous set, to reveal significant sex differences in the retina. Results: We observed that the peripapillary area is darker in males compared to females ($p=.023, d=.243$). We also observed that males have richer retinal vasculature networks by showing a higher number of branches ($p=.016, d=.272$) and nodes ($p=.014, d=.299$) and a larger total length of branches ($p=.045, d=.206$) in the vessel graph. Also, vessels cover a greater area in the superior temporal quadrant of the retina in males compared to females ($p=0.048, d=.194$). Conclusions: Our methodology reveals retinal features in fundus photographs that allow CNNs to predict traits currently unknown, but meaningful to experts.


title: Explainable, Interpretable & Trustworthy AI for Intelligent Digital Twin: Case Study on Remaining Useful Life

abstract: Artificial intelligence (AI) and Machine learning (ML) are increasingly used in energy and engineering systems, but these models must be fair, unbiased, and explainable. It is critical to have confidence in AI's trustworthiness. ML techniques have been useful in predicting important parameters and in improving model performance. However, for these AI techniques to be useful for making decisions, they need to be audited, accounted for, and easy to understand. Therefore, the use of explainable AI (XAI) and interpretable machine learning (IML) is crucial for the accurate prediction of prognostics, such as remaining useful life (RUL), in a digital twin system, to make it intelligent while ensuring that the AI model is transparent in its decision-making processes and that the predictions it generates can be understood and trusted by users. By using AI that is explainable, interpretable, and trustworthy, intelligent digital twin systems can make more accurate predictions of RUL, leading to better maintenance and repair planning, and ultimately, improved system performance. The objective of this paper is to explain the ideas of XAI and IML and to justify the important role of AI/ML in the digital twin framework and components, which requires XAI to understand the prediction better. This paper explains the importance of XAI and IML in both local and global aspects to ensure the use of trustworthy AI/ML applications for RUL prediction. We used the RUL prediction for the XAI and IML studies and leveraged the integrated Python toolbox for interpretable machine learning~(PiML).


title: Towards $21$-cm intensity mapping at $z=2.28$ with uGMRT using the tapered gridded estimator II: Cross-polarization power spectrum

abstract: Neutral hydrogen ($\rm{HI}$) $21$-cm intensity mapping (IM) offers an efficient technique for mapping the large-scale structures in the universe. We introduce the 'Cross' Tapered Gridded Estimator (Cross TGE), which cross-correlates two cross-polarizations (RR and LL) to estimate the multi-frequency angular power spectrum (MAPS) $C_{\ell}(\Delta\nu)$. We expect this to mitigate several effects like noise bias, calibration errors etc., which affect the 'Total' TGE which combines the two polarizations. Here we apply the Cross TGE on a $24.4 \,\rm{MHz}$ bandwidth uGMRT Band $3$ data centred at $432.8 \,\rm{MHz}$ aiming $\rm{HI}$ IM at $z=2.28$. The measured $C_{\ell}(\Delta\nu)$ is modelled to yield maximum likelihood estimates of the foregrounds and the spherical power spectrum $P(k)$ in several $k$ bins. Considering the mean squared brightness temperature fluctuations, we report a $2\sigma$ upper limit $\Delta_{UL}^{2}(k) \le (58.67)^{2} \, {\rm mK}^{2}$ at $k=0.804 \, {\rm Mpc}^{-1}$ which is a factor of $5.2$ improvement on our previous estimate based on the Total TGE. Assuming that the $\rm{HI}$ traces the underlying matter distribution, we have modelled $C_{\ell}(\Delta\nu)$ to simultaneously estimate the foregrounds and $[\Omega_{\rm{HI}} b_{\rm{HI}}] $ where $\Omega_{\rm{HI}}$ and $b_{\rm{HI}}$ are the $\rm{HI}$ density and linear bias parameters respectively. We obtain a best fit value of $[\Omega_{\rm{HI}}b_{\rm{HI}}]^2 = 7.51\times 10^{-4} \pm 1.47\times 10^{-3}$ which is consistent with noise. Although the $2\sigma$ upper limit $[\Omega_{\rm{HI}}b_{\rm{HI}}]_{UL} \leq 0.061$ is $\sim 50$ times larger than the expected value, this is a considerable improvement over earlier works at this redshift.


title: Feature-based Image Matching for Identifying Individual K\=ak\=a

abstract: This report investigates an unsupervised, feature-based image matching pipeline for the novel application of identifying individual k\=ak\=a. Applied with a similarity network for clustering, this addresses a weakness of current supervised approaches to identifying individual birds which struggle to handle the introduction of new individuals to the population. Our approach uses object localisation to locate k\=ak\=a within images and then extracts local features that are invariant to rotation and scale. These features are matched between images with nearest neighbour matching techniques and mismatch removal to produce a similarity score for image match comparison. The results show that matches obtained via the image matching pipeline achieve high accuracy of true matches. We conclude that feature-based image matching could be used with a similarity network to provide a viable alternative to existing supervised approaches.


title: Rethinking Lightweight Salient Object Detection via Network Depth-Width Tradeoff

abstract: Existing salient object detection methods often adopt deeper and wider networks for better performance, resulting in heavy computational burden and slow inference speed. This inspires us to rethink saliency detection to achieve a favorable balance between efficiency and accuracy. To this end, we design a lightweight framework while maintaining satisfying competitive accuracy. Specifically, we propose a novel trilateral decoder framework by decoupling the U-shape structure into three complementary branches, which are devised to confront the dilution of semantic context, loss of spatial structure and absence of boundary detail, respectively. Along with the fusion of three branches, the coarse segmentation results are gradually refined in structure details and boundary quality. Without adding additional learnable parameters, we further propose Scale-Adaptive Pooling Module to obtain multi-scale receptive filed. In particular, on the premise of inheriting this framework, we rethink the relationship among accuracy, parameters and speed via network depth-width tradeoff. With these insightful considerations, we comprehensively design shallower and narrower models to explore the maximum potential of lightweight SOD. Our models are purposed for different application environments: 1) a tiny version CTD-S (1.7M, 125FPS) for resource constrained devices, 2) a fast version CTD-M (12.6M, 158FPS) for speed-demanding scenarios, 3) a standard version CTD-L (26.5M, 84FPS) for high-performance platforms. Extensive experiments validate the superiority of our method, which achieves better efficiency-accuracy balance across five benchmarks.


title: DIGITOUR: Automatic Digital Tours for Real-Estate Properties

abstract: A virtual or digital tour is a form of virtual reality technology which allows a user to experience a specific location remotely. Currently, these virtual tours are created by following a 2-step strategy. First, a photographer clicks a 360 degree equirectangular image; then, a team of annotators manually links these images for the "walkthrough" user experience. The major challenge in the mass adoption of virtual tours is the time and cost involved in manual annotation/linking of images. Therefore, this paper presents an end-to-end pipeline to automate the generation of 3D virtual tours using equirectangular images for real-estate properties. We propose a novel HSV-based coloring scheme for paper tags that need to be placed at different locations before clicking the equirectangular images using 360 degree cameras. These tags have two characteristics: i) they are numbered to help the photographer for placement of tags in sequence and; ii) bi-colored, which allows better learning of tag detection (using YOLOv5 architecture) in an image and digit recognition (using custom MobileNet architecture) tasks. Finally, we link/connect all the equirectangular images based on detected tags. We show the efficiency of the proposed pipeline on a real-world equirectangular image dataset collected from the Housing.com database.


title: Cross-domain Self-supervised Framework for Photoacoustic Computed Tomography Image Reconstruction

abstract: Accurate image reconstruction is crucial for photoacoustic (PA) computed tomography (PACT). Recently, deep learning has been used to reconstruct the PA image with a supervised scheme, which requires high-quality images as ground truth labels. In practice, there are inevitable trade-offs between cost and performance since the use of more channels is an expensive strategy to access more measurements. Here, we propose a cross-domain unsupervised reconstruction (CDUR) strategy with a pure transformer model, which overcomes the lack of ground truth labels from limited PA measurements. The proposed approach exploits the equivariance of PACT to achieve high performance with a smaller number of channels. We implement a self-supervised reconstruction in a model-based form. Meanwhile, we also leverage the self-supervision to enforce the measurement and image consistency on three partitions of measured PA data, by randomly masking different channels. We find that dynamically masking a high proportion of the channels, e.g., 80%, yields nontrivial self-supervisors in both image and signal domains, which decrease the multiplicity of the pseudo solution to efficiently reconstruct the image from fewer PA measurements with minimum error of the image. Experimental results on in-vivo PACT dataset of mice demonstrate the potential of our unsupervised framework. In addition, our method shows a high performance (0.83 structural similarity index (SSIM) in the extreme sparse case with 13 channels), which is close to that of supervised scheme (0.77 SSIM with 16 channels). On top of all the advantages, our method may be deployed on different trainable models in an end-to-end manner.


title: The structure of heavily doped impurity band in crystalline host

abstract: We study the properties of the impurity band in heavily-doped non-magnetic semiconductors using the Jacobi-Davidson algorithm and the supervised deep learning method. The disorder averaged inverse participation ratio (IPR) and thouless number calculation show us the rich structure inside the impurity band. A Convolutional Neural Network(CNN) model, which is trained to distinguish the extended/localized phase of the Anderson model with high accuracy, shows us the results in good agreement with the conventional approach. Together, we find that there are three mobility edges in the impurity band for a specific on-site impurity potential, which means the presence of the extended states while filling the impurity band.


title: Surgical Aggregation: Federated Class-Heterogeneous Learning

abstract: The release of numerous chest x-ray datasets has spearheaded the development of deep learning models with expert-level performance. However, they have limited interoperability due to class-heterogeneity -- a result of inconsistent labeling schemes and partial annotations. Therefore, it is challenging to leverage these datasets in aggregate to train models with a complete representation of abnormalities that may occur within the thorax. In this work, we propose surgical aggregation, a federated learning framework for aggregating knowledge from class-heterogeneous datasets and learn a model that can simultaneously predict the presence of all disease labels present across the datasets. We evaluate our method using simulated and real-world class-heterogeneous datasets across both independent and identically distributed (iid) and non-iid settings. Our results show that surgical aggregation outperforms current methods, has better generalizability, and is a crucial first step towards tackling class-heterogeneity in federated learning to facilitate the development of clinically-useful models using previously non-interoperable chest x-ray datasets.


title: Co-analytic Counterexamples to Marstrand's Projection Theorem

abstract: Assuming $V=L$, we construct a plane set $E$ of Hausdorff dimension $1$ whose every orthogonal projection onto straight lines through the origin has Hausdorff dimension $0$. This is a counterexample to J. M. Marstrand's seminal projection theorem. While counterexamples had already been constructed decades ago, initially by R. O. Davies, the novelty of our result lies in the fact that $E$ is co-analytic. Following Marstrand's original proof (and R. Kaufman's newer, and now standard, approach based on capacities), a counterexample to the projection theorem cannot be analytic, hence our counterexample is optimal. We then extend the result in a strong way: we show that for each $\epsilon \in (0, 1)$ there exists a co-analytic set $E_{\epsilon}$ of dimension $1 + \epsilon$, each of whose orthogonal projections onto straight lines through the origin has Hausdorff dimension $\epsilon$. The constructions of $E$ and $E_{\epsilon}$ are by induction on the countable ordinals, applying a theorem by Z. Vidnyanszky.


title: Distribution Aligned Feature Clustering for Zero-Shot Sketch-Based Image Retrieval

abstract: Zero-Shot Sketch-Based Image Retrieval (ZS-SBIR) is a challenging cross-modal retrieval task. In prior arts, the retrieval is conducted by sorting the distance between the query sketch and each image in the gallery. However, the domain gap and the zero-shot setting make neural networks hard to generalize. This paper tackles the challenges from a new perspective: utilizing gallery image features. We propose a Cluster-then-Retrieve (ClusterRetri) method that performs clustering on the gallery images and uses the cluster centroids as proxies for retrieval. Furthermore, a distribution alignment loss is proposed to align the image and sketch features with a common Gaussian distribution, reducing the domain gap. Despite its simplicity, our proposed method outperforms the state-of-the-art methods by a large margin on popular datasets, e.g., up to 31% and 39% relative improvement of mAP@all on the Sketchy and TU-Berlin datasets.


title: Well-posedness and convergence analysis of PML method for time-dependent acoustic scattering problems over a locally rough surface

abstract: We aim to analyze and calculate time-dependent acoustic wave scattering by a bounded obstacle and a locally perturbed non-selfintersecting curve. The scattering problem is equivalently reformulated as an initial-boundary value problem of the wave equation in a truncated bounded domain through a well-defined transparent boundary condition. Well-posedness and stability of the reduced problem are established. Numerically, we adopt the perfect matched layer (PML) scheme for simulating the propagation of perturbed waves. By designing a special absorbing medium in a semi-circular PML, we show well-posedness and stability of the truncated initial-boundary value problem. Finally, we prove that the PML solution converges exponentially to the exact solution in the physical domain. Numerical results are reported to verify the exponential convergence with respect to absorbing medium parameters and thickness of the PML.


title: DQNAS: Neural Architecture Search using Reinforcement Learning

abstract: Convolutional Neural Networks have been used in a variety of image related applications after their rise in popularity due to ImageNet competition. Convolutional Neural Networks have shown remarkable results in applications including face recognition, moving target detection and tracking, classification of food based on the calorie content and many more. Designing of Convolutional Neural Networks requires experts having a cross domain knowledge and it is laborious, which requires a lot of time for testing different values for different hyperparameter along with the consideration of different configurations of existing architectures. Neural Architecture Search is an automated way of generating Neural Network architectures which saves researchers from all the brute-force testing trouble, but with the drawback of consuming a lot of computational resources for a prolonged period. In this paper, we propose an automated Neural Architecture Search framework DQNAS, guided by the principles of Reinforcement Learning along with One-shot Training which aims to generate neural network architectures that show superior performance and have minimum scalability problem.


title: Ergodicity in the dynamics of holomorphic correspondences

abstract: This paper studies ergodic properties of certain measures arising in the dynamics of holomorphic correspondences. These measures, in general, are not invariant in the classical sense of ergodic theory. We define a notion of ergodicity, and prove a version of Birkhoff's ergodic theorem in this setting. In fact, we strengthen this classical result in the setting of rational maps on the Riemann sphere with the Lyubich measure. We also show the existence of ergodic measures when a holomorphic correspondence is defined on a compact complex manifold. Lastly, we give an explicit class of dynamically interesting measures that are ergodic as in our definition.


title: Ember-IO: Effective Firmware Fuzzing with Model-Free Memory Mapped IO

abstract: Exponential growth in embedded systems is driving the research imperative to develop fuzzers to automate firmware testing to uncover software bugs and security vulnerabilities. But, employing fuzzing techniques in this context present a uniquely challenging proposition; a key problem is the need to deal with the diverse and large number of peripheral communications in an automated testing framework. Recent fuzzing approaches: i) employ re-hosting methods by executing code in an emulator because fuzzing on resource limited embedded systems is slow and unscalable; and ii) integrate models of hardware behaviour to overcome the challenges faced by the massive input-space to be explored created by peripheral devices and to generate inputs that are effective in aiding a fuzzer to make progress. Our efforts expounds upon program execution behaviours unique to firmware to address the resulting input-space search problem. The techniques we propose improve the fuzzer's ability to generate values likely to progress execution and avoids time consumed on mutating inputs that are functionally equivalent to other test cases. We demonstrate the methods are highly efficient and effective at overcoming the input-space search problem. Our emulation-based implementation, Ember-IO, when compared to the existing state-of-the-art fuzzing framework across 21 firmware binaries, demonstrates up to 255% improvement in blocks covered. Further Ember-IO discovered 6 new bugs in the real-world firmware, previously not identified by state-of-the-art fuzzing frameworks. Importantly, Ember-IO integrated with the state-of-the-art fuzzer, Fuzzware, demonstrates similar or improved coverage across all firmware binaries whilst reproducing 3 of the 6 new bugs discovered by Ember-IO.


title: Audio2Gestures: Generating Diverse Gestures from Audio

abstract: People may perform diverse gestures affected by various mental and physical factors when speaking the same sentences. This inherent one-to-many relationship makes co-speech gesture generation from audio particularly challenging. Conventional CNNs/RNNs assume one-to-one mapping, and thus tend to predict the average of all possible target motions, easily resulting in plain/boring motions during inference. So we propose to explicitly model the one-to-many audio-to-motion mapping by splitting the cross-modal latent code into shared code and motion-specific code. The shared code is expected to be responsible for the motion component that is more correlated to the audio while the motion-specific code is expected to capture diverse motion information that is more independent of the audio. However, splitting the latent code into two parts poses extra training difficulties. Several crucial training losses/strategies, including relaxed motion loss, bicycle constraint, and diversity loss, are designed to better train the VAE. Experiments on both 3D and 2D motion datasets verify that our method generates more realistic and diverse motions than previous state-of-the-art methods, quantitatively and qualitatively. Besides, our formulation is compatible with discrete cosine transformation (DCT) modeling and other popular backbones (\textit{i.e.} RNN, Transformer). As for motion losses and quantitative motion evaluation, we find structured losses/metrics (\textit{e.g.} STFT) that consider temporal and/or spatial context complement the most commonly used point-wise losses (\textit{e.g.} PCK), resulting in better motion dynamics and more nuanced motion details. Finally, we demonstrate that our method can be readily used to generate motion sequences with user-specified motion clips on the timeline.


title: Circumventing the Phonon Bottleneck by Multiphonon-Mediated Hot Exciton Cooling at the Nanoscale

abstract: In semiconductor materials, hot exciton cooling is the process by which highly excited carriers nonradiatively relax to form a band edge exciton. While cooling plays an important role in determining the thermal losses and quantum yield of a system, the timescales and mechanism of cooling are not well understood in confined semiconductor nanocrystals (NCs). A mismatch between electronic energy gaps and phonon frequencies in NCs has led to the hypothesis of a phonon bottleneck, in which cooling would be extremely slow, while enhanced electron-hole interactions in NCs have been used to explain cooling that would occur on ultrafast timescales. Here, we develop an atomistic approach for describing phonon-mediated exciton dynamics, and we use it to simulate hot exciton cooling in NCs of experimentally relevant sizes. Our framework includes electron-hole correlations as well as multiphonon processes, both of which are necessary to accurately describe the cooling process. We find that cooling occurs on timescales of tens of femtoseconds in CdSe cores, in agreement with experimental measurements, through a cascade of relaxation events that are mediated by efficient multiphonon emission. Cooling timescales increase with increasing NC size due to decreased exciton-phonon coupling (EXPC), and they are an order of magnitude larger in CdSe-CdS core-shell NCs because of reduced EXPC to low- and mid-frequency acoustic modes.


title: Using Fano factors to determine certain types of gene autoregulation

abstract: The expression of one gene might be regulated by its corresponding protein, which is called autoregulation. Although gene regulation is a central topic in biology, autoregulation is much less studied. In general, it is extremely difficult to determine the existence of autoregulation with direct biochemical approaches. Nevertheless, some papers have observed that certain types of autoregulations are linked to noise levels in gene expression. We generalize these results by two propositions on discrete-state continuous-time Markov chains. These two propositions form a simple but robust method to infer the existence of autoregulation in certain scenarios from gene expression data. This method only depends on the Fano factor, namely the ratio of variance and mean of the gene expression level. Compared to other methods for inferring autoregulation, our method only requires non-interventional one-time data, and does not need to estimate parameters. Besides, our method has few restrictions on the model. We apply this method to four groups of experimental data and find some genes that might have autoregulation. Some inferred autoregulations have been verified by experiments or other theoretical works.


title: Representations via differential algebras and equationally Noetherian algebras

abstract: We show that free algebras of the variety of algebras generated by the Witt algebra $W_n$, the left-symmetric Witt algebra $L_n$, and the symplectic Poisson algebra $P_n$ can be described as subalgebras of differential polynomial algebras with respect to appropriately defined products. Using these representations, we prove that $W_n$, $L_n$, $P_n$, and the free algebras of the varieties of algebras generated by these algebras are equationally Noetherian.


title: Photonic Characterization of Oxygen and Air-Annealed Zn3N2 Thin Films

abstract: Zinc nitride films were synthesized by reactive radio frequency (rf) magnetron sputtering of a zinc target using an Ar+N2 mixture. The as-deposited films were annealed in the air and O2 at 300 {\deg}C for 1 hr. The XRD measurements indicated that the films had a polycrystalline structure with a preferred (400) Zn3N2 orientation. The annealing process enhanced the crystallinity. After annealing, the AFM and SEM morphology revealed no significant change in the surface morphology and surface roughness. The direct bandgap of Zn3N2 was estimated to be in the range of 1.15 -1.35 eV where annealing resulted in a reduction of the bandgap. The films were confirmed to be p-type conduction and the resistivity was slightly increased by annealing. The photoconductivity measurements indicated that the as-deposited films did not have any photoresponse, whereas the annealed films exhibited photoconductivity.


title: Quantifying and Managing Impacts of Concept Drifts on IoT Traffic Inference in Residential ISP Networks

abstract: Millions of vulnerable consumer IoT devices in home networks are the enabler for cyber crimes putting user privacy and Internet security at risk. Internet service providers (ISPs) are best poised to play key roles in mitigating risks by automatically inferring active IoT devices per household and notifying users of vulnerable ones. Developing a scalable inference method that can perform robustly across thousands of home networks is a non-trivial task. This paper focuses on the challenges of developing and applying data-driven inference models when labeled data of device behaviors is limited and the distribution of data changes (concept drift) across time and space domains. Our contributions are three-fold: (1) We collect and analyze network traffic of 24 types of consumer IoT devices from 12 real homes over six weeks to highlight the challenge of temporal and spatial concept drifts in network behavior of IoT devices; (2) We analyze the performance of two inference strategies, namely "global inference" (a model trained on a combined set of all labeled data from training homes) and "contextualized inference" (several models each trained on the labeled data from a training home) in the presence of concept drifts; and (3) To manage concept drifts, we develop a method that dynamically applies the ``closest'' model (from a set) to network traffic of unseen homes during the testing phase, yielding better performance in 20% of scenarios.


title: Excitation and Ionization Properties of Star-forming Galaxies at z=2.0-9.3 with JWST/NIRSpec

abstract: We utilize medium-resolution JWST/NIRSpec observations of 164 galaxies at $z=2.0-9.3$ from the Cosmic Evolution Early Release Science (CEERS) survey to investigate the evolution of the excitation and ionization properties of galaxies at high redshifts. Our results represent the first statistical constraints on the evolution of the [OIII]/H$\beta$ vs. [NII]/H$\alpha$, [SII]/H$\alpha$, and [OI]/H$\alpha$ ``BPT'' diagrams at $z>2.7$, and the first analysis of the O32 vs. R23 diagram at $z>4$ with a large sample. We divide the sample into five redshift bins containing 30-40 galaxies each. The subsamples at $z\sim2.3$, $z\sim3.3$, and $z\sim4.5$ are representative of the main-sequence star-forming galaxy population at these redshifts, while the $z\sim5.6$ and $z\sim7.5$ samples are likely biased toward high specific star-formation rate due to selection effects. Using composite spectra, we find that each subsample at $z=2.0-6.5$ falls on the same excitation sequence in the [NII] and [SII] BPT diagrams and the O32-R23 diagram on average, offset from the sequences followed by $z=0$ HII regions in the same diagrams. The direction of these offsets are consistent with high-redshift star-forming galaxies uniformly having harder ionizing spectra than typical local galaxies at fixed nebular metallicity. The similarity of the average line ratios suggests that the ionization conditions of the interstellar medium do not strongly evolve between $z\sim2$ and $z\sim6$. Overall, the rest-optical line ratios suggest the $z=2.7-9.3$ CEERS/NIRSpec galaxies at log($M_*/M_{\odot})\sim7.5-10$ have high degrees of ionization and moderately low oxygen abundances ($\sim0.1-0.3~Z_{\odot}$), but are not extremely metal poor ($<0.1~Z_{\odot}$) even at $z>6.5$.


title: Estimation of Policy-Relevant Causal Effects in the Presence of Interference with an Application to the Philadelphia Beverage Tax

abstract: To comprehensively evaluate a public policy intervention, researchers must consider the effects of the policy not just on the implementing region, but also nearby, indirectly-affected regions. For example, an excise tax on sweetened beverages in Philadelphia was shown to not only be associated with a decrease in volume sales of taxed beverages in Philadelphia, but also an increase in sales in bordering counties not subject to the tax. The latter association may be explained by cross-border shopping behaviors of Philadelphia residents and indicate a causal effect of the tax on nearby regions, which may offset the total effect of the intervention. To estimate causal effects in this setting, we extend difference-in-differences methodology to account for such interference between regions and adjust for potential confounding present in quasi-experimental evaluations. Our doubly robust estimators for the average treatment effect on the treated and neighboring control relax standard assumptions on interference and model specification. We apply these methods to evaluate the change in volume sales of taxed beverages in 231 Philadelphia and bordering county stores due to the Philadelphia beverage tax. We also use our methods to explore the heterogeneity of effects across geographic features.


title: Tactile Tool Manipulation

abstract: Humans can effortlessly perform very complex, dexterous manipulation tasks by reacting to sensor observations. In contrast, robots can not perform reactive manipulation and they mostly operate in open-loop while interacting with their environment. Consequently, the current manipulation algorithms either are inefficient in performance or can only work in highly structured environments. In this paper, we present closed-loop control of a complex manipulation task where a robot uses a tool to interact with objects. Manipulation using a tool leads to complex kinematics and contact constraints that need to be satisfied for generating feasible manipulation trajectories. We first present an open-loop controller design using Non-Linear Programming (NLP) that satisfies these constraints. In order to design a closed-loop controller, we present a pose estimator of objects and tools using tactile sensors. Using our tactile estimator, we design a closed-loop controller based on Model Predictive Control (MPC). The proposed algorithm is verified using a 6 DoF manipulator on tasks using a variety of objects and tools. We verify that our closed-loop controller can successfully perform tool manipulation under several unexpected contacts. Video summarizing this work and hardware experiments are found https://youtu.be/VsClK04qDhk.


title: Self-Tuning Network Control Architectures

abstract: We formulate a general mathematical framework for self-tuning network control architecture design. This problem involves jointly adapting the locations of active sensors and actuators in the network and the feedback control policy to all available information about the time-varying network state and dynamics to optimize a performance criterion. We propose a general solution structure analogous to the classical self-tuning regulator from adaptive control. We show that a special case with full-state feedback can be solved in principle with dynamic programming, and in the linear quadratic setting the optimal cost functions and policies are piecewise quadratic and piecewise linear, respectively. For large networks where exhaustive architecture search is prohibitive, we describe a greedy heuristic for joint architecture-policy design. We demonstrate in numerical experiments that self-tuning architectures can provide dramatically improved performance over fixed architectures. Our general formulation provides an extremely rich and challenging problem space with opportunities to apply a wide variety of approximation methods from stochastic control, system identification, reinforcement learning, and static architecture design.


title: Conformal flatness of compact three-dimensional Cotton-parallel manifolds

abstract: A three-dimensional pseudo-Riemannian manifold is called essentially conformally symmetric (ECS) if its Cotton tensor is parallel but nowhere-vanishing. In this note we prove that three-dimensional ECS manifolds must be noncompact or, equivalently, that every compact three-dimensional Cotton-parallel pseudo-Riemannian manifold must be conformally flat.


title: Improved generalization with deep neural operators for engineering systems: Path towards digital twin

abstract: Neural Operator Networks (ONets) represent a novel advancement in machine learning algorithms, offering a robust and generalizable alternative for approximating partial differential equations (PDEs) solutions. Unlike traditional Neural Networks (NN), which directly approximate functions, ONets specialize in approximating mathematical operators, enhancing their efficacy in addressing complex PDEs. In this work, we evaluate the capabilities of Deep Operator Networks (DeepONets), an ONets implementation using a branch/trunk architecture. Three test cases are studied: a system of ODEs, a general diffusion system, and the convection/diffusion Burgers equation. It is demonstrated that DeepONets can accurately learn the solution operators, achieving prediction accuracy scores above 0.96 for the ODE and diffusion problems over the observed domain while achieving zero shot (without retraining) capability. More importantly, when evaluated on unseen scenarios (zero shot feature), the trained models exhibit excellent generalization ability. This underscores ONets vital niche for surrogate modeling and digital twin development across physical systems. While convection-diffusion poses a greater challenge, the results confirm the promise of ONets and motivate further enhancements to the DeepONet algorithm. This work represents an important step towards unlocking the potential of digital twins through robust and generalizable surrogates.


title: Shackled: a 3D Rendering Engine Programmed Entirely in Ethereum Smart Contracts

abstract: The Ethereum blockchain permits the development and deployment of smart contracts which can store and execute code 'on-chain' - that is, entirely on nodes in the blockchain's network. Smart contracts have traditionally been used for financial purposes, but since smart contracts are Turing-complete, their algorithmic scope is broader than any single domain. To that end, we design, develop, and deploy a comprehensive 3D rendering engine programmed entirely in Ethereum smart contracts, called Shackled. Shackled computes a 2D image from a 3D scene, executing every single computation on-chain, on Ethereum. To our knowledge, Shackled is the first and only fully on-chain 3D rendering engine for Ethereum. In this work, we 1) provide three unique datasets for the purpose of using and benchmarking Shackled, 2) execute said benchmarks and provide results, 3) demonstrate a potential use case of Shackled in the domain of tokenised generative art, 4) provide a no-code user interface to Shackled, 5) enumerate the challenges associated with programming complex algorithms in Solidity smart contracts, and 6) outline potential directions for improving the Shackled platform. It is our hope that this work increases the Ethereum blockchain's native graphics processing capabilities, and that it enables increased use of smart contracts for more complex algorithms, thus increasing the overall richness of the Ethereum ecosystem.


title: Non-Gaussianity from Schwinger-Keldysh Effective Field Theory

abstract: We present a systematic treatment of non-Gaussianity in stochastic systems using the Schwinger-Keldysh effective field theory framework, in which the non-Gaussianity is realized as nonlinear terms in the fluctuation field. We establish two stochastic formulations of the Schwinger-Keldysh effective field theory, with those nonlinear terms manifested as multiple non-Gaussian noises in the Langevin equation and as higher order diffusive terms in the Fokker-Planck equation. The equivalence of the stochastic formulations with the original Schwinger-Keldysh effective field theory is demonstrated with non-trivial examples for arbitrary non-Gaussian parameters. The stochastic formulations will be more flexible and effective in studying non-equilibrium dynamics. We also reveal an ambiguity when coarse-graining time scale and non-Gaussian parameters vanish simultaneously, which may be responsible for the unphysical divergence found in perturbative analysis.


title: Legendrian contact instanton cohomology and its spectral invariants on the one-jet bundle

abstract: In the present paper, we develop the Floer-style elliptic Morse theory for the Hamiltonian-perturbed contact action functional attached to the Legendrian links. Motivated by the present authors' construction [OY2] of the a perturbed action functional defined on the Carnot path space introduced in [OY2] as the canonical generating function, we apply a Floer-type theory to the aforementioned functional and associate the Legendrian contact instanton cohomology, denote by $HI^*(J^1B,H;R)$, to each Legendrian submanifold contact isotopic to the zero section of one-jet bundle. Then we give a Floer theoretic construction of Legendrian spectral invariants and establish their basic properties. This theory subsumes the Lagrangian intersection theory and spectral invariants on the cotangent bundle previously developed by the first-named author in [Oh1,Oh2]. The main ingredient for the study is the interplay between the geometric analysis of the Hamiltonian-perturbed contact instantons and the calculus of contact Hamiltonian geometry.


title: Black hole spectroscopy by mode cleaning

abstract: We formulate a Bayesian framework to analyze ringdown gravitational waves from colliding binary black holes and test the no-hair theorem. The idea hinges on mode cleaning -- revealing subdominant oscillation modes by removing dominant ones using newly proposed ${\it rational~filters}$. By incorporating the filter into Bayesian inference, we construct a likelihood function that depends only on the mass and spin of the remnant black hole (no dependence on mode amplitudes and phases) and implement an efficient pipeline to constrain the remnant mass and spin without Markov chain Monte Carlo (MCMC). We test ringdown models by cleaning combinations of different modes and evaluating the consistency between the residual data and pure noise. The model evidence and Bayes factor are used to demonstrate the presence of a particular mode and to infer the mode starting time. In addition, we design a hybrid approach to estimate the remnant black hole properties exclusively from a single mode using MCMC after mode cleaning. We apply the framework to GW150914 and demonstrate more definitive evidence of the first overtone by cleaning the fundamental mode. This new framework provides a powerful tool for black hole spectroscopy in future gravitational-wave events.


title: Feasibility Analysis of Grover-meets-Simon Algorithm

abstract: Quantum algorithm is a key tool for cryptanalysis. At present, people are committed to building powerful quantum algorithms and tapping the potential of quantum algorithms, so as to further analyze the security of cryptographic algorithms under quantum computing. Recombining classical quantum algorithms is one of the current ideas to construct quantum algorithms. However, they cannot be easily combined, the feasibility of quantum algorithms needs further analysis in quantum environment. This paper reanalyzes the existing combined algorithm Grover-meets-Simon algorithm in terms of the principle of deferred measurement. First of all, due to the collapse problem caused by the measurement, we negate the measurement process of Simon's algorithm during the process of the Grover-meets-Simon algorithm. Second, since the output of the unmeasured Simon algorithm is quantum linear systems of equations, we discuss the solution of quantum linear systems of equations and find it feasible to consider the deferred measurement of the parallel Simon algorithm alone. Finally, since the Grover-meets-Simon algorithm involves an iterative problem, we reconsider the feasibility of the algorithm when placing multiple measurements at the end. According to the maximum probability of success and query times, we get that the Grover-meets-Simon algorithm is not an effective attack algorithm when putting the measurement process of the Simon algorithm in the iterative process at the end of Grover-meets-Simon algorithm.


title: Progression of Digital-Receiver Architecture: From MWA to SKA1-Low,and beyond

abstract: Backed by advances in digital electronics, signal processing, computation, and storage technologies, aperture arrays, which had strongly influenced the design of telescopes in the early years of radio astronomy, have made a comeback. Amid all these developments, an international effort to design and build the world's largest radio telescope, the Square Kilometre Array (SKA), is ongoing. With its vast collecting area of 1 sq-km, the SKA is envisaged to provide unsurpassed sensitivity and leverage technological advances to implement a complex receiver to provide a large field of view through multiple beams on the sky. Many pathfinders and precursor aperture array telescopes for the SKA, operating in the frequency range of 10-300 MHz, have been constructed and operationalized to obtain valuable feedback on scientific, instrumental, and functional aspects. This review article looks explicitly into the progression of digital-receiver architecture from the Murchison Widefield Array (precursor) to the SKA1-Low. It highlights the technological advances in analog-to-digital converters (ADCs),field-programmable gate arrays (FPGAs), and central processing unit-graphics processing unit (CPU-GPU) hybrid platforms around which complex digital signal processing systems implement efficient channelizers, beamformers, and correlators. The article concludes with a preview of the design of a new generation signal processing platform based on radio frequency system-on-chip (RFSoC).


title: Antennas for low-frequency radio telescope of SKA

abstract: The low-frequency radio telescope of the Square Kilometre Array (SKA) is being built by the international radio astronomical community to (i) have orders of magnitude higher sensitivity and (ii) be able to map the sky several hundred times faster, than any other existing facilities over the frequency range of 50 - 350 MHz. The sensitivity of a radio telescope array is in general, dependent upon the number of electromagnetic sensors used to receive the sky signal. The total number of them is further constrained by the effects of mutual coupling between the sensor elements, allowable grating lobes in their radiation patterns, etc. The operating frequency band is governed by the desired spatial and spectral responses, acceptable sidelobe and backlobe levels, radiation efficiency, polarization purity and calibratability of sensors' response. This paper presents a brief review of several broadband antennas considered as potential candidates by various engineering groups across the globe, for the low-frequency radio telescope of SKA covering the frequency range of 50 - 350 MHz, on the basis of their suitability for conducting primary scientific objectives.


title: Identification of the weak-to-strong transition in Alfv\'enic turbulence from space plasma

abstract: Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described phenomena in nature, its core prediction of an Alfvenic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report the first observational evidence for the Alfvenic weak-to-strong transition in MHD turbulence in the terrestrial magnetosheath using the four Cluster spacecraft. The observed transition indicates the universal existence of strong turbulence regardless of the initial level of MHD fluctuations. Moreover, the observations demonstrate that the nonlinear interactions of MHD turbulence play a crucial role in the energy cascade, widening the directions of the energy cascade and broadening the fluctuating frequencies. Our work takes a critical step toward understanding the complete picture of turbulence cascade, connecting the weak and strong MHD turbulence systems. It will have broad implications in star formation, energetic particle transport, turbulent dynamo, and solar corona or solar wind heating.

categories: astro-ph.SR astro-ph.GA physics.plasm-ph physics.space-ph


title: Word Embeddings as Statistical Estimators

abstract: Word embeddings are a fundamental tool in natural language processing. Currently, word embedding methods are evaluated on the basis of empirical performance on benchmark data sets, and there is a lack of rigorous understanding of their theoretical properties. This paper studies word embeddings from a statistical theoretical perspective, which is essential for formal inference and uncertainty quantification. We propose a copula-based statistical model for text data and show that under this model, the now-classical Word2Vec method can be interpreted as a statistical estimation method for estimating the theoretical pointwise mutual information (PMI). Next, by building on the work of Levy and Goldberg (2014), we develop a missing value-based estimator as a statistically tractable and interpretable alternative to the Word2Vec approach. The estimation error of this estimator is comparable to Word2Vec and improves upon the truncation-based method proposed by Levy and Goldberg (2014). The proposed estimator also performs comparably to Word2Vec in a benchmark sentiment analysis task on the IMDb Movie Reviews data set.


title: Numerical experiments on coefficients of instanton partition functions

abstract: We analyze the coefficients of partition functions of Vafa-Witten theory for the complex projective plane $\mathbb{CP}^2$. We experimentally study the growth of the coefficients for gauge group $SU(2)$ and $SU(3)$, which are examples of mock modular forms of depth $1$ and 2 respectively. We also introduce the notion of ``mock cusp form'', and study an example of weight 3 related to the $SU(3)$ partition function. Numerical experiments on the first 200 coefficients suggest that the coefficients of a mock modular form of weight $k$ grow as the coefficients of a modular form of weight $k$, that is to say as $n^{k-1}$. On the other hand the coefficients of the mock cusp form appear to grow as $n^{3/2}$, which exceeds the growth of classical cusp forms of weight 3. We provide bounds using saddle point analysis, which however largely exceed the experimental observation.


title: $(k,q)$-core decomposition of hypergraphs

abstract: In complex networks, many elements interact with each other in different ways. A hypergraph is a network in which group interactions occur among more than two elements. In this study, first, we propose a method to identify influential subgroups in hypergraphs, named $(k,q)$-core decomposition. The $(k,q)$-core is defined as the maximal subgraph in which each vertex has at least $k$ hypergraph degrees \textit{and} each hyperedge contains at least $q$ vertices. The method contains a repeated pruning process until reaching the $(k,q)$-core, which shares similarities with a widely used $k$-core decomposition technique in a graph. Second, we analyze the pruning dynamics and the percolation transition with theoretical and numerical methods in random hypergraphs. We set up evolution equations for the pruning process, and self-consistency equations for the percolation properties. Based on our theory, we find that the pruning process generates a hybrid percolation transition for either $k\ge 3$ \textit{or} $q\ge 3$. The critical exponents obtained theoretically are confirmed with finite-size scaling analysis. Next, when $k=q=2$, we obtain a unconventional degree-dependent critical relaxation dynamics analytically and numerically. Finally, we apply the $(k,q)$-core decomposition to a real coauthorship dataset and recognize the leading groups at an early stage.


title: Open RAN: Evolution of Architecture, Deployment Aspects, and Future Directions

abstract: The Open Radio Access Network (Open RAN) aims to enable disaggregated, virtualized, programmable, and data-driven intelligent network with open interfaces to support various real-time and non-real-time applications for different classes of users and multiple industry verticals in beyond 5G and 6G networks while providing interoperability among multi-vendor network functions and components. In this article, we first discuss the evolution of RAN and then the O-RAN Alliance standardization activities and objectives to provide a comprehensive overview of O-RAN from a standardization point of view. Then, we discuss the O-RAN security aspects, use cases, deployment aspects, and open source projects and related activities in other forums. Finally, we summarize the open issues, challenges, and future research directions to explore further for in-depth study and analysis.


title: Testing Linearly Coasting Cosmology by Strong Lensing System and Pantheon+ Data

abstract: The standard model of cosmology ($\Lambda$CDM) is facing a serious crisis caused by the inconsistencies in the measurements of some fundamental cosmological parameters (Hubble constant $H_{0}$ and cosmic curvature parameter $\Omega_{k}$ for example). On the other hand, a strictly linear evolution of the cosmological scale factor is found to be an excellent fit to a host of observations. Any model that can support such a coasting presents itself as a falsifiable model as far as the cosmological tests are concerned. In this article the observational data of strong gravitational lensing (SGL) systems from SLACS, BELLS, LSD and SL2S surveys has been used to test the viability of linearly coasting cosmology. Assuming the spherically symmetric mass distribution in lensing galaxies, the ratio of angular diameter distance from lens to source and angular diameter distance of the source is evaluated and is used to constrain the power law cosmology. Further, updated type Ia supernovae dataset (Pantheon+) with covariance matrix incorporating all statistical and systematic uncertainties is used to constrain the power law cosmology. It is found that the linear coasting is consistent with the SGL data within 1-$\sigma$ uncertainties but Pantheon+ sample does not support linear coasting.


title: SwinDepth: Unsupervised Depth Estimation using Monocular Sequences via Swin Transformer and Densely Cascaded Network

abstract: Monocular depth estimation plays a critical role in various computer vision and robotics applications such as localization, mapping, and 3D object detection. Recently, learning-based algorithms achieve huge success in depth estimation by training models with a large amount of data in a supervised manner. However, it is challenging to acquire dense ground truth depth labels for supervised training, and the unsupervised depth estimation using monocular sequences emerges as a promising alternative. Unfortunately, most studies on unsupervised depth estimation explore loss functions or occlusion masks, and there is little change in model architecture in that ConvNet-based encoder-decoder structure becomes a de-facto standard for depth estimation. In this paper, we employ a convolution-free Swin Transformer as an image feature extractor so that the network can capture both local geometric features and global semantic features for depth estimation. Also, we propose a Densely Cascaded Multi-scale Network (DCMNet) that connects every feature map directly with another from different scales via a top-down cascade pathway. This densely cascaded connectivity reinforces the interconnection between decoding layers and produces high-quality multi-scale depth outputs. The experiments on two different datasets, KITTI and Make3D, demonstrate that our proposed method outperforms existing state-of-the-art unsupervised algorithms.


title: Effective Shielding of $\lesssim$ 10 GeV Cosmic Rays from Dense Molecular Clumps

abstract: The density of cosmic rays inside molecular clouds determines the ionization rate in the dense cores where stars form. It is also one of the drivers of astrochemistry leading to the creation of complex molecules. Through Fermi Large Area Telescope observations of nearby giant molecular clouds, we observed deficits (holes) in the gamma-ray residual map when modelling with the expected gamma-ray diffuse emission from uniform cosmic rays interacting with the molecular content. We propose that the deficit is due to the lack of penetration of the low-energy (sub-GeV to GeV) cosmic rays into denser regions or clumps. This differs from the prevailing view of fast cosmic ray transport in giant molecular clouds where the magnetic turbulence is suppressed by neutral-ion damping, as our results require a slow diffusion inside dense molecular clumps. Through modelling we find that while the shielding is negligible on the cloud scale, it becomes important in the denser, parsec-sized regions where the gravitational collapse is already at play, changing the initial condition of star formation and astrochemistry.


title: Hypersurface Curvatures of Geological Features

abstract: Reflector-normal angles and reflector-curvature parameters are the principal geometric attributes used in seismic interpretation for characterizing the orientations and shapes, respectively, of geological reflecting surfaces. Commonly, the input dataset for their computation consists of fine 3D grids of scalar fields representing either the seismic-driven reflectivities (e.g., amplitudes of 3D seismic migrated volumes) or model-driven reflectivities, computed, for example, from the derived elastic impedance parameters. Conventionally, the computation of curvature parameters at each grid point is based on analyzing the local change in the inline/crossline dips, considering the potential existence of a local quadratic reflecting surface in the vicinity of that point. This assumption breaks down for subsurface points in the vicinity of either complex reflecting surfaces (e.g., brittle/rough/tilted synclines/anticlines, ridges/troughs, and saddles) and/or sharp discontinuous geological features (e.g., fault edges/tips, pinch-outs, fracture systems, channels, small geobodies), where the values of the computed curvature become extremely high. While these high values can indicate the existence of non-reflecting objects, they do not deliver their specific geometric characteristics. We present a novel method that better characterizes the shapes of these complex geological features by extending the assumption of local surfaces (2D surfaces in 3D space) into hypersurfaces (3D hypersurfaces in 4D space), with their corresponding principal and effective curvature parameters. We demonstrate the advantages of our method by comparing the conventional dip-based surface curvature parameters with the hypersurface curvature parameters, using a synthetic model/image with different types and shapes of geological features and two seismic images of real data containing complex faults and hidden buried channels.


title: Sparse and Integrative Principal Component Analysis for Multiview Data

abstract: We consider dimension reduction of multiview data, which are emerging in scientific studies. Formulating multiview data as multi-variate data with block structures corresponding to the different views, or views of data, we estimate top eigenvectors from multiview data that have two-fold sparsity, elementwise sparsity and blockwise sparsity. We propose a Fantope-based optimization criterion with multiple penalties to enforce the desired sparsity patterns and a denoising step is employed to handle potential presence of heteroskedastic noise across different data views. An alternating direction method of multipliers (ADMM) algorithm is used for optimization. We derive the l2 convergence of the estimated top eigenvectors and establish their sparsity and support recovery properties. Numerical studies are used to illustrate the proposed method.


title: FemtoDet: An Object Detection Baseline for Energy Versus Performance Tradeoffs

abstract: Efficient detectors for edge devices are often optimized for parameters or speed count metrics, which remain in weak correlation with the energy of detectors. However, some vision applications of convolutional neural networks, such as always-on surveillance cameras, are critical for energy constraints. This paper aims to serve as a baseline by designing detectors to reach tradeoffs between energy and performance from two perspectives: 1) We extensively analyze various CNNs to identify low-energy architectures, including selecting activation functions, convolutions operators, and feature fusion structures on necks. These underappreciated details in past work seriously affect the energy consumption of detectors; 2) To break through the dilemmatic energy-performance problem, we propose a balanced detector driven by energy using discovered low-energy components named \textit{FemtoDet}. In addition to the novel construction, we improve FemtoDet by considering convolutions and training strategy optimizations. Specifically, we develop a new instance boundary enhancement (IBE) module for convolution optimization to overcome the contradiction between the limited capacity of CNNs and detection tasks in diverse spatial representations, and propose a recursive warm-restart (RecWR) for optimizing training strategy to escape the sub-optimization of light-weight detectors by considering the data shift produced in popular augmentations. As a result, FemtoDet with only 68.77k parameters achieves a competitive score of 46.3 AP50 on PASCAL VOC and 1.11 W $\&$ 64.47 FPS on Qualcomm Snapdragon 865 CPU platforms. Extensive experiments on COCO and TJU-DHD datasets indicate that the proposed method achieves competitive results in diverse scenes.


title: Testing Firm Conduct

abstract: Evaluating policy in imperfectly competitive markets requires understanding firm behavior. While researchers test conduct via model selection and assessment, we present advantages of Rivers and Vuong (2002) (RV) model selection under misspecification. However, degeneracy of RV invalidates inference. With a novel definition of weak instruments for testing, we connect degeneracy to instrument strength, derive weak instrument properties of RV, and provide a diagnostic for weak instruments by extending the framework of Stock and Yogo (2005) to model selection. We test vertical conduct (Villas-Boas, 2007) using common instrument sets. Some are weak, providing no power. Strong instruments support manufacturers setting retail prices.


title: On Delay-Doppler Plane Orthogonal Pulse

abstract: In this paper, we analyze the recently discovered delay-Doppler plane orthogonal pulse (DDOP), which is essential for delay-Doppler plane multi-carrier modulation waveform. In particular, we introduce a local orthogonality property of pulses corresponding to Weyl-Heisenberg (WH) subset and justify the DDOP's existence, in contrast to global orthogonality corresponding to WH set governed by the WH frame theory. Then, sufficient conditions for locally-orthogonal pulses are presented and discussed. Based on the analysis, we propose a general DDOP design. We also derive the frequency domain representation of the DDOP, and compare the DDOP-based orthogonal delay-Doppler division multiplexing (ODDM) modulation with other modulation schemes, in terms of TF signal localization. Interestingly, we show perfect local orthogonality property of the DDOP with respect to delay-Doppler resolutions using its ambiguity function.


title: Joint constraints on cosmological parameters using future multi-band gravitational wave standard siren observations

abstract: Gravitational waves (GWs) from the compact binary coalescences can be used as standard sirens to explore the cosmic expansion history. In the next decades, it is anticipated that we could obtain the multi-band GW standard siren data (from nanohertz to a few hundred hertz), which are expected to play an important role in cosmological parameter estimation. In this work, we give for the first time the joint constraints on cosmological parameters using the future multi-band GW standard siren observations. We simulate the multi-band GW standard sirens based on the SKA-era pulsar timing array (PTA), the Taiji observatory, and the Cosmic Explorer (CE) to perform cosmological analysis. In the $\Lambda$CDM model, we find that the joint PTA+Taiji+CE data could provide a tight constraint on the Hubble constant with a $0.5\%$ precision. Moreover, PTA+Taiji+CE could break the cosmological parameter degeneracies generated by CMB, especially in the dynamical dark energy models. When combining the PTA+Taiji+CE data with the CMB data, the constraint precisions of $\Omega_{\rm m}$ and $H_0$ are $1.0\%$ and $0.3\%$, meeting the standard of precision cosmology. The joint CMB+PTA+Taiji+CE data give $\sigma(w)=0.028$ in the $w$CDM model and $\sigma(w_0)=0.11$ and $\sigma(w_a)=0.32$ in the $w_0w_a$CDM model, which are comparable with or close to the latest constraint results by CMB+BAO+SN. In conclusion, it is worth expecting to use the future multi-band GW observations to explore the nature of dark energy and measure the Hubble constant.


title: Center to limb variation of transition region Doppler shift in active regions

abstract: Studying Doppler shifts provides deeper insights into the flow of mass and energy in the solar atmosphere. We perform a comprehensive measurement of Doppler shifts in the transition region and its center-to-limb variation (CLV) in the strong field regions ($|\textbf{B}| \geq$ 50 G) of 50 active regions (ARs), using the \ion{Si}{4} 1394~{\AA} line recorded by the Interface Region Imaging Spectrometer(IRIS). To locate the ARs and identify strong field regions, we have used the magnetograms obtained by the Helioseismic and Magnetic Imager (HMI). We find that in strong field regions, on average, all the ARs show mean redshifts ranging between 4{--}11~ km/s, which varies with ARs. These flows show a mild CLV, with sizable magnitudes at the limb and substantial scatter at the mid-longitude range. Our observations do not support the idea that redshifts in the lower transition region (T $<\sim$ 0.1 MK) are produced by field-aligned downflows as a result of impulsive heating and warrant alternative interpretation, such as downflow of type-\rm{II} spicules in the presence of a chromospheric wall created by cooler type-\rm{I} spicules.


title: Calculating lexicase selection probabilities is NP-Hard

abstract: Calculating the probability of an individual solution being selected under lexicase selection is an important problem in attempts to develop a deeper theoretical understanding of lexicase selection, a state-of-the art parent selection algorithm in evolutionary computation. Discovering a fast solution to this problem would also have implications for efforts to develop practical improvements to lexicase selection. Here, I prove that this problem, which I name lex-prob, is NP-Hard. I achieve this proof by reducing SAT, a well-known NP-Complete problem, to lex-prob in polynomial time. This reduction involves an intermediate step in which a popular variant of lexicase selection, epsilon-lexicase selection, is reduced to standard lexicase selection. This proof has important practical implications for anyone needing a fast way of calculating the probabilities of individual solutions being selected under lexicase selection. Doing so in polynomial time would be incredibly challenging, if not all-together impossible. Thus, finding approximation algorithms or practical optimizations for speeding up the brute-force solution is likely more worthwhile. This result also has deeper theoretical implications about the relationship between epsilon-lexicase selection and lexicase selection and the relationship between lex-prob and other NP-Hard problems.


title: Optimal Placement of Active and Passive Elements in Hybrid RIS-assisted Communication Systems

abstract: Hybrid reconfigurable intelligent surfaces (HRIS) are RIS architectures having both active and passive elements. The received signal-to-noise ratio (SNR) in HRIS-assisted communication systems depends on the placement of active elements. In this paper, we show that received SNR can be improved with a channel-aware placement of the active elements. We jointly design the transmit precoder, the RIS coefficients, and the location of active and passive elements of the HRIS to maximize the SNR. We solve the underlying combinatorial nonconvex optimization problem using alternating optimization and propose a low-complexity solver, which is provably nearly optimal. Through numerical simulations, we demonstrate that the proposed method offers significantly improved performance compared to communication systems having a fully passive RIS array or a hybrid RIS array with channel agnostic active element placement and performance comparable to that of communication systems with a fully active RIS array.


title: Self-exciting negative binomial distribution process and critical properties of intensity distribution

abstract: We study the continuous time limit of a self-exciting negative binomial process and discuss the critical properties of its intensity distribution. In this limit, the process transforms into a marked Hawkes process. The probability mass function of the marks has a parameter $\omega$, and the process reduces to a "pure" Hawkes process in the limit $\omega\to 0$. We investigate the Lagrange--Charpit equations for the master equations of the marked Hawkes process in the Laplace representation close to its critical point and extend the previous findings on the power-law scaling of the probability density function (PDF) of intensities in the intermediate asymptotic regime to the case where the memory kernel is the superposition of an arbitrary finite number of exponentials. We develop an efficient sampling method for the marked Hawkes process based on the time-rescaling theorem and verify the power-law exponents.


title: Roadmap for Unconventional Computing with Nanotechnology

abstract: In the "Beyond Moore's Law" era, with increasing edge intelligence, domain-specific computing embracing unconventional approaches will become increasingly prevalent. At the same time, adopting a variety of nanotechnologies will offer benefits in energy cost, computational speed, reduced footprint, cyber resilience, and processing power. The time is ripe for a roadmap for unconventional computing with nanotechnologies to guide future research, and this collection aims to fill that need. The authors provide a comprehensive roadmap for neuromorphic computing using electron spins, memristive devices, two-dimensional nanomaterials, nanomagnets, and various dynamical systems. They also address other paradigms such as Ising machines, Bayesian inference engines, probabilistic computing with p-bits, processing in memory, quantum memories and algorithms, computing with skyrmions and spin waves, and brain-inspired computing for incremental learning and problem-solving in severely resource-constrained environments. These approaches have advantages over traditional Boolean computing based on von Neumann architecture. As the computational requirements for artificial intelligence grow 50 times faster than Moore's Law for electronics, more unconventional approaches to computing and signal processing will appear on the horizon, and this roadmap will help identify future needs and challenges. In a very fertile field, experts in the field aim to present some of the dominant and most promising technologies for unconventional computing that will be around for some time to come. Within a holistic approach, the goal is to provide pathways for solidifying the field and guiding future impactful discoveries.


title: Some Appell-type orthogonal polynomials on lattices

abstract: We investigate on some Appel-type orthogonal polynomial sequences on q-quadratic lattices and we provide some entire new characterizations of the Al-Salam Chihara polynomials (including the Rogers q-Hermite polynomials). The corresponding forms are well described. The proposed method can be applied to similar and to more general problems involving the Askey-Wilson and the Averaging operators, in order to obtain new characterization theorems for classical and semiclassical orthogonal polynomials on lattices.


title: Existence Theorems on Quasi-variational Inequalities over Banach Spaces and its Applications to Time-dependent Pure Exchange Economy

abstract: We study a class of quasi-variational inequality problems defined over infinite dimensional Banach space and deduce sufficient conditions for ensuring solutions to such problems under the upper semi-continuity and pseudomonotonicity assumptions on the map defining the inequalities. The special structure of the quasi-variational inequality enables us to show the occurrence of solutions for such inequalities based on the classical existence theorem for variational inequalities. This special type of quasi-variational inequalities is motivated by the pure exchange economic problems and Radner equilibrium problems for sequential trading game. Further, we study the solvability of the specific class of quasi-variational inequalities on Banach spaces in which the constraint map may admit unbounded values. Finally, we demonstrate the occurrence of dynamic competitive equilibrium for a time-dependent pure exchange economy as an application.


title: Bag of States: A Non-sequential Approach to Video-based Engagement Measurement

abstract: Automatic measurement of student engagement provides helpful information for instructors to meet learning program objectives and individualize program delivery. Students' behavioral and emotional states need to be analyzed at fine-grained time scales in order to measure their level of engagement. Many existing approaches have developed sequential and spatiotemporal models, such as recurrent neural networks, temporal convolutional networks, and three-dimensional convolutional neural networks, for measuring student engagement from videos. These models are trained to incorporate the order of behavioral and emotional states of students into video analysis and output their level of engagement. In this paper, backed by educational psychology, we question the necessity of modeling the order of behavioral and emotional states of students in measuring their engagement. We develop bag-of-words-based models in which only the occurrence of behavioral and emotional states of students is modeled and analyzed and not the order in which they occur. Behavioral and affective features are extracted from videos and analyzed by the proposed models to determine the level of engagement in an ordinal-output classification setting. Compared to the existing sequential and spatiotemporal approaches for engagement measurement, the proposed non-sequential approach improves the state-of-the-art results. According to experimental results, our method significantly improved engagement level classification accuracy on the IIITB Online SE dataset by 26% compared to sequential models and achieved engagement level classification accuracy as high as 66.58% on the DAiSEE student engagement dataset.


title: On discrete-time dissipative port-Hamiltonian (descriptor) systems

abstract: Port-Hamiltonian (pH) systems have been studied extensively for linear continuous-time dynamical systems. This manuscript presents a discrete-time pH descriptor formulation for linear, completely causal, scattering passive dynamical systems based on the system coefficients. The relation of this formulation to positive and bounded real systems and the characterization via positive semidefinite solutions of Kalman-Yakubovich-Popov inequalities is also studied.


title: Predictive Modeling of Coronal Hole Areas Using Long Short-Term Memory Networks

abstract: In the era of space exploration, the implications of space weather have become increasingly evident. Central to this is the phenomenon of coronal holes, which can significantly influence the functioning of satellites and aircraft. These coronal holes, present on the sun, are distinguished by their open magnetic field lines and comparatively cooler temperatures, leading to the emission of solar winds at heightened rates. To anticipate the effects of these coronal holes on Earth, our study harnesses computer vision to pinpoint the coronal hole regions and estimate their dimensions using imagery from the Solar Dynamics Observatory (SDO). Further, we deploy deep learning methodologies, specifically the Long Short-Term Memory (LSTM) approach, to analyze the trends in the data related to the area of the coronal holes and predict their dimensions across various solar regions over a span of seven days. By evaluating the time series data concerning the area of the coronal holes, our research seeks to uncover patterns in the behavior of coronal holes and comprehend their potential influence on space weather occurrences. This investigation marks a pivotal stride towards bolstering our capacity to anticipate and brace for space weather events that could have ramifications for Earth and its technological apparatuses.

categories: astro-ph.SR astro-ph.EP cs.LG physics.space-ph


title: Face Inverse Rendering via Hierarchical Decoupling

abstract: Previous face inverse rendering methods often require synthetic data with ground truth and/or professional equipment like a lighting stage. However, a model trained on synthetic data or using pre-defined lighting priors is typically unable to generalize well for real-world situations, due to the gap between synthetic data/lighting priors and real data. Furthermore, for common users, the professional equipment and skill make the task expensive and complex. In this paper, we propose a deep learning framework to disentangle face images in the wild into their corresponding albedo, normal, and lighting components. Specifically, a decomposition network is built with a hierarchical subdivision strategy, which takes image pairs captured from arbitrary viewpoints as input. In this way, our approach can greatly mitigate the pressure from data preparation, and significantly broaden the applicability of face inverse rendering. Extensive experiments are conducted to demonstrate the efficacy of our design, and show its superior performance in face relighting over other state-of-the-art alternatives. {Our code is available at \url{https://github.com/AutoHDR/HD-Net.git}}


title: Towards real-time reconstruction of velocity fluctuations in turbulent channel flow

abstract: We develop a framework for efficient streaming reconstructions of turbulent velocity fluctuations from limited sensor measurements with the goal of enabling real-time applications. The reconstruction process is simplified by computing linear estimators using flow statistics from an initial training period and evaluating their performance during a subsequent testing period with data obtained from direct numerical simulation. We address cases where (i) no, (ii) limited, and (iii) full-field training data are available using estimators based on (i) resolvent modes, (ii) resolvent-based estimation, and (iii) spectral proper orthogonal decomposition modes. During training, we introduce blockwise inversion to accurately and efficiently compute the resolvent operator in an interpretable manner. During testing, we enable efficient streaming reconstructions by using a temporal sliding discrete Fourier transform to recursively update Fourier coefficients using incoming measurements. We use this framework to reconstruct with minimal time delay the turbulent velocity fluctuations in a minimal channel at ${\rm Re}_\tau \approx 186$ from sparse planar measurements. We evaluate reconstruction accuracy in the context of the extent of data required and thereby identify potential use cases for each estimator. The reconstructions capture large portions of the dynamics from relatively few measurement planes when the linear estimators are computed with sufficient fidelity. We also evaluate the efficiency of our reconstructions and show that the present framework has the potential to help enable real-time reconstructions of turbulent velocity fluctuations in an analogous experimental setting.


title: Two Stage Contextual Word Filtering for Context bias in Unified Streaming and Non-streaming Transducer

abstract: It is difficult for an E2E ASR system to recognize words such as entities appearing infrequently in the training data. A widely used method to mitigate this issue is feeding contextual information into the acoustic model. Previous works have proven that a compact and accurate contextual list can boost the performance significantly. In this paper, we propose an efficient approach to obtain a high quality contextual list for a unified streaming/non-streaming based E2E model. Specifically, we make use of the phone-level streaming output to first filter the predefined contextual word list then fuse it into non-casual encoder and decoder to generate the final recognition results. Our approach improve the accuracy of the contextual ASR system and speed up the inference process. Experiments on two datasets demonstrates over 20% CER reduction comparing to the baseline system. Meanwhile, the RTF of our system can be stabilized within 0.15 when the size of the contextual word list grows over 6,000.


title: Syllable Subword Tokens for Open Vocabulary Speech Recognition in Malayalam

abstract: In a hybrid automatic speech recognition (ASR) system, a pronunciation lexicon (PL) and a language model (LM) are essential to correctly retrieve spoken word sequences. Being a morphologically complex language, the vocabulary of Malayalam is so huge and it is impossible to build a PL and an LM that cover all diverse word forms. Usage of subword tokens to build PL and LM, and combining them to form words after decoding, enables the recovery of many out of vocabulary words. In this work we investigate the impact of using syllables as subword tokens instead of words in Malayalam ASR, and evaluate the relative improvement in lexicon size, model memory requirement and word error rate.


title: The Chow ring of a sequence of point blow-ups

abstract: Given a sequence of point blow-ups of smooth n-dimensional projective varieties $Z_{i}$ defined over an algebraically closed field k, $Z_{s}\rightarrow Z_{s-1}\rightarrow ...\rightarrow Z_{1}\rightarrow Z_{0}$, we give two presentations of the Chow ring of its sky $A^{\bullet}(Z_{s})$. The first one using the classes of the total transforms of the exceptional components as generators and the second one using the classes of the strict transforms ones. We prove that the skies of two sequences of point blow-ups of the same length have isomorphic Chow rings. Finally we give a characterization of final divisor of a sequence of point blow-ups in terms of some relations defined over the Chow group of zero-cycles of its sky $A_{0}(Z_{s})$.


title: HUBO and QUBO models for Prime factorization

abstract: The security of the RSA cryptosystem is based on the difficulty of factoring a large number N into prime numbers p and q satisfying N=p*q . This paper presents a prime factoriaation method using D-Wave quantum computer that can threaten the RSA cryptosystem. The starting point for this method is very simple, representing two prime numbers as qubits. Then, set the difference between the product of two prime numbers expressed in qubits and N as a cost function, and find the solution when the cost function becomes the minimum. D-Wave's quantum annealer can find the minimum value of any quadratic problem. However, the cost function is to be a higher-order unconstrained optimiaation (HUBO) model because it contains the second or higher order terms. We used a hybrid solver and dimod package provided by -Wave Ocean software development kit (SDK) to solve the HUBO problem. We also successfully factoriaed 102,454,763 with 26 logical qubits. In addition, we factoriaed 1,000,070,001,221 using the range dependent Hamiltonian algorithm.


title: Recoverability and estimation of causal effects under typical multivariable missingness mechanisms

abstract: In the context of missing data, the identifiability or "recoverability" of the average causal effect (ACE) depends on causal and missingness assumptions. The latter can be depicted by adding variable-specific missingness indicators to causal diagrams, creating "missingness-directed acyclic graphs" (m-DAGs). Previous research described ten canonical m-DAGs, representing typical multivariable missingness mechanisms in epidemiological studies, and determined the recoverability of the ACE in the absence of effect modification. We extend the research by determining the recoverability of the ACE in settings with effect modification and conducting a simulation study evaluating the performance of widely used missing data methods when estimating the ACE using correctly specified g-computation, which has not been previously studied. Methods assessed were complete case analysis (CCA) and various multiple imputation (MI) implementations regarding the degree of compatibility with the outcome model used in g-computation. Simulations were based on an example from the Victorian Adolescent Health Cohort Study (VAHCS), where interest was in estimating the ACE of adolescent cannabis use on mental health in young adulthood. In the canonical m-DAGs that excluded unmeasured common causes of missingness indicators, we derived the recoverable ACE if no incomplete variable causes its missingness, and non-recoverable otherwise. Besides, the simulation showed that compatible MI approaches may enable approximately unbiased ACE estimation, unless the outcome causes its missingness or it causes the missingness of a variable that causes its missingness. Researchers must consider sensitivity analysis methods incorporating external information in the latter setting. The VAHCS case study illustrates the practical implications of these findings.


title: Contributions of negative-energy states to the E2-M1 polarizability of the Sr clock

abstract: With the improvement of high-precision optical clock, the higher-order multipolar interaction between atoms and light needs quantitative evaluation. However for the Sr clock, the differential dynamic E2-M1 polarizability at the magic wavelength has contradictions among available theoretical and experimental results. Recently, the new experimental measurement of S. D\"{o}rscher {\em et al.} [arXiv: 2210. 14727] is consistent with measurement of Ushijima {\em et al.}, which poses new challenges to theory and urgently calls for theoretical explanations. In present work, we investigate contributions of negative-energy states to the E2 and M1 polarizabilities. We find that for the M1 polarizability, the contribution from negative-energy states is crucial and dominant. Our new theoretical result for E2-M1 polarizability difference is $-7.74(3.92)\times 10^{-5}$ a.u., which is in good agreement with the recent experiment of S. D\"{o}rscher et al., so the inconsistency problem of E2-M1 polarizability in the Sr clock between theory and experiment is eliminated.


title: Faster Sinkhorn's Algorithm with Small Treewidth

abstract: Computing optimal transport (OT) distances such as the earth mover's distance is a fundamental problem in machine learning, statistics, and computer vision. In this paper, we study the problem of approximating the general OT distance between two discrete distributions of size $n$. Given the cost matrix $C=AA^\top$ where $A \in \mathbb{R}^{n \times d}$, we proposed a faster Sinkhorn's Algorithm to approximate the OT distance when matrix $A$ has treewidth $\tau$. To approximate the OT distance, our algorithm improves the state-of-the-art results [Dvurechensky, Gasnikov, and Kroshnin ICML 2018] from $\widetilde{O}(\epsilon^{-2} n^2)$ time to $\widetilde{O}(\epsilon^{-2} n \tau)$ time.


title: Robust Realized Integrated Beta Estimator with Application to Dynamic Analysis of Integrated Beta

abstract: In this paper, we develop a robust non-parametric realized integrated beta estimator using high-frequency financial data contaminated by microstructure noises, which is robust to the stylized features, such as the time-varying beta and the dependence structure of microstructure noises. With this robust realized integrated beta estimator, we investigate dynamic structures of integrated betas and find an auto-regressive--moving-average (ARMA) structure. To model this dynamic structure, we utilize the ARMA model for daily integrated market betas. We call this the dynamic realized beta (DR Beta). We further introduce a high-frequency data generating process by filling the gap between the high-frequency-based non-parametric estimator and low-frequency dynamic structure. Then, we propose a quasi-likelihood procedure for estimating the model parameters with the robust realized integrated beta estimator as the proxy. We also establish asymptotic theorems for the proposed estimator and conduct a simulation study to check the performance of finite samples of the estimator. The empirical study with the S&P 500 index and the top 50 large trading volume stocks from the S&P 500 illustrates that the proposed DR Beta model with the robust realized beta estimator effectively accounts for dynamics in the market beta of individual stocks and better predicts future market betas.


title: On the Gamma-Ray Emission of the Andromeda Galaxy M31

abstract: Using the $\gamma$-ray data obtained with the Large Area Telescope (LAT) onboard {\it the Fermi Gamma-ray Space Telescope (Fermi)} for $\sim$14 years, we examine the high energy emission emanating from the center of the Andromeda Galaxy M31. Different from previously reported results, which show a seemingly extended source, we instead find two individual point sources, one consistent with being at the center and one 0\fdg4 south-east of the center. The emission of the former is well described using a Log-Parabola model, similar to those of previous studies, and that of the latter can be fitted with a power law. We discuss the possible origins for the two sources. M31's central source, now consistent with being a point source, necessitates a revisit of its previously discussed originations with this new property taken into consideration, in particular those cosmic rays or dark matter scenarios involving extended source distributions. The SE source appears to have a projected distance of $\sim$6\,kpc from M31's center, and the investigation is required as to whether it is a source locally associated with M31, or is instead a background extra-galactic one.


title: Two-sided heat kernel estimates for Schr\"{o}dinger operators with unbounded potentials

abstract: Consider the Schr\"odinger operator $ \mathcal L^V=-\Delta+V $ on $\R^d$, where $V:\R^d\to [0,\infty)$ is a nonnegative and locally bounded potential on $\R^d$ so that for all $x\in \R^d$ with $|x|\ge 1$, $c_1g(|x|)\le V(x)\le c_2g(|x|)$ with some constants $c_1,c_2>0$ and a nondecreasing and strictly positive function $g:[0,\infty)\to [1,+\infty)$ that satisfies $g(2r)\le c_0 g(r)$ for all $r>0$ and $\lim_{r\to \infty} g(r)=\infty.$ We establish global in time and qualitatively sharp bounds for the heat kernel of the associated Schr\"{o}dinger semigroup by the probabilistic method. In particular, we can present global in space and time two-sided bounds of heat kernel even when the Schr\"{o}dinger semigroup is not intrinsically ultracontractive. Furthermore, two-sided estimates for the corresponding Green's functions are also obtained.