While hybrid approaches of relativistic hydrodynamics+transport have been well established for the dynamical description of heavy-ion collisions at high beam energies, moving to lower beam energies is challenging. In this work, we propose dynamic initial conditions for the viscous hydrodynamic evolution in heavy-ion collisions at low to intermediate beam energies. They are comprised of core hadrons based on the local energy density during the pre-equilibrium hadronic evolution. The SMASH-vHLLE hybrid approach is then applied to lower beam energies, achieving good agreement with measured bulk observables between $\sqrt{s_\mathrm{NN}} = 3$ and $9.1\ \mathrm{GeV}$, thus providing guidance for measurements in STAR-BES and CBM at FAIR. This framework has the advantage that the equation of state of dense QCD matter is a direct input to the calculation, so the current results can be regarded as a baseline with respect to future studies of the phase transition.

We present optical + near-infrared (NIR) + mid-infrared (MIR) observations of the normal Type Ia supernovae (SN Ia) 2022aaiq and 2024gy in the nebular phase, continuously spanning 0.35-28 microns. Medium-resolution JWST spectroscopy reveals novel narrow (v_{\mathrm{FWHM}}<1500 km s$^{-1}$) [Ni II] 1.94 and 6.64 micron cores in both events. The MIR [Ni II] 6.64 micron line exhibits a distinct narrow core atop a broader base, indicating a central enhancement of stable Ni. This structure points to high central densities consistent with a near-Chandrasekhar-mass ($M_{Ch}$) progenitor or a high-metallicity sub-$M_{Ch}$ progenitor. From detailed line-profile inversions of SN 2024gy, we derive emissivity profiles for stable iron-group elements (IGEs), radioactive material, and intermediate-mass elements (IMEs), revealing spatially distinct ejecta zones. The [Ni III] 7.35 micron line shows a shallow-to-steep slope transition -- a "broken-slope" morphology -- that matches predictions for delayed detonation explosions with separated deflagration and detonation ashes. We also reanalyze and compare to archival JWST spectra of SN 2021aefx and the subluminous SN 2022xkq. We estimate a stable $^{58}$Ni mass of $\sim0.1$ M$_\odot$ for SN 2024gy, consistent with delayed detonation models, and $\sim0.01$ M$_\odot$ for SN 2022xkq, favoring sub-$M_{Ch}$ scenarios. These results demonstrate that resolved line profiles, now accessible with JWST, provide powerful diagnostics of explosion geometry, central density, and progenitor mass in SN Ia.

We present first-principles numerical calculations of the depolarization rate of spin-polarized deuterium and tritium nuclei in realistic tokamak plasmas, driven by resonant interactions with plasma waves. Backed up by first-of-a-kind linear and nonlinear simulations, we find that alpha particle-driven Alfvénic modes cause only negligible depolarization, which is contrary to expectations in prior literature. Other Alfvénic instabilities can in principle degrade polarization, but only under conditions unlikely to be realized on transport timescales. By combining full-orbit particle tracing with a dedicated depolarization solver, we demonstrate that wave-driven depolarization is surprisingly weak in SPARC and ITER-scale devices. These results provide strong evidence that spin-polarized fuel can maintain its polarization long enough to boost fusion reactivity, opening a viable path toward substantially enhanced performance in magnetic confinement fusion power plants.

Heavy-ion collisions at varying collision energies provide access to different regions of the QCD phase diagram. In particular collisions at intermediate energies are promising candidates to experimentally identify the postulated first order phase transition and critical end point. While heavy-ion collisions at low and high collision energies are theoretically well described by transport approaches and hydrodynamics+transport hybrid approaches, respectively, intermediate energy collisions remain a challenge. In this work, a modular hybrid approach, the SMASH-vHLLE-hybrid coupling 3+1D viscous hydrodynamics (vHLLE) to hadronic transport (SMASH), is introduced. It is validated and subsequently applied in Au+Au/Pb+Pb collisions between $\sqrt{s_\mathrm{NN}}$ = 4.3 GeV and $\sqrt{s_\mathrm{NN}}$ = 200.0 GeV to study the rapidity and transverse mass distributions of identified particles as well as excitation functions for $\mathrm{dN}/\mathrm{d}y|_{y = 0}$ and $\langle p_\mathrm{T} \rangle$. A good agreement with experimental measurements is obtained, including the baryon stopping dynamics. The transition from a Gaussian rapidity spectrum of protons at lower energies to the double-hump structure at high energies is reproduced. The centrality and energy dependence of charged particle $v_2$ is also described reasonably well. This work serves as a basis for further studies, e.g. systematic investigations of different equations of state or transport coefficients.

