Search

We present the first complete high-precision QCD corrections to the inclusive decay width $\mathrm{\Gamma}_t$, the $W$-helicity fractions $f_{\mathrm{L,R,0}}$ and semi-inclusive distributions for the top-quark decay process $t \rightarrow b + W^+ + X_{\mathrm{\tiny QCD}}$ at NNNLO in the strong coupling constant $\alpha_s$. In particular, the pure NNNLO QCD correction decreases the $\mathrm{\Gamma}_t$ by about $0.8\%$ of the previous NNLO result at the top-quark pole mass scale, exceeding the error estimated by the usual scale-variation prescription. After taking into account all sources of errors, we get $\mathrm{\Gamma}_t = 1.3148^{+0.003}_{-0.005} + 0.027\,(m_t - 172.69)\,\text{GeV} $, the error of which meets the request by future colliders. On the other hand, the NNNLO QCD effects on $f_{\mathrm{L,R,0}}$ are found to be much smaller, at the level of one per-mille for the dominating $f_{0}$, predestining them to act as precision observables for the top-quark decay process., Comment: 7 pages, 3 figures

We present the results of an unbiased $^{12}$CO/$^{13}$CO/C$^{18}$O ($J$ = 1-0) survey in a portion of the third Galactic quadrant (TGQ): $l$ = [219.75, 229.75]$^\circ$ and $b$ = [-5.25, 5.25]$^\circ$. The high-resolution and high-sensitivity data sets help to unravel the distributions and physical properties of the molecular clouds (MCs) in the mapped area. In the LSR velocity range from -1 to 85 km/s, the molecular material successfully traces the Local, Perseus, and Outer arms. In the TGQ, the Outer arm appears to be more prominent than that in the second Galactic quadrant (SGQ), but the Perseus arm is not as conspicuous as that in the SGQ. A total of 1,502 $^{12}$CO, 570 $^{13}$CO, and 53 C$^{18}$O molecular structures are identified, spanning over $\sim2$ and $\sim6$ orders of magnitude in size and mass, respectively. Tight mass-radius correlations and virial parameter-mass anticorrelations are observable. Yet, it seems that no clear correlations between velocity dispersion and effective radius can be found over the full dynamic range. The vertical distribution of the MCs renders evident pictures of the Galactic warp and flare., Comment: 22 pages, 13 figures, 7 tables (with machine-readable versions), published in ApJS

We present the successful realization of four-dimensional (4D) semimetal bands featuring tensor monopoles, achieved using superconducting quantum circuits. Our experiment involves the creation of a highly tunable diamond energy diagram with four coupled transmons, and the parametric modulation of their tunable couplers, effectively mapping momentum space to parameter space. This approach enables us to establish a 4D Dirac-like Hamiltonian with fourfold degenerate points. Moreover, we manipulate the energy of tensor monopoles by introducing an additional pump microwave field, generating effective magnetic and pseudo-electric fields and simulating topological parity magnetic effects emerging from the parity anomaly. Utilizing non-adiabatic response methods, we measure the fractional second Chern number for a Dirac valley with a varying mass term, signifying a nontrivial topological phase transition connected to a 5D Yang monopole. Our work lays the foundation for further investigations into higher-dimensional topological states of matter and enriches our comprehension of topological phenomena.

