Cluster radioactivity has been successfully described as a super-asymmetric fission mode within the microscopic self-consistent Gogny Hartree-Fock-Bogoliubov approximation [Phys. Rev. C 84, 044608 (2011)]. For nuclei preserving the neutron-to-proton $N/Z$ ratio of the doubly magic $^{208}$Pb, a cluster radioactivity fission valley has been identified. Such a valley can also be found both in actinides and super-heavy nuclei. In this paper, chains of isotopes and isotones are examined to determine the limits of existence of the cluster radioactivity fission mode. It is shown that the super-asymmetric valley can be found in a wide range of the nuclear chart. Nevertheless, the valley flattens more and more when diverging from the isospin asymmetry of $^{208}$Pb. For neutron-deficient nuclei with $N/Z <$ 1.41, it is found that the valley diminishes before reaching the scission point, and cluster radioactivity can not be observed.

Neutron skins encode essential information about nuclear geometry, surface structure, and isovector response, yet a compact description across the nuclear chart remains elusive. We present an empirical analysis of neutron-excess surface systematics using a mass-normalized, charge-radius-derived proxy ("neutron-skin curvature") built from evaluated experimental charge radii. By normalizing radii to the reduced Compton length $r_B=\hbar/(mc)$, we form a dimensionless curvature ratio that enables comparison across isotopic chains of widely varying mass. When expressed versus normalized neutron excess, data for more than 800 nuclei spanning 88 elements collapse onto a single empirical curve without element-specific rescaling or interaction-model tuning; the curve is used only as a fixed baseline for residual analysis. The collapse accounts for approximately 88% of the variance and is substantially tighter than droplet-style baselines fit to the same dataset. Residuals show structured deviations: three finite-size regimes (skin formation, relaxation toward bulk geometry, and saturation) and a distinct few-body domain for very light nuclei ($Z\le 4$). Stratifying residuals by periodic-table families reveals tighter submanifolds for several groups, suggesting additional geometric constraints layered on the global trend. These results are obtained directly from evaluated experimental data and physical constants, without introducing new interaction terms, and motivate further study of geometric correlations with other nuclear and atomic observables.

The Pauli exclusion principle imposes an important structural constraint in cluster descriptions of light nuclei and is commonly taken into account using methods such as the Resonating Group Method (RGM), the Orthogonality Condition Model (OCM), and the Orthogonalizing Pseudopotential (OPP) approach. The latter provides a practical implementation for suppressing Pauli-forbidden states in few-body calculations through the introduction of a large auxiliary coupling constant $\lambda_0$ and an associated limiting procedure. Exact $\lambda_0$-eliminated formulations have appeared in the literature. It is shown that the OPP method may be interpreted as the singular $\lambda_0\to\infty$ limit of the Feshbach-Schur projection, and that the Schur complement provides a natural operator-level framework for understanding this connection. In contrast to earlier approaches, the elimination is derived explicitly as a closed Schur-complement operator identity. When formulated in terms of operators and Green's functions, the Feshbach-Schur projection eliminates Pauli-forbidden components algebraically and avoids the explicit introduction of large pseudopotential parameters. This reformulation clarifies the implementation of Pauli projection in cluster models and provides a convenient framework for few-body calculations.

The properties of nuclear matter and the structures of neutron stars are analyzed with a baryonic extended linear sigma model in mean-field approximation, where the masses of baryons and mesons are generated via the spontaneous chiral symmetry breaking. The couplings between the iso-scalar scalar meson and nucleons, $g_{\sigma NN}$, the iso-vector scalar meson and nucleons, $g_{a_0 NN}$, and the four-vector meson couplings play an important role in the properties of nuclear matter and neutron stars. The introduction of the $\delta$ meson leads to a plateau structure of the symmetry energy, $E_{\rm sym}(n)$, at intermediate densities, which is crucial to the consistency of neutron skin thickness of $^{208}$Pb and the tidal deformability of a canonical neutron star. The explicit chiral symmetry breaking term is then introduced with a constant background field, $\xi$, which can be related to the current quark mass and thus the pion-nucleon sigma term, $\sigma_{\pi N}$. A negative $\sigma_{\pi N}$ leads to a stiffer EOS of neutron star matter and thus a larger maximum mass of neutron stars, but the value of $\sigma_{\pi N}$ needed to satisfy the astrophysical constraints is negative, not positive as the vacuum value. The study may provide insights into the running behaviors of the parameters in the low-energy effective model to give the density-dependent description for the EOS of neutron star matter.

