Search

Let $p$ be an odd prime, Jianqiang Zhao has established a curious congruence $$ \sum_{i+j+k=p \atop i,j,k > 0} \frac{1}{ijk} \equiv -2B_{p-3}\pmod p , $$ where $B_{n}$ denotes the $n-$th Bernoulli numbers. In this paper, we will generalize this problem by using congruent theory and combinatorial methods, and we get some curious congruences.

Visualizing high-dimensional data is an important routine for understanding biomedical data and interpreting deep learning models. Neighbor embedding methods, such as t-SNE, UMAP, and LargeVis, among others, are a family of popular visualization methods which reduce high-dimensional data to two dimensions. However, recent studies suggest that these methods often produce visual artifacts, potentially leading to incorrect scientific conclusions. Recognizing that the current limitation stems from a lack of data-independent notions of embedding maps, we introduce a novel conceptual and computational framework, LOO-map, that learns the embedding maps based on a classical statistical idea known as the leave-one-out. LOO-map extends the embedding over a discrete set of input points to the entire input space, enabling a systematic assessment of map continuity, and thus the reliability of the visualizations. We find for many neighbor embedding methods, their embedding maps can be intrinsically discontinuous. The discontinuity induces two types of observed map distortion: ``overconfidence-inducing discontinuity," which exaggerates cluster separation, and ``fracture-inducing discontinuity," which creates spurious local structures. Building upon LOO-map, we propose two diagnostic point-wise scores -- perturbation score and singularity score -- to address these limitations. These scores can help identify unreliable embedding points, detect out-of-distribution data, and guide hyperparameter selection. Our approach is flexible and works as a wrapper around many neighbor embedding algorithms. We test our methods across multiple real-world datasets from computer vision and single-cell omics to demonstrate their effectiveness in enhancing the interpretability and accuracy of visualizations., Comment: 43 pages, 15 figures

Deep supervised models possess significant capability to assimilate extensive training data, thereby presenting an opportunity to enhance model performance through training on multiple datasets. However, conflicts arising from different label spaces among datasets may adversely affect model performance. In this paper, we propose a novel approach to automatically construct a unified label space across multiple datasets using graph neural networks. This enables semantic segmentation models to be trained simultaneously on multiple datasets, resulting in performance improvements. Unlike existing methods, our approach facilitates seamless training without the need for additional manual reannotation or taxonomy reconciliation. This significantly enhances the efficiency and effectiveness of multi-dataset segmentation model training. The results demonstrate that our method significantly outperforms other multi-dataset training methods when trained on seven datasets simultaneously, and achieves state-of-the-art performance on the WildDash 2 benchmark.

Embedding high-dimensional data into a low-dimensional space is an indispensable component of data analysis. In numerous applications, it is necessary to align and jointly embed multiple datasets from different studies or experimental conditions. Such datasets may share underlying structures of interest but exhibit individual distortions, resulting in misaligned embeddings using traditional techniques. In this work, we propose \textit{Entropic Optimal Transport (EOT) eigenmaps}, a principled approach for aligning and jointly embedding a pair of datasets with theoretical guarantees. Our approach leverages the leading singular vectors of the EOT plan matrix between two datasets to extract their shared underlying structure and align the datasets accordingly in a common embedding space. We interpret our approach as an inter-data variant of the classical Laplacian eigenmaps and diffusion maps embeddings, showing that it enjoys many favorable analogous properties. We then analyze a data-generative model where two observed high-dimensional datasets share latent variables on a common low-dimensional manifold, but each dataset is subject to data-specific translation, scaling, nuisance structures, and noise. We show that in a high-dimensional asymptotic regime, the EOT plan recovers the shared manifold structure by approximating a kernel function evaluated at the locations of the latent variables. Subsequently, we provide a geometric interpretation of our embedding by relating it to the eigenfunctions of population-level operators encoding the density and geometry of the shared manifold. Finally, we showcase the performance of our approach for data integration and embedding through simulations and analyses of real-world biological data, demonstrating its advantages over alternative methods in challenging scenarios.

Low-dimensional embeddings (LDEs) of high-dimensional data are ubiquitous in science and engineering. They allow us to quickly understand the main properties of the data, identify outliers and processing errors, and inform the next steps of data analysis. As such, LDEs have to be faithful to the original high-dimensional data, i.e., they should represent the relationships that are encoded in the data, both at a local as well as global scale. The current generation of LDE approaches focus on reconstructing local distances between any pair of samples correctly, often out-performing traditional approaches aiming at all distances. For these approaches, global relationships are, however, usually strongly distorted, often argued to be an inherent trade-off between local and global structure learning for embeddings. We suggest a new perspective on LDE learning, reconstructing angles between data points. We show that this approach, Mercat, yields good reconstruction across a diverse set of experiments and metrics, and preserve structures well across all scales. Compared to existing work, our approach also has a simple formulation, facilitating future theoretical analysis and algorithmic improvements.

Sulfur chemistry in the formation process of low-mass stars and planets remains poorly understood. The protoplanetary disks (PPDs) are the birthplace of planets and its distinctive environment provides an intriguing platform for investigating models of sulfur chemistry. We analyzed the ALMA observations of CS 7-6 transitions in the HD 163296 disk and perform astrochemical modeling to explore its sulfur chemistry. We simulated the distribution of sulfur-containing molecules and compared it with observationally deduced fractional column densities. We have found that the simulated column density of CS is consistent with the observationally deduced fractional column densities, while the simulated column density of C$_2$S is lower than the observationally deduced upper limits on column densities. This results indicate that we have a good understanding of the chemical properties of CS and C$_2$S in the disk. We also investigated the influence of the C/O ratio on sulfur-containing molecules and found that the column densities of SO, SO$_2$, and H$_2$S near the central star are dependent on the C/O ratio. Additionally, we found that the $N$[CS]/$N$[SO] ratio can serve as a promising indicator of the disk's C/O ratio in the HD 163296. Overall, the disk of HD 163296 provides a favorable environment for the detection of sulfur-containing molecules., Comment: 16 pages, 7 figures

