An intrinsic antiferromagnetic topological insulator ${\mathrm{MnBi}}_2{\mathrm{Te}}_4$ arises when intercalating a Mn-Te bilayer chain in a topological insulator, ${\mathrm{Bi}}_2{\mathrm{Te}}_3$. We present observations on the inter- and intralayer phonon modes of the generalized $\mathrm{Mn}{\mathrm{Bi}}_{2n}{\mathrm{Te}}_{3n+1}\left(n=1,2,3,4\right)$ family using cryogenic low-frequency Raman spectroscopy with various polarization configurations. Two peaks at 66 and $112\mathrm{c}{\mathrm{m}}^{\text{–}1}$ show abnormal perturbation in Raman linewidths below magnetic transition temperature due to spin-phonon coupling. In ${\mathrm{MnBi}}_4{\mathrm{Te}}_7$, ${\mathrm{Bi}}_2{\mathrm{Te}}_3$ layers induce Davydov splitting of the $A_{1g}$ mode around $137\mathrm{c}{\mathrm{m}}^{\text{–}1}$ at 5 K. The out-of-plane interlayer force constant estimated using the linear chain model was $(3.98\pm 0.14)\times {10}^{19}\mathrm{N}/{\mathrm{m}}^3$, three times weaker than that of ${\mathrm{Bi}}_2{\mathrm{Te}}_3$. Adding more ${\mathrm{Bi}}_2{\mathrm{Te}}_3$ layers, such as ${\mathrm{MnBi}}_6{\mathrm{Te}}_{10}$ and ${\mathrm{MnBi}}_8{\mathrm{Te}}_{13}$, makes ${\mathrm{Bi}}_2{\mathrm{Te}}_3$ properties more dominant than magnetic properties. Our work experimentally and theoretically discovers the dynamics of phonon modes of $\mathrm{Mn}{\mathrm{Bi}}_{2n}{\mathrm{Te}}_{3n+1}$ family, facilitating utilization of magnetic topological heterostructures.

Dipolar magnon-magnon coupling has long been predicted in nanopatterned artificial spin systems. However, observation of such phenomena and related collective spin-wave signatures have until recently proved elusive or been limited to low-power edge modes which are difficult to measure experimentally. Here we describe the requisite conditions for dipolar mode-hybridization, how it may be controlled, why it was not observed earlier, and how strong coupling may occur between nanomagnet bulk modes. We experimentally investigate four nanopatterned artificial spin system geometries: chevron arrays, square, staircase, and brickwork artificial spin ices. We observe significant dynamic dipolar-coupling in all systems with relative coupling strengths and avoided-crossing gaps supported by micromagnetic-simulation results. We demonstrate reconfigurable mode-hybridization regimes in each system via microstate control, and in doing so elucidate the underlying dynamics governing dynamic dipolar-coupling with implications across reconfigurable magnonics. We demonstrate that confinement of the bulk modes via edge effects plays a critical role in dipolar hybridized modes, and treating each nanoisland as a coherently precessing macro-spin or a standing spin-wave is insufficient to capture experimentally observed coupling phenomena. Finally, we present a parameter-space search detailing how coupling strength may be tuned via nanofabrication dimensions and material properties.

Density functional theory is the standard theory for computing the electronic structure of materials, which is based on a functional that maps the electron density to the energy. However, a rigorous form of the functional is not known and has been heuristically constructed by interpolating asymptotic constraints known for extreme situations, such as isolated atoms and uniform electron gas. Recent studies have demonstrated that the functional can be effectively approximated using machine learning (ML) approaches. However, most ML models do not satisfy asymptotic constraints. In this paper, by applying a ML model architecture, we demonstrate a neural network-based exchange-correlation functional satisfying physical asymptotic constraints. Calculations reveal that the trained functional is applicable to various materials with an accuracy higher than that of existing functionals, even for materials whose electronic properties are not included in training dataset. Our proposed method thus improves the accuracy and generalization performance of the ML-based functional by combining the advantages of ML and analytical modeling.

