We explore the recently introduced $\eta$-ensemble approach to compute the free energy directly from ab initio path integral Monte Carlo (PIMC) simulations [T. Dornheim et al., Phys. Rev. B 111, L041114 (2025)] and apply it to the archetypal uniform electron gas model both in the warm dense matter and strongly coupled regimes. Specifically, we present an in-depth study of the relevant algorithmic details such as the choice of the free weighting parameter and the choice of the optimum number of intermediate $\eta$-steps to connect the real, nonideal system ($\eta =1$) with the ideal limit ($\eta =0$). Moreover, we explore the inherent decomposition of the full free energy into its ideal bosonic, ideal-to-interacting, and bosonic-to-fermionic contributions for different parameter regimes. Finally, we compare our new free energy data with an existing free energy parametrization [Groth et al., Phys. Rev. Lett. 119, 135001 (2017)] obtained via adiabatic connection formula evaluations, and we find very good agreement in its range of applicability, i.e., for density parameters $r_s\le 20$; in addition, we present the first PIMC results for the free energy in the low density regime of $20<r_s\le 100$. We expect our results to be of interest both for the study of matter under extreme conditions, as well as to the more general field of PIMC simulations of interacting quantum many-body systems.

Space-time structured laser pulses feature an intensity peak that can travel at an arbitrary velocity while maintaining a near-constant profile. These pulses can propagate in uniform media, where their frequencies are correlated with continuous transverse wave vectors, or in structured media, such as a waveguide, where their frequencies are correlated with discrete mode numbers. Here, we demonstrate the formation and propagation of arbitrary-velocity laser pulses in a plasma waveguide where the intensity can be orders of magnitude higher than in a solid-state waveguide. The flexibility to control the velocity of the peak intensity in a plasma waveguide enables new configurations for plasma-based sources of radiation and energetic particles, including THz generation, laser wakefield acceleration, and direct laser acceleration.

We simulate interaction quenches crossing from a superfluid to a supersolid state in a dipolar quantum gas of ${}^{164}\mathrm{Dy}$ atoms, trapped in an elongated tube with periodic boundary conditions, via the extended Gross-Pitaevskii equation. A freeze-out time is observed through a delay in supersolid formation after crossing the critical point. We compute the density-density correlations at the freeze-out time and extract the frozen correlation length for the solid order. An analysis of the freeze-out time and correlation length versus the interaction quench rate allows us to extract universal exponents corresponding to the relaxation time and correlation length based on predictions of the Kibble-Zurek mechanism. Over several orders of magnitude, clear power-law scaling is observed for both the freeze-out time and the correlation length, and the corresponding exponents are compatible with predictions based on the excitation spectrum calculated via Bogoliubov theory. Defects due to independent local breaking of translational symmetry, contributing to globally incommensurate supersolid order, are identified, and their number at the freeze-out time is found to also scale as a power law. Our results support the hypothesis of a continuous transition whose universality class remains to be determined but appears to differ from that of the (1 + 1)-dimensional XY model.

We establish a toolbox for studying and applying spin-adapted generalized Pauli constraints (GPCs) in few-electron quantum systems. By exploiting the spin symmetry of realistic $N$-electron wave functions, the underlying one-body pure $N$-representability problem simplifies, allowing us to calculate the GPCs for larger system sizes than previously accessible. We then uncover and rigorously prove a superselection rule that highlights the significance of GPCs: whenever a spin-adapted GPC is (approximately) saturated—referred to as (quasi)pinning—the corresponding $N$-electron wave function assumes a simplified structure. Specifically, in a configuration interaction expansion based on natural orbitals only very specific spin configuration state functions may contribute. To assess the nontriviality of (quasi)pinning, we introduce a geometric measure that contrasts it with the (quasi)pinning induced by simple (spin-adapted) Pauli constraints. Applications to few-electron systems suggest that previously observed quasipinning largely stems from spin symmetries.

We present a novel protocol for implementing quantum gates between distant atomic qubits connected by an array of neutral atoms that play the role of a quantum bus. The protocol is based on adiabatically transferring the atoms in the array to an antiferromagnetic-like state of Rydberg excitations using chirped laser pulses. Upon exciting and deexciting the atoms in the array under the blockage of nearest neighbors, depending on the state of the two qubits, the system acquires a conditional geometric $\pi$-phase, while the dynamical phase cancels exactly, even when the atomic positions are disordered but nearly frozen in time, which requires sufficiently low temperatures. With experimentally relevant parameters, using smooth pulses minimizing the Rydberg-state decay and nonadiabatic errors, we obtain the gate times of 2–$3\mu \mathrm{s}$ and gate fidelities of 0.99–0.98 for a pair of atoms separated by $L=20$–$30\mu \mathrm{m}$ and connected by a quantum bus of several (three to six) atoms. Optimizing the pulses, we can obtain faster gates exhibiting even better fidelities than those with smooth adiabatic pulses.