We introduce a wall model (WM) for large-eddy simulation (LES) applicable to rough surfaces with Gaussian and non-Gaussian distributions for both transitionally and fully rough regimes. The model is applicable to arbitrary complex geometries where roughness elements are assumed to be underresolved. The wall model is implemented using a feedforward neural network, with the geometric properties of the roughness topology and near-wall flow quantities serving as input. The optimal set of non-dimensional input features is identified using information theory, selecting variables that maximize information about the output while minimizing redundancy among inputs. The model incorporates a confidence score based on Gaussian process modeling, enabling the detection of low model performance for unseen rough surfaces. The model is trained using a direct numerical simulation roughness database comprising approximately 200 cases. The roughness geometries for the database are selected from a large repository through active learning. This approach ensures that the rough surfaces incorporated into the database are the most informative. The model performance is evaluated both a-priori and a-posteriori in WMLES of turbulent channel flows with rough walls. Over 120 channel flow cases are considered, including untrained roughness geometries, roughness Reynolds numbers, and grid resolutions for both transitionally- and fully-rough regimes. The results show that the rough-wall model predicts the wall shear stress within 15% accuracy. The model is also assessed on a high-pressure turbine blade with two different rough surfaces. The WM predicts the skin friction and the mean velocity deficit within 10% accuracy except the region with shock waves.

Integrative analysis of multiple heterogeneous datasets has become standard practice in many research fields, especially in single-cell genomics and medical informatics. Existing approaches oftentimes suffer from limited power in capturing nonlinear structures, insufficient account of noisiness and effects of high-dimensionality, lack of adaptivity to signals and sample sizes imbalance, and their results are sometimes difficult to interpret. To address these limitations, we propose a novel kernel spectral method that achieves joint embeddings of two independently observed high-dimensional noisy datasets. The proposed method automatically captures and leverages possibly shared low-dimensional structures across datasets to enhance embedding quality. The obtained low-dimensional embeddings can be utilized for many downstream tasks such as simultaneous clustering, data visualization, and denoising. The proposed method is justified by rigorous theoretical analysis. Specifically, we show the consistency of our method in recovering the low-dimensional noiseless signals, and characterize the effects of the signal-to-noise ratios on the rates of convergence. Under a joint manifolds model framework, we establish the convergence of ultimate embeddings to the eigenfunctions of some newly introduced integral operators. These operators, referred to as duo-landmark integral operators, are defined by the convolutional kernel maps of some reproducing kernel Hilbert spaces (RKHSs). These RKHSs capture the either partially or entirely shared underlying low-dimensional nonlinear signal structures of the two datasets. Our numerical experiments and analyses of two single-cell omics datasets demonstrate the empirical advantages of the proposed method over existing methods in both embeddings and several downstream tasks., Comment: 32 pages, 5 figures; comments are welcome

To realize large-scale quantum information processes, an ideal scheme for two-qubit operations should enable diverse operations with given hardware and physical interaction. However, for spin qubits in semiconductor quantum dots, the common two-qubit operations, including CPhase gates, SWAP gates, and CROT gates, are realized with distinct parameter regions and control waveforms, posing challenges for their simultaneous implementation. Here, taking advantage of the inherent Heisenberg interaction between spin qubits, we propose and verify a fast composite two-qubit gate scheme to extend the available two-qubit gate types as well as reduce the requirements for device properties. Apart from the formerly proposed CPhase (controlled-phase) gates and SWAP gates, theoretical results indicate that the iSWAP-family gate and Fermionic simulation (fSim) gate set are additionally available for spin qubits. Meanwhile, our gate scheme limits the parameter requirements of all essential two-qubit gates to a common J~{\Delta}E_Z region, facilitate the simultaneous realization of them. Furthermore, we present the preliminary experimental demonstration of the composite gate scheme, observing excellent match between the measured and simulated results. With this versatile composite gate scheme, broad-spectrum two-qubit operations allow us to efficiently utilize the hardware and the underlying physics resources, helping accelerate and broaden the scope of the upcoming noise intermediate-scale quantum (NISQ) computing., Comment: 23 pages, 6 figures

We present an unbiased $\lambda$ 3 mm spectral line survey (between 84.5 and 115.8 GHz), conducted by the Purple Mountain Observatory 13.7 meter radio telescope, together with updated modeling results, towards the carbon-rich Asymptotic Giant Branch star, IRC+10216 (CW Leo). A total of 75 spectral lines (96 transitions) are detected, and identified to arise from 19 molecules: C$_2$H, $l$-C$_3$H, C$_4$H, CN, C$_3$N, HC$_3$N, HC$_5$N, HCN, HNC, CH$_3$CN, MgNC, CO, $c$-C$_3$H$_2$, SiC$_2$, SiO, SiS, CS, C$_2$S, C$_3$S, and their isotopologues. Among them, one molecular emission line (H$^{13}$CCCN $J=13-12$) is discovered in IRC+10216 for the first time. The excitation temperature, column density, and fractional abundance of the detected species are deduced by assuming they are in local thermodynamic equilibrium. In addition, the isotopic ratios of [$^{12}$C]/[$^{13}$C], [$^{32}$S]/[$^{34}$S], [$^{28}$Si]/[$^{29}$Si], and [$^{12}$C$^{34}$S]/[$^{13}$C$^{32}$S] are obtained and found to be consistent with previous studies. Finally, we summarize all of the 106 species detected in IRC+10216 to date with their observed and modeled column densities for the convenience of future studies., Comment: 71 pages, 39 figures, 10 tables. Accepted for publication in ApJS

Nonreciprocal transport in uniform systems has attracted great research interest recently and the existing theories mainly focus on the diffusive regime. In this study, we uncover a novel scenario for nonreciprocal charge transport in the ballistic regime enabled by asymmetric band structures of the system. The asymmetry of the bands induces unequal Coulomb potentials within the system as the bias voltage imposed by the electrodes inverts its sign. As a result, the bands undergo different energy shifts as the current flows in opposite directions, giving rise to the nonreciprocity. Utilizing the gauge-invariant nonlinear transport theory, we show that the nonreciprocal transport predominantly originates from the second-order conductance, which violates the Onsager reciprocal relation but fulfills a generalized reciprocal relation similar to that of unidirectional magnetoresistance. The ballistic nonreciprocal transport phenomena differ from the diffusive ones by considering the internal asymmetric Coulomb potential, a factor not accounted for in diffusive cases but undeniably crucial in ballistic scenarios. Our work opens a avenue for implementing nonreciprocal transport in the ballistic regime and provides an alternative perspective for further experimental explorations for nonreciprocal transport., Comment: 7 pages, 4 figures

