The three-dimensional (3D) quantum XYZ product can construct a class of non-Calderbank-Shor-Stearne (non-CSS) quantum codes by using three classical codes. However, there has been limited study on their error-correcting performance so far and whether this code construction can be generalized to higher dimension is an open question. In this paper, we first study the error-correcting performance of the 3D Chamon code, which is an instance of the 3D XYZ product of three repetition codes. Second, we show that the 3D XYZ product can be generalized to four dimensions and propose four-dimensional (4D) XYZ product code construction, which constructs a class of non-CSS quantum codes by using either four classical codes or two CSS quantum codes. Compared with the 4D homological product, we show that the 4D XYZ product can construct non-CSS codes with higher code dimension or code distance. Third, we consider two instances of the 4D XYZ product, which we refer to as the 4D Chamon code and the 4D XYZ product concatenated code, respectively. Our simulation results show that the 4D XYZ product can construct non-CSS codes with better error-correcting performance against Pauli-$Z$-biased noise than CSS codes constructed by the 4D homological product. Finally, we present the geometric arrangement of the 4D Chamon code within a 4D cubic lattice, demonstrating that it possesses two key characteristics of fracton models, which strongly suggest that it is a novel 4D fracton model.

This study provides a comprehensive analysis of $S$-wave exotic hydrogen-like three-body systems $(pp{\mu }^{-},$ $pp{\tau }^{-},$ ${\mu }^{-}{\mu }^{-}p,$ ${\tau }^{-}{\tau }^{-}p,$ $p{\mu }^{-}{\tau }^{-})$ with spin parity $J^P=1/2^{+}$ and $3/2^{+}$, and four-body systems $(pp{\mu }^{-}{\mu }^{-},$ $pp{\tau }^{-}{\tau }^{-})$ with $J^P=0^{+}$, $1^{+}$, and $2^{+}$. The wave functions are constructed using complete spin configurations and spatial components with various Jacobi coordinate structures. We use complex scaling and Gaussian expansion methods to solve the complex-scaled Schrödinger equation and obtain possible bound and quasibound states. The resulting binding energies range from $-33.8$ keV to $-340\mathrm{eV}$, with the $p{\mu }^{-}{\tau }^{-}$ system being the only case where no bound state is found. Furthermore, we present theoretical estimation of the bound-state energy levels of $pp{\mu }^{-}{\mu }^{-}$ and $pp{\tau }^{-}{\tau }^{-}$, which is of significant importance for understanding exotic few-body Coulomb systems. We also analyze spin configurations and root-mean-square radii to elucidate the spatial structure of these bound and quasibound states. Our results reveal that $K$-type spatial configurations play a crucial role in accurately describing bound and quasibound states in the hydrogen-molecule-like systems $pp{\mu }^{-}{\mu }^{-}$ and $pp{\tau }^{-}{\tau }^{-}$. Incorporating $K$-type configurations significantly alters the mass spectra of these states. Future muon colliders and muon facilities may offer promising platforms for the possible copious production of such heavy flavored hydrogen molecules and molecular ions. For instance, scattering processes such as $2{\mu }^{-}+{\mathrm{H}}_2\to {\mathrm{H}}_{2\mu }+2e^{-}$, ${\mu }^{-}+{\mathrm{H}}_2\to {\mathrm{H}}_{\mu \mathrm{e}}+e^{-}$, and ${\mu }^{-}+{\mathrm{H}}_2^{+}\to {\mathrm{H}}_{2\mu }^{+}+e^{-}$ could be utilized, facilitating detailed studies of intriguing states such as ${\mathrm{H}}_{2\mu }$, ${\mathrm{H}}_{\mu \mathrm{e}}$, and ${\mathrm{H}}_{2\mu }^{+}$.