From ultracold atoms to quantum chromodynamics, reliable ab initio studies of strongly interacting fermions require numerical methods, typically in some form of quantum Monte Carlo calculation. Unfortunately, (non)relativistic systems at finite density (spin polarization) generally have a sign problem, such that those ab initio calculations are impractical. It is well-known, however, that in the relativistic case imaginary chemical potentials solve this problem, assuming the data can be analytically continued to the real axis. Is this feasible for nonrelativistic systems? Are the interesting features of the phase diagram accessible in this manner? By introducing complex chemical potentials, for real total particle number and imaginary polarization, the sign problem is avoided in the nonrelativistic case. To give a first answer to the above questions, we perform a mean-field study of the finite-temperature phase diagram of spin-1/2 fermions with imaginary polarization.

We employ the SMASH transport model to provide event-by-event initial conditions for the energy-momentum tensor and conserved charge currents in hydrodynamic simulations of relativistic heavy-ion collisions. We study the fluctuations and dynamical evolution of three conserved charge currents (net baryon, net electric charges, and net strangeness) with a 4D lattice-QCD-based equation of state, NEOS-4D, in the hydrodynamic phase. Out-of-equilibrium corrections at the particlization are generalized to finite densities to ensure the conservation of energy, momentum, and the three types of charges. These theoretical developments are integrated within X-SCAPE as a unified framework for studying the nuclear matter properties in the Beam Energy Scan program.

The doubly charm tetraquark with flavor $cc\bar u\bar d$ and isospin $I\!=\!0$ is investigated by calculating the $DD^*$ scattering amplitude with lattice QCD. The simulation is done on CLS ensembles with dynamical $u/d,s$ quarks and $m_\pi\simeq 280~$MeV for two charm quark masses, one slightly larger and one slightly lower than the physical value. The scattering amplitudes for partial waves $l=0,1$ are extracted near-threshold via the Lüscher's method by considering systems with total momenta $PL/(2\pi)=0,1,\sqrt{2},2$ on two spatial volumes. A virtual bound state pole in the $DD^*$ scattering amplitude with $l=0$ is found $9.9_{-7.1}^{+3.6}~$MeV below $DD^*$ threshold for the charm quark mass closer to the physical value. This pole is likely related to the doubly charm tetraquark discovered by LHCb less than $1~$MeV below $D^0D^{*+}$ threshold. Future lattice simulations closer to the continuum limit and physical quark masses would be valuable to establish this connection systematically.

We investigate the angular distribution of electrons that are emitted in the ionization of hydrogen-like ions by twisted photons. Analysis is performed based on the first-order perturbation theory and the non-relativistic Schrödinger equation. Special attention is paid to the dependence of the electron emission pattern on the impact parameter b of the ion with respect to the centre of the twisted wave front. In order to explore such a dependence, detailed calculations were carried out for the photoionization of the 1s ground and 2 py excited states of neutral hydrogen atoms. Based on these calculations, we argue that for relatively small impact parameters the electron angular distributions may be strongly affected by altering the position of the atom within the wave front. In contrast, if the atom is placed far from the front centre, the emission pattern of the electrons is independent on the impact parameter b and resembles that observed in the photoionization by plane wave photons.

We compute net-proton probability distributions in heavy ion collisions within the hadron resonance gas model. The model results are compared with data taken by the STAR Collaboration in Au-Au collisions at sqrt(s_{NN})= 200 GeV for different centralities. We show that in peripheral Au-Au collisions the measured distributions, and the resulting first four moments of net-proton fluctuations, are consistent with results obtained from the hadron resonance gas model. However, data taken in central Au-Au collisions differ from the predictions of the model. The observed deviations can not be attributed to uncertainties in model parameters. We discuss possible interpretations of the observed deviations.