Renormalization group methods generate low-resolution Hamiltonians that are more diagonal and easier to solve. This chapter reviews the similarity renormalization group for nuclear Hamiltonians, which is a popular method for generating low-resolution nuclear forces. It presents the similarity renormalization group flow equations, analyzes how the similarity renormalization group drives the Hamiltonian towards the diagonal, and studies the effect of induced many-body interactions. It concludes by highlighting the progress in first-principles calculations of nuclei driven by low-resolution nuclear Hamiltonians.

We investigate the $T_{cc}$ tetraquark, treating it as a bound state of a heavy diquark and a light antidiquark. Using the Silvestre-Brac potential and solving the Schrödinger equation via the Gaussian Expansion Method, we find that the excitation energy between the heavy diquark and light antidiquark is unexpectedly larger than that between the two light anti-quarks within the anti-diquark -- contrary to the naive expectation where the former is smaller than the latter. We trace this inversion of the mass hierarchy to the centrifugal force acting on the light degree of freedom. Applying the same framework to other systems ($T_{bb}, \Lambda_b, \Lambda_c$) yields qualitatively identical behavior, demonstrating the robustness of the mechanism. These results provide new insights into diquark dynamics and the mass structure of exotic hadrons.

The Collins-Soper (CS) kernel may be obtained through the TMD soft function by formulating the Wilson line in terms of 1-dimensional auxiliary fermion fields on the lattice. Our computation takes place in the region of the lattice that corresponds to the "spacelike" region in Minkowski space, i.e., Collins' scheme. We explore two methods for obtaining the CS kernel. The "ratio method"; which would allow us to obtain the soft function as well as the CS kernel. And the "double ratio"; which allows us to achieve a high degree of statistical precision, but only produces the CS kernel. The matching of our result to Minkowski space is achieved through the mapping of the complex auxiliary field directional vector to the Wilson line rapidity. We present a preliminary extraction of the CS kernel using the "double ratio", and discuss the methodology employed.

Within the framework of SU(2) chiral perturbation theory, we derive the general solution of the QCD $\theta$-vacuum for an arbitrary vacuum phase, explicitly incorporating isospin-breaking effects from the light quark mass difference, and compute the temperature dependence of the topological susceptibility, higher-order cumulants, and the domain wall tension up to next-to-leading order. We find that the topological susceptibility agrees with lattice data at low temperatures but deviates at higher temperatures as expected from the breakdown of the chiral expansion; moreover, we demonstrate that the normalized fourth-order cumulant and the domain wall tension decrease monotonically with increasing temperature, while the normalized sixth-order cumulant exhibits the opposite behavior. These results extend earlier analyses by showing how isospin breaking reshapes the full hierarchy of topological charge cumulants and the dynamics of $\theta$-vacuum domain walls, thereby offering new theoretical input on the $\theta$-vacuum properties, which are relevant for axion-related effective theories in hot QCD matter.

We develop a quantum kinetic theory for QCD, which incorporates all leading order collision terms. At lowest order in gradient expansion, it reproduces the spin-averaged Boltzmann equation with both elastic and inelastic collisions. At next order in gradient expansion, the solution to the quantum kinetic equations give spin polarization of on-shell quarks and gluons in quark-gluon plasma when the gradients are of hydrodynamic ones. A power counting in the coupling shows the spin polarization behaves differently in vortical and non-vortical gradients: the former is free of collisional contribution to leading order, while the latter contains a collisional contribution at parametrically the same order as the free theory counterpart. We also find the inelastic collision in a spin basis provides a possible mechanism for conversion between spin and orbital angular momentum.

We present theoretical and experimental evidence for a new phase of matter in neutron-rich systems that we call multimodal superfluidity. Using ab initio lattice calculations, we show that the condensate consists of coexisting s-wave pairs, p-wave pairs in entangled double pair combinations, and quartets composed of bound states of two s-wave pairs. We identify multimodal superfluidity as a general feature of single-flavor spin-1/2 fermionic systems with attractive s-wave and p-wave interactions, provided the system is stable against collapse into a dense droplet. Beyond neutrons at sub-saturation densities, we demonstrate that this phase appears in generalized attractive extended Hubbard models in one, two, and three dimensions. We elucidate the mechanism for this coexistence using self-consistent few-body Cooper models and compare with Bardeen-Cooper-Schrieffer theory. We also derive the form of the effective action and show that spin, rotational, and parity symmetries remain unbroken. Finally, we analyze experimental data to show that p-wave pair gaps and quartet gaps are present in atomic nuclei, and we discuss the consequences of this new phase for the structure and dynamics of neutron star crusts.