We experimentally compare nanoplasma ignition in a helium (He) nanodroplet doped with either oxygen molecules $\left({\mathrm{O}}_2\right)$ or xenon atoms (Xe) under intense near-infrared femtosecond laser pulses. At low dopant concentrations, Xe demonstrates higher ignition efficiency than ${\mathrm{O}}_2$. However, as dopant concentration increases, ${\mathrm{O}}_2$ surpasses Xe due to the weaker intermolecular interactions of ${\mathrm{O}}_2$ with He atoms, which reduce He evaporation losses. The yields and relative abundances of ${\mathrm{He}}_n·\mathrm{H}{\mathrm{e}}^{+}$ ions $(n=0, 1, 2)$ exhibit distinct trends for ${\mathrm{O}}_2$ and Xe, reflecting a competition between dopant-assisted ionization and neutral He depletion. At higher laser intensities, ${\mathrm{O}}_2$ exhibits reduced efficiency compared to Xe, underscoring the critical influence of molecular electronic structures on nanoplasma-ignition dynamics. We also examine the role of droplet size in these processes. Our results reveal the potential of molecular dopants like ${\mathrm{O}}_2$ to enable tailored ignition strategies for nanoplasma applications.

We present a Bethe-ansatz-solvable dissipative model of one-dimensional hard-core anyons under dephasing, where the particles obey fractional statistics. We establish an exact mapping between the Liouvillian and a non-Hermitian Hubbard model with special twisted boundary conditions. With the help of string hypothesis, we rigorously calculate the finite-size corrections and the Liouvillian gap in terms of the statistical parameter, total particle number, lattice size, and dissipation strength. The statistical parameter breaks the $\eta$-pairing symmetry of the system, partially lifts out the degeneracy of the spectrum, and contributes to the Liouvillian gap in finite-size corrections. We also derive the general form of the equation of motion for expectation values of the $2w$-point correlator of hard-core anyons and present a proposal for experimentally measuring the correlation function. Our results shed light on the Liouvillians of exactly solvable dissipative models of anyonic particles that obey the fractional statistics.

In traditional research, chaos is frequently accompanied by nonlinearity, which typically stems from nonlinear interactions or external driving forces. However, in this paper, we present the chaotic behavior that is completely attributed to the nonlinear back-reaction of the non-Markovian environment. To be specific, we derive the dynamical equations of an optomechanical system and demonstrate that the nonlinearity (cause of chaos) in the equations arises entirely from the time-domain convolutions (TDCs) induced by non-Markovian corrections. Under Markovian conditions, these TDCs are reduced into constants, thereby losing the nonlinearity and ultimately leading to the disappearance of chaos. Furthermore, we also observe chaos generation in the absence of optomechanical couplings, which further confirms that the non-Markovian effect is the sole inducement of chaos and the environmental parameters play important roles in the generation of chaos. We hope these results may open a new direction to investigate chaotic dynamics purely caused by non-Markovian environments.

The cutoff energy of high-order harmonics generation (HHG) in doped systems is investigated theoretically with the different laser wavelengths. Our study reveals that in the doped system, when driven by lasers with wavelengths ranging from $1.4$ to $6.6\mu \mathrm{m}$, the cutoff energy exhibits three distinct regimes: linear dependence on wavelength, rapid enhancement, and wavelength independence. In the regime where the cutoff energy exhibits linear dependence on the wavelength, interband transitions between localized electron states and the conduction bands dominate the cutoff of HHG in the doped system, with the cutoff energy being significantly higher than that of HHG in the pristine crystal. When the cutoff energy increases rapidly with wavelength, interband transitions involving valence band states and the second and third conduction bands govern the HHG cutoff energy. In the wavelength-independent regime, the intensity of intraband HHG increases rapidly: electrons tend to move within the same band rather than transiting between different bands, causing LES to lose the advantage in modulating interband processes. Consequently, the doped system no longer exhibits the higher HHG cutoff energy relative to the pristine crystal. This work explores the correlation between HHG in doped systems and the laser wavelength, providing a theoretical basis for the modulation of HHG through doping.