We compare the observational properties between $^{12}$CO, $^{13}$CO, and C$^{18}$O and summarize the observational parameters based on 7069 clouds sample from the Milky Way Imaging Scroll Painting (MWISP) CO survey in a section of the third Galactic quadrant. We find that the $^{13}$CO angular area ($A_{\rm ^{13}CO}$) generally increases with that of $^{12}$CO ($A_{\rm ^{12}CO}$), and the ratio of $A_{\rm ^{13}CO}$ to $A_{\rm ^{12}CO}$ is 0.38 by linear fitting. We find that the $^{12}$CO and $^{13}$CO flux are tightly correlated as $F_{\rm ^{13}CO}~=~0.17~ F_{\rm ^{12}CO}$ with both fluxes calculated within the $^{13}$CO-bright region. This indicates that the abundance $X_{\rm ^{13}CO}$ is a constant to be 6.5$^{+0.1}_{-0.5}$ $\times 10^{-7}$ for all samples under assumption of local thermodynamic equilibrium (LTE). Additionally, we observed that the X-factor is approximately constant in large sample molecular clouds. Similarly, we find $F_{\rm C^{18}O}~=~0.11~F_{\rm ^{13}CO}$ with both fluxes calculated within C$^{18}$O-bright region, which indicates that the abundance ratios ${X_{\rm ^{13}CO}/X_{\rm C^{18}O}}$ stays the same value 9.7$^{+0.6}_{-0.8}$ across the molecular clouds under LTE assumption. The linear relationships of $F_{\rm ^{12}CO}$ vs. $F_{\rm ^{13}CO}$ and $F_{\rm ^{13}CO}$ vs. $F_{\rm C^{18}O}$ hold not only for the $^{13}$CO-bright region or C$^{18}$O-bright region, but also for the entire molecular cloud scale with lower flux ratio. The abundance ratio ${X_{\rm ^{13}CO}/X_{\rm C^{18}O}}$ inside clouds shows a strong correlation with column density and temperature. This indicates that the ${X_{\rm ^{13}CO}/X_{\rm C^{18}O}}$ is dominated by a combination of chemical fractionation, selectively dissociation, and self-shielding effect inside clouds., Comment: 11 pages, 16 figures, 1 table, accepted by AJ

Quantum computing offers potential solutions for finding ground states in condensed-matter physics and chemistry. However, achieving effective ground state preparation is also computationally hard for arbitrary Hamiltonians. It is necessary to propose certain assumptions to make this problem efficiently solvable, including preparing a trial state of a non-trivial overlap with the genuine ground state. Here, we propose a classical-assisted quantum ground state preparation method for quantum many-body systems, combining Tensor Network States (TNS) and Monte Carlo (MC) sampling as a heuristic method to prepare a trial state with a non-trivial overlap with the genuine ground state. We extract a sparse trial state by sampling from TNS, which can be efficiently prepared by a quantum algorithm on early fault-tolerant quantum computers. Our method demonstrates a polynomial improvement in scaling of overlap between the trial state and genuine ground state compared to random trial states, as evidenced by numerical tests on the spin-$1/2$ $J_1$-$J_2$ Heisenberg model. Furthermore, our method is a novel approach to hybridize a classical numerical method and a quantum algorithm and brings inspiration to ground state preparation in other fields.

In quantum logic operations, information is carried by the wavefunction rather than the energy distribution. Therefore, the relative phase is essential. Abelian and non-Abelian phases can be emulated in classical waves using passive coupled waveguides with geometric modulation. However, the dynamic phases interference induced by waveguide structure variation is inevitable.To overcome the challenges, we introduce an electroacoustic coupled system that enables the precise control of phase distribution through dynamic modulation of hopping. Such effective hopping is electronically controlled and is utilized to construct various paths in parameter space. These paths lead to state evolution with matrix-valued geometric phases, which correspond to logic operations. We report experimental realizations of several logic operations, including $Y$ gate, $Z$ gate, Hadamard gate and non-Abelian braiding. Our work introduces a temporal process to manipulate transient modes in a compact structure, providing a versatile experimental testbed for exploring other logic gates and exotic topological phenomena.