Interacting fermionic chains exhibit extended regions of topological degeneracy of their ground states as a result of the presence of Majorana or parafermionic zero modes localized at the edges. In the opposite limit of infinite temperature, the corresponding nonintegrable spin chains, obtained via generalized Jordan-Wigner mapping, are known to host so-called almost strong zero modes, which are long-lived with respect to any bulk excitations. Here we study the fairly unexplored territory that bridges these two extreme cases of zero and infinite temperature. We blend two established techniques for states, the Lanczos series expansion and a tensor network ansatz, uplifting them to the level of operator algebra. This allows us to efficiently simulate large system sizes for arbitrarily long timescales and to extract the temperature-dependent decay rates. We observe that for the Kitaev-Hubbard model, the decay rate of the edge mode depends exponentially on the inverse temperature $\beta$, and on an effective energy scale ${\mathrm{\Delta }}_{\mathrm{eff}}$ that is greater than the thermodynamic gap of the system $\mathrm{\Delta }$.

We explore how the spectral phase of attosecond pulse trains influences the optical cross section in transient absorption (TA) spectroscopy. The interaction of extreme ultraviolet (XUV) and time-delayed near-infrared (NIR) fields with an atomic or molecular system governs the dynamics. As already shown in RABBITT experiments (Reconstruction of Attosecond Beating by Interference of Two-Photon Transitions), the spectral phase of the XUV pulses can be extracted from the photoionization spectrum as a function of the time delay. Similarly, this XUV phase imprints itself on delay-dependent optical cross-section oscillations. With a perturbative analytical approach and by simulating the quantum dynamics both in a few-level model and via solving the time-dependent Schrödinger equation for atomic hydrogen, we reveal the similarity between the spectral phase in RABBITT and TA spectroscopy.

This study provides a detailed theoretical description of periodic $\gamma$-ray intensity signals. These signals are generated by the tabletop time coincidence experimental setups which utilize nuclear resonant interaction of $\gamma$-photons from a radioactive source with an absorber exhibiting ultrasonic vibrations. Transmission of the $\gamma$-radiation field through the absorber is described by semiclassical theory considering arbitrary periodic motion of the absorber. A key aspect of the presented work is the application of the Fourier analysis. A fundamental relation between the Fourier coefficients of the $\gamma$-ray intensity signal and the Fourier coefficients of a complex exponential of the absorber motion is derived. Application of this relation is shown on the example of the trapezoid motion, which demonstrates its usability for systematic inspection of the discussed signals. The results of the analysis expand the knowledge of the coherent manipulation of $\gamma$-photons and may find usage in a broad range of research fields utilizing $\gamma$-ray sources.

Quantum magnonics based on yttrium iron garnet (YIG) spheres provides a new arena for observing macroscopic quantum states. Here, we propose to prepare two kinds of non-Gaussian magnonic states by adding a single magnon to two Gaussian states, namely, coherent and thermal states. We adopt an optomagnonic system of a YIG sphere and use fast optical pulses to weakly activate the magnon-induced Stokes scattering. Conditioned on the detection of a polarized single photon, a single magnon can be added to an initial Gaussian state. We use a microwave cavity to prepare an initial magnon coherent state and finally read out the generated single-magnon added coherent or thermal state. Both the non-Gaussian states of a magnon mode in a large YIG sphere are macroscopic quantum states, which exhibit many nonclassical properties, such as sub-Poissonian statistics, quadrature squeezing, and a negative Wigner function. Further, both states show a smooth transition from a quantum state to a classical state with the change of the coherent or thermal magnon excitations and thus can be used for the fundamental study of the quantum-to-classical transition.

