2026-02-10 | | Total: 4
We propose a forward-mode automatic differentiation (AD) framework for tensor renormalization group (TRG) methods. In this approach, evaluating the derivatives of the partition function up to order $k$ increases the matrix-multiplication cost by a factor of $(k+1)(k+2)/2$ compared to computing the free energy alone, while the memory footprint is only $k$ times that of the original calculation. In the limit where the derivatives of the SVD are neglected, we establish a theoretical correspondence between our forward-mode AD and conventional impurity methods. Numerically, we find that the proposed AD algorithm can calculate internal energy and specific heat significantly higher accuracy than the impurity method at comparable computational cost. We also provide a practical procedure to extract critical exponents from derivatives of the renormalized tensor in TRG calculations in both two and three dimensions.
We present the first systematic study of the radiative decays of charmed mesons using $2+1$-flavor clover fermion gauge ensembles generated by the CLQCD collaboration. One of the ensembles is at the physical pion mass, and one has a fine lattice spacing $a\sim 0.05 ~\text{fm}$. We determine the coupling constants to be $g_{D^{\ast+} D^+ γ} = -0.205(35)$ GeV$^{-1}$, $g_{D^{\ast0} D^0 γ} = 2.20(27)$ GeV$^{-1}$, and $g_{D_s^{\ast+} D_s^+ γ} = -0.125(21)$ GeV$^{-1}$, respectively. Compared with previous studies, our results demonstrate significant improvements in precision. In particular, we carefully estimate the systematic uncertainty arising from matrix element fits, momentum transfer extrapolations, and chiral and continuum limit extrapolations, which are included in the reported total uncertainties. These couplings yield the following predictions of decay widths: $Γ_{D^{\ast+} \rightarrow D^+ γ} = 0.255(87)$ keV, $Γ_{D^{\ast0} \rightarrow D^0 γ} = 29.4(7.2)$ keV, and $Γ_{D_s^{\ast+}\rightarrow D_s^+ γ} = 0.102(34)$ keV. This work establishes first-principles results of the charmed meson radiative transitions and provides inputs for understanding the structure and properties of heavy-light mesons.
Neutrinoless double beta ($0νββ$) decays of spin-1/2 hyperons are investigated in a covariant baryon chiral perturbation theory framework, extended by a $ΔL=2$ operator proportional to the Majorana neutrino mass, where $L$ denotes the lepton number. Within the light Majorana neutrino exchange mechanism, the decay amplitudes are found to emerge at the one-loop level, representing the long-range contribution. The extended-on-mass-shell scheme is employed to renormalize the one-loop amplitudes and restore consistent chiral power counting. Consequently, the differential decay rates for all accessible hyperon $0νββ$ channels are predicted and the corresponding branching ratios are more than 20 orders of magnitude smaller than the current experimental upper bounds. Interestingly, it is found that the leading contribution to hyperon $0νββ$ decay is actually from short-range counterterm operators, as required by the renormalization argument. Neutrinoless transition form factors are proposed to determine this leading contribution through future lattice QCD simulations.
We review recent results of net-proton multiplicity fluctuations from STAR experiment, aiming to locate the QCD critical point in high-energy nuclear collisions at RHIC. We show net-proton number cumulant and proton number factorial cumulant ratios up to fourth order using experimental data from RHIC BES-II Au+Au collisions in collider mode and fixed-target mode. The comparison is made between experimental data and non-critical model calculations from Lattice QCD, HRG, hydrodynamic simulations and transport model UrQMD. In addition, we discuss initial volume fluctuation effect, which plays significant role in fixed-target energies. Finally, an outlook on experimental research on the QCD critical point in future experiments will be presented.