With one- and two-qubit gate fidelities approaching the fault-tolerance threshold for spin qubits in silicon, how to scale up the architecture and make large arrays of spin qubits become the more pressing challenges. In a scaled-up structure, qubit-to-qubit connectivity has crucial impact on gate counts of quantum error correction and general quantum algorithms. In our toolbox of quantum gates for spin qubits, SWAP gate is quite versatile: it can help solve the connectivity problem by realizing both short- and long-range spin state transfer, and act as a basic two-qubit gate, which can reduce quantum circuit depth when combined with other two-qubit gates. However, for spin qubits in silicon quantum dots, high fidelity SWAP gates have not been demonstrated due to the requirements of large circuit bandwidth and a highly adjustable ratio between the strength of the exchange coupling J and the Zeeman energy difference Delta E_z. Here we demonstrate a fast SWAP gate with a duration of ~25 ns based on quantum dots in isotopically enriched silicon, with a highly adjustable ratio between J and Delta E_z, for over two orders of magnitude in our device. We are also able to calibrate the single-qubit local phases during the SWAP gate by incorporating single-qubit gates in our circuit. By independently reading out the qubits, we probe the anti-correlations between the two spins, estimate the operation fidelity and analyze the dominant error sources for our SWAP gate. These results pave the way for high fidelity SWAP gates, and processes based on them, such as quantum communication on chip and quantum simulation by engineering the Heisenberg Hamiltonian in silicon., Comment: 25 pages, 5 figures

Preserving qubit coherence and maintaining high-fidelity qubit control under complex noise environment is an enduring challenge for scalable quantum computing. Here we demonstrate an addressable fault-tolerant single spin qubit with an average control fidelity of 99.12% via randomized benchmarking on a silicon quantum dot device with an integrated micromagnet. Its dephasing time T2* is 1.025 us and can be enlarged to 264 us by using the Hahn echo technique, reflecting strong low-frequency noise in our system. To break through the noise limitation, we introduce geometric quantum computing to obtain high control fidelity by exploiting its noise-resilient feature. However, the control fidelities of the geometric quantum gates are lower than 99%. According to our simulation, the noise-resilient feature of geometric quantum gates is masked by the heating effect. With further optimization to alleviate the heating effect, geometric quantum computing can be a potential approach to reproducibly achieving high-fidelity qubit control in a complex noise environment., Comment: 10 pages, 8 figures

High-fidelity singlet-triplet state readout is essential for large-scale quantum computing. However, the widely used threshold method of comparing a mean value with the fixed threshold will limit the judgment accuracy, especially for the relaxed triplet state, under the restriction of relaxation time and signal-to-noise ratio. Here, we achieve an enhanced latching readout based on Pauli spin blockade in a Si-MOS double quantum dot device and demonstrate an average singlet-triplet state readout fidelity of 97.59% by the threshold method. We reveal the inherent deficiency of the threshold method for the relaxed triplet state classification and introduce machine learning as a relaxation-independent readout method to reduce the misjudgment. The readout fidelity for classifying the simulated single-shot traces can be improved to 99.67% by machine learning method, better than the threshold method of 97.54% which is consistent with the experimental result. This work indicates that machine learning method can be a strong potential candidate for alleviating the restrictions of stably achieving high-fidelity and high-accuracy singlet-triplet state readout in large-scale quantum computing., Comment: 11 pages,11 figures

Pulse distortion, as one of the coherent error sources, hinders the characterization and control of qubits. In the semiconductor quantum dot system, the distortions on measurement pulses and control pulses disturb the experimental results, while no effective calibration procedure has yet been reported. Here, we demonstrate two different calibration methods to calibrate and correct the distortion using the two-qubit system as a detector. The two calibration methods have different correction accuracy and complexity. One is the coarse predistortion (CPD) method, with which the distortion is partly relieved. The other method is the all predistortion (APD) method, with which we measure the transfer function and significantly improve the exchange oscillation homogeneity. The two methods use the exchange oscillation homogeneity as the metric and are appropriate for any qubit that oscillates with a diabatic pulse. With the APD procedure, an arbitrary control waveform can be accurately delivered to the device, which is essential for characterizing qubits and improving gate fidelity., Comment: 19 pages, 5 figures

Single-cell data integration can provide a comprehensive molecular view of cells, and many algorithms have been developed to remove unwanted technical or biological variations and integrate heterogeneous single-cell datasets. Despite their wide usage, existing methods suffer from several fundamental limitations. In particular, we lack a rigorous statistical test for whether two high-dimensional single-cell datasets are alignable (and therefore should even be aligned). Moreover, popular methods can substantially distort the data during alignment, making the aligned data and downstream analysis difficult to interpret. To overcome these limitations, we present a spectral manifold alignment and inference (SMAI) framework, which enables principled and interpretable alignability testing and structure-preserving integration of single-cell data with the same type of features. SMAI provides a statistical test to robustly assess the alignability between datasets to avoid misleading inference, and is justified by high-dimensional statistical theory. On a diverse range of real and simulated benchmark datasets, it outperforms commonly used alignment methods. Moreover, we show that SMAI improves various downstream analyses such as identification of differentially expressed genes and imputation of single-cell spatial transcriptomics, providing further biological insights. SMAI's interpretability also enables quantification and a deeper understanding of the sources of technical confounders in single-cell data.

Whether SiC$_2$ is a parent species, that is formed in the photosphere or as a by-product of high-temperature dust formation, or a daughter species, formed in a chemistry driven by the photodestruction of parent species in the outer envelope, has been debated for a long time. Here, we analyze the ALMA observations of four SiC$_2$ transitions in the CSEs of three C-rich AGB stars (AI Vol, II Lup, and RAFGL 4211), and found that SiC$_2$ exhibits an annular, shell-like distribution in these targets, suggesting that SiC$_2$ can be a daughter species in the CSEs of carbon-rich AGB stars. The results can provide important references for future chemical models., Comment: Accepted in Frontiers in Astronomy and Space Sciences