In this paper, we presented the 23.3 years of pulsar timing results of PSR J1456-6413 based on the observation of Parkes 64m radio telescope. We detected two new glitches at MJD 57093(3) and 59060(12) and confirmed its first glitch at MJD 54554(10). Using the "Cholesky" timing analysis method, we have determined its position, proper motion, and two-dimensional transverse velocities from the data segments before and after the second glitch, respectively. Furthermore, we detected exponential recovery behavior after the first glitch, with a recovery time scale of approximately 200 days and a corresponding exponential recovery factor Q of approximately 0.15(2), while no exponential recovery was detected for the other two glitches. More interestingly, we found that the leading component of the integral pulse profile after the second glitch became stronger, while the main component became weaker. Our results will expand the sample of pulsars with magnetosphere fluctuation triggered by the glitch event.

We propose a universal spin superconducting diode effect (SDE) induced by spin-orbit coupling (SOC), where the critical spin supercurrents in opposite directions are unequal. By analysis from both the Ginzburg-Landau theory and energy band analysis, we show that the spin-$\uparrow \uparrow$ and spin-$\downarrow \downarrow$ Cooper pairs possess opposite phase gradients and opposite momenta from the SOC, which leads to the spin SDE. Two superconductors with SOC, a $p$-wave superconductor as a toy model and a practical superconducting nanowire, are numerically studied and they both exhibit spin SDE. In addition, our theory also provides a unified picture for both spin and charge SDEs. Besides, we propose spin-polarized detection and nonlocal spin transport, as mature experimental technologies, to confirm the spin SDE in superconducting nanowires., Comment: 6 pages, 4 figures (+ supplementary materials 4 pages, 1 figure)

We study the dissociation effect of $J/\Psi$ in magnetized, rotating QGP matter at finite temperature and chemical potential using gauge/gravity duality. By incorporating angular velocity into the holographic magnetic catalysis model, we analyze the influence of temperature, chemical potential, magnetic field, and angular velocity on the properties of $J/\Psi$ meson. The results reveal that temperature, chemical potential, and rotation enhance the dissociation effect and increase the effective mass in the QGP phase. However, the magnetic field suppresses dissociation, and its effect on the effective mass is non-trivial. Additionally, we explore the interplay between magnetic field and rotation, identifying a critical angular velocity that determines the dominant effect. As a parallel study, we also examine the rotation effect in the holographic inverse magnetic catalysis model, although the magnetic field exhibits distinctly different behaviors in these two models, the impact of rotation on the dissociation effect of $J/\Psi$ is similar. Finally, we investigate the influence of electric field and demonstrate that it also speeds up the $J/\Psi$ dissociation., Comment: 24 pages, 13 figures

Deriving a comprehensive set of reduction rules for Feynman integrals has been a longstanding challenge. In this paper, we present a proposed solution to this problem utilizing generating functions of Feynman integrals. By establishing and solving differential equations of these generating functions, we are able to derive a system of reduction rules that effectively reduce any associated Feynman integrals to their bases. We illustrate this method through various examples and observe its potential value in numerous scenarios., Comment: 11 pages, 4 figures, references added

To better understand the effects of strong coupling and QCD at high temperature in QGP, by using holographic model, we investigate the dissociation effect of bottomonium under the higher-order curvature corrections to the supergravity action corresponding to the corrections of large N expansion of boundary CFT in the side of field theory. The results show that effective potential is not a good physical quantity to estimate the dissociation strength of bottomonium in the case of finite wave number and considering the higher-order curvature corrections. Therefore, we calculate the quasinormal spectra(QNMs) and the differential configuration entropy(DCE). It is found that the dissociation effect is stronger for the stronger coupling., Comment: 8 pages, 8 figures

It was found that, using nonrelativistic QCD factorization, the predicted $\chi_{cJ}$ hadroproduction cross section at large $p_T$ can be negative. The negative cross sections originate from terms proportional to plus function in ${^{3}\hspace{-0.6mm}P_{J}^{[1]}}$ channels, which are remnants of the infrared subtraction in matching the ${^{3}\hspace{-0.6mm}P_{J}^{[1]}}$ short-distance coefficients. In this article, we find that the above terms can be factorized into the nonperturbative ${^{3}\hspace{-0.6mm}S_{1}^{[8]}}$ soft gluon distribution function in the soft gluon factorization (SGF) framework. Therefore, the problem can be naturally resolved in SGF. With an appropriate choice of nonperturbative parameters, the SGF can indeed give positive predictions for $\chi_{cJ}$ production rates within the whole $p_T$ region. The production of $\psi(2S)$ is also discussed, and there is no negative cross section problem., Comment: 15 pages, 8 figures

