We observe a wealth of multimagnon bound states in the strongly anisotropic spin-1 triangular antiferromagnet ${\mathrm{FeI}}_2$ using time-domain terahertz spectroscopy. These unconventional excitations can arise in ordered magnets due to attractive magnon-magnon interactions and alter their properties. We analyze the energy-magnetic field spectrum via an exact diagonalization method for a dilute gas of interacting magnons and detect up to 4- and 6-magnon bound states, hybridized with single magnons. This zoo of tunable interacting quasiparticles provides a unique platform to study decay and renormalization, reminiscent of the few-body problems found in cold-atom, nuclear, and particle physics.

The low-temperature properties of interacting quantum systems are believed to require exponential resources to compute in the general case. Quantifying the extent to which such properties can be approximated using efficient algorithms remains a significant open challenge. Here, we consider the task of approximating the ground state energy of two-local quantum Hamiltonians with bounded-degree interaction graphs. Most existing algorithms optimize the energy over the set of product states. We propose and analyze a family of shallow quantum circuits that can be used to improve the approximation ratio achieved by a given product state. The algorithm takes as input an $n$-qubit product state with variance Var and improves its energy by an amount proportional to ${\mathrm{Var}}^2/n$. In a typical case, this results in an extensive improvement in the estimated energy. We extend our results to $k$-local Hamiltonians and entangled initial states.

The Gibbs state is widely taken to be the equilibrium state of a system in contact with an environment at temperature $T$. However, non-negligible interactions between system and environment can give rise to an altered state. Here, we derive general expressions for this mean force Gibbs state, valid for any system that interacts with a bosonic reservoir. First, we derive the state in the weak coupling limit and find that, in general, it maintains coherences with respect to the bare system Hamiltonian. Second, we develop a new expansion method suited to investigate the ultrastrong coupling regime. This allows us to derive the explicit form for the mean force Gibbs state, and we find that it becomes diagonal in the basis set by the system-reservoir interaction instead of the system Hamiltonian. Several examples are discussed including a single qubit, a three-level $V$-system, and two coupled qubits all interacting with bosonic reservoirs. The results shed light on the presence of coherences in the strong coupling regime, and provide key tools for nanoscale thermodynamics investigations.

Systems evolving through aggregation and fragmentation may possess an intriguing supercluster state (SCS). Clusters constituting this state are mostly very large, so the SCS resembles a gelling state, but the formation of the SCS is controlled by fluctuations and in this aspect, it is similar to a critical state. The SCS is nonextensive, that is, the number of clusters varies sublinearly with the system size. In the parameter space, the SCS separates equilibrium and jamming (extensive) states. The conventional methods, such as, e.g., the van Kampen expansion, fail to describe the SCS. To characterize the SCS we propose a scaling approach with a set of critical exponents. Our theoretical findings are in good agreement with numerical results.

We study the elliptic double-box integral, which contributes to generic massless QFTs and is the only contribution to a particular 10-point scattering amplitude in $\mathcal{N}=4$ SYM theory. Based on a Feynman parametrization, we express this integral in terms of elliptic polylogarithms. We then study its symbol, finding a rich structure and remarkable similarity with the nonelliptic case. In particular, the first entry of the symbol is expressible in terms of logarithms of dual-conformal cross ratios, and elliptic letters only occur in the last two entries.

We propose a world sheet formula for tree-level correlation functions describing a scalar field with arbitrary mass and quartic self-interaction in de Sitter space, which is a simple model for inflationary cosmology. The correlation functions are located on the future boundary of the spacetime and are Fourier-transformed to momentum space. Our formula is supported on mass-deformed scattering equations involving conformal generators in momentum space and reduces to the CHY formula for ${\varphi }^4$ amplitudes in the flat space limit. Using the global residue theorem, we verify that it reproduces the Witten diagram expansion at four and six points, and sketch the extension to $n$ points.

