Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
62works
0followers
31topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

62 published item(s)

preprint2026arXiv

Revisiting The Cosmological Time Dilation of Distant Quasars: Influence of Source Properties and Evolution

After decades of searching, cosmological time dilation was recently identified in the timescale of variability seen in distant quasars. Here, we expand on the previous analysis to disentangle this cosmological signal from the influence of the properties of the source population, specifically the quasar bolometric luminosity and the rest-frame emission wavelength at which the variability was observed. Furthermore, we consider the potential influence of the evolution of the quasar population over cosmic time. We find that a significant intrinsic scatter of 0.288 +- 0.021 dex in the variability timescales, which was not considered in the previous analysis, is favoured by the data. This slightly increases the uncertainty in the results. However, the expected cosmological dependence of the variability timescales is confirmed to be robust to changes in the underlying assumptions. We find that the variability timescales increase smoothly with both wavelength and bolometric luminosity, and that black hole mass has no effect on the variability timescale once rest wavelength and bolometric luminosity are accounted for. Moreover, if the standard cosmological model is correct, governed by relativistic expansion, we also find very little cosmological evolution in the intrinsic variability timescales of distant quasars.

preprint2026arXiv

Robust Multimodal Recommendation via Graph Retrieval-Enhanced Modality Completion

Multimodal data plays a critical role in web-based recommendation systems, where information from diverse modalities such as vision and text enhances representation learning. However, real-world multimodal datasets often suffer from modality incompleteness due to sensor failures, annotation scarcity, or privacy constraints, which substantially degrade model performance and reliability. One effective solution to address this issue is modality completion, which reconstructs missing features to provide modality-complete graphs for downstream tasks. Given a query node with missing multimodal features, existing modality completion methods typically infer information from the node itself or its neighbors to reconstruct the missing modality. However, these methods may overlook semantically relevant context in the graph, which contains valuable cues that are non-trivial to capture through simple methods like neighborhood aggregation. In this work, we propose GRE-MC, a Graph Retrieval-Enhanced Modality Completion framework, to overcome these limitations. By introducing a modality-aware subgraph retrieval mechanism, GRE-MC selects semantically relevant subgraphs from the entire graph, providing richer contextual information for completing missing modalities. Subsequently, a graph transformer jointly encodes the query node and the retrieved subgraph via global attention to complete the missing features, while a learnable sparse-routing codebook regularizes latent embeddings into compact bases for improved robustness. Extensive experiments on multimodal recommendation benchmarks demonstrate that GRE-MC consistently outperforms state-of-the-art methods, validating the effectiveness of subgraph retrieval and joint-encoding graph transformer for robust modality completion.

preprint2026arXiv

Rotation of the Globular Cluster Population of the Dark Matter Deficient Galaxy NGC 1052-DF4: Implication for the Total Mass

We explore the globular cluster population of NGC 1052-DF4, a dark matter deficient galaxy, using Bayesian inference to search for the presence of rotation. The existence of such a rotating component is relevant to the estimation of the mass of the galaxy, and therefore the question of whether NGC 1052-DF4 is truly deficient of dark matter, similar to NGC 1052-DF2 another galaxy in the same group. The rotational characteristics of seven globular clusters in NGC 1052-DF4 were investigated, finding that a non-rotating kinematic model has a higher Bayesian evidence than a rotating model, by a factor of approximately 2.5. In addition, we find that under the assumption of rotation, its amplitude must be small. This distinct lack of rotation strengthens the case that, based on its intrinsic velocity dispersion, NGC 1052-DF4 is a truly dark matter deficient galaxy.

preprint2026arXiv

Rotational Kinematics in the Globular Cluster System of M31: Insights from Bayesian Inference

As ancient stellar systems, globular clusters (GCs) offer valuable insights into the dynamical histories of large galaxies. Previous studies of GC populations in the inner and outer regions of the Andromeda Galaxy (M31) have revealed intriguing subpopulations with distinct kinematic properties. Here, we build upon earlier studies by employing Bayesian modelling to investigate the kinematics of the combined inner and outer GC populations of M31. Given the heterogeneous nature of the data, we examine subpopulations defined by GCs' metallicity and by associations with substructure, in order to characterise possible relationships between the inner and outer GC populations. We find that lower-metallicity GCs and those linked to substructures exhibit a common, more rapid rotation, whose alignment is distinct from that of higher-metallicity and non-substructure GCs. Furthermore, the higher-metallicity GCs rotate in alignment with Andromeda's stellar disk. These pronounced kinematic differences reinforce the idea that different subgroups of GCs were accreted to M31 at distinct epochs, shedding light on the complex assembly history of the galaxy.

preprint2026arXiv

The Hilbert-Schmidt norms of quantum channels and matrix integrals over the unit sphere

The dynamics of quantum systems are generally described by a family of quantum channels (linear, completely positive and trace preserving maps). In this note, we mainly study the range of all possible values of $\|\mathcal{E}\|_2^2+\|\widetilde{\mathcal{E}}\|_2^2$ for quantum channels $\mathcal{E}$ and give the equivalent characterizations for quantum channels that achieve these maximum and minimum values, respectively, where $\|\mathcal{E}\|_2$ is the Hilbert-Schmidt norm of $\mathcal{E}$ and $\widetilde{\mathcal{E}}$ is a complementary channel of $\mathcal{E}.$ Also, we get a concrete description of completely positive maps on infinite dimensional systems preserving pure states. Moreover, the equivalency of several matrix integrals over the unit sphere is demonstrated and some extensions of these matrix integrals are obtained.