It is significant to study congruences involving multiple harmonic sums. Let $p$ be an odd prime, in recent years, the following curious congruence $$\sum_{\substack{i+j+k=p \\ i, j, k>0}} \frac{1}{i j k} \equiv-2 B_{p-3}\pmod p$$ has been generalized along different directions, where $B_n$ denote the $n$th Bernoulli number. In this paper, we obtain several new generalizations of the above congruence by applying congruences involving multiple harmonic sums. For example, we have $$\sum_{\substack{k_1+k_2+\cdots+k_n=p \\ k_i> 0, 1 \le i \le n}} \dfrac{(-1)^{k_1}\left(\dfrac{k_1}{3}\right)}{k_1 \cdots k_n} \equiv \dfrac{(n-1)!}{n}\dfrac{2^{n-1}+1}{3\cdot6^{n-1}}B_{p-n}\left(\dfrac{1}{3}\right)\pmod p,$$ where $n$ is even, $B_n(x)$ denote the Bernoulli polynomials., Comment: 12 pages

The influence of roughness spacing on boundary layer transition over distributed roughness elements is studied using direct numerical simulation (DNS) and global stability analysis, and compared to isolated roughness elements at the same Reh. Small spanwise spacing ($\lambda_z = 2.5h$) inhibits the formation of counter-rotating vortices (CVP) and as a result, hairpin vortices are not generated and the downstream shear layer is steady. For $\lambda_z = 5h$, the CVP and hairpin vortices are induced by the first row of roughness, perturbing the downstream shear layer and causing transition. The temporal periodicity of the primary hairpin vortices appears to be independent of the streamwise spacing. Distributed roughness leads to a lower critical Reh for instability to occur and a more significant breakdown of the boundary layer compared to isolated roughness. When the streamwise spacing is comparable to the region of flow separation ($\lambda_x = 5h$), the high-momentum fluid barely moves downward into the cavities and the wake flow has little impact on the following roughness elements. The leading unstable varicose mode is associated with the central low-speed streaks along the aligned roughness elements, and its frequency is close to the shedding frequency of hairpin vortices in the isolated case. For larger streamwise spacing ($\lambda_x = 10h$), two distinct modes are obtained from global stability analysis. The first mode shows varicose symmetry, corresponding to the primary hairpin vortex shedding induced by the first-row roughness. The high-speed streaks formed in the longitudinal grooves are also found to be unstable and interacting with the varicose mode. The second mode is a sinuous mode with lower frequency, induced as the wake flow of the first-row roughness runs into the second row. It extracts most energy from the spanwise shear between the high- and low-speed streaks.

Let $[x]$ be the integral part of $x$, $n>1$ be a positive integer and $\chi_n$ denote the trivial Dirichlet character modulo $n$. In this paper, we use an identity established by Z. H. Sun to get congruences of $T_{m,k}(n)=\sum_{x=1}^{[n/m]}\frac{\chi_n(x)}{x^k}\left(\bmod n^{r+1}\right)$ for $r\in \{1,2\}$, any positive integer $m $ with $n \equiv \pm 1 \left(\bmod m \right)$ in terms of Bernoulli polynomials. As its an application, we also obtain some new congruences involving binomial coefficients modulo $n^4$ in terms of generalized Bernoulli numbers., Comment: 21pages

Dimension reduction and data visualization aim to project a high-dimensional dataset to a low-dimensional space while capturing the intrinsic structures in the data. It is an indispensable part of modern data science, and many dimensional reduction and visualization algorithms have been developed. However, different algorithms have their own strengths and weaknesses, making it critically important to evaluate their relative performance for a given dataset, and to leverage and combine their individual strengths. In this paper, we propose an efficient spectral method for assessing and combining multiple visualizations of a given dataset produced by diverse algorithms. The proposed method provides a quantitative measure -- the visualization eigenscore -- of the relative performance of the visualizations for preserving the structure around each data point. Then it leverages the eigenscores to obtain a consensus visualization, which has much improved { quality over the individual visualizations in capturing the underlying true data structure.} Our approach is flexible and works as a wrapper around any visualizations. We analyze multiple simulated and real-world datasets from diverse applications to demonstrate the effectiveness of the eigenscores for evaluating visualizations and the superiority of the proposed consensus visualization. Furthermore, we establish rigorous theoretical justification of our method based on a general statistical framework, yielding fundamental principles behind the empirical success of consensus visualization along with practical guidance., Comment: Under revision of Nature Communications

Spin qubits based on silicon metal-oxide semiconductor (Si-MOS) quantum dots (QDs) are promising platforms for large-scale quantum computers. To control spin qubits in QDs, electric dipole spin resonance (EDSR) has been most commonly used in recent years. By delocalizing an electron across a double quantum dots charge state, flopping-mode EDSR has been realized in Si/SiGe QDs. Here, we demonstrate a flopping-mode spin qubit in a Si-MOS QD via Elzerman single-shot readout. When changing the detuning with a fixed drive power, we achieve s-shape spin resonance frequencies, an order of magnitude improvement in the spin Rabi frequencies, and virtually constant spin dephasing times. Our results offer a route to large-scale spin qubit systems with higher control fidelity in Si-MOS QDs., Comment: 5 pages, 4 figures