The strong suppression of high-$p_T$ jets in heavy-ion collisions is a result of elastic and inelastic energy loss suffered by the jet multiprong collection of color charges that are resolved by medium interactions. Hence, quenching effects depend on the fluctuations of the jet substructure that are probed by the cone-size dependence of the spectrum. In this Letter, we present the first complete, analytic calculation of the inclusive $R$-dependent jet spectrum in PbPb collisions at LHC energies, including resummation of energy loss effects from hard, vacuumlike emissions occurring in the medium and modeling of soft energy flow and recovery at the jet cone. Both the geometry of the collision and the local medium properties, such as the temperature and fluid velocity, are given by a hydrodynamic evolution of the medium, leaving only the coupling constant in the medium as a free parameter. The calculation yields a good description of the centrality and $p_T$ dependence of jet suppression for $R=0.4$ together with a mild cone-size dependence, which is in agreement with recent experimental results. Gauging the theoretical uncertainties, we find that the largest sensitivity resides in the leading logarithmic approximation of the phase space of resolved splittings, which can be improved systematically, while nonperturbative modeling of the soft-gluon sector is of relatively minor importance up to large cone sizes.

Determination of nuclear moments for many nuclei relies on the computation of hyperfine constants, with theoretical uncertainties directly affecting the resulting uncertainties of the nuclear moments. In this work, we improve the precision of such a method by including for the first time an iterative solution of equations for the core triple cluster amplitudes into the relativistic coupled-cluster method, with large-scale complete basis sets. We carried out calculations of the energies and magnetic dipole and electric quadrupole hyperfine structure constants for the low-lying states of ${{}^{229}\mathrm{Th}}^{3+}$ in the framework of such a relativistic coupled-cluster single double triple method. We present a detailed study of various corrections to all calculated properties. Using the theory results and experimental data, we found the nuclear magnetic dipole and electric quadrupole moments to be $\mu =0.366(6){\mu }_N$ and $Q=3.11(2)e\mathrm{b}$, respectively, and reduce the uncertainty of the quadrupole moment by a factor of 3. The Bohr-Weisskopf effect of the finite nuclear magnetization is investigated, with bounds placed on the deviation of the magnetization distribution from the uniform one.

Because of quantum noise fluctuations, the rate of error achievable in decision problems involving several possible configurations of a scattering system is subject to a fundamental limit known as the Helstrom bound. Here, we present a general framework to calculate and minimize this bound using coherent probe fields with tailored spatial distributions. As an example, we experimentally study a target located in between two disordered scattering media. We first show that the optimal field distribution can be directly identified using a general approach based on scattering matrix measurements. We then demonstrate that this optimal light field successfully probes the presence of the target with a number of photons that is reduced by more than 2 orders of magnitude as compared to unoptimized fields.

Topological properties of physical systems play a crucial role in our understanding of nature, yet their experimental determination remains elusive. We show that the mean helicity, a dynamical invariant in ideal flows, quantitatively affects trajectories of fluid elements: the linking number of Lagrangian trajectories depends on the mean helicity. Thus, a global topological invariant and a topological number of fluid trajectories become related, and we provide an empirical expression linking them. The relation shows the existence of long-term memory in the trajectories: the links can be made of the trajectory up to a given time, with particles positions in the past. This property also allows experimental measurements of mean helicity.

The discovery of topologically protected boundary states in topological insulators opens a new avenue toward exploring novel transport phenomena. The one-way feature of boundary states against disorders and impurities prospects great potential in applications of electronic and classical wave devices. Particularly, for the 3D higher-order topological insulators, it can host hinge states, which allow the energy to transport along the hinge channels. However, the hinge states have only been observed along a single hinge, and a natural question arises: whether the hinge states can exist simultaneously on all the three independent directions of one sample? Here we theoretically predict, numerically simulate, and experimentally observe the hinge states on three different directions of a higher-order topological phononic crystal, and demonstrate their robust one-way transport from hinge to hinge. Therefore, 3D topological hinge transport is successfully achieved. The novel sound transport may serve as the basis for acoustic devices of unconventional functions.