We investigate the previously unexplored quantum dynamics of nonrelativistic, spinless particles propagating in curved spaces with torsion. Our findings demonstrate that while torsion has been predominantly associated with spin, it can also influence the quantum behavior of spinless particles by inducing a logarithmic nonlinearity in the Schrödinger equation through quantum fluctuations, even in flat space. To facilitate quantization in curved spaces, we introduce an extended stochastic variational method. Unlike canonical quantization, this approach is naturally suited to general coordinate systems, with quantum fluctuations arising from a noise term in the stochastic process that is directly influenced by torsion. By requiring consistency with quantum dynamics, we derive an upper bound on the magnitude of torsion. Our results reveal a previously unrecognized mechanism by which torsion, as predicted in certain extensions of general relativity, can influence quantum systems, with potential implications for early-universe physics and dark matter or energy models.

We study the structure of rotation sequences composed of piecewise displacements by an angle $\beta$ about arbitrary axes in three dimensions. We find that when $\beta =2\pi /m$, repeated toggling-frame transformations return the sequence to its original form after $m$ applications. We develop this cyclicity property to unify several families of error-tuned control sequences used across quantum technologies, including NMR, MRI, atomic clocks, and atom-scale computing. For the $m=2$ case, we show a duality between narrowband and broadband $\pi$-element sequences. For $m>2$, we relate the cyclic structure to polyhedral symmetry, leading to original sequence designs based on $m$-fold rotational axes. The result may guide the construction of robust control protocols grounded in geometric and group-theoretic principles.

In the field of ab initio polaritonic chemistry, the use of quantum computers may offer, in the long term, a promising computational pathway for simulating and exploring novel cavity-induced light-matter phenomena. Addressing such problems on emerging platforms, however, already presents significant challenges that need to be considered for near-term applications. Within the framework of the variational quantum eigensolver (VQE), a key question concerns the development of quantum circuit Ansätze capable of encoding strongly correlated electron-photon states at a reasonable resource cost. Motivated by this question, in this work we develop different Ansatz strategies covering not only conventional qubit-based approaches but also higher-dimensional information units. Specifically, we introduce a set of circuits tailored to three different quantum computing architectures: (i) qubits, (ii) qudits, and (iii) hybrid qubit-qumode units of information. All circuits are designed to be both compact and physically motivated, with a central element being the development of spin-adapted electron-photon entangling fragment circuits tailored to the native capabilities of each platform. Integrating the different Ansätze within the state-averaged VQE, we numerically benchmark the three strategies on a cavity-embedded ${\mathrm{H}}_2$ molecule using state-vector simulations. Our results show that each proposed Ansatz achieves high accuracy in encoding the ground, first, and second polaritonic eigenstates of the system, even reproducing characteristic light-matter phenomena such as light-induced avoided crossings and conical intersections. This prospective work, beyond introducing multiple computational strategies adapted to different platforms, also highlights the potential of higher-dimensional quantum architectures to address hybrid fermion-boson systems.

An experimental setup for accurately determining the frequency of alkali-metal-atom two-photon transitions is proposed and demonstrated, with which we avoid laser modulation by an offset-locking technique; keep away collision shift via a high-vacuum pumping; compensate for Zeeman shift directly from the information of the spectral linewidth; remove residual Doppler background by the two-photon crossover spectra; observe the frequency reproducibility between normal days and aftershocks; and compare different approaches of counting laser frequencies, that is, different time bases from either a cesium atomic beam clock or a precalibrated rubidium clock. The frequency of cesium $6S_{1/2}\to 6D_{3/2}$ hyperfine transitions and the corresponding hyperfine constants are hence updated, which are both significant for improving the calculation precision of atomic parity nonconservation and for enhancing the accuracy of the two-photon transition-based optical clocks.