A central challenge in every field of biology is to use existing measurements to predict the outcomes of future experiments. In this work, we consider the wealth of antibody inhibition data against variants of the influenza virus. Due to this viru's genetic diversity and evolvability, the variants examined in one study will often have little-to-no overlap with other studies, making it difficult to discern common patterns or unify datasets for further analysis. To that end, we develop a computational framework that predicts how an antibody or serum would inhibit any variant from any other study. We use this framework to greatly expand seven influenza datasets utilizing hemagglutination inhibition, validating our method upon 200,000 existing measurements and predicting 2,000,000 new values along with their uncertainties. With these new values, we quantify the transferability between seven vaccination and infection studies in humans and ferrets, show that the serum potency is negatively correlated with breadth, and present a tool for pandemic preparedness. This data-driven approach does not require any information beyond each virus's name and measurements, and even datasets with as few as 5 viruses can be expanded, making this approach widely applicable. Future influenza studies using hemagglutination inhibition can directly utilize our curated datasets to predict newly measured antibody responses against ~80 H3N2 influenza viruses from 1968-2011, whereas immunological studies utilizing other viruses or a different assay only need a single partially-overlapping dataset to extend their work. In essence, this approach enables a shift in perspective when analyzing data from "what you see is what you get" into "what anyone sees is what everyone gets."

The single-shot readout data process is essential for the realization of high-fidelity qubits and fault-tolerant quantum algorithms in semiconductor quantum dots. However, the fidelity and visibility of the readout process is sensitive to the choice of the thresholds and limited by the experimental hardware. By demonstrating the linear dependence between the measured spin state probabilities and readout visibilities along with dark counts, we describe an alternative threshold-independent method for the single-shot readout of spin qubits in semiconductor quantum dots. We can obtain the extrapolated spin state probabilities of the prepared probabilities of the excited spin state through the threshold-independent method. Then, we analyze the corresponding errors of the method, finding that errors of the extrapolated probabilities cannot be neglected with no constraints on the readout time and threshold voltage. Therefore, by limiting the readout time and threshold voltage we ensure the accuracy of the extrapolated probability. Then, we prove that the efficiency and robustness of this method is 60 times larger than that of the most commonly used method. Moreover, we discuss the influence of the electron temperature on the effective area with a fixed external magnetic field and provide a preliminary demonstration for a single-shot readout up to 0.7 K/1.5T in the future., Comment: 18 pages, 6 figures

The past two decades have witnessed the rapid development of personalized recommendation techniques. Despite significant progress made in both research and practice of recommender systems, to date, there is a lack of a widely-recognized benchmarking standard in this field. Many existing studies perform model evaluations and comparisons in an ad-hoc manner, for example, by employing their own private data splits or using different experimental settings. Such conventions not only increase the difficulty in reproducing existing studies, but also lead to inconsistent experimental results among them. This largely limits the credibility and practical value of research results in this field. To tackle these issues, we present an initiative project (namely BARS) aiming for open benchmarking for recommender systems. In comparison to some earlier attempts towards this goal, we take a further step by setting up a standardized benchmarking pipeline for reproducible research, which integrates all the details about datasets, source code, hyper-parameter settings, running logs, and evaluation results. The benchmark is designed with comprehensiveness and sustainability in mind. It covers both matching and ranking tasks, and also enables researchers to easily follow and contribute to the research in this field. This project will not only reduce the redundant efforts of researchers to re-implement or re-run existing baselines, but also drive more solid and reproducible research on recommender systems. We would like to call upon everyone to use the BARS benchmark for future evaluation, and contribute to the project through the portal at: https://openbenchmark.github.io/BARS., Comment: Accepted by SIGIR 2022. Note that version v5 is updated to keep consistency with the ACM camera-ready version

Orbital angular momentum (OAM) multiplexed entangled state is an important quantum resource for high dimensional quantum information processing. In this paper, we experimentally quantify quantum coherence of OAM multiplexed continuous-variable (CV) entangled state and characterize its evolution in a noisy environment. We show that the quantum coherence of the OAM multiplexed CV entangled state carrying topological charges $l=1$ and $l=2$ are the same as that of the Gaussian mode with $l=0$ in a noisy channel. Furthermore, we show that the quantum coherence of OAM multiplexed entangled state is robust to noise, even though the sudden death of entanglement is observed. Our results provide reference for applying quantum coherence of OAM multiplexed CV entangled state in a noisy environment., Comment: 8 pages, 3 figures. Comments are welcome. arXiv admin note: substantial text overlap with arXiv:2203.15976

Orbital angular momentum (OAM) multiplexing provides an efficient method to improve data-carrying capacity in various quantum communication protocols. It is a precondition to distribute OAM multiplexed quantum resources in quantum channels for implementing quantum communication. However, quantum steering of OAM multiplexed optical fields and the effect of channel noise on OAM multiplexed quantum resources remain unclear. Here, we generate OAM multiplexed continuous-variable (CV) entangled states and distribute them in lossy or noisy channels. We show that the decoherence property of entanglement and quantum steering of the OAM multiplexed states carrying topological charges $l=1$ and $l=2$ are the same as that of the Gaussian mode with $l=0$ in lossy and noisy channels. The sudden death of entanglement and quantum steering of high-order OAM multiplexed states is observed in the presence of excess noise. Our results demonstrate the feasibility to realize high data-carrying capacity quantum information processing by utilizing OAM multiplexed CV entangled states., Comment: 9 pages, 7 figures. Comments are welcome

We propose a kernel-spectral embedding algorithm for learning low-dimensional nonlinear structures from high-dimensional and noisy observations, where the datasets are assumed to be sampled from an intrinsically low-dimensional manifold and corrupted by high-dimensional noise. The algorithm employs an adaptive bandwidth selection procedure which does not rely on prior knowledge of the underlying manifold. The obtained low-dimensional embeddings can be further utilized for downstream purposes such as data visualization, clustering and prediction. Our method is theoretically justified and practically interpretable. Specifically, we establish the convergence of the final embeddings to their noiseless counterparts when the dimension and size of the samples are comparably large, and characterize the effect of the signal-to-noise ratio on the rate of convergence and phase transition. We also prove convergence of the embeddings to the eigenfunctions of an integral operator defined by the kernel map of some reproducing kernel Hilbert space capturing the underlying nonlinear structures. Numerical simulations and analysis of three real datasets show the superior empirical performance of the proposed method, compared to many existing methods, on learning various manifolds in diverse applications., Comment: Accepted to the Annals of Statistics