The manipulation of topological states in quantum matter is an essential pursuit of fundamental physics and next-generation quantum technology. Here we report the magnetic manipulation of Weyl fermions in the kagome spin-orbit semimetal ${\mathrm{Co}}_3{\text{Sn}}_2{\mathrm{S}}_2$, observed by high-resolution photoemission spectroscopy. We demonstrate the exchange collapse of spin-orbit-gapped ferromagnetic Weyl loops into paramagnetic Dirac loops under suppression of the magnetic order. We further observe that topological Fermi arcs disappear in the paramagnetic phase, suggesting the annihilation of exchange-split Weyl points. Our findings indicate that magnetic exchange collapse naturally drives Weyl fermion annihilation, opening new opportunities for engineering topology under correlated order parameters.

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, characterized by Cooper pairs condensed at finite momentum, has been a long-sought state that remains unresolved in many classes of fermionic systems, including superconductors and ultracold atoms. A fascinating aspect of the FFLO state is the emergence of periodic nodal planes in real space, but its observation is still lacking. Here we investigate the superconducting order parameter at high magnetic fields $\mathit{H}$ applied perpendicular to the $ab$ plane in a high-purity single crystal of FeSe. The heat capacity and magnetic torque provide thermodynamic evidence for a distinct superconducting phase at the low-temperature/high-field corner of the phase diagram. Despite the bulk superconductivity, spectroscopic-imaging scanning tunneling microscopy performed on the same crystal demonstrates that the order parameter vanishes at the surface upon entering the high-field phase. These results provide the first demonstration of a pinned planar node perpendicular to $\mathit{H}$, which is consistent with a putative FFLO state.

We argue that a superconducting state with a Fermi surface of Bogoliubov quasiparticles, a Bogoliubov Fermi surface (BG-FS), can be identified by the dependence of physical quantities on disorder. In particular, we show that a linear dependence of the residual density of states at weak disorder distinguishes a BG-FS state from other nodal superconducting states. We further demonstrate the stability of supercurrent against impurities and a characteristic Drude-like behavior of the optical conductivity. Our results can be directly applied to electron irradiation experiments on candidate materials of BG-FSs, including ${\mathrm{Sr}}_2{\mathrm{RuO}}_4$, ${\mathrm{FeSe}}_{1-x}{\mathrm{S}}_x$, and ${\mathrm{UBe}}_{13}$.

We present micromagnetic simulations on resonant spin wave modes of magnetic Hopfions up to 15 GHz driven by external magnetic fields. A sharp transition is found around 66 mT coinciding with a transition from Hopfions to magnetic torons. The modes exhibit characteristic amplitudes in frequency space accompanied by unique localization patterns in real space and are found to be robust to damping around topological features, particularly vortex lines in Hopfions and Bloch points in torons. The marked differences in spin wave spectra between Hopfions, torons, and target skyrmions can serve as fingerprints in future experimental validation studies of these novel 3D topological spin textures.

By using constrained density functional theory modeling, we demonstrate that ultrafast optical pumping unveils hidden charge orders in group VI monolayer transition metal ditellurides. We show that irradiation of the insulating $2H$ phases stabilizes multiple transient charge density wave orders with light-tunable distortion, periodicity, electronic structure, and band gap. Moreover, optical pumping of the semimetallic $1T^{\prime }$ phases generates a transient charge ordered metallic phase composed of 2D diamond clusters. For each transient phase we identify the critical fluence at which it is observed and the specific optical and Raman fingerprints to directly compare with future ultrafast pump-probe experiments. Our work demonstrates that it is possible to stabilize charge density waves even in insulating 2D transition metal dichalcogenides by ultrafast irradiation.

Nature’s most powerful high-energy sources are capable of accelerating particles to high energy and radiating it away on extremely short timescales, even shorter than the light crossing time of the system. It is yet unclear what physical processes can produce such an efficient acceleration, despite the copious radiative losses. By means of radiative particle-in-cell simulations, we show that magnetically dominated turbulence in pair plasmas subject to strong synchrotron cooling generates a nonthermal particle spectrum with a hard power-law range (slope $p\sim 1$) within a few eddy turnover times. Low pitch-angle particles can significantly exceed the nominal radiation-reaction limit, before abruptly cooling down. The particle spectrum becomes even harder ($p<1$) over time owing to particle cooling with an energy-dependent pitch-angle anisotropy. The resulting synchrotron spectrum is hard ($\nu F_{\nu }\propto {\nu }^s$ with $s\sim 1$). Our findings have important implications for understanding the nonthermal emission from high-energy astrophysical sources, most notably the prompt phase of gamma-ray bursts and gamma-ray flares from the Crab nebula.