preprint2025arXiv

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Black holes (BH), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to tera-electronvolt (TeV) gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectrum extending to the ultrahigh-energy (UHE; photon energy $E>100$ TeV) band and one microquasar (Cygnus X-1) of spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory (LHAASO). Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles out of 1 peta-electronvolt (PeV), suggesting microquasars to be important contributors to Galactic cosmic rays especially around the `knee' region.

preprint2023arXiv

The Covering Radius of the Third-Order Reed-Muller Code RM(3,7) is 20

We prove the covering radius of the third-order Reed-Muller code RM(3,7) is 20, which was previously known to be between 20 and 23 (inclusive). The covering radius of RM(3, 7) is the maximum third-order nonlinearity among all 7-variable Boolean functions. It was known that there exist 7-variable Boolean functions with third-order nonlinearity 20. We prove the third-order nonlinearity cannot achieve 21. According to the classification of the quotient space of RM(6,6)/RM(3,6), we classify all 7-variable Boolean functions into 66 types. Firstly, we prove 62 types (among 66) cannot have third-order nonlinearity 21; Secondly, we prove function of the remaining 4 types can be transformed into a type (6, 10) function, if its third-order nonlinearity is 21; Finally, we transform type (6, 10) functions into a specific form, and prove the functions in that form cannot achieve third-order nonlinearity 21 (with the assistance of computers). By the way, we prove that the affine transformation group over any finite field can be generated by two elements.

preprint2022arXiv

BBA-net: A bi-branch attention network for crowd counting

In the field of crowd counting, the current mainstream CNN-based regression methods simply extract the density information of pedestrians without finding the position of each person. This makes the output of the network often found to contain incorrect responses, which may erroneously estimate the total number and not conducive to the interpretation of the algorithm. To this end, we propose a Bi-Branch Attention Network (BBA-NET) for crowd counting, which has three innovation points. i) A two-branch architecture is used to estimate the density information and location information separately. ii) Attention mechanism is used to facilitate feature extraction, which can reduce false responses. iii) A new density map generation method combining geometric adaptation and Voronoi split is introduced. Our method can integrate the pedestrian's head and body information to enhance the feature expression ability of the density map. Extensive experiments performed on two public datasets show that our method achieves a lower crowd counting error compared to other state-of-the-art methods.

preprint2022arXiv

Blowup dynamics for inhomogeneous mass critical half-wave equation

We consider the focusing inhomogeneous mass critical half-wave equation in one dimension. Under the mild conditions of the inhomogeneous factor, we show that the existence of the radial blowup solutions with ground state mass $\|u_0\|_2=\|Q\|_2$, where $Q$ is the unique positive ground state solution of equation $DQ+Q=Q^3$ and obtain the blowup rate $\|D^{\frac{1}{2}}u(t)\|_{L^2}\sim\frac{1}{|t|}$ as $t\nearrow 0^-$.

preprint2022arXiv

Classification of solutions to equations involving Higher-order fractional Laplacian

In this paper, we are concerned with the following equation involving higher-order fractional Lapalacian \begin{equation*} \left\{\begin{aligned} &(-Δ)^{p+{\fracα{2}}}u(x)=u_+^γ~~ \mbox{ in }\mathbb{R}^n,\\ &\int_{\mathbb{R}^n}u_+^γdx<+\infty, \end{aligned}\right. \end{equation*} where $p\geq 1$ is an integer, $0<\alp<2$, $n> 2p+α$ and $γ\in (1,\frac{n}{n-2p-\alp})$. We establish an integral representation formula for any nonconstant classical solution satisfying certain growth at infinity. From this we prove that these solutions are radially symmetric about some point in $\R^n$ and monotone decreasing in the radial direction via method of moving planes in integral forms.

preprint2022arXiv

Classification of solutions to several semi-linear polyharmonic equations and fractional equations

We are concerned with the following semi-linear polyharmonic equation with integral constraint \begin{align} \left\{\begin{array}{rl} &(-Δ)^pu=u^γ_+ ~~ \mbox{ in }{\mathbb{R}^n},\\ \nonumber &\int_{\mathbb{R}^n}u_+^γdx<+\infty, \end{array}\right. \end{align} where $n>2p$, $p\geq2$ and $p\in\mathbb{Z}$. We obtain for $γ\in(1,\frac{n}{n-2p})$ that any nonconstant solution satisfying certain growth at infinity is radial symmetric about some point in $\mathbb{R}^{n}$ and monotone decreasing in the radial direction. In the case $p=2$, the same results are established for more general exponent $γ\in(1,\frac{n+4}{n-4})$. For the following fractional equation with integral constraint \begin{equation*} \left\{\begin{array}{rl} &(-Δ)^sv=v^γ_+ ~~ \mbox{ in }{\mathbb{R}^n},~~~~\\ &\int_{\mathbb{R}^n}v_+^{\frac{n(γ-1)}{2s}}dx<+\infty,~~~~~ \end{array}\right. \end{equation*} where $s\in(0,1)$, $γ\in (1, \frac{n+2s}{n-2s})$ and $n\geq 2$, we also complete the classification of solutions with certain growth at infinity. In addition, observe that the assumptions of the maximum principle named decay at infinity in \cite{chen} can be weakened slightly. Based on this observation, we classify all positive solutions of two semi-linear fractional equations without integral constraint.

preprint2022arXiv

Connection-oriented and Connectionless Quantum Internet Considering Quantum Repeaters

With the rapid development of quantum information and technology in recent years, the construction of quantum internet for interconnecting all kinds of quantum devices, such as quantum processors and sensors, will be the next trend for practical quantum applications. In this paper, we propose the protocols for construction of connection-oriented and connectionless quantum networks by considering the concrete quantum repeater (QR) nodes. Four classes of QRs networks are considered first and designed with two types of protocols in link layer, i.e. simultaneous and one-by-one link. Based on those two link models, the connection-oriented protocol is presented for all classes of QRs networks and the connectionless protocol is proposed for the first, second and third classes QRs networks by only one-by-one link. Furthermore, we introduce a new hybrid connection model of quantum networks combined with connection-oriented and connectionless for practical uses. Our work is a new attempt to study the model of the network layer for different kinds of QR networks and paves the way for developing the protocol stack of universal large-scale quantum internet.

preprint2022arXiv

Deterministic Identification over Channels without CSI

Identification capacities of randomized and deterministic identification were proved to exceed channel capacity for Gaussian channels \emph{with} channel side information (CSI). In this work, we extend deterministic identification to the block fading channels without CSI by applying identification codes for both channel estimation and user identification. We prove that identification capacity is asymptotically higher than transmission capacity even in the absence of CSI. And we also analyze the finite-length performance theoretically and numerically. The simulation results verify the feasibility of the proposed blind deterministic identification in finite blocklength regime.

preprint2022arXiv

Enhancing and Dissecting Crowd Counting By Synthetic Data

In this article, we propose a simulated crowd counting dataset CrowdX, which has a large scale, accurate labeling, parameterized realization, and high fidelity. The experimental results of using this dataset as data enhancement show that the performance of the proposed streamlined and efficient benchmark network ESA-Net can be improved by 8.4\%. The other two classic heterogeneous architectures MCNN and CSRNet pre-trained on CrowdX also show significant performance improvements. Considering many influencing factors determine performance, such as background, camera angle, human density, and resolution. Although these factors are important, there is still a lack of research on how they affect crowd counting. Thanks to the CrowdX dataset with rich annotation information, we conduct a large number of data-driven comparative experiments to analyze these factors. Our research provides a reference for a deeper understanding of the crowd counting problem and puts forward some useful suggestions in the actual deployment of the algorithm.

preprint2022arXiv

Large-scale protein-protein post-translational modification extraction with distant supervision and confidence calibrated BioBERT

Protein-protein interactions (PPIs) are critical to normal cellular function and are related to many disease pathways. However, only 4% of PPIs are annotated with PTMs in biological knowledge databases such as IntAct, mainly performed through manual curation, which is neither time nor cost-effective. We use the IntAct PPI database to create a distant supervised dataset annotated with interacting protein pairs, their corresponding PTM type, and associated abstracts from the PubMed database. We train an ensemble of BioBERT models - dubbed PPI-BioBERT-x10 to improve confidence calibration. We extend the use of ensemble average confidence approach with confidence variation to counteract the effects of class imbalance to extract high confidence predictions. The PPI-BioBERT-x10 model evaluated on the test set resulted in a modest F1-micro 41.3 (P =5 8.1, R = 32.1). However, by combining high confidence and low variation to identify high quality predictions, tuning the predictions for precision, we retained 19% of the test predictions with 100% precision. We evaluated PPI-BioBERT-x10 on 18 million PubMed abstracts and extracted 1.6 million (546507 unique PTM-PPI triplets) PTM-PPI predictions, and filter ~ 5700 (4584 unique) high confidence predictions. Of the 5700, human evaluation on a small randomly sampled subset shows that the precision drops to 33.7% despite confidence calibration and highlights the challenges of generalisability beyond the test set even with confidence calibration. We circumvent the problem by only including predictions associated with multiple papers, improving the precision to 58.8%. In this work, we highlight the benefits and challenges of deep learning-based text mining in practice, and the need for increased emphasis on confidence calibration to facilitate human curation efforts.

preprint2022arXiv

Lattice QCD Study of Transverse-Momentum Dependent Soft Function

In this work, we perform a lattice QCD study of the intrinsic, rapidity-independent soft function within the framework of large momentum effective theory. The computation is carried out using a gauge ensemble of $N_f=2+1+1$ clover-improved twisted mass fermion. After applying an appropriate renormalization procedure and the removal of significant higher-twist contamination, we obtain the intrinsic soft function that is comparable to the one-loop perturbative result at large external momentum. The determination of the nonperturbative soft function from first principles is crucial to sharpen our understanding of the processes with small transverse momentum such as the Drell-Yan production and the semi-inclusive deep inelastic scattering. Additionally, we calculate the Collins-Soper evolution kernel using the quasi-transverse-momentum-dependent wave function as input.

preprint2022arXiv

Limit theorems for the realised semicovariances of multivariate Brownian semistationary processes

In this article we will introduce the realised semicovariance for Brownian semistationary (BSS) processes, which is obtained from the decomposition of the realised covariance matrix into components based on the signs of the returns, and study its in-fill asymptotic properties. More precisely, a weak convergence in the space of càdlàg functions endowed with the Skorohod topology for the realised semicovariance of a general Gaussian process with stationary increments is proved first. The methods are based on Breuer-Major theorems and on a moment bound for sums of products of Gaussian vector&#39;s functions. Furthermore, we establish a corresponding stable convergence. Finally, a weak law of large numbers and a central limit theorem for the realised semicovariance of multivariate BSS processes are established. These results extend the limit theorems for the realised covariation to a result for non-linear functionals.

preprint2022arXiv

Magnetic molecular orbitals in MnSi

A large body of knowledge about magnetism is attained from models of interacting spins, which usually reside on magnetic ions. Proposals beyond the ionic picture are uncommon and seldom verified by direct observations in conjunction with microscopic theory. Here, using inelastic neutron scattering to study the itinerant near-ferromagnet MnSi, we find that the system&#39;s fundamental magnetic units are interconnected, extended molecular orbitals consisting of three Mn atoms each, rather than individual Mn atoms. This result is further corroborated by magnetic Wannier orbitals obtained by ab initio calculations. It contrasts the ionic picture with a concrete example, and presents a novel regime of the spin waves where the wavelength is comparable to the spatial extent of the molecular orbitals. Our discovery brings important insights into not only the magnetism of MnSi, but also a broad range of magnetic quantum materials where structural symmetry, electron itinerancy and correlations act in concert.

preprint2022arXiv

Numerical Investigation of Dynamical and Morphological Trends in Relativistic Jets

Active galactic nuclei (AGN) show a range of morphologies and dynamical properties, which are determined not only by parameters intrinsic to the central engine but also their interaction with the surrounding environment. We investigate the connection of kiloparsec scale AGN jet properties to their intrinsic parameters and surroundings. This is done using a suite of 40 relativistic hydrodynamic simulations spanning a wide range of engine luminosities and opening angles. We explore AGN jet propagation with different ambient density profiles, including $r^{-2}$ (self-similar solution) and $r^{-1}$, which is more relevant for AGN host environments. The Fanaroff-Riley (FR) morphological dichotomy arises naturally in our models. Jets with low energy density compared to the ambient medium produce a center-brightened emissivity distribution, while emissivity from relatively higher energy density jets is dominated by a terminal bright spot. We observe recollimation shocks in our simulations that can generate bright spots along the spine of the jet, providing a possible explanation for &#34;knots&#34; observed in AGN jets. We additionally find a scaling relation between the number of knots and the jet-head-to-surroundings energy density ratio. This scaling relation is generally consistent with the observations of the jets in M87 and Cygnus A. Our model also correctly predicts M87 as FR I and Cygnus A as FR II. Our model can be used to relate jet dynamical parameters such as jet head velocity, jet opening angle, and external pressure to jet power and ambient density estimates.

preprint2022arXiv

On mass - critical NLS with local and non-local nonlinearities

We consider the following nonlinear Schrödinger equation with the double $L^2$-critical nonlinearities \begin{align*} iu_t+Δu+|u|^\frac{4}{3}u+μ\left(|x|^{-2}*|u|^2\right)u=0\ \ \ \text{in $\mathbb{R}^3$,} \end{align*} where $μ>0$ is small enough. Our first goal is to prove the existence and the non-degeneracy of the ground state $Q_μ$. In particular, we develop an appropriate perturbation approach to prove the radial non-degeneracy property and then obtain the general non-degeneracy of the ground state $Q_μ$. We then show the existence of finite time blowup solution with minimal mass $\|u_0\|_{L^2}=\|Q_μ\|_{L^2}$. More precisely, we construct the minimal mass blowup solutions that are parametrized by the energy $E_μ(u_0)>0$ and the momentum $P_μ(u_0)$. In addition, the non-degeneracy property plays crucial role in this construction.

preprint2022arXiv

PACSan: Enforcing Memory Safety Based on ARM PA

Memory safety is a key security property that stops memory corruption vulnerabilities. Existing sanitizers enforce checks and catch such bugs during development and testing. However, they either provide partial memory safety or have overwhelmingly high performance overheads. Our novel sanitizer PACSan enforces spatial and temporal memory safety with no false positives at low performance overheads. PACSan removes the majority of the overheads involved in pointer tracking by sealing metadata in pointers through ARM PA (Pointer Authentication), and performing the memory safety checks when pointers are dereferenced. We have developed a prototype of PACSan and systematically evaluated its security and performance on the Magma, Juliet, Nginx, and SPEC CPU2017 test suites, respectively. In our evaluation, PACSan shows no false positives together with negligible false negatives, while introducing stronger security guarantees and lower performance overheads than state-of-the-art sanitizers, including HWASan, ASan, SoftBound+CETS, Memcheck, LowFat, and PTAuth. Specifically, PACSan has 0.84x runtime overhead and 1.92x memory overhead on average. Compared to the widely deployed ASan, PACSan has no false positives and much fewer false negatives and reduces 7.172% runtime overheads and 89.063%memory overheads.

preprint2022arXiv

Performance Analysis for Covert Communications Under Faster-than-Nyquist Signaling

In this letter, we analyze the performance of covert communications under faster-than-Nyquist (FTN) signaling in the Rayleigh block fading channel. Both Bayesian criterion- and Kullback-Leibler (KL) divergence-based covertness constraints are considered. Especially, for KL divergence-based one, we prove that both the maximum transmit power and covert rate under FTN signaling are higher than those under Nyquist signaling. Numerical results coincide with our analysis and validate the advantages of FTN signaling to realize covert data transmission.

preprint2022arXiv

Preferred Spin Excitations in the Bilayer Iron-Based Superconductor CaK(Fe$_{0.96}$Ni$_{0.04}$)$_4$As$_4$ with Spin-Vortex Crystal Order

The spin-orbit coupling (SOC) is a key to understand the magnetically driven superconductivity in iron-based superconductors, where both local and itinerant electrons are present and the orbital angular momentum is not completely quenched. Here, we report a neutron scattering study on the bilayer compound CaK(Fe$_{0.96}$Ni$_{0.04}$)$_4$As$_4$ with superconductivity coexisting with a non-collinear spin-vortex crystal magnetic order that preserves the tetragonal symmetry of Fe-Fe plane. In the superconducting state, two spin resonance modes with odd and even $L$ symmetries due to the bilayer coupling are found similar to the undoped compound CaKFe$_4$As$_4$ but at lower energies. Polarization analysis reveals that the odd mode is $c-$axis polarized, and the low-energy spin anisotropy can persist to the paramagnetic phase at high temperature, which closely resembles other systems with in-plane collinear and $c-$axis biaxial magnetic orders. These results provide the missing piece of the puzzle on the SOC effect in iron-pnictide superconductors, and also establish a common picture of $c-$axis preferred magnetic excitations below $T_c$ regardless of the details of magnetic pattern or lattice symmetry.

preprint2022arXiv

Switching modulation of spin transport in ferromagnetic tetragonal silicene

We study the band structure and transport properties of ferromagnetic tetragonal silicene nanoribbons by using the non-equilibrium Green&#39;s function method. The band structure and spin-dependent conductance are discussed under the combined effect of the external electric field, potential energy, exchange field and the spin-orbit coupling. One can easily realize a phase transition from a semimetallic to a semiconducting state by changing the transverse width of the nanoribbon. Separation of spin-dependent conductances arises from the effect of exchange field and the spin-orbit coupling, while zero-conductance behaviors exhibit spin-dependent band gaps induced by the electric field. We propose a device configuration of four-terminal tetragonal silicene nanoribbon with two central channels. It is found that spin current can be controlled by utilizing two switches. The switch with a high potential barrier can block electrons flowing from the central scattering region into other terminals. Interestingly, applying only one switch can realize spin-dependent zero conductance and large spin polarization. Two switches can provide multiple operations for controlling spin-dependent transport properties. The two-channel ferromagnetic tetragonal silicene nanoribbon can realize an effective separation of spin current, which may be a potential candidate for spintronic devices.

preprint2022arXiv

The Complete SC-invariant Affine Automorphisms of Polar Codes

Automorphism ensemble (AE) decoding for polar codes was proposed by decoding permuted codewords with successive cancellation (SC) decoders in parallel and hence has lower latency compared to that of successive cancellation list (SCL) decoding. However, some automorphisms are SC-invariant, thus are redundant in AE decoding. In this paper, we find a necessary and sufficient condition related to the block lower-triangular structure of transformation matrices to identify SC-invariant automorphisms. Furthermore, we provide an algorithm to determine the complete SC-invariant affine automorphisms under a specific polar code construction.

preprint2022arXiv

The OARF Benchmark Suite: Characterization and Implications for Federated Learning Systems

This paper presents and characterizes an Open Application Repository for Federated Learning (OARF), a benchmark suite for federated machine learning systems. Previously available benchmarks for federated learning have focused mainly on synthetic datasets and use a limited number of applications. OARF mimics more realistic application scenarios with publicly available data sets as different data silos in image, text and structured data. Our characterization shows that the benchmark suite is diverse in data size, distribution, feature distribution and learning task complexity. The extensive evaluations with reference implementations show the future research opportunities for important aspects of federated learning systems. We have developed reference implementations, and evaluated the important aspects of federated learning, including model accuracy, communication cost, throughput and convergence time. Through these evaluations, we discovered some interesting findings such as federated learning can effectively increase end-to-end throughput.

preprint2022arXiv

Turbulence in Milky Way Star-Forming Regions Traced by Young Stars and Gas

The interstellar medium (ISM) is turbulent on all scales and in all phases. In this paper, we study turbulence with different tracers in four nearby star-forming regions: Orion, Ophiuchus, Perseus, and Taurus. We combine the APOGEE-2 and Gaia surveys to obtain the full 6-dimensional measurements of positions and velocities of young stars in these regions. The velocity structure functions (VSFs) of the stars show a universal scaling of turbulence. We also obtain Hα gas kinematics in these four regions from the Wisconsin H-Alpha Mapper. The VSFs of the Hα are more diverse compared to the stars. In regions with recent supernova activities, they show characteristics of local energy injections and higher amplitudes compared to the VSFs of stars and of CO from the literature. Such difference in amplitude of the VSFs can be explained by the different energy and momentum transport from supernovae into different phases of the ISM, thus resulting in higher levels of turbulence in the warm ionized phase traced by Hα. In regions without recent supernova activities, the VSFs of young stars, Hα, and CO are generally consistent, indicating well-coupled turbulence between different phases. Within individual regions, the brighter parts of the Hα gas tend to have a higher level of turbulence than the low-emission parts. Our findings support a complex picture of the Milky Way ISM, where turbulence can be driven at different scales and inject energy unevenly into different phases.

preprint2021arXiv

Certificate complexity and symmetry of nested canalizing functions

Boolean nested canalizing functions (NCFs) have important applications in molecular regulatory networks, engineering and computer science. In this paper, we study their certificate complexity. For both Boolean values $b\in\{0,1\}$, we obtain a formula for $b$-certificate complexity and consequently, we develop a direct proof of the certificate complexity formula of an NCF. Symmetry is another interesting property of Boolean functions and we significantly simplify the proofs of some recent theorems about partial symmetry of NCFs. We also describe the algebraic normal form of $s$-symmetric NCFs. We obtain the general formula of the cardinality of the set of $n$-variable $s$-symmetric Boolean NCFs for $s=1,\dots,n$. In particular, we enumerate the strongly asymmetric Boolean NCFs.

preprint2021arXiv

Extremal solution and Liouville theorem for anisotropic elliptic equations

We study the quasilinear Dirichlet boundary problem \begin{equation}\nonumber \left\{ \begin{aligned} -Qu&=λe^{u} \quad \mbox{in}\quadΩ\\ u&=0 \quad \mbox{on}\quad\partialΩ,\\ \end{aligned} \right. \end{equation} where $λ>0$ is a parameter, $Ω\subset\mathbb{R}^{N}$ with $N\geq2$ be a bounded domain, and the operator $Q$, known as Finsler-Laplacian or anisotropic Laplacian, is defined by $$Qu:=\sum_{i=1}^{N}\frac{\partial}{\partial x_{i}}(F(\nabla u)F_{ξ_{i}}(\nabla u)). $$ Here, $F_{ξ_{i}}=\frac{\partial F}{\partialξ_{i}}$ and $F: \mathbb{R}^{N}\rightarrow[0,+\infty)$ is a convex function of $ C^{2}(\mathbb{R}^{N}\setminus\{0\})$, that satisfies certain assumptions. We derive the existence of extremal solution and obtain that it&#39;s regular, if $N\leq9$. We also concern the Hénon type anisotropic Liouville equation, namely, $$-Qu=(F^{0}(x))^αe^{u}\quad\mbox{in}\quad\mathbb{R}^{N}$$ where $α>-2$, $N\geq2$ and $F^{0}$ is the support function of $K:=\{x\in\mathbb{R}^{N}:F(x)<1\}$ which is defined by $$F^{0}(x):=\sup_{ξ\in K}\langle x,ξ\rangle.$$ We obtain the Liouville theorem for stable solutions and the finite Morse index solutions for $2\leq N<10+4α$ and $3\leq N<10+4α^{-}$ respectively, where $α^{-}=\min\{α,0\}$.

preprint2021arXiv

Geometrical and electrical modulation on the transport property of silicene constrictions

We study the electrical modulation of the transport properties of silicene constrictions with different geometrical structures by adopting the tight-binding model and non-equilibrium Green&#39;s function method. The band structure and transmission properties are discussed under the influence of the external electric field and potential energy. Especially, we investigate the effects of the position and width of the central scattering region on the conductance with increasing of Fermi energy. We find that the conductance significantly depends on the position and the width. Interestingly, the symmetrical structure of the central region can induce a resonance effect and significantly enlarge the system&#39;s conductance. Obviously, we obtain an effective method to adjust the transport property of the silicene heterojunctions. Correspondingly, we propose a novel two-channel structure with an excellent performance on the conductance compared to the one-channel structure with the same total width.

preprint2021arXiv

Giant anisotropic photocurrent modulated by strain in type-II Weyl semimetal Td-MoTe2

We build a Cu-MoTe2-Cu device model and use first-principles density functional theory to study the transport properties of single-layer Td-MoTe2. We obtained the effect of strain on the energy band structure, transport properties, and photocurrent. The strain-induced photocurrent shows an anisotropy that reflects the modulation of the energy bands, including the Weyl point, by strain. The photocurrent can be suppressed to almost zero when the strain is applied along the vacuum direction. In contrast, the photocurrent can be significantly increased when the strain is applied along the transport direction. The transport properties and magnitude of the photocurrent in the MoTe2-based device can be effectively modulated by adjusting the strength and direction of the strain.

preprint2021arXiv

Linear Stability of an Impulsively Accelerated Density Interface in an Ideal Two-Fluid Plasma

We investigate the linear evolution of Richtmyer-Meshkov (RM) instability in the framework of an ideal two-fluid plasma model. The two-fluid plasma equations of motion are separated into a base state and a set of linearized equations governing the evolution of the perturbations. Different coupling regimes between the charged species are distinguished based on a non-dimensional Debye length parameter $d_{D,0}$. When $d_{D,0}$ is large, the coupling between ions and electrons is sufficiently small that the induced Lorentz force is very weak and the two species evolve as two separate fluids. When $d_{D,0}$ is small, the coupling is strong and the induced Lorentz force is strong enough that the difference between state of ions and electrons is rapidly decreased by the force. As a consequence, the ions and electrons are tightly coupled and evolve like one fluid. The temporal dynamics is divided into two phases: an early phase wherein electron precursor waves are prevalent, and a post ion shock-interface interaction phase during which the RM instability manifests itself. We also examine the effect of an initially applied magnetic field in the streamwise direction characterized by the non-dimensional parameter $β_0$. For a short duration after the ion shock-interface interaction, the growth rate is similar for different initial magnetic field strengths. As time progresses the suppression of the instability due to the magnetic field is observed. The growth rate shows oscillations with a frequency that is related to the ion or electron cyclotron frequency. The instability is suppressed due to the vorticity being transported away from the interface.

preprint2021arXiv

Measuring Turbulence with Young Stars in the Orion Complex

Stars form in molecular clouds in the interstellar medium (ISM) with a turbulent kinematic state. Newborn stars therefore should retain the turbulent kinematics of their natal clouds. Gaia DR2 and APOGEE-2 surveys in combination provide three-dimensional (3D) positions and 3D velocities of young stars in the Orion Molecular Cloud Complex. Using the full 6D measurements, we compute the velocity structure functions (VSFs) of the stars in six different groups within the Orion Complex. We find that the motions of stars in all diffuse groups exhibit strong characteristics of turbulence. Their first-order VSFs have a power-law exponent ranging from $\sim0.2-0.5$ on scales of a few to a few tens of pc, generally consistent with Larson&#39;s relation. On the other hand, dense star clusters, such as the Orion Nebula Cluster (ONC), have experienced rapid dynamical relaxation, and have lost the memory of the initial turbulent kinematics. The VSFs of several individual groups and the whole Complex all show features supporting local energy injection from supernovae. The measured strength of turbulence depends on the location relative to the supernova epicenters and the formation history of the groups. Our detection of turbulence traced by young stars introduces a new method of probing the turbulent kinematics of the ISM. Unlike previous gas-based studies with only projected measurements accessible to observations, we utilize the full 6D information of stars, presenting a more complete picture of the 3D interstellar turbulence.

preprint2021arXiv

Unusual heat transport of the Kitaev material Na$_2$Co$_2$TeO$_6$: putative quantum spin liquid and low-energy spin excitations

We studied the field dependent thermal conductivity ($κ$) of Na$_2$Co$_2$TeO$_6$, a compound considered as the manifestation of the Kitaev model based on the high-spin $d^7$ Co$^{2+}$ ions. We found that in-plane magnetic fields beyond a critical value $B_c \approx$~10 T are able to drastically enhance $κ$ at low temperatures, resulting in a double-peak structure of $κ(T)$ that closely resembles the behavior of $α$-RuCl$_3$. This result suggests that heat transport in Na$_2$Co$_2$TeO$_6$ is primarily phononic, and it is strongly affected by scattering from magnetic excitations that are highly tunable by external fields. Interestingly, for magnetic fields $B // a$ (i.e., along the zigzag direction of the Co-Co bonds), there is an extended field range which separates the long-range magnetic order for $B\leq B_c\approx10$ T and the partially spin-polarized gapped high-field phase for $B\gtrsim 12$ T. The low-energy phonon scattering is particularly strong in this field range, consistent with the notion that the system becomes a quantum spin liquid with prominent spin fluctuations down to energies of no more than 2 meV.

preprint2021arXiv

Valley-dependent transport in strain engineering graphene heterojunctions

We study the effect of the strain on the band structure and the valley-dependent transport property of graphene heterojunctions. It is found that valley-dependent separation of electrons can be achieved by utilizing the strain and on-site energies. In the presence of the strain, the values of the transmission can be effectively adjusted by changing the strengths of the strain, while the transport angle basically keeps unchanged. When an extra on-site energy is simultaneously applied to the central scattering region, not only are the electrons of valleys K and K&#39; separated into two distinct transmission lobes in opposite transverse directions, but the transport angles of two valleys can be significantly changed. Therefore, one can realize an effective modulation of valley-dependent transport by changing the strength and stretch angle of the strain and on-site energies, which can be exploited for graphene-based valleytronics devices.

preprint2020arXiv

AFDet: Anchor Free One Stage 3D Object Detection

High-efficiency point cloud 3D object detection operated on embedded systems is important for many robotics applications including autonomous driving. Most previous works try to solve it using anchor-based detection methods which come with two drawbacks: post-processing is relatively complex and computationally expensive; tuning anchor parameters is tricky. We are the first to address these drawbacks with an anchor free and Non-Maximum Suppression free one stage detector called AFDet. The entire AFDet can be processed efficiently on a CNN accelerator or a GPU with the simplified post-processing. Without bells and whistles, our proposed AFDet performs competitively with other one stage anchor-based methods on KITTI validation set and Waymo Open Dataset validation set.

preprint2020arXiv

Anomalous spectral weight transfer in the nematic state of iron-selenide superconductor

Nematic phase intertwines closely with high-Tc superconductivity in iron-based superconductors. Its mechanism, which is closely related to the pairing mechanism of superconductivity, still remains controversial. Comprehensive characterization of how the electronic state reconstructs in the nematic phase is thus crucial. However, most experiments focus only on the reconstruction of band dispersions. Another important characteristic of electronic state, the spectral weight, has not been studied in details so far. Here, we studied the spectral weight transfer in the nematic phase of FeSe$_{0.9}$S$_{0.1}$ using angle-resolved photoemission spectroscopy and in-situ detwinning technique. There are two elliptical electron pockets overlapping with each other orthogonally at the Brillouin zone corner. We found that, upon cooling, one electron pocket loses spectral weight and fades away, while the other electron pocket gains spectral weight and becomes pronounced. Our results show that the symmetry breaking of electronic state is manifested by not only the anisotropic band dispersion but also the band-selective modulation of spectral weight. Our observation completes our understanding of the nematic electronic state, and put strong constraints on the theoretical models. It further provide crucial clues to understand the gap anisotropy and orbital-selective pairing in iron-selenide superconductors.

preprint2020arXiv

AutoPose: Searching Multi-Scale Branch Aggregation for Pose Estimation

We present AutoPose, a novel neural architecture search(NAS) framework that is capable of automatically discovering multiple parallel branches of cross-scale connections towards accurate and high-resolution 2D human pose estimation. Recently, high-performance hand-crafted convolutional networks for pose estimation show growing demands on multi-scale fusion and high-resolution representations. However, current NAS works exhibit limited flexibility on scale searching, they dominantly adopt simplified search spaces of single-branch architectures. Such simplification limits the fusion of information at different scales and fails to maintain high-resolution representations. The presentedAutoPose framework is able to search for multi-branch scales and network depth, in addition to the cell-level microstructure. Motivated by the search space, a novel bi-level optimization method is presented, where the network-level architecture is searched via reinforcement learning, and the cell-level search is conducted by the gradient-based method. Within 2.5 GPU days, AutoPose is able to find very competitive architectures on the MS COCO dataset, that are also transferable to the MPII dataset. Our code is available at https://github.com/VITA-Group/AutoPose.

preprint2020arXiv

Band insulator to Mott insulator transition in 1T-TaS$_2$

1T-TaS$_2$ undergoes successive phase transitions upon cooling and eventually enters an insulating state of mysterious origin. Some consider this state to be a band insulator with interlayer stacking order, yet others attribute it to Mott physics that support a quantum spin liquid state.Here, we determine the electronic and structural properties of 1T-TaS$_2$ using angle-resolved photoemission spectroscopy and X-Ray diffraction. At low temperatures, the 2$π$/2c-periodic band dispersion, along with half-integer-indexed diffraction peaks along the c axis, unambiguously indicates that the ground state of 1T-TaS$_2$ is a band insulator with interlayer dimerization. Upon heating, however, the system undergoes a transition into a Mott insulating state, which only exists in a narrow temperature window. Our results refute the idea of searching for quantum magnetism in 1T-TaS$_2$ only at low temperatures, and highlight the competition between on-site Coulomb repulsion and interlayer hopping as a crucial aspect for understanding the material&#39;s electronic properties.

preprint2020arXiv

Densely Connected Search Space for More Flexible Neural Architecture Search

Neural architecture search (NAS) has dramatically advanced the development of neural network design. We revisit the search space design in most previous NAS methods and find the number and widths of blocks are set manually. However, block counts and block widths determine the network scale (depth and width) and make a great influence on both the accuracy and the model cost (FLOPs/latency). In this paper, we propose to search block counts and block widths by designing a densely connected search space, i.e., DenseNAS. The new search space is represented as a dense super network, which is built upon our designed routing blocks. In the super network, routing blocks are densely connected and we search for the best path between them to derive the final architecture. We further propose a chained cost estimation algorithm to approximate the model cost during the search. Both the accuracy and model cost are optimized in DenseNAS. For experiments on the MobileNetV2-based search space, DenseNAS achieves 75.3% top-1 accuracy on ImageNet with only 361MB FLOPs and 17.9ms latency on a single TITAN-XP. The larger model searched by DenseNAS achieves 76.1% accuracy with only 479M FLOPs. DenseNAS further promotes the ImageNet classification accuracies of ResNet-18, -34 and -50-B by 1.5%, 0.5% and 0.3% with 200M, 600M and 680M FLOPs reduction respectively. The related code is available at https://github.com/JaminFong/DenseNAS.

preprint2020arXiv

Fast Neural Network Adaptation via Parameter Remapping and Architecture Search

Deep neural networks achieve remarkable performance in many computer vision tasks. Most state-of-the-art (SOTA) semantic segmentation and object detection approaches reuse neural network architectures designed for image classification as the backbone, commonly pre-trained on ImageNet. However, performance gains can be achieved by designing network architectures specifically for detection and segmentation, as shown by recent neural architecture search (NAS) research for detection and segmentation. One major challenge though, is that ImageNet pre-training of the search space representation (a.k.a. super network) or the searched networks incurs huge computational cost. In this paper, we propose a Fast Neural Network Adaptation (FNA) method, which can adapt both the architecture and parameters of a seed network (e.g. a high performing manually designed backbone) to become a network with different depth, width, or kernels via a Parameter Remapping technique, making it possible to utilize NAS for detection/segmentation tasks a lot more efficiently. In our experiments, we conduct FNA on MobileNetV2 to obtain new networks for both segmentation and detection that clearly out-perform existing networks designed both manually and by NAS. The total computation cost of FNA is significantly less than SOTA segmentation/detection NAS approaches: 1737$\times$ less than DPC, 6.8$\times$ less than Auto-DeepLab and 7.4$\times$ less than DetNAS. The code is available at https://github.com/JaminFong/FNA.

preprint2020arXiv

FasterSeg: Searching for Faster Real-time Semantic Segmentation

We present FasterSeg, an automatically designed semantic segmentation network with not only state-of-the-art performance but also faster speed than current methods. Utilizing neural architecture search (NAS), FasterSeg is discovered from a novel and broader search space integrating multi-resolution branches, that has been recently found to be vital in manually designed segmentation models. To better calibrate the balance between the goals of high accuracy and low latency, we propose a decoupled and fine-grained latency regularization, that effectively overcomes our observed phenomenons that the searched networks are prone to &#34;collapsing&#34; to low-latency yet poor-accuracy models. Moreover, we seamlessly extend FasterSeg to a new collaborative search (co-searching) framework, simultaneously searching for a teacher and a student network in the same single run. The teacher-student distillation further boosts the student model&#39;s accuracy. Experiments on popular segmentation benchmarks demonstrate the competency of FasterSeg. For example, FasterSeg can run over 30% faster than the closest manually designed competitor on Cityscapes, while maintaining comparable accuracy.

preprint2020arXiv

Field sparsening for the construction of the correlation functions in lattice QCD

Two field-sparsening methods, namely the sparse-grid method and the random field selection method, are used in this paper for the construction of the 2-point and 3-point correlation functions in lattice QCD. We argue that, due to the high correlation among the lattice correlators at different field points associated with source, current, and sink locations, one can save a lot of computational time by performing the summation over a subset of the lattice sites. Furthermore, with this strategy, one only needs to store a small fraction of the full quark propagators. It is found that the number of field points can be reduced by a factor of $\sim$100 for the point-source operator and a factor of $\sim$1000 for the Gaussian-smeared operator, while the uncertainties of the correlators only increase by $\sim$15\%. Therefore, with a modest cost of the computational resources, one can approach the precision of the all-to-all correlators using the field-sparsening methods.

preprint2020arXiv

Giant pressure-enhancement of multiferroicity in CuBr2

Type-II multiferroic materials, in which ferroelectric polarization is induced by inversion non-symmetric magnetic order, promise new and highly efficient multifunctional applications based on the mutual control of magnetic and electric properties. Although this phenomenon has to date been limited to low temperatures, here we report a giant pressure-dependence of the multiferroic critical temperature in CuBr$_2$. At 4.5 GPa, $T_\mathrm{C}$ is enhanced from 73.5 to 162 K, to our knowledge the highest value yet reported for a non-oxide type-II multiferroic. This growth shows no sign of saturating and the dielectric loss remains small under these high pressures. We establish the structure under pressure and demonstrate a 60\% increase in the two-magnon Raman energy scale up to 3.6 GPa. First-principles structural and magnetic energy calculations provide a quantitative explanation in terms of dramatically pressure-enhanced interactions between CuBr$_2$ chains. These large, pressure-tuned magnetic interactions motivate structural control in cuprous halides as a route to applied high-temperature multiferroicity.

preprint2020arXiv

Momentum-resolved measurement of electronic nematic susceptibility in the FeSe$_{0.9}$S$_{0.1}$ superconductor

Unveiling the driving force for a phase transition is normally difficult when multiple degrees of freedom are strongly coupled. One example is the nematic phase transition in iron-based superconductors. Its mechanism remains controversial due to a complex intertwining among different degrees of freedom. In this paper, we report a method for measuring the nematic susceptibly of FeSe$_{0.9}$S$_{0.1}$ using angle-resolved photoemission spectroscopy (ARPES) and an $in$-$situ$ strain-tuning device. The nematic susceptibility is characterized as an energy shift of band induced by a tunable uniaxial strain. We found that the temperature-dependence of the nematic susceptibility is strongly momentum dependent. As the temperature approaches the nematic transition temperature from the high temperature side, the nematic susceptibility remains weak at the Brillouin zone center while showing divergent behavior at the Brillouin zone corner. Our results highlight the complexity of the nematic order parameter in the momentum space, which provides crucial clues to the driving mechanism of the nematic phase transition. Our experimental method which can directly probe the electronic susceptibly in the momentum space provides a new way to study the complex phase transitions in various materials.

preprint2020arXiv

Nondispersive solutions to the mass critical half-wave equation in two dimensions

We consider the half-wave equation with mass critical in two dimension \begin{eqnarray*} \begin{cases} iu_t=Du-|u|u,\,\,\, \\ u(0,x)=u_0(x), \end{cases} \end{eqnarray*} First, we prove the existence of a family of traveling solitary waves. We then show the existence of finite-time blowup solutions with minimal mass $\|u_0\|_2=\|Q\|_2$, where $Q$ is the ground state solution of equation $DQ+Q=Q^2$.

preprint2020arXiv

Paramagnons and high-temperature superconductivity in mercury-based cuprates

We present a comparative study of magnetic excitations in the first two Ruddlesden-Popper members of the Hg-family of high-temperature superconducting cuprates, which are chemically nearly identical and have the highest critical temperature ($T_\mathrm{c}$) among all cuprate families. Our inelastic photon scattering experiments reveal that the two compounds&#39; paramagnon spectra are nearly identical apart from an energy scale factor of $\sim130\%$ that matches the ratio of $T_\mathrm{c}$&#39;s, as expected in magnetic Cooper pairing theories. By relating our observations to other cuprates, we infer that the strength of magnetic interactions determines how high $T_\mathrm{c}$ can reach. Our finding can be viewed as a magnetic analogue of the isotope effect, thus firmly supporting models of magnetically mediated high-temperature superconductivity.

preprint2020arXiv

Partial regularity and Liouville theorems for stable solutions of anisotropic elliptic equations

We study the quasilinear elliptic equation \begin{equation*} -Qu=e^u \ \ \text{in} \ \ Ω\subset \mathbb{R}^{N} \end{equation*} where the operator $Q$, known as Finsler-Laplacian (or anisotropic Laplacian), is defined by $$Qu:=\sum_{i=1}^{N}\frac{\partial}{\partial x_{i}}(F(\nabla u)F_{ξ_{i}}(\nabla u)),$$ where $F_{ξ_{i}}=\frac{\partial F}{\partialξ_{i}}$ and $F: \mathbb{R}^{N}\rightarrow[0,+\infty)$ is a convex function of $ C^{2}(\mathbb{R}^{N}\setminus\{0\})$, that satisfies certain assumptions. For bounded domain $Ω$ and for a stable weak solution of the above equation, we prove that the Hausdorff dimension of singular set does not exceed $N-10$. For the entire space, we apply Moser iteration arguments, established by Dancer-Farina and Crandall-Rabinowitz in the context, to prove Liouville theorems for stable solutions and for finite Morse index solutions in dimensions $N<10$ and $2<N<10$, respectively. We also provide an explicit solution that is stable outside a compact set in $N=2$. In addition, we provide similar Liouville theorems for the power-type nonlinearities.

preprint2020arXiv

Parton distribution functions of $Δ^+$ on the lattice

We perform a first calculation for the unpolarized parton distribution function of the $Δ^+$ baryon using lattice QCD simulations within the framework of Large Momentum Effective Theory. Two ensembles of $N_f=2+1+1$ twisted mass fermions are utilized with a pion mass of 270 MeV and 360 MeV, respectively. The baryon, which is treated as a stable single-particle state, is boosted with momentum $P_3$ with values $\{0.42,0.83,1.25\}$ GeV, and we utilize momentum smearing to improve the signal. The unpolarized parton distribution function of $Δ^+$ is obtained using a non-perturbative renormalization and a one-loop formula for the matching, with encouraging precision. In particular, we compute the $\overline{d}(x)-\overline{u}(x)$ asymmetry and compare it with the same quantity in the nucleon, in a first attempt towards resolving the physical mechanism responsible for generating such asymmetry.

preprint2020arXiv

Properties of the Hot Ambient Medium of Early-type Galaxies Hosting Powerful Radio Sources

We present an archival analysis of Chandra X-ray observations for twelve nearby early-type galaxies hosting radio sources with radio power $>10^{23} \, \rm{W}~\rm{Hz}^{-1}$ at 1.4 GHz, similar to the radio power of the radio source in NGC 4261. Previously, in a similar analysis of eight nearby X-ray and optically-bright elliptical galaxies, Werner et al. 2012, found that NGC 4261 exhibited unusually low central gas entropy compared to the full sample. In the central 0.3 kpc of NGC 4261, the ratio of cooling time to freefall time ($t_{\rm{cool}}/t_{\rm ff}$) is less than $10$, indicating that cold clouds may be precipitating out of the hot ambient medium and providing fuel for accretion in the central region. NGC 4261 also hosts the most powerful radio source in the original sample. Because NGC 4261 may represent an important phase during which powerful feedback from a central active galactic nucleus (AGN) is fueled by multiphase condensation in the central kpc, we searched the Chandra archive for analogs to NGC 4261. We present entropy profiles of those galaxies as well as profiles of $t_{\rm{cool}}/t_{\rm ff}$. We find that one of them, IC 4296, exhibits properties similar to NGC 4261, including the presence of only single phase gas outside of $r \sim 2$ kpc and a similar central velocity dispersion. We compare the properties of NGC 4261 and IC 4296 to hydrodynamic simulations of AGN feedback fueled by precipitation. Over the course of those simulations, the single phase galaxy has an entropy gradient that remains similar to the entropy profiles inferred from our observations.

preprint2020arXiv

Properties of the Simulated Circumgalactic Medium

The circumgalactic medium (CGM) is closely linked to galaxy formation and evolution, but difficult to characterize observationally and typically poorly resolved in cosmological simulations. We use spherically-symmetric, idealized, high-resolution simulations of the CGM in $10^{12}M_\odot$ and $10^{11}M_\odot$ dark matter halos to characterize the gas pressure, turbulent and radial velocities, and degree of thermal and effective dynamic pressure support in the overall CGM as well as in its high- and low-temperature phases. We find that the $10^{12}M_\odot$ halo contains a CGM mostly formed of a hot gas halo in hydrostatic equilibrium out of which cold gas condenses and falls onto the central galaxy, while the $10^{11}M_\odot$ halo&#39;s CGM is not in hydrostatic equilibrium, has a wider spread of properties at a given galactocentric radius, does not have a clear separation of hot and cold phases, and is dominated by bulk motions. We also find that the degree of pressure support in the $10^{11}M_\odot$ halo is strongly dependent on the parameters of the galactic winds of the central galaxy. These results promote the idea that there is no &#34;average&#34; CGM and care must be taken when setting the initial conditions for a small-box simulation of a patch of the CGM.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: I. Formation of the Multiphase Interstellar medium

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: II. Energetics and Turbulence

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.

preprint2020arXiv

The Interplay of Kinetic and Radiative Feedback in Galaxy Clusters

Recent observations provide evidence that some cool-core clusters (CCCs) host quasars in their brightest cluster galaxies (BCGs). Motivated by these findings we use 3D radiation-hydrodynamic simulations with the code Enzo to explore the joint role of the kinetic and radiative feedback from supermassive black holes (SMBHs) in BCGs. We implement kinetic feedback as sub-relativistic plasma outflows and model radiative feedback using the ray-tracing radiative transfer or thermal energy injection. In our simulations the central SMBH transitions between the radiatively efficient and radiatively inefficient states on timescales of a few Gyr, as a function of its accretion rate. The timescale for this transition depends primarily on the fraction of power allocated to each feedback mode, and to a lesser degree on the overall feedback luminosity of the active galactic nucleus (AGN). Specifically, we find that (a) kinetic feedback must be present at both low and high accretion rates in order to prevent the cooling catastrophe, and (b) its contribution likely accounts for > 10% of the total AGN feedback power, since below this threshold simulated BCGs tend to host radio-loud quasars most of the time, in apparent contrast with observations. We also find a positive correlation between the AGN feedback power and the mass of the cold gas filaments in the cluster core, indicating that observations of H$α$ filaments can be used as a measure of AGN feedback.

preprint2019arXiv

Correlations Between Black Holes and Host Galaxies in the Illustris and IllustrisTNG Simulations

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-σ$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-σ$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-σ$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $σ>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

preprint2019arXiv

Direct Detection of Black Hole-Driven Turbulence in the Centers of Galaxy Clusters

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this &#34;feedback&#34; process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

preprint2019arXiv

Distinct fingerprints of charge density waves and electronic standing waves in ZrTe$_3$

Experimental signatures of charge density waves (CDW) in high-temperature superconductors have evoked much recent interest, yet an alternative interpretation has been theoretically raised based on electronic standing waves resulting from quasiparticles scattering off impurities or defects, also known as Friedel oscillations (FO). Indeed the two phenomena are similar and related, posing a challenge to their experimental differentiation. Here we report a resonant X-ray diffraction study of ZrTe$_3$, a model CDW material. Near the CDW transition, we observe two independent diffraction signatures that arise concomitantly, only to become clearly separated in momentum while developing very different correlation lengths in the well-ordered state. Anomalously slow dynamics of mesoscopic ordered nanoregions are further found near the transition temperature, in spite of the expected strong thermal fluctuations. These observations reveal that a spatially-modulated CDW phase emerges out of a uniform electronic fluid via a process that is promoted by self-amplifying FO, and identify a viable experimental route to distinguish CDW and FO.

preprint2019arXiv

Doping-dependent phonon anomaly and charge-order phenomena in the HgBa$_{2}$CuO$_{4+δ}$ and HgBa$_{2}$CaCu$_{2}$O$_{6+δ}$

Using resonant X-ray diffraction and Raman spectroscopy, we study charge correlations and lattice dynamics in two model cuprates, HgBa$_{2}$CuO$_{4+δ}$ and HgBa$_{2}$CaCu$_{2}$O$_{6+δ}$. We observe a maximum of the characteristic charge order temperature around the same hole concentration ($p \approx 0.09$) in both compounds, and concomitant pronounced anomalies in the lattice dynamics that involve the motion of atoms in and/or adjacent to the CuO$_2$ layers. These anomalies are already present at room temperature, and therefore precede the formation of the static charge correlations, and we attribute them to an instability of the CuO$_2$ layers. Our finding implies that the charge order in the cuprates is an emergent phenomenon, driven by a fundamental variation in both lattice and electronic properties as a function of doping.

preprint2019arXiv

Ferrimagnetism and anisotropic phase tunability by magnetic fields in Na$_2$Co$_2$TeO$_6$

Na$_2$Co$_2$TeO$_6$ has recently been proposed to be a Kitaev-like honeycomb magnet. To assess how close it is to realizing Kitaev quantum spin liquids, we have measured magnetization and specific heat on high-quality single crystals in magnetic fields applied along high-symmetry directions. Small training fields reveal a weak but canonical ferrimagnetic behavior below 27 K, which cannot be explained by the zigzag antiferromagnetic order alone and suggests coexisting Néel-type order of moments canted away from the zigzag chains. Moderate fields in the honeycomb plane suppress the thermal transition at 27 K, and seem to partly reverse the moment-canting when applied perpendicular to the zigzag chains. In contrast, out-of-plane fields leave the transition largely unaffected, but promotes another transition below 10 K, possibly also related to canting reversal. The magnetism in Na$_2$Co$_2$TeO$_6$ is highly anisotropic and close to tipping points between competing phases.

preprint2019arXiv

The Fate of AGB Wind in Massive Galaxies and the ICM

Asymptotic Giant Branch (AGB) winds from evolved stars not only provide a non-trivial amount of mass and energy return, but also produce dust grains in massive elliptical galaxies. Due to the fast stellar velocity and the high ambient temperature, the wind is thought to form a comet-like tail, similar to Mira in the Local Bubble. Many massive elliptical galaxies and cluster central galaxies host extended dusty cold filaments. The fate of the cold dusty stellar wind and its relation to cold filaments are not well understood. In this work, we carry out both analytical and numerical studies of the interaction between an AGB wind and the surrounding hot gas. We find that the cooling time of the tail is inversely proportional to the ambient pressure. In the absence of cooling, or in low pressure environments (e.g., the outskirts of elliptical galaxies), AGB winds are quickly mixed into the hot gas, and all the AGB winds have similar appearance and head-to-tail ratio. In high pressure environments, such as the Local Bubble and the central regions of massive elliptical galaxies, some of the gas in the mixing layer between the stellar wind and the surrounding hot gas can cool efficiently and cause the tail to become longer. Our simulated tail of Mira itself has similar length and velocity to that observed, and appears similar to the simulated AGB tail in the central regions of massive galaxies. We speculate that instead of thermal instability, the induced condensation at the mixing layer of AGB winds may be the origin of cold filaments in massive galaxies and galaxy clusters. This naturally explains the existence of dust and PAH in the filaments.

preprint2019arXiv

The relationship between black hole mass and galaxy properties: Examining the black hole feedback model in IllustrisTNG

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$.