This paper studies the problem of statistical inference for genetic relatedness between binary traits based on individual-level genome-wide association data. Specifically, under the high-dimensional logistic regression models, we define parameters characterizing the cross-trait genetic correlation, the genetic covariance and the trait-specific genetic variance. A novel weighted debiasing method is developed for the logistic Lasso estimator and computationally efficient debiased estimators are proposed. The rates of convergence for these estimators are studied and their asymptotic normality is established under mild conditions. Moreover, we construct confidence intervals and statistical tests for these parameters, and provide theoretical justifications for the methods, including the coverage probability and expected length of the confidence intervals, as well as the size and power of the proposed tests. Numerical studies are conducted under both model generated data and simulated genetic data to show the superiority of the proposed methods. By analyzing a real data set on autoimmune diseases, we demonstrate its ability to obtain novel insights about the shared genetic architecture between ten pediatric autoimmune diseases.

The dynamic factors in the environment will lead to the decline of camera localization accuracy due to the violation of the static environment assumption of SLAM algorithm. Recently, some related works generally use the combination of semantic constraints and geometric constraints to deal with dynamic objects, but problems can still be raised, such as poor real-time performance, easy to treat people as rigid bodies, and poor performance in low dynamic scenes. In this paper, a dynamic scene-oriented visual SLAM algorithm based on object detection and coarse-to-fine static probability named CFP-SLAM is proposed. The algorithm combines semantic constraints and geometric constraints to calculate the static probability of objects, keypoints and map points, and takes them as weights to participate in camera pose estimation. Extensive evaluations show that our approach can achieve almost the best results in high dynamic and low dynamic scenarios compared to the state-of-the-art dynamic SLAM methods, and shows quite high real-time ability.

Motivated by applications in single-cell biology and metagenomics, we investigate the problem of matrix reordering based on a noisy disordered monotone Toeplitz matrix model. We establish the fundamental statistical limit for this problem in a decision-theoretic framework and demonstrate that a constrained least squares estimator achieves the optimal rate. However, due to its computational complexity, we analyze a popular polynomial-time algorithm, spectral seriation, and show that it is suboptimal. To address this, we propose a novel polynomial-time adaptive sorting algorithm with guaranteed performance improvement. Simulations and analyses of two real single-cell RNA sequencing datasets demonstrate the superiority of our algorithm over existing methods., Comment: accepted by IEEE Transactions on Information Theory

Let [t] be the integral part of the real number t and let 1 P be the characteristic function of the primes. Denote by $\pi$ G (x) the number of primes in the floor function set G(x) := {[ x n ] : 1 n x} and by S 1 P (x) the number of primes in the sequence {[ x n ]} n 1. Very recently, Heyman proves