We describe a new class of nonequilibrium quantum many-body phenomena in the form of networks of caustics that dominate the many-body wave function in the semiclassical regime following a sudden quench. It includes the light cone-like propagation of correlations as a particular case. Caustics are singularities formed by the birth and death of waves and form a hierarchy of universal patterns whose natural mathematical description is via catastrophe theory. Examples in classical waves range from rainbows and gravitational lensing in optics to tidal bores and rogue waves in hydrodynamics. Quantum many-body caustics are discretized by second quantization (“quantum catastrophes”) and live in Fock space, which can potentially have many dimensions. We illustrate these ideas using the Bose Hubbard dimer and trimer models, which are simple enough that the caustic structure can be elucidated from first principles and yet run the full range from integrable to nonintegrable dynamics. The dimer gives rise to discretized versions of fold and cusp catastrophes whereas the trimer allows for higher catastrophes including the codimension-3 hyperbolic and elliptic umbilics, which are organized by, and projections of, an 8-dimensional corank-2 catastrophe known as $X_9$. These results describe a hitherto unrecognized form of universality in quantum dynamics organized by singularities that manifest as strong fluctuations in mode population probabilities.

Ultrafast magnetization dynamics are governed by energy flow between electronic, magnetic, and lattice degrees of freedom. A quantitative understanding of these dynamics must be based on a model that agrees with experimental results for all three subsystems. However, ultrafast dynamics of the lattice remain largely unexplored experimentally. Here we combine femtosecond electron diffraction experiments of the lattice dynamics with energy-conserving atomistic spin dynamics (ASD) simulations and ab initio calculations to study the intrinsic energy flow in the $3d$ ferromagnets cobalt (Co) and iron (Fe). The simulations yield a good description of experimental data, in particular an excellent description of our experimental results for the lattice dynamics. We find that the lattice dynamics are influenced significantly by the magnetization dynamics due to the energy cost of demagnetization. Our results highlight the role of the spin system as the dominant heat sink in the first hundreds of femtoseconds. Together with previous findings for nickel [Zahn et al., Phys. Rev. Research 3, 023032 (2021)], our work demonstrates that energy-conserving ASD simulations provide a general and consistent description of the laser-induced dynamics in all three elemental $3d$ ferromagnets.

To provide a physical example of quantum scars, we study the many-body scars in the spin-1 Kitaev chain where the so-called PXP Hamiltonian is exactly embedded in the spectra. Regarding the conserved quantities, the Hilbert space is fragmented into disconnected subspaces and we explore the associated constrained dynamics. The continuous revivals of the fidelity and the entanglement entropy when the initial state is prepared in $|{\mathbb{Z}}_k\rangle$ ($k=2,3$) state illustrate the essential physics of the PXP model. We study the quantum phase transitions in the one-dimensional spin-1 Kitaev-Heisenberg model using the density-matrix renormalization group and Lanczos exact diagonalization methods, and determine the phase diagram. We parametrize the two terms in the Hamiltonian by the angle $\varphi$, where the Kitaev term is $K\equiv sin\left(\varphi \right)$ and competes with the Heisenberg $J\equiv cos\left(\varphi \right)$ term. One finds a rich ground-state phase diagram as a function of the angle $\varphi$. Depending on the ratio $K/J\equiv tan\left(\varphi \right)$, the system either breaks the symmetry to one of distinct symmetry broken phases, or preserves the symmetry in a quantum spin-liquid phase with frustrated interactions. We find that the scarred state is stable for perturbations which obey ${\mathbb{Z}}_2$-symmetry, while it becomes unstable against Heisenberg-type perturbations.