Photoelectron interferometry in strong-field photoionization serves as a powerful tool for probing atomic and molecular structure, and ultrafast electron dynamics. In the photoelectron momentum spectrum, different types of interference patterns are intertwined, which complicates the extraction of information from specific interference patterns. Enhancing the desired interference features is thus essential for their applications. Here, we propose a scheme to control the interference in the photoelectron momentum spectrum through ultraviolet (XUV)-assisted photoionization. Our scheme employs an XUV pulse to govern the initial ionization or excitation process, and an infrared (IR) field to steer the subsequent electron dynamics, offering more degrees of freedom to tailor photoelectron interference. By numerically solving the three-dimensional time-dependent Schrödinger equation, we demonstrate that the intracycle interference and the near-forward rescattering photoelectron holographic interference are selectively enhanced or suppressed by tuning the XUV-IR time delay and the XUV pulse parameters. Our findings pave the way for attosecond-scale tracking of electron motion via photoelectron interferometry.

We analyze the many-body Hamiltonian describing a mobile impurity immersed in a Bose-Einstein condensate (BEC). Using exact unitary transformations and rigorous error estimates, we show the validity of the Bogoliubov-Fröhlich Hamiltonian for the Bose polaron in the regime of moderately strong, repulsive interactions with a dilute BEC. Moreover, we calculate analytically the universal logarithmic correction to the ground state energy that arises from an impurity mediated phonon-phonon interaction.

We numerically study the dynamical properties of a mixture consisting of a dipolar condensate and a degenerate Fermi gas in a quasi-one-dimensional geometry. In particular, we focus on the system's response to a temporal variation in the interaction strength between bosons and fermions. When the interspecies attraction becomes sufficiently strong, we observe a phase transition to a supersolid state. This conclusion is supported by the emergence of an out-of-phase Goldstone mode in the excitation spectrum.

We consider excitations of a spin-1 Bose-Einstein-condensate in the vicinity of different mean-field configurations and derive mappings to emergent relativistic quantum field theories minimally coupled to curved acoustic spacetimes. The quantum fields are typically identified with Nambu-Goldstone bosons, such that the structure of the analog quantum field theories on curved spacetimes depends on the (spontaneous) symmetry breaking pattern of the respective ground state. The emergent spacetime geometries are independent of each other and exhibit bimetricity in the polar and antiferromagnetic phase, whereas one has trimetricity in the ferromagnetic phase. Compared to scalar BECs, the spinor degrees of freedom allow us to investigate massive vector and scalar fields where the former is a spin-nematic rotation mode in the polar phase, which can be cast into a Proca field that is minimally coupled to a curved spacetime that emerges on length scales larger than the spin-healing length. Finally, we specify the Zeeman couplings and the condensate trap to be spacetime-dependent such that a cosmological Friedmann-Lemaître-Robertson-Walker metric can be achieved. This work enables a pathway towards quantum-simulating cosmological particle production of Proca quanta via quenching the quadratic Zeeman coefficient or via magnetic field ramps, which both result in the creation of spin-nematic squeezed states.

Recently discovered antichiral edge states, which are pinned to a pair of offset Dirac points with different degenerate energies, exhibit the same propagating direction along parallel boundaries. However, few studies focus on accompanying opposite bulk propagation in terms of energy conservation in nondissipative systems. Here, we fabricate a two-dimensional honeycomb-lattice antiferromagnetic photonic crystal to construct a photonic Dirac semimetal. Both nonreciprocal bulk states and antichiral edge states are observed in the vicinity of the Dirac point frequency in microwave experiments. Using a rhombus-shaped sample, we realize one-way bulk to boundary and boundary to bulk couplings. Furthermore, we demonstrate a nonreciprocal tunnelinglike phenomenon from the upper to the lower boundaries across the bulk region. Our results highlight the underlying relation between the nonreciprocal bulk and antichiral edge states, which may provide an ideal strategy to design perfect bulk-boundary couplers, splitters, and various other applications.