${\gamma}$-GeSe is a new type of layered bulk material that was recently successfully synthesized. By means of density functional theory first-principles calculations, we systematically studied the physical properties of two-dimensional (2D) few-layer ${\gamma}$-GeSe. It is found that few-layer ${\gamma}$-GeSe are semiconductors with band gaps decreasing with increasing layer number; and 2D ${\gamma}$-GeSe with layer number $n \geq 2$ are ferroelectric with rather low transition barriers, consistent with the sliding ferroelectric mechanism. Particularly, spin-orbit coupling induced spin splitting is observed at the top of the valence band, which can be switched by the ferroelectric reversal; furthermore, their negative piezoelectricity also enables the regulation of spin splitting by strain. Finally, excellent optical absorption was also revealed. These intriguing properties make 2D few-layer ${\gamma}$-GeSe promising in spintronic and optoelectric applications., Comment: 18 pages,6 figures

Drone-based quantum link has the potential to realize mobile quantum network, and entanglement distribution has been demonstrated using one and two drones. Here we report the first drone-based quantum key distribution (QKD), with average secure key rate larger than 8 kHz using decoy-state BB84 protocol with polarization coding. Compact acquisition, pointing, and tracking (APT) system and QKD modules are developed and loaded on a home-made octocopter, within takeoff weight of 30 kg. A robust link is established between the flying octocopter and a ground station separated 200 meters away and real-time QKD is performed for 400 seconds. This work shows potential of drone-based quantum communication for the future mobile quantum networks.