Recently significant attention has been paid to magnetic-field-dependent photoluminescence (PL) features of the negatively charged nitrogen-vacancy (NV) centers in diamond. These features are used for microwave-free sensing and are indicative of the spin-bath properties in the diamond sample. Examinating the temperature dependence of the PL features allows to identify both temperature dependent and independent features, and to utilize them in diamond-based quantum sensing and dynamic nuclear polarization applications. Here, we study the thermal variability of many different features visible in a wide range of magnetic fields. To this end, we first discuss the origin of the features and tentatively assign the previously unidentified features to cross relaxation of NV center containing multi-spin systems. The experimental results are compared with theoretically predicted temperature shifts deduced from a combination of thermal expansion and electron-phonon interactions. A deeper insight into the thermal behavior of a wide array of the features may come with important consequences for various applications in high-precision NV thermometry, gyroscopes, solid-state clocks, and biomagnetic measurements.

The fragmentation cross sections of carbon ion beams with kinetic energies of $115 - 353~\text{MeV/u}$ impinging on thin targets of graphite (C), polyvinyl-toluene (C$_9$H$_{10}$) and PMMA (C$_2$O$_5$H$_8$) have been measured at 90$^{\text{o}}$ and 60$^{\text{o}}$ at the CNAO particle therapy center (Pavia, Italy). The presented measurements are a complete reanalysis by the FOOT collaboration of already published elemental cross section on composite targets, in order to refine the analysis, improve the systematic uncertainties and show the comparison with the FLUKA Monte Carlo code calculations. In this work, the kinetic energy at production of measured fragments has been completely redefined, together with the efficiencies computation. The new analysis strategy has been successfully validated against the Monte Carlo cross sections. Two detection arms were positioned at two different angles to perform the measurement at 90$^{\text{o}}$ and 60$^{\text{o}}$. The fragment species have been identified in charge (Z$_{id}$ = H) and mass (M$_{id}$ = $^1$H, $^2$H, $^3$H) combining the information of the deposited energy in thin plastic scintillators, of the deposited energy in a thick LYSO crystal and of the fragments Time of Flight (ToF) measurement. The ToF was also used to compute the fragments measured kinetic energy. The cross sections are presented as a function of the fragments kinetic energy at production thanks to an unfolding technique applied to the data.

We propose a resolution of the discrepancy between the proton yield predicted by the statistical hadronization approach and data on hadron production in ultra-relativistic nuclear collisions at the LHC. Applying the S-matrix formulation of statistical mechanics to include pion-nucleon interactions, we reexamine their contribution to the proton yield, taking resonance widths and the presence of nonresonant correlations into account. The effect of multi-pion-nucleon interactions is estimated using lattice QCD results on the baryon-charge susceptibility. We show that a consistent implementation of these features in the statistical hadronization model, leads to a reduction of the predicted proton yield, which then quantitatively matches data of the ALICE collaboration for Pb-Pb collisions at the LHC.

We develop a (3+1)-dimensional hybrid evolution model for heavy-ion collisions with dynamical sources for the energy-momentum tensor and baryon current. During an initial pre-equilibrium stage based on UrQMD, the four-momenta and baryon numbers carried by secondary particles created within UrQMD are fed continuously, after a short thermalization time, into a (3+1)-dimensional viscous hydrodynamic evolution module including baryon transport. The sensitivity of the initial conditions to model parameters and the effect of baryon diffusion on the hydrodynamic evolution are studied.

Precision optical spectroscopy of exotic ions reveals accurate information about nuclear properties such as charge radii and magnetic and quadrupole moments. Thorium ions exhibit unique nuclear properties with high relevance for testing symmetries of nature. We report loading and trapping of single $^{232}$Th$^+$ ions in a linear Paul trap, embedded into and sympathetically cooled by small crystals of trapped $^{40}$Ca$^+$ ions. Trapped Th ions are identified in a non-destructive manner from the voids in the laser-induced Ca fluorescence pattern emitted by the crystal, and alternatively, by means of a time-of-flight signal when extracting ions from the Paul trap and steering them into an external detector. We have loaded and handled a total of 231 individual Th ions. We reach a time-of-flight detection efficiency of $\gtrsim 95\, \%$, consistent with the quantum efficiency of the detector. The sympathetic cooling technique is expected to be applicable for other isotopes and various charge states of Th e.g., for future studies of $^{229m}$Th.