Quantum correlations in the spatial domain hold great promise for applications in quantum imaging, quantum cryptography and quantum information processing, owing to the infinite dimensionality of the associated Hilbert space. Here, we present a theoretical investigation, complemented by experimental measurements, of the propagation dynamics of the spatial quantum correlations in bright structured twin beams generated via a four-wave mixing process in a double-$\mathrm{\Lambda }$ configuration in atomic vapor. We derive an analytical expression describing the evolution of the spatial quantum correlation distribution from the near field to the far field. To qualitatively support the theoretical predictions, we perform experiments measuring intensity-difference noise between different spatial subregions of the twin beams as they propagate from the near field to the far field. The presence of quantum correlations is manifested as squeezing in the intensity-difference noise measurement. With a Gaussian pump, we observe localized correlations in the near field and localized anticorrelations in the far field. In contrast, with a structured Laguerre-Gaussian pump, there is a transition from localized correlations in the near field to delocalized correlations in the far field. The present results offer valuable insights into the fundamental behavior of spatial quantum correlations and open possibilities for potential applications in quantum information, quantum imaging, and sensing.

We investigate the steady-state dynamical regimes of boundary-driven, dissipative bosonic chains subjected to $n$-photon drives. Using the truncated Wigner approximation, we explore how multiphoton drives shape the interplay between quantum fluctuations, nonlinear interactions, and dissipative processes in such quantum systems. We identify two main regimes: a chaotic hydrodynamic regime characterized by the restoration of a local $\text{U}\left(1\right)$ symmetry, photon saturation due to Kerr nonlinearity, and spatial prethermalization effects and a nonchaotic resonant nonlinear wave (RNW) regime exhibiting sub-Poissonian photon statistics, persistent ${\mathbb{Z}}_n$ symmetry, and quantum-driven phase decoherence. Our findings reveal the universal nature of the hydrodynamic regime and highlight the RNW regime's sensitivity to boundary driving conditions, suggesting novel routes for quantum state engineering in driven-dissipative quantum devices. These results are experimentally relevant for state-of-the-art circuit quantum electrodynamics platforms.

Conservation laws are of fundamental importance in science and engineering, and so are ubiquitous in applications. For continuum systems such as classical fluids and plasmas, they give rise to coupled systems of nonlinear partial differential equations (PDEs) whose solutions can develop discontinuities such as shocks. Finite-volume (FV) methods are an important tool for finding approximate solutions to such PDEs as they are based on an integral formulation of the conservation laws that can handle the presence of discontinuities. Weighted essentially nonoscillatory (WENO) methods are a class of high-resolution methods that can (i) produce solutions of arbitrarily high-order accuracy in regions where the solution is smooth, while (ii) simultaneously quashing the development of nonphysical oscillations in the vicinity of discontinuities that often cause numerical simulations to fail. In this paper we show how FV and WENO methods can be incorporated into a quantum algorithm for solving systems of nonlinear PDEs (QPDE) introduced by one of the authors. We numerically simulate application of the modified QPDE algorithm to three verification problems. The first has a smooth solution, while the other two are well-known one-dimensional Riemann problems whose exact solutions contain moving shocks, contact discontinuities, and rarefaction waves. For all three problems excellent agreement is found between the results of the QPDE simulations and the exact solutions.

Majorana zero mode (MZM) exhibits inherent resilience to local parametric fluctuations, due to the topological protection mechanism in the non-Abelian braiding statistics of the anyonic quasiparticles. In this paper, we construct the braiding operations between an arbitrary pair in three MZMs under the theoretical framework of universal quantum control. Largely detuned driving fields on the mediator, a local defect of lattice, enable indirect and tunable exchange interaction between arbitrary two MZMs. The braiding operations can then emerge through the time evolution along the universal nonadiabatic passages, whose robustness against driving-field errors and quasiparticle poisoning can be substantially enhanced by the rapid modulation over the passage-dependent global phase. Moreover, a chiral transfer for population on MZMs along the universal passage can be perfectly demonstrated in both clockwise and counterclockwise manners without eliminating the mediator defect. Our protocol is linear scalable and provides an avenue towards the universal quantum control over MZMs, which is fundamental and essential for topological quantum computation.