Dynamic and steady-state aspects of wave propagation are deeply connected in lossless open systems in which the scattering matrix is unitary. There is then an equivalence among the energy excited within the medium through all channels, the Wigner time delay, which is the sum of dwell times in all channels coupled to the medium, and the density of states. But these equivalences fall away in the presence of material loss or gain. In this paper, we use microwave measurements, numerical simulations, and theoretical analysis to discover the changing relationships among fundamental wave properties with loss and gain, and their dependence upon dimensionality and spectral overlap. We begin with the demonstrations that the transmission time in random 1D media is equal to the density of states even in the presence of ultrastrong absorption and that its ensemble average is independent of the strengths of scattering and absorption. In contrast, the Wigner time becomes imaginary in the presence of loss, with real and imaginary parts that fall with absorption. In multichannel media, the transmission time remains equal to the density of states and is independent of the scattering strength in unitary systems but falls with absorption to a degree that increases with the strengths of absorption and scattering, and the number of channels coupled to the medium. We show that the relationships between key propagation variables in non-Hermitian systems can be understood in terms of the singularities of the phase of the determinant of the transmission matrix. The poles of the transmission matrix are the same as those of the scattering matrix, but the transmission zeros are fundamentally different. Whereas the zeros of the scattering matrix are the complex conjugates of the poles, the transmission zeros are topological: in unitary systems they occur only singly on the real axis or as conjugate pairs. We follow the evolution and statistics of zeros in the complex plane as random samples are deformed. The sensitivity of the spacing of zeros in the complex plane with deformation of the sample has a square-root singularity at a zero point at which two single zeros and a complex pair interconvert. The transmission time is a sum of Lorentzian functions associated with poles and zeros. The sum over poles is the density of states with an average that is independent of scattering and dissipation. But the sum over zeros changes with loss, gain, scattering strength and the number of channels in ways that make it possible to control ultranarrow spectral features in transmission and transmission time. We show that the field, including the contribution of the still coherent incident wave, is a sum over modal partial fractions with amplitudes that are independent of loss and gain. The energy excited may be expressed in terms of the resonances of the medium and is equal to the dwell time even in the presence of loss or gain.

We propose an enantiodetection method of chiral molecules in a cavity with an external drive. The chiral molecules are coupled with a quantized cavity field and two classical light fields to form the cyclic three-level systems. The chirality-dependent cavity-assisted three-photon process in the three-level systems leads to the generation of intracavity photons. Simultaneously, the drive field also results in a chirality-independent process of the generation of intracavity photons. Based on the interference between the intracavity photons generated from these two processes, one can detect the enantiomeric excess of chiral mixture via monitoring the transmission rate of the drive field. It is also shown that our method can work in the single-molecule weak-coupling region and in the absence of additional enantiopure samples. These advantages make it offer promising applications in enantiodetection of chiral molecules.

Continuous-variable (CV) quantum key distribution (QKD) employs the quadratures of a bosonic mode to establish a secret key between two remote parties, and this is usually achieved via a Gaussian modulation of coherent states. The resulting secret key rate depends not only on the loss and noise in the communication channel, but also on a series of data processing steps that are needed for transforming shared correlations into a final string of secret bits. Here we consider a Gaussian-modulated coherent-state protocol with homodyne detection in the general setting of composable finite-size security. After simulating the process of quantum communication, the output classical data are postprocessed via procedures of parameter estimation, error correction, and privacy amplification. In particular, we analyze the high signal-to-noise regime, which requires the use of high-rate (nonbinary) low-density parity check codes. We implement all these steps in a python-based library that allows one to investigate and optimize the protocol parameters to be used in practical experimental implementations of short-range CV-QKD.

Shallow nitrogen-vacancy (NV) centers in diamond are promising for nanomagnetometry, for they can be placed proximate to targets. To study the intrinsic magnetic properties, zero-field magnetometry is desirable. However, for shallow NV centers under zero field, the strain near diamond surfaces would cause level anticrossing between the spin states, leading to clock transitions whose frequencies are insensitive to magnetic signals. Furthermore, the charge noises from the surfaces would induce extra spin decoherence and hence reduce the magnetic sensitivity. Here, we demonstrate that the relatively strong hyperfine coupling (130 MHz) from a first-shell ${}^{13}\mathrm{C}$ nuclear spin can provide an effective bias field to an NV center spin so that the clock-transition condition is broken and the charge noises are suppressed. The hyperfine bias enhances the dc magnetic sensitivity by a factor of 22 in our setup. With the charge noises suppressed by the strong hyperfine field, the ac magnetometry under zero field also reaches the limit set by decoherence due to the nuclear spin bath. In addition, the 130 MHz splitting of the NV center spin transitions allows relaxometry of magnetic noises simultaneously at two well-separated frequencies (∼2.870 ± 0.065 GHz), providing (low-resolution) spectral information of high-frequency noises under zero field. The hyperfine-bias-enhanced zero-field magnetometry can be combined with dynamical decoupling to enhance single-molecule magnetic resonance spectroscopy and to improve the frequency resolution in nanoscale magnetic resonance imaging.