The advent of memristive devices offers a promising avenue for efficient and scalable analog computing, particularly for linear algebra operations essential in various scientific and engineering applications. This paper investigates the potential of memristive crossbars in implementing matrix inversion algorithms. We explore both static and dynamic approaches, emphasizing the advantages of analog and in-memory computing for matrix operations beyond multiplication. In particular, we demonstrate that the electrical properties of memristive crossbars uniquely suit them for the evolution of a family of matrix exponentials, which can be exploited for the efficient computation of matrix inverses and online solutions for linear problems. Our results demonstrate that memristive arrays can reduce computational complexity. We also study power consumption and show a tradeoff between precision and energy. Furthermore, we address the challenges of device variability, precision, and scalability, providing insights into the practical implementation of these algorithms.

The quantum approximate optimization algorithm (QAOA) is a hybrid quantum-classical algorithm that seeks to achieve approximate solutions to optimization problems by iteratively alternating between intervals of controlled quantum evolution. Here, we examine the effect of analog precision errors on QAOA performance from the perspective of both algorithmic training and performance guarantees. Leveraging cumulant expansions, we recast the faulty QAOA as a control problem in which precision errors are expressed as multiplicative control noise and derive bounds on the performance of QAOA. We show using both analytical techniques and numerical simulations that fixed precision implementations of QAOA circuits are subject to an exponential degradation in performance dependent upon the number of optimal QAOA layers and magnitude of the precision error. Despite this significant reduction, we show that it is possible to mitigate precision errors in QAOA via digitization of the variational parameters at the cost of increasing circuit depth.

When the strong-field ionization of atoms is induced by an ultrashort pulse instead of a long pulse with a flat envelope, many symmetry properties of the photoelecton momentum distribution are broken. The induced asymmetry is measured using the asymmetry parameter which depends on the values of the driving-pulse parameters and the type of the target. We investigate the driving pulses with two carrier frequencies because in this case the dependence of the asymmetry on the characteristics of the target is more robust. Particular attention is devoted to the pulse which consists of two circularly polarized few-cycle waves and the pulse which has two linearly polarized components with mutually orthogonal polarizations. In the former case, we show that the asymmetry parameter is highly sensitive to the ionization potential and to the structure of the ground state. This is particularly the case for the photoelectron energy just above the value for which the contribution of the electrons which do not interact with the core after the ionization becomes negligible. We explain this sensitivity by investigating the dependence of the short-travel-time saddle-point solutions on the characteristics of the target. On the other hand, for the driving pulse with linearly polarized components, the dependence of the asymmetry parameter on the ionization potential is significant, while the dependence on the structure of the ground state is relatively small. In conclusion, we show that the characteristics of the target are imprinted in the asymmetry parameter and this signature is more pronounced for two-component pulses than for the linearly polarized driving pulse with one carrier frequency.

Quasiperiodic potentials and dipolar interactions each impose long-range order in quantum systems, but their interplay unlocks a rich landscape of unexplored quantum phases. In this work, we investigate how dipolar bosonic crystals respond to correlated disorder in the form of quasiperiodic potentials. Using exact numerical simulations and a suite of observables—including order parameters, energy, density distributions, and two-body coherence measures—we explore one-dimensional dipolar bosons in quasiperiodic lattices at both commensurate and incommensurate fillings. Our results reveal a complex competition between superfluid, Mott-insulator, density-wave, and crystalline phases, governed by the intricate balance of dipolar interactions, kinetic energy, and disorder strength. Crucially, we identify mechanisms that influence dipolar crystals, showing their surprising robustness even in the presence of strong quasiperiodic disorder. Strikingly, we challenge previous claims by demonstrating that a kinetic crystal phase—expected to precede full crystallization—does not emerge in the ground state. Instead, its traits appear only under moderate disorder, but never fully develop, giving way to a direct transition from a charge density wave to a crystal state. These findings provide new insights into the resilience of many-body quantum phases in complex environments and pave the way for engineering exotic quantum states in ultracold atomic systems.

The evolution of a quantum system interacting with an environment can be described as a unitary process acting on both the system and the environment. In this framework, the system's evolution can be predicted by tracing out the environmental degrees of freedom. Here, we establish a precise mapping between the global unitary dynamics and the quantum operation involving the system, wherein the system is a single qubit, and the environment is modeled as a discrete lattice space. This approach enables the implementation of arbitrary noise operations on single-polarization qubits using a minimal set of three liquid-crystal metasurfaces, whose transverse distribution of the optic axes can be patterned to reproduce the target process. We experimentally validate this method by simulating common noise processes, such as phase errors and depolarization. Besides providing a practical solution for quantum state purification, this work demonstrates a versatile approach for the simulation of open qubit dynamics, with potential implications for quantum error correction and environment-induced quantum phase transitions.