Global stability analysis and direct numerical simulation (DNS) are performed to study boundary layer flows with an isolated roughness element. Wall-attached cuboids with aspect ratios $\eta=1$ and $\eta=0.5$ are investigated for fixed ratio of roughness height to displacement boundary layer thickness $h/\delta^*=2.86$. Global stability analysis is able to capture the frequency of the primary vortical structures. For $\eta=1$, only varicose instability is seen. For the thinner roughness element ($\eta=0.5$), the varicose instability dominates the sinuous instability, and the sinuous instability becomes more pronounced as $Re_h$ increases, due to increased spanwise shear in the near-wake region. The unstable modes mainly extract energy from the central streak, although the lateral streaks also contribute. The DNS results show that different instability features lead to different behavior and development of vortical structures in the nonlinear transition process. For $\eta=1$, the varicose mode is associated with the shedding of hairpin vortices. As $Re_h$ increases, the breakdown of hairpin vortices occurs closer to the roughness and sinuous breakdown behavior promoting transition to turbulence is seen in the farther wake. A fully-developed turbulent flow is established in both the inner and outer layers farther downstream when $Re_h$ is sufficiently high. For $\eta=0.5$, the sinuous wiggling of hairpin vortices is prominent at higher $Re_h$, leading to stronger interactions in the near wake, as a result of combined varicose and sinuous instabilities. A sinuous mode captured by dynamic mode decomposition (DMD) analysis, and associated with the `wiggling' of streaks persists far downstream.

Let $q\ge3$ be an integer, $\chi$ be a Dirichlet character modulo $q$, and $L(s,\chi)$ denote the Dirichlet $L$-functions corresponding to $\chi$. In this paper, we show some special function series, and give some new identities for the Dirichlet $L$-functions involving Gauss sums. Specially, we give specific identities for $L(2,\chi)$., Comment: 10 pages

Blockwise missing data occurs frequently when we integrate multisource or multimodality data where different sources or modalities contain complementary information. In this paper, we consider a high-dimensional linear regression model with blockwise missing covariates and a partially observed response variable. Under this framework, we propose a computationally efficient estimator for the regression coefficient vector based on carefully constructed unbiased estimating equations and a blockwise imputation procedure, and obtain its rate of convergence. Furthermore, building upon an innovative projected estimating equation technique that intrinsically achieves bias-correction of the initial estimator, we propose a nearly unbiased estimator for each individual regression coefficient, which is asymptotically normally distributed under mild conditions. Based on these debiased estimators, asymptotically valid confidence intervals and statistical tests about each regression coefficient are constructed. Numerical studies and application analysis of the Alzheimer's Disease Neuroimaging Initiative data show that the proposed method performs better and benefits more from unsupervised samples than existing methods., Comment: V2: 40 pages, 2 figures. Accepted at Statistica Sinica

This paper investigates the theoretical foundations of the t-distributed stochastic neighbor embedding (t-SNE) algorithm, a popular nonlinear dimension reduction and data visualization method. A novel theoretical framework for the analysis of t-SNE based on the gradient descent approach is presented. For the early exaggeration stage of t-SNE, we show its asymptotic equivalence to power iterations based on the underlying graph Laplacian, characterize its limiting behavior, and uncover its deep connection to Laplacian spectral clustering, and fundamental principles including early stopping as implicit regularization. The results explain the intrinsic mechanism and the empirical benefits of such a computational strategy. For the embedding stage of t-SNE, we characterize the kinematics of the low-dimensional map throughout the iterations, and identify an amplification phase, featuring the intercluster repulsion and the expansive behavior of the low-dimensional map, and a stabilization phase. The general theory explains the fast convergence rate and the exceptional empirical performance of t-SNE for visualizing clustered data, brings forth interpretations of the t-SNE visualizations, and provides theoretical guidance for applying t-SNE and selecting its tuning parameters in various applications., Comment: Accepted by Journal of Machine Learning Research

Let $q>2$ be an odd integer. For each integer $x$ with $0

Tunable synthetic spin-orbit coupling (s-SOC) is one of the key challenges in various quantum systems, such as ultracold atomic gases, topological superconductors, and semiconductor quantum dots. Here we experimentally demonstrate controlling the s-SOC by investigating the anisotropy of spin-valley resonance in a silicon quantum dot. As we rotate the applied magnetic field in-plane, we find a striking nonsinusoidal behavior of resonance amplitude that distinguishes s-SOC from the intrinsic spin-orbit coupling (i-SOC), and associate this behavior with the previously overlooked in-plane transverse magnetic field gradient. Moreover, by theoretically analyzing the experimentally measured s-SOC field, we predict the quality factor of the spin qubit could be optimized if the orientation of the in-plane magnetic field is rotated away from the traditional working point., Comment: 26 pages, 10 figures

This paper studies the high-dimensional mixed linear regression (MLR) where the output variable comes from one of the two linear regression models with an unknown mixing proportion and an unknown covariance structure of the random covariates. Building upon a high-dimensional EM algorithm, we propose an iterative procedure for estimating the two regression vectors and establish their rates of convergence. Based on the iterative estimators, we further construct debiased estimators and establish their asymptotic normality. For individual coordinates, confidence intervals centered at the debiased estimators are constructed. Furthermore, a large-scale multiple testing procedure is proposed for testing the regression coefficients and is shown to control the false discovery rate (FDR) asymptotically. Simulation studies are carried out to examine the numerical performance of the proposed methods and their superiority over existing methods. The proposed methods are further illustrated through an analysis of a dataset of multiplex image cytometry, which investigates the interaction networks among the cellular phenotypes that include the expression levels of 20 epitopes or combinations of markers.

Operation speed and coherence time are two core measures for the viability of a qubit. Strong spin-orbit interaction (SOI) and relatively weak hyperfine interaction make holes in germanium (Ge) intriguing candidates for spin qubits with rapid, all-electrical coherent control. Here we report ultrafast single-spin manipulation in a hole-based double quantum dot in a germanium hut wire (GHW). Mediated by the strong SOI, a Rabi frequency exceeding 540 MHz is observed at a magnetic field of 100 mT, setting a record for ultrafast spin qubit control in semiconductor systems. We demonstrate that the strong SOI of heavy holes (HHs) in our GHW, characterized by a very short spin-orbit length of 1.5 nm, enables the rapid gate operations we accomplish. Our results demonstrate the potential of ultrafast coherent control of hole spin qubits to meet the requirement of DiVincenzo's criteria for a scalable quantum information processor., Comment: 19 pages, 4 figures

Modularity is a popular metric for quantifying the degree of community structure within a network. The distribution of the largest eigenvalue of a network's edge weight or adjacency matrix is well studied and is frequently used as a substitute for modularity when performing statistical inference. However, we show that the largest eigenvalue and modularity are asymptotically uncorrelated, which suggests the need for inference directly on modularity itself when the network size is large. To this end, we derive the asymptotic distributions of modularity in the case where the network's edge weight matrix belongs to the Gaussian Orthogonal Ensemble, and study the statistical power of the corresponding test for community structure under some alternative model. We empirically explore universality extensions of the limiting distribution and demonstrate the accuracy of these asymptotic distributions through type I error simulations. We also compare the empirical powers of the modularity based tests with some existing methods. Our method is then used to test for the presence of community structure in two real data applications.

According to the World Health Organization, the number of mental disorder patients, especially depression patients, has grown rapidly and become a leading contributor to the global burden of disease. However, the present common practice of depression diagnosis is based on interviews and clinical scales carried out by doctors, which is not only labor-consuming but also time-consuming. One important reason is due to the lack of physiological indicators for mental disorders. With the rising of tools such as data mining and artificial intelligence, using physiological data to explore new possible physiological indicators of mental disorder and creating new applications for mental disorder diagnosis has become a new research hot topic. However, good quality physiological data for mental disorder patients are hard to acquire. We present a multi-modal open dataset for mental-disorder analysis. The dataset includes EEG and audio data from clinically depressed patients and matching normal controls. All our patients were carefully diagnosed and selected by professional psychiatrists in hospitals. The EEG dataset includes not only data collected using traditional 128-electrodes mounted elastic cap, but also a novel wearable 3-electrode EEG collector for pervasive applications. The 128-electrodes EEG signals of 53 subjects were recorded as both in resting state and under stimulation; the 3-electrode EEG signals of 55 subjects were recorded in resting state; the audio data of 52 subjects were recorded during interviewing, reading, and picture description. We encourage other researchers in the field to use it for testing their methods of mental-disorder analysis.

Driven by a wide range of applications, many principal subspace estimation problems have been studied individually under different structural constraints. This paper presents a unified framework for the statistical analysis of a general structured principal subspace estimation problem which includes as special cases non-negative PCA/SVD, sparse PCA/SVD, subspace constrained PCA/SVD, and spectral clustering. General minimax lower and upper bounds are established to characterize the interplay between the information-geometric complexity of the structural set for the principal subspaces, the signal-to-noise ratio (SNR), and the dimensionality. The results yield interesting phase transition phenomena concerning the rates of convergence as a function of the SNRs and the fundamental limit for consistent estimation. Applying the general results to the specific settings yields the minimax rates of convergence for those problems, including the previous unknown optimal rates for non-negative PCA/SVD, sparse SVD and subspace constrained PCA/SVD.

Let $q\ge3$ be an integer, $\chi$ denote a Dirichlet character modulo $q$, for any real number $a\ge 0$, we define the generalized Dirichlet $L$-functions $$ L(s,\chi,a)=\sum_{n=1}^{\infty}\frac{\chi(n)}{(n+a)^s}, $$ where $s=\sigma+it$ with $\sigma>1$ and $t$ both real. It can be extended to all $s$ by analytic continuation. For any integer $m$, the famous Gauss sum $G(m,\chi)$ is defined as follows: $$G(m,\chi)=\sum_{a=1}^{q}\chi(a)e\left(\frac{am}{q}\right), $$ where $e(y)=e^{2\pi iy}$. The main purpose of this paper is to use the analytic method to study the mean value properties of the generalized Dirichlet $L$-functions with the weight of the Gauss Sums, and obtain a sharp asymptotic formula., Comment: 14 pages

In silicon quantum dots (QDs), at a certain magnetic field commonly referred to as the "hot spot", the electron spin relaxation rate (T_1^(-1)) can be drastically enhanced due to strong spin-valley mixing. Here, we experimentally find that with a valley splitting of 78.2 ${\pm}$ 1.6 ${\mu}$eV, this hot spot in spin relaxation can be suppressed by more than 2 orders of magnitude when the in-plane magnetic field is oriented at an optimal angle, about 9{\deg} from the [100] sample plane. This directional anisotropy exhibits a sinusoidal modulation with a 180{\deg} periodicity. We explain the magnitude and phase of this modulation using a model that accounts for both spin-valley mixing and intravalley spin-orbit mixing. The generality of this phenomenon is also confirmed by tuning the electric field and the valley splitting up to 268.2 ${\pm}$ 0.7 ${\mu}$eV., Comment: 13 pages and 4 figures for Manuscript, 17 pages and 7 figures for Supplementary material

Motivated by the problem of estimating the bacterial growth rates for genome assemblies from shotgun metagenomic data, we consider the permuted monotone matrix model $Y=\Theta\Pi+Z$, where $Y\in \mathbb{R}^{n\times p}$ is observed, $\Theta\in \mathbb{R}^{n\times p}$ is an unknown approximately rank-one signal matrix with monotone rows, $\Pi \in \mathbb{R}^{p\times p}$ is an unknown permutation matrix, and $Z\in \mathbb{R}^{n\times p}$ is the noise matrix. This paper studies the estimation of the extreme values associated to the signal matrix $\Theta$, including its first and last columns, as well as their difference. Treating these estimation problems as compound decision problems, minimax rate-optimal estimators are constructed using the spectral column sorting method. Numerical experiments through simulated and synthetic microbiome metagenomic data are presented, showing the superiority of the proposed methods over the alternatives. The methods are illustrated by comparing the growth rates of gut bacteria between inflammatory bowel disease patients and normal controls.

Motivated by recent research on quantifying bacterial growth dynamics based on genome assemblies, we consider a permuted monotone matrix model $Y=\Theta\Pi+Z$, where the rows represent different samples, the columns represent contigs in genome assemblies and the elements represent log-read counts after preprocessing steps and Guanine-Cytosine (GC) adjustment. In this model, $\Theta$ is an unknown mean matrix with monotone entries for each row, $\Pi$ is a permutation matrix that permutes the columns of $\Theta$, and $Z$ is a noise matrix. This paper studies the problem of estimation/recovery of $\Pi$ given the observed noisy matrix $Y$. We propose an estimator based on the best linear projection, which is shown to be minimax rate-optimal for both exact recovery, as measured by the 0-1 loss, and partial recovery, as quantified by the normalized Kendall's tau distance. Simulation studies demonstrate the superior empirical performance of the proposed estimator over alternative methods. We demonstrate the methods using a synthetic metagenomics dataset of 45 closely related bacterial species and a real metagenomic dataset to compare the bacterial growth dynamics between the responders and the non-responders of the IBD patients after 8 weeks of treatment.

Fresh water scarcity is one of the critical challenges for global sustainable development. Several novel water resources such as passive seawater solar desalination and atmospheric water harvesting have made some progress in recent years. However, no investigation has referred to harvesting water from shallow subsurface soils, which are potential huge water reservoirs. Here, we introduce a method of solar-driven water harvesting from soils, which can provide cheap fresh water in impoverished, arid and decentralized areas. The concentrated solar energy is used to heat the soils to evaporate the soil moisture. Then vapors flow to the condenser through tubes and condense as freshwater. Sustainable water harvesting is realized by water migration due to capillary pumping effect within soils. In the laboratory condition, an experimental setup is designed and its water-harvesting ability from soils is investigated. The maximum water mass harvesting rate was 99.8 g h-1. In about 12 h, the total harvesting water could be as high as about 900 ml. The whole process is solar-driven and spontaneous without other mechanical or electrical ancillaries. The water harvesting rate under one sun energy flux (1 kW m-2) is estimated to be about 360 g h-1 with a 1 m2 solar concentrator. Our proposal provides a potential onsite and sustainable fresh water supply solution to deal with the water scarcity problem., Comment: 26pages,7 figures

Nondegenerate four-wave mixing (FWM) process based on a double-$\Lambda$ scheme in hot alkali metal vapor is a versatile tool in quantum state engineering, quantum imaging, and quantum precision measurements. In this Letter, we investigate the generation of quantum correlated twin beams which carry nonzero orbital angular momentums (OAMs) based on the FWM process in hot cesium vapor. The amplified probe beam and the newly generated conjugate beam in the FWM process have the same and opposite topological charge as the seed beam, respectively. We also explore the FWM process operated in a nonamplifying regime where quantum correlated twin beams carrying OAMs can still be generated. In this regime, the FWM process plays the role of quantum beam splitter for the OAM of light, that is, a device that can split a coherent light beam carrying OAM into quantum-correlated twin beams carrying OAMs. More generally, our setup can be used as a quantum beam splitter of images., Comment: 5 pages, 5 figures

Let $q\ge3$ be an integer, $\chi$ denote a Dirichlet character modulo $q$, for any real number $a\ge 0$, we define the generalized Dirichlet $L$-functions $$ L(s,\chi,a)=\sum_{n=1}^{\infty}\frac{\chi(n)}{(n+a)^s}, $$ where $s=\sigma+it$ with $\sigma>1$ and $t$ both real. It can be extended to all $s$ by analytic continuation. In this paper, we study the mean value properties of the generalized Dirichlet $L$-functions, and obtain several sharp asymptotic formulae by using analytic method., Comment: 15 pages,accepted by IJNT