Deep quantum neural networks may provide a promising way to achieve a quantum learning advantage with noisy intermediate-scale quantum devices. Here, we use deep quantum feed-forward neural networks capable of universal quantum computation to represent the mixed states for open quantum many-body systems and introduce a variational method with quantum derivatives to solve the master equation for dynamics and stationary states. Owning to the special structure of the quantum networks, this approach enjoys a number of notable features, including an efficient quantum analog of the back-propagation algorithm, resource-saving reuse of hidden qubits, general applicability independent of dimensionality and entanglement properties, as well as the convenient implementation of symmetries. As proof-of-principle demonstrations, we apply this approach to both one-dimensional transverse field Ising and two-dimensional $J_1\text{−}{J}_2$ models with dissipation, and show that it can efficiently capture their dynamics and stationary states with a desired accuracy.

Polynomial systems over the binary field have important applications, especially in symmetric and asymmetric cryptanalysis, multivariate-based postquantum cryptography, coding theory, and computer algebra. In this paper, we study the quantum annealing model for solving Boolean systems of multivariate equations of degree 2, usually referred to as the multivariate quadratic problem. We present different methodologies to embed the problem into a Hamiltonian that can be solved by available quantum annealing platforms. In particular, we provide three embedding options, and we highlight their differences in terms of quantum resources. Moreover, we design a machine-agnostic algorithm that adopts an iterative approach to better solve the problem Hamiltonian by repeatedly reducing the search space. Finally, we use D-Wave devices to successfully implement our methodologies on several instances of the multivariate quadratic problem.

Using first-principles calculations and group-theory-based models, we study the stabilization of ferrielectricity (FiE) in ${\mathrm{CuInP}}_2{\mathrm{Se}}_6$. We find that the FiE ground state is stabilized by a large anharmonic coupling between the polar mode and a fully symmetric Raman-active mode. Our results open possibilities for controlling the single-step switching barrier for polarization by tuning the Raman-active mode. We discuss the implications of our findings in the context of designing next-generation optoelectronic devices that can overcome the voltage-time dilemma.

We investigate the classical chiral Ashkin-Teller model on a square lattice with the corner transfer matrix renormalization group algorithm. We show that the melting of the period-4 phase in the presence of a chiral perturbation takes different forms depending on the coefficient of the four-spin term in the Ashkin-Teller model. Close to the clock limit of two decoupled Ising models, the system undergoes a two-step commensurate-incommensurate transition as soon as the chirality is introduced, with an intermediate critical floating phase bounded by a Kosterlitz-Thouless transition at high temperature and a Pokrovsky-Talapov transition at low temperature. By contrast, close to the four-states Potts model, we argue for the existence of a unique commensurate-incommensurate transition that belongs to the chiral universality class, and for the presence of a Lifshitz point where the ordered, disordered, and floating phases meet. Finally, we map the whole phase diagram, which turns out to be in qualitative agreement with the 40-year-old prediction by Huse and Fisher.

The quantum logic gates used in the design of a quantum computer should be both universal, meaning arbitrary quantum computations can be performed, and fault-tolerant, meaning the gates keep errors from cascading out of control. A number of no-go theorems constrain the ways in which a set of fault-tolerant logic gates can be universal. These theorems are very restrictive, and conventional wisdom holds that a universal fault-tolerant logic gate set cannot be implemented natively, requiring us to use costly distillation procedures for quantum computation. Here, we present a general framework for universal fault-tolerant logic with stabilizer codes, together with a no-go theorem that reveals the very broad conditions constraining such gate sets. Our theorem applies to a wide range of stabilizer code families, including concatenated codes and conventional topological stabilizer codes such as the surface code. The broad applicability of our no-go theorem provides a new perspective on how the constraints on universal fault-tolerant gate sets can be overcome. In particular, we show how nonunitary implementations of logic gates provide a general approach to circumvent the no-go theorem, and we present a rich landscape of constructions for logic gate sets that are both universal and fault-tolerant. That is, rather than restricting what is possible, our no-go theorem provides a signpost to guide us to new, efficient architectures for fault-tolerant quantum computing.