We study some interesting aspects of the spectral properties of SU(3) gauge theory, both with and without dynamical quarks (QCD) at thermal equilibrium using lattice gauge theory techniques. By calculating the eigenstates of a massless overlap Dirac operator on the gauge configurations, we implement a gauge-invariant method to study spectral properties of non-Abelian gauge theories. We have unambiguously categorized Dirac eigenvalues into different regimes based on a quantity defined in terms of the ratios of nearest neighbor spacings. While majority of these eigenstates below the magnetic scale are similar to those of random matrices belonging to the Gaussian Unitary ensemble at temperatures much higher than the chiral crossover transition in QCD, a few among them start to become prominent only near the crossover. These form fractal-like clusters with the median value for their fractal dimensions hinting at the universality class of the chiral transition in QCD. We further demonstrate that momentum modes below the magnetic scale in a particular non-equilibrium state of QCD are classically chaotic and estimate an upper bound on the thermalization time $\sim 1.44$ fm/c by matching this magnetic scale with that of a thermal state at $\sim 600$ MeV.

Dark matter may accumulate in neutron stars given its gravitational interaction and abundance. We investigate the influence of strongly-interacting dark matter, described by a QCD-like one-flavor $G_2$ gauge theory, on neutron stars. This choice allows to test, for the first time, a first-principles-determined non-Abelian dark matter equation of state, which supports composite fermionic dark matter and thus a Fermi-pressure-stabilized dark matter component. The ordinary matter part of the mixed star is described by available model-agnostic equations of state that interpolate between the low-density regime and high-density regime. We find that strongly-interacting dark matter has a similar impact on neutron stars as other model equation of states and confirm that strongly-interacting dark matter can be accommodated by constraints from neutron star observations within our uncertainties.

In this paper we include kinematic power corrections up to twist-four to the deeply virtual Compton scattering dispersion relation. We demonstrate that, both for (pseudo-)scalar and spin-$1/2$ targets, the formal expression of the $n$-subtracted leading-twist dispersion relations is preserved. However, the expression of the subtracted constants is modified by the kinematic powers. Importantly, the minimal-subtracted dispersion relation for the helicity-conserving amplitude, previously thought to depend only on the Polyakov-Weiss $D$-term, now also depends on the double distributions $F$ and $K$. These results are consistent with the ones obtained previously in the literature. Such a mixing may be critical for the Jefferson Lab kinematic range, as it is not suppressed for typical values of $t$ and $Q^{2}$ in the valence region. We therefore expect a strong impact on attempts to extract pressure forces from DVCS data.

In this letter, we analyse and interpret the kinematic power corrections to deeply virtual Compton scattering dispersion relation. We show that the kinematic corrections at twist-4 can be connected to other form factors of the Energy-Momentum Tensor beyond the pressure distribution involved at leading-power, namely the ones related to Momentum and total Angular Momentum distributions. In the nucleon case, these corrections are not negligible at presently accessible virtualities. The DVCS subtraction constant becomes an experimental constraint on momentum distributions, pressure forces distributions, and total angular momentum distributions. Finally, we use continuum and lattice-QCD results to predict the expected size of the DVCS subtraction constant, and conclude that momentum distributions are responsible of roughly one-third of the experimental signal at $Q^2 = 2\textrm{GeV}^2$.

We study the spatial distributions of pressure, energy density, and shear forces inside the nucleon within the two-flavor Skyrme model including vector mesons. This framework has the advantage that nucleon configurations can be stabilized without the Skyrme term. In contrast to the model without vector mesons, however, we realize that the energy-momentum tensor (EMT) becomes pseudo-gauge dependent. We explicitly demonstrate that all these distributions differ between the canonical and Belinfante forms of the EMTs. We identify the pseudo-gauge ambiguity as originating from nonvanishing surface terms associated with spin currents generated by the vector-meson field strength tensors. Furthermore, we show that the pressure and shear-force distributions in the canonical EMT develop singularities at the nucleon center, whereas the corresponding Belinfante distributions remain finite. Finally, we discuss the implications of pseudo-gauge dependence for extracting the confining force and for constructing the equation of state inside the nucleon.