Inspired by recent experimental measurements [Guo \textit{et al.}, Phys. Rev. Lett.~\textbf{124}, 206602 (2020); Jim\'enez \textit{et al.}, Nature \textbf{592}, 370 (2021)] on frustrated quantum magnet SrCu$_2$(BO$_3$)$_2$ under combined pressure and magnetic fields, we study the related spin-$1/2$ Shastry-Sutherland (SS) model using state-of-the-art tensor network methods. By calculating thermodynamics, correlations and susceptibilities, we find, in zero magnetic field, not only a line of first-order plaquette-singlet (PS) to dimer-singlet phase transition ending with a critical point, but also signatures of the ordered PS transition with its critical endpoint terminating on this first-order line. Moreover, we uncover prominent magnetic barocaloric responses, a novel type of quantum correlation induced cooling effect, in the strongly fluctuating supercritical regime. Under finite fields, we identify a quantum phase transition from the PS phase to the spin supersolid phase that breaks simultaneously lattice translational and spin rotational symmetries. The present findings on the SS model are accessible in current experiments and would shed new light on exotic critical and supercritical phenomena in archetypal frustrated quantum magnets., Comment: Close to the published version. 7 pages, 4 figures (SM 9 pages, 12 figures)

Ultrarelativistic gamma-ray burst (GRB) jets are strong gravitational wave (GW) sources with memory-type signals. The plateau (or shallow decay) phases driven by the energy injection might appear in the early X-ray afterglows of GRBs. In this paper, we investigate the GW signal as well as X-ray afterglow emission in the framework of GRB jets with energy injection, and both short- and long-duration GRBs are considered. We find that, regardless of the case, because of the antibeaming and time delay effects, a rising slope emerging in the waveform of GW signal due to the energy injection lags far behind the energy ejection, and the typical frequency of the characteristic amplitudes falls within a low-frequency region of $\sim10^{-4}-10^{-6} \,{\rm Hz}$; and we consider that the GW memory triggered by GRB jets with energy injection are previously unaware and the nearby GRBs with strong energy injection might disturb the measurement of the stochastic GW background. Such GW memory detection would provide a direct test for models of energy injection in the scenario of GRB jets., Comment: 9 pages, 4 figures, accepted for publication in ApJ

Two-dimensional Euler insulators are novel kind of systems that host multi-gap topological phases, quantified by a quantised first Euler number in their bulk. Recently, these phases have been experimentally realised in suitable two-dimensional synthetic matter setups. Here we introduce the second Euler invariant, a familiar invariant in both differential topology (Chern-Gauss-Bonnet theorem) and in four-dimensional Euclidean gravity, whose existence has not been explored in condensed matter systems. Specifically, we firstly define two specific novel models in four dimensions that support a non-zero second Euler number in the bulk together with peculiar gapless boundary states. Secondly, we discuss its robustness in general spacetime-inversion invariant phases and its role in the classification of topological degenerate real bands through real Grassmannians. In particular, we derive from homotopy arguments the minimal Bloch Hamiltonian form from which the tight-binding models of any second Euler phase can be generated. Considering more concretely the gapped Euler phase associated with the tangent bundle of the four-sphere, we show that the bulk band structure of the nontrivial 4D Euler phase necessarily exhibits triplets of linked nodal surfaces (where the three types of nodal surfaces are formed by the crossing of the three successive pairs of bands within one four-band subspace). Finally, we show how to engineer these new topological phases in a four-dimensional ultracold atom setup. Our results naturally generalize the second Chern and spin Chern numbers to the case of four-dimensional phases that are characterised by real Hamiltonians and open doors for implementing such unexplored higher-dimensional phases in artificial engineered systems, ranging from ultracold atoms to photonics and electric circuits., Comment: 15 pages, 5 figures

Reciprocal spin-orbit coupling (SOC) via geometric phase with flat optics provides a promising platform for shaping and controlling paraxial structured light. Current devices, from the pioneering q-plates to the recent J-plates, provide only spin-dependent wavefront modulation without amplitude control. However, achieving control over all the spatial dimensions of paraxial SOC states requires spin-dependent control of corresponding complex amplitude, which remains challenging for flat optics. Here, to address this issue, we present a new type of flat-optics elements termed structured geometric phase gratings that is capable of conjugated complex-amplitude control for orthogonal input circular polarizations. By using a microstructured liquid crystal photoalignment technique, we engineered a series of flat-optics elements and experimentally showed their excellent precision in arbitrary SOC control. This principle unlocks the full-field control of paraxial structured light via flat optics, providing a promising way to develop an information exchange and processing units for general photonic SOC states, as well as extra-/intracavity mode convertors for high-precision laser beam shaping.

Kilonovae are generally believed to originate from the ejecta of binary neutron stars (NSs) or black hole and NS mergers. Free neutrons might be retained in the outermost layer of the ejecta to produce a precursor via $\beta$-decay. During the propagation of kilonovae to observers, a small percentage of them might be gravitationally lensed by foreground objects. In this paper, three lens models, i.e., the point-mass model, the singular isothermal sphere (SIS) model, and the Chang-Refsdal model, were taken into consideration to explore the light curves and polarizations of gravitationally lensed kilonovae. We found that if the time delay between two images exceeds the ejecta heating timescale for the lens mass $\sim 10^{10}~M_\odot$ in the SIS model, a tiny bump-like signal will be generated in the light curve, and the total luminosity will be magnified in all cases. The polarization of lensed kilonovae is significantly enhanced in most cases. Future detections of lensed kilonovae will impose constraints on the morphology of the ejecta and aid in the determination of the nature of compact object mergers and the search for strong gravitational lenses., Comment: 11 pages, 4 figures, accepted for publication in ApJ

We study the rotation effects of the hot and dense QCD matter in a non-perturbative regime by the gauge/gravity duality. We use the gravitational model that is designated to match the state-of-the-art lattice data on the thermal properties of (2+1)-flavor QCD and predict the location of the critical endpoint and the first-order phase transition line at large baryon chemical potential without rotation. After introducing the angular velocity via a local Lorentz boost, we investigate the thermodynamic quantities for the system under rotation in a self-consistent way. We find that the critical temperature and baryon chemical potential associated with the QCD phase transition decrease as the angular velocity increases. Moreover, some interesting phenomena are observed near the critical endpoint. We then construct the 3-dimensional phase diagram of the QCD matter in terms of temperature, baryon chemical potential, and angular velocity. As a parallel investigation, we also consider the gravitational model of $SU(3)$ pure gluon system, for which the 2-dimensional phase diagram associated with temperature and angular velocity has been predicted. The corresponding thermodynamic quantities with rotation are investigated., Comment: 22 pages, 24 figures

We look into the 2851 $^{12}$CO molecular clouds harboring $^{13}$CO structures to reveal the distribution of the projected angular separations and radial velocity separations between their internal $^{13}$CO structures. The projected angular separations are determined using the minimal spanning tree algorithm. We find that $\sim$ 50$\%$ of the angular separations fall in a narrow range of $\sim$ 3 - 7 arcmin with a median of $\sim$ 5 arcmin, and the corresponding radial velocity separations mainly range from $\sim$ 0.3 km s$^{-1}$ to 2.5 km s$^{-1}$. The mean and standard deviation of the angular separations of the internal $^{13}$CO structures within $^{12}$CO clouds appear to be universal, independent of the $^{12}$CO cloud angular areas and the counts of their internal $^{13}$CO structures. We also reveal a scaling relation between the $^{12}$CO cloud angular area and its harbored $^{13}$CO structure count. These results suggest there is a preferred angular separation between $^{13}$CO structures in these $^{12}$CO clouds, considering the distance effects. According to that, we propose an alternative picture for the assembly and destruction of molecular clouds: there is a fundamental separation for the internal structures of molecular clouds, the build-up and destruction of molecular clouds proceeds under this fundamental unit., Comment: 14 pages, 12 figures, accepted for publication in ApJ

We study the collective motion of self-propelled particles affected by the spatial-dependent noise based on the Vicsek rules. Only the particles inside the special region will affected by noise. The consideration of the spatial-dependent noise is closer to reality because of the complexity of the environment. Interestingly, we find that there exists an optimal amplitude of noise to adjust the average motional direction of the system. Particular orientation of the noisy region makes the motional direction of the system parallel to the orientation of the noisy region. The adjustment of the motional direction of the system also depends on the shape, the proportion and the spatial distribution of the noisy region. Our findings may inspire the capture of the key features of collective motion underlying various phenomena., Comment: 12 pages, 6 figures, 39 references

In the coalescence events of binary neutron star (NS) or a black hole (BH) and an NS, a BH hyperaccretion disk might be eventually formed. At very high mass accretion rates, MeV neutrinos will be emitted from this disk, which is called a neutrino-dominated accretion flow (NDAF). Neutrino annihilation in the space out of the disk is energetic enough to launch ultrarelativistic jets to power gamma-ray bursts. Moreover, vertical advection might exist in NDAFs, which can generate the magnetic buoyancy bubbles to release gamma-ray photons. In this paper, we visit the effects of the vertical advection in NDAFs on the disk structure and gamma-ray and neutrino luminosities for different accretion rates. Then we study the anisotropic emission of kilonovae and the following gravitational waves (GWs) driven by the gamma-ray photons and neutrinos from NDAFs. Comparing NDAFs without vertical advection, the neutrino luminosity and GW strains slightly decrease for the case with vertical advection, and the kilonovae will be brightened by the injected gamma-ray photons. The future joint multimessenger observations might distinguish whether the vertical advection exists in NDAFs or not after compact binary coalescences., Comment: 14 pages, 7 figures, accepted for publication in ApJ

We present wide-field and high-sensitivity CO(1-0) molecular line observations toward the Cassiopeia region, using the 13.7m millimeter telescope of the Purple Mountain Observatory (PMO). The CO observations reveal a large-scale highly filamentary molecular cloud within the Galactic region of 132\fdg0\,$\geq$\,$l$\,$\geq$\,122\fdg0 and -1\fdg0\,$\leq$\,$b$\,$\leq$\,3\fdg0 and the velocity range from approximately +1 to +4 km/s. The measured length of the large-scale filament, referred to as the Cassiopeia Filament, is about 390 pc. The observed properties of the Cassiopeia Filament, such as length, column density, and velocity gradient, are consistent with those synthetic large-scale filaments in the inter-arm regions. Based on its observed properties and location on the Galactic plane, we suggest that the Cassiopeia Filament is a spur of the Local arm, which is formed due to the galactic shear. The western end of the Cassiopeia Filament shows a giant arc-like molecular gas shell, which is extending in the velocity range from roughly -1 to +7 km/s. Finger-like structures, with systematic velocity gradients, are detected in the shell. The CO kinematics suggest that the large shell is expanding at a velocity of ~6.5 km/s. Both the shell and finger-like structures outline a giant bubble with a radius of ~16 pc, which is likely produced by stellar wind from the progenitor star of a supernova remnant. The observed spectral linewidths suggest that the whole Cassiopeia Filament was quiescent initially until its west part was blown by stellar wind and became supersonically turbulent., Comment: 46 pages, 19 figures, to be published by the AJ

We propose to uncover the topology of a pseudo-Hermitian Chern insulator by quantum quench dynamics. The Bloch Hamiltonian of the pseudo-Hermitian Chern insulator is defined in the basis of the q-deformed Pauli matrices, which are related to the representation of the deformed algebras. We show the bulk-surface duality of the pseudo-Hermitian phases, then further build a concrete relation between the static band topology and quench dynamics, in terms of the time-averaged spin textures. The results are also generalized into a fully nonequilibrium case where the postquench evolution is governed by a Floquet pseudo-Hermitian Hamiltonian. Furthermore, we propose a possible scheme to realize the seemingly challenging model in a bilayer lattice and detect the dynamics with a double-quench protocol., Comment: 10 pages, 6 figures

Applying the system of linear partial differential equations derived from the Mellin-Barnes representation and the Miller transformation, we present the GKZ-system of the Feynman integral of the 2-loop self energy diagram with 4 propagators. The codimension of the derived GKZ-system equals the number of independent dimensionless ratios among the external momentum squared and virtual mass squared. In total 536 hypergeometric functions are obtained in the neighborhoods of the origin and infinity, in which 30 linearly independent hypergeometric functions whose convergent regions have nonempty intersection constitute a fundamental solution system in a proper subset of the whole parameter space., Comment: latex, 306 pages, including 1 figure + 254 pages of supplementary material, accepted for publication in the European Physical Journal C. arXiv admin note: text overlap with arXiv:2206.04224

We compute the complete two-loop electroweak corrections to the Higgsstralung process $e^+e^-\rightarrow HZ$ at the future Higgs factory. The Feynman integrals involved in the computation are decomposed into linear combinations of a minimal set of master integrals taking advantage of the recent developments of integral reduction techniques. The master integrals are then evaluated by differential equations with boundary conditions provided by the auxiliary mass flow method. Our final result for given $\sqrt{s}$ is expressed as a piecewise function defined by several deeply expanded power series, which has high precision and can be further manipulated efficiently. Our calculation presents the first complete two-loop electroweak corrections for processes with four external particles., Comment: 5 pages, 2 figures. arXiv admin note: text overlap with arXiv:2209.14259