The doubly charm tetraquark with flavor $cc\bar u\bar d$ and isospin $I\!=\!0$ is investigated by calculating the $DD^*$ scattering amplitude with lattice QCD. The simulation is done on CLS ensembles with dynamical $u/d,s$ quarks and $m_\pi\simeq 280~$MeV for two charm quark masses, one slightly larger and one slightly lower than the physical value. The scattering amplitudes for partial waves $l=0,1$ are extracted near-threshold via the Lüscher's method by considering systems with total momenta $PL/(2\pi)=0,1,\sqrt{2},2$ on two spatial volumes. A virtual bound state pole in the $DD^*$ scattering amplitude with $l=0$ is found $9.9_{-7.1}^{+3.6}~$MeV below $DD^*$ threshold for the charm quark mass closer to the physical value. This pole is likely related to the doubly charm tetraquark discovered by LHCb less than $1~$MeV below $D^0D^{*+}$ threshold. Future lattice simulations closer to the continuum limit and physical quark masses would be valuable to establish this connection systematically.

We report on the first application of the stochastic Laplacian Heaviside method for computing multi-particle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a positive definite set of two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance in order to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the $SU(3)$-flavor symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of $a\sim0.086$~fm, lattice volume of $L=48a\simeq4.1$~fm and pion mass $m_\pi\simeq714$ MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

The relativistic method of moments is one of the most successful approaches to extract second order viscous hydrodynamics from a kinetic underlying background. The equations can be systematically improved to higher order, and they have already shown a fast convergence to the kinetic results. In order to generalize that method, we introduced long range effects in the form of effective (medium dependent) masses and gauge (coherent) fields. The most straightforward generalization of the hydrodynamic expansion is problematic, or simply ill-defined, at higher order. Instead of introducing an additional set of approximations, we propose to rewrite the series in terms of moments resumming the contributions of infinite non-hydrodynamics modes. The resulting equations are consistent with hydrodynamics and well defined at all order. We tested the new approximation against the exact solutions of the Boltzmann-Maxwell-Vlasov equations in $(0+1)$-dimensions, finding a fast and stable convergence to the exact results.

Anisotropic hydrodynamics improves upon standard dissipative fluid dynamics by treating certain large dissipative corrections non-perturbatively. Relativistic heavy-ion collisions feature two such large dissipative effects: (i) Strongly anisotropic expansion generates a large shear stress component which manifests itself in very different longitudinal and transverse pressures, especially at early times. (ii) Critical fluctuations near the quark-hadron phase transition lead to a large bulk viscous pressure on the conversion surface between hydrodynamics and a microscopic hadronic cascade description of the final collision stage. We present a new dissipative hydrodynamic formulation for non-conformal fluids where both of these effects are treated nonperturbatively. The evolution equations are derived from the Boltzmann equation in the 14-moment approximation, using an expansion around an anisotropic leading-order distribution function with two momentum-space deformation parameters, accounting for the longitudinal and transverse pressures. To obtain their evolution we impose generalized Landau matching conditions for the longitudinal and transverse pressures. We describe an approximate anisotropic equation of state that relates the anisotropy parameters with the macroscopic pressures. Residual shear stresses are smaller and are treated perturbatively, as in standard second-order dissipative fluid dynamics. The resulting optimized viscous anisotropic hydrodynamic evolution equations are derived in 3+1 dimensions and tested in a (0+1)-dimensional Bjorken expansion, using a state-of-the-art lattice equation of state. Comparisons with other viscous hydrodynamical frameworks are presented.

The microscopic description of heavy-ion reactions at low beam energies is achieved within hadronic transport approaches. In this article a new approach SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) is introduced and applied to study the production of non-strange particles in heavy-ion reactions at $E_{\rm kin}=0.4-2A$ GeV. First, the model is described including details about the collision criterion, the initial conditions and the resonance formation and decays. To validate the approach, equilibrium properties such as detailed balance are presented and the results are compared to experimental data for elementary cross sections. Finally results for pion and proton production in C+C and Au+Au collisions is confronted with HADES and FOPI data. Predictions for particle production in $\pi+A$ collisions are made.

We argue that features of hadron production in relativistic nuclear collisions, mainly at CERN-SPS energies, may be explained by the existence of three forms of matter: Hadronic Matter, Quarkyonic Matter, and a Quark-Gluon Plasma. We suggest that these meet at a triple point in the QCD phase diagram. Some of the features explained, both qualitatively and semi-quantitatively, include the curve for the decoupling of chemical equilibrium, along with the non-monotonic behavior of strange particle multiplicity ratios at center of mass energies near 10 GeV. If the transition(s) between the three phases are merely crossover(s), the triple point is only approximate.