Jet quenching - the modification of high-energy jets in the quark-gluon plasma - has been extensively studied through weakly coupled scattering amplitudes embedded in parton-shower frameworks. These models, often combined with bulk hydrodynamic evolution, successfully describe a wide range of observables, though they typically rely on assumptions of rapid thermalization and simplified treatments of medium response. Parallel to these developments, jet thermalization has been investigated within the finite-temperature QCD effective kinetic theory, which provides our best microscopic understanding of equilibration in heavy-ion collisions. Early studies of linearized perturbations have highlighted both the promise and the limitations of current approaches, as existing MC implementations face challenges - particularly in the treatment of recoils and particle merging. Building on this foundation, we introduce a new parton-shower algorithm that exactly reproduces the dynamics of the linearized EKT, enabling a first-principles description of jet thermalization with proper inclusion of recoils, holes, quantum statistics, and merging processes.

We present the first next-to-leading-logarithmic QCD analysis of the electromagnetic corrections to the semileptonic weak Hamiltonian, including the mixed $\mathcal{O}(\alpha\,\alpha_s^2)$ corrections to the vector coupling $g_V$. The analysis combines the evaluation of three-loop anomalous dimensions and two-loop matching corrections with a consistent factorization of short-distance QCD effects. The latter is implemented through a scheme change based on a $d$-dimensional operator product expansion performed inside the loop integrals. The resulting renormalization-group--improved expression for the radiative correction $\Delta^V_R = 2.436(16)\%$ can be systematically refined using input from lattice QCD and perturbation theory and improves the consistency of first-row CKM unitarity tests.

Hybrid qubit-qumode quantum computing platforms provide a natural setting for simulating interacting bosonic quantum field theories. However, existing continuous-variable gate constructions rely predominantly on polynomial functions of canonical quadratures. In this work, we introduce a complementary universality paradigm based on trigonometric continuous-variable gates, which enable a Fourier-like representation of bosonic operators and are particularly well suited for periodic and non-perturbative interactions. We present a deterministic ancilla-based method for implementing unitary and non-unitary trigonometric gates whose arguments are arbitrary Hermitian functions of qumode quadratures. As a concrete application, we develop a hybrid qubit-qumode quantum simulation of the lattice sine-Gordon model. Using these gates, we prepare ground states via quantum imaginary-time evolution, simulate real-time dynamics, compute time-dependent vertex two-point correlation functions, and extract quantum kink profiles under topological boundary conditions. Our results demonstrate that trigonometric continuous-variable gates provide a physically natural framework for simulating interacting field theories on near-term hybrid quantum hardware, while establishing a parallel route to universality beyond polynomial gate constructions. We expect that the trigonometric gates introduced here to find broader applications, including quantum simulations of condensed matter systems, quantum chemistry, and biological models.

A size-extensive, converging, black-box, ab initio coupled-cluster ($\Delta$CC) ansatz is introduced that computes the energies and wave functions of stationary states from any degenerate or nondegenerate Slater-determinant references with any numbers of $\alpha$- and $\beta$-spin electrons, any patterns of orbital occupancy, any spin multiplicities, and any spatial symmetries. For a nondegenerate reference, it reduces to the single-reference coupled-cluster ansatz. For a degenerate multireference, it is a natural coupled-cluster extension of degenerate Rayleigh-Schrödinger perturbation ($\Delta$MP) theory. For ionized and electron-attached references, it can be viewed as a coupled-cluster Green's function, although the present theory is convergent toward the full-configuration-interaction (FCI) limits, while Feynman-Dyson many-body Green's function (MBGF) theory generally is not. Additionally, a new state-universal multireference coupled-cluster theory for general model spaces is developed by slightly modifying the $\Delta$CC ansatz. This quasidegenerate coupled-cluster (QCC) theory is size-extensive, converging, but not black-box, which is expected to be well suited for strong correlation. Determinant-based, general-order algorithms of $\Delta$CC and QCC theories are implemented, which are compared with configuration-interaction (CI) and equation-of-motion coupled-cluster (EOM-CC) theories through octuple excitations and with $\Delta$MP and MBGF theories up to the nineteenth order. For transition energies, the order of performance is: QCC $\approx$ $\Delta$CC $>$ EOM-CC $>$ CI at the same excitation order or QCC $\approx$ $\Delta$CC $>$ $\Delta$MP $>$ MBGF at the same cost scaling.