Boson Hamiltonian Research Articles

Efficient construction of matrix-product representations of many-body Gaussian states

We present a procedure to construct tensor-network representations of many-body Gaussian states efficiently and with a controllable error. These states include the ground and thermal states of bosonic and fermionic quadratic Hamiltonians, which are essential in the study of quantum many-body systems. The procedure improves computational time requirements for constructing many-body Gaussian states by up to five orders of magnitude for reasonable parameter values, thus allowing simulations beyond the range of what was hitherto feasible. Our procedure combines ideas from the theory of Gaussian quantum information with tensor-network based numerical methods thereby opening the possibility of exploiting the rich tool-kit of Gaussian methods in tensor-network simulations.

Simple approach to characterizing band topology in bosonic pairing Hamiltonians

We revisit the problem of characterizing band topology in dynamically-stable quadratic bosonic Hamiltonians that do not conserve particle number. We show this problem can be rigorously addressed by a smooth and local adiabatic mapping procedure to a particle number conserving Hamiltonian. In contrast to a generic fermionic pairing Hamiltonian, such a mapping can always be constructed for bosons. Our approach shows that particle non-conserving bosonic Hamiltonians can be classified using known approaches for fermionic models. It also provides a simple means for identifying and calculating appropriate topological invariants. We also explicitly study dynamically stable but non-positive definite Hamiltonians (as arise frequently in driven photonic systems). We show that in this case, each band gap is characterized by two distinct invariants.

Gaussian continuous tensor network states for simple bosonic field theories

Tensor networks states allow one to find the low-energy states of local lattice Hamiltonians through variational optimization. Recently, a construction of such states in the continuum was put forward, providing a first step towards the goal of solving quantum field theories (QFTs) variationally. However, the proposed manifold of continuous tensor network states (CTNSs) is difficult to study in full generality, because the expectation values of local observables cannot be computed analytically. In this paper we study a tractable subclass of CTNSs, the Gaussian CTNSs (GCTNSs), and benchmark them on simple quadratic and quartic bosonic QFT Hamiltonians. We show that GCTNSs provide arbitrarily accurate approximations to the ground states of quadratic Hamiltonians and decent estimates for quartic ones at weak coupling. Since they capture the short distance behavior of the theories we consider exactly, GCTNSs even allow one to renormalize away simple divergences variationally. In the end our study makes it plausible that CTNSs are indeed a good manifold to approximate the low-energy states of QFTs.Received 18 December 2020Accepted 9 March 2021DOI:https://doi.org/10.1103/PhysRevResearch.3.023059Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Open access publication funded by the Max Planck Society.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasQuantum field theory (low energy)TechniquesTensor network methodsVariational approachQuantum Information

Quadrupole-octupole coupling and the onset of octupole deformation in actinides

The evolution of quadrupole and octupole collectivity and their coupling is investigated in a series of even-even isotopes of the actinide Ra, Th, U, Pu, Cm, and Cf with neutron number in the interval $130\leqslant N\leqslant 150$. The Hartree-Fock-Bogoliubov approximation, based on the parametrization D1M of the Gogny energy density functional, is employed to generate potential energy surfaces depending upon the axially-symmetric quadrupole and octupole shape degrees of freedom. The mean-field energy surface is then mapped onto the expectation value of the $sdf$ interacting-boson-model Hamiltonian in the boson condensate state as to determine the strength parameters of the boson Hamiltonian. Spectroscopic properties related to the octupole degree of freedom are produced by diagonalizing the mapped Hamiltonian. Calculated low-energy negative-parity spectra, $B(E3;3^{-}_{1}\to 0^{+}_{1})$ reduced transition rates, and effective octupole deformation suggest that the transition from nearly spherical to stable octupole-deformed, and to octupole vibrational states occurs systematically in the actinide region.

Scalable W -type entanglement resource in neutral-atom arrays with Rydberg-dressed resonant dipole-dipole interaction

While the Rydberg-blockade regime provides the natural setting for creating $W$-type entanglement with cold neutral atoms, it is demonstrated here that a scalable entanglement resource of this type can even be obtained under completely different physical circumstances. To be more precise, a special instance of twisted $W$ states -- namely, $\pi$-twisted ones -- can be engineered in one-dimensional arrays of cold neutral atoms with Rydberg-dressed resonant dipole-dipole interaction. In particular, it is shown here that this is possible even when a (dressed) Rydberg excitation is coupled to the motional degrees of freedom of atoms in their respective, nearly-harmonic optical-dipole microtraps, which are quantized into dispersionless (zero-dimensional) bosons. For a specially chosen ("sweet-spot") detuning of the off-resonant dressing lasers from the relevant internal atomic transitions, the desired $\pi$-twisted $W$ state of Rydberg-dressed qubits is the ground state of the effective excitation -- boson Hamiltonian of the system in a broad window of the relevant parameters. Being at the same time separated from the other eigenstates by a gap equal to the single-boson energy, this $W$ state can be prepared using a Rabi-type driving protocol. The corresponding preparation times are independent of the system size and several orders of magnitude shorter than the effective lifetimes of the relevant atomic states.

Rotation-time symmetry in bosonic systems and the existence of exceptional points in the absence of {mathscr{PT}} symmetry

We study symmetries of open bosonic systems in the presence of laser pumping. Non-Hermitian Hamiltonians describing these systems can be parity-time ({{mathscr{PT}}}) symmetric in special cases only. Systems exhibiting this symmetry are characterised by real-valued energy spectra and can display exceptional points, where a symmetry-breaking transition occurs. We demonstrate that there is a more general type of symmetry, i.e., rotation-time ({mathscr{RT}}) symmetry. We observe that {mathscr{RT}}-symmetric non-Hermitian Hamiltonians exhibit real-valued energy spectra which can be made singular by symmetry breaking. To calculate the spectra of the studied bosonic non-diagonalisable Hamiltonians we apply diagonalisation methods based on bosonic algebra. Finally, we list a versatile set rules allowing to immediately identifying or constructing {mathscr{RT}}-symmetric Hamiltonians. We believe that our results on the {mathscr{RT}}-symmetric class of bosonic systems and their spectral singularities can lead to new applications inspired by those of the {{mathscr{PT}}}-symmetric systems.

Pairing vibrations in the interacting boson model based on density functional theory

We propose a method to incorporate the coupling between shape and pairing collective degrees of freedom in the framework of the interacting boson model (IBM), based on the nuclear density functional theory. To account for pairing vibrations, a boson-number non-conserving IBM Hamiltonian is introduced. The Hamiltonian is constructed by using solutions of self-consistent mean-field calculations based on a universal energy density functional and pairing force, with constraints on the axially-symmetric quadrupole and pairing intrinsic deformations. By mapping the resulting quadrupole-pairing potential energy surface onto the expectation value of the bosonic Hamiltonian in the boson condensate state, the strength parameters of the boson Hamiltonian are determined. An illustrative calculation is performed for $^{122}$Xe, and the method is further explored in a more systematic study of rare-earth $N=92$ isotones. The inclusion of the dynamical pairing degree of freedom significantly lowers the energies of bands based on excited $0^+$ states. The results are in quantitative agreement with spectroscopic data, and are consistent with those obtained using the collective Hamiltonian approach.

Restoring number conservation in quadratic bosonic Hamiltonians with dualities

Number-non-conserving terms in quadratic bosonic Hamiltonians can induce unwanted dynamical instabilities. By exploiting the pseudo-Hermitian structure built in to these Hamiltonians, we show that as long as dynamical stability holds, one may always construct a non-trivial dual (unitarily equivalent) number-conserving quadratic bosonic Hamiltonian. We exemplify this construction for a gapped harmonic chain and a bosonic analogue to Kitaev's Majorana chain. Our duality may be used to identify local number-conserving models that approximate stable bosonic Hamiltonians without the need for parametric amplification and to implement non-Hermitian -symmetric dynamics in non-dissipative number-conserving bosonic systems. Implications for computing topological invariants are addressed.

Deconstructing effective non-Hermitian dynamics in quadratic bosonic Hamiltonians

Unlike their fermionic counterparts, the dynamics of Hermitian quadratic bosonic Hamiltonians are governed by a generally non-Hermitian Bogoliubov-de Gennes effective Hamiltonian. This underlying non-Hermiticity gives rise to a dynamically stable regime, whereby all observables undergo bounded evolution in time, and a dynamically unstable one, whereby evolution is unbounded for at least some observables. We show that stability-to-instability transitions may be classified in terms of a suitably generalized symmetry, which can be broken when diagonalizability is lost at exceptional points in parameter space, but also when degenerate real eigenvalues split off the real axis while the system remains diagonalizable. By leveraging tools from Krein stability theory in indefinite inner-product spaces, we introduce an indicator of stability phase transitions, which naturally extends the notion of phase rigidity from non-Hermitian quantum mechanics to the bosonic setting. As a paradigmatic example, we fully characterize the stability phase diagram of a bosonic analogue to the Kitaev–Majorana chain under a wide class of boundary conditions. In particular, we establish a connection between phase-dependent transport properties and the onset of instability, and argue that stable regions in parameter space become of measure zero in the thermodynamic limit. Our analysis also reveals that boundary conditions that support Majorana zero modes in the fermionic Kitaev chain are precisely the same that support stability in the bosonic chain.

AN APPLICATION OF THE DEMOCRATIC MAPPING IN REALISTIC SYSTEMS

The democratic mapping is used for the calculation of low lying states of nuclei in the sd and fp shells. In addition to demonstrating the applicability of the method in realistic cases where many non-degenerate levels are present, the method allows for the ranking of the various bosons according to their importance as building blocks of low lying states. It is proven that the s and d bosons are the most important building blocks, followed by the d' and g bosons. Thus one of the basic assumptions of the Interacting Boson Model (IBM) is proven to be correct. Very good agreement between the boson calculation and the shell model results is obtained for A = 20 nuclei when 12 bosons are taken into account, while an even larger number of bosons is required to reproduce the low-lying states of the A = 44 nuclei. In order to obtain equally good results with a smaller number of bosons one needs to introduce effective boson hamiltonians which correspond to truncated fermion spaces.

An equation of motion approach for the vibrational transition energies in the effective harmonic oscillator formalism: the Random phase approximation

A theory for calculating vibrational energy levels and infrared intensities is developed in the equation of motion framework at the random phase approximation level. The vibrational Hamiltonian is expanded in the harmonic oscillator ladder operators making a Hamiltonian a bosonic Hamiltonian. The excitation operator is expanded to include at most two creations and two annihilation operators making it equivalent to the random phase approximation. The method is applied for the calculation of vibrational spectral properties of two molecules. The results are found to be satisfactory, making this approach a viable option for large molecular systems. A theory for calculating vibrational energy levels and infrared intensities is developed in the equation of motion framework at the random phase approximation level. The vibrational Hamiltonian is expressed in bosonic representation. The excitation operator is expanded to include at most two creation and annihilation operators.

Boson Description of Band Crossing in Even Barium Isotopes

Properties of the states of yrast bands in even Ba isotopes are studied on the basis of the boson representation of fermion operators. The parameters of the boson Hamiltonian and the interaction between collective quadrupole bosons, positive parity bosons, and spins $${{J}^{\pi }}$$ from $${{0}^{ + }}$$ to $${{10}^{ + }}$$ are calculated microscopically. This approach describes band crossing in the yrast bands and reproduces their energies and $$B(E2)$$ within these bands up to $${{I}^{\pi }} = {{18}^{ + }}.$$

All-creation and all-annihilation time-dependent PT -symmetric bosonic Hamiltonians: An infinite squeezing degree at a finite time

All-creation and all-annihilation time-dependent <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi mathvariant="script">PT</mml:mi></mml:math> -symmetric bosonic Hamiltonians: An infinite squeezing degree at a finite time

Optimal Upper Bound for the Correlation Energy of a Fermi Gas in the Mean-Field Regime

While Hartree–Fock theory is well established as a fundamental approximation for interacting fermions, it has been unclear how to describe corrections to it due to many-body correlations. In this paper we start from the Hartree–Fock state given by plane waves and introduce collective particle–hole pair excitations. These pairs can be approximately described by a bosonic quadratic Hamiltonian. We use Bogoliubov theory to construct a trial state yielding a rigorous Gell-Mann–Brueckner–type upper bound to the ground state energy. Our result justifies the random-phase approximation in the mean-field scaling regime, for repulsive, regular interaction potentials.

Non-Hermitian dynamics without dissipation in quantum systems

Models based on non-Hermitian Hamiltonians can exhibit a range of surprising and potentially useful phenomena. Physical realizations typically involve couplings to sources of incoherent gain and loss; this is problematic in quantum settings, because of the unavoidable fluctuations associated with this dissipation. Here, we present several routes for obtaining unconditional non-Hermitian dynamics in non-dissipative quantum systems. We exploit the fact that quadratic bosonic Hamiltonians that do not conserve particle number give rise to non-Hermitian dynamical matrices. We discuss the nature of these mappings from non-Hermitian to Hermitian Hamiltonians, and explore applications to quantum sensing, entanglement dynamics and topological band theory. The systems we discuss could be realized in a variety of photonic and phononic platforms using the ubiquitous resource of parametric driving.

Density–driven superfluid transition of the constrained bosons on a lattice

We formulate the Landau theory and study its consequences for lattice bosons, which are stabilized by a three-body on-site constraint. This is relevant for experiments with cold atoms in optical lattices, where the presence of three-body dissipative loss processes renders the occupation number n for bosons on short timescales effectively constrained to n≤2. We elaborate on the Landau free energy expansion resulting from the Hamiltonian of bosons in the restricted Hilbert space, and present the phase diagram of the constrained Bose lattice gas.

Floquet time crystals in clock models

We construct a class of period-$n$-tupling discrete time crystals based on $\mathbb{Z}_n$ clock variables, for all the integers $n$. We consider two classes of systems where this phenomenology occurs, disordered models with short-range interactions and fully connected models. In the case of short-range models we provide a complete classification of time-crystal phases for generic $n$. For the specific cases of $n=3$ and $n=4$ we study in details the dynamics by means of exact diagonalisation. In both cases, through an extensive analysis of the Floquet spectrum, we are able to fully map the phase diagram. In the case of infinite-range models, the mapping onto an effective bosonic Hamiltonian allows us to investigate the scaling to the thermodynamic limit. After a general discussion of the problem, we focus on $n=3$ and $n=4$, representative examples of the generic behaviour. Remarkably, for $n=4$ we find clear evidence of a new crystal-to-crystal transition between period $n$-tupling and period $n/2$-tupling.

Analyticity of resonances and eigenvalues and spectral properties of the massless Spin–Boson model

We extend the method of Pizzo multiscale analysis for resonances introduced in [5] in order to infer analytic properties of resonances and eigenvalues (and their eigenprojections) as well as estimates for the localization of the spectrum of dilated Hamiltonians and norm-bounds for the corresponding resolvent operators, in neighborhoods of resonances and eigenvalues. We apply our method to the massless Spin–Boson model assuming a slight infrared regularization. We prove that the resonance and the ground-state eigenvalue (and their eigenprojections) are analytic with respect to the dilation parameter and the coupling constant. Moreover, we prove that the spectrum of the dilated Spin–Boson Hamiltonian in the neighborhood of the resonance and the ground-state eigenvalue is localized in two cones in the complex plane with vertices at the location of the resonance and the ground-state eigenvalue, respectively. Additionally, we provide norm-estimates for the resolvent of the dilated Spin–Boson Hamiltonian near the resonance and the ground-state eigenvalue. The topic of analyticity of eigenvalues and resonances has let to several studies and advances in the past. However, to the best of our knowledge, this is the first time that it is addressed from the perspective of Pizzo multiscale analysis. Once the multiscale analysis is set up our method gives easy access to analyticity: Essentially, it amounts to proving it for isolated eigenvalues only and use that uniform limits of analytic functions are analytic. The type of spectral and resolvent estimates that we prove are needed to control the time evolution including the scattering regime. The latter will be demonstrated in a forthcoming publication. The introduced multiscale method to study spectral and resolvent estimates follows its own inductive scheme and is independent (and different) from the method we apply to construct resonances.

Broken-Hermiticity phase transition in the Bose-Hubbard model

A new version of the change of the "phase" (i.e., of the set of observable characteristics) of a quantum system is proposed. In a general scenario the evolution is assumed generated, before the phase transition, by some standard Hermitian Hamiltonian $H^{(before)}$, and, after the phase transition, by one of the recently very popular non-standard, non-Hermitian (but hiddenly Hermitian, i.e., still unitarity-guaranteeing) Hamiltonians $H^{(after)}$. For consistency, a smoothness of matching between the two operators as well as between the related physical Hilbert spaces must be guaranteed. The feasibility of the idea is illustrated via the two-mode $(N-1)-$bosonic Bose-Hubbard Hamiltonian. In $H^{(before)}=H^{(BH)}(\varepsilon)$ we use the decreasing real $\varepsilon^{(before)} \to 0$. In the hiddenly Hermitian continuation $H^{(after)}=H^{(BH)}(\tilde{\varepsilon})$ the imaginary part of the purely imaginary $\tilde{\varepsilon}^{(after)}$ grows. The smoothness of the transition occurring at the interface $\varepsilon=\tilde{\varepsilon}=0$ is then guaranteed by an {\it ad hoc\,} amendment of the inner product in Hilbert space "after". The trivial Hilbert-space metric $\Theta^{(before)}=I$ must match $\Theta^{(after)} \neq I$ smoothly. This is confirmed and illustrated by the explicit constructions of a few $\Theta^{(after)}$s in closed form.

Stochastic unravelings of non-Markovian completely positive and trace-preserving maps

We consider open quantum systems with factorized initial states, providing the structure of the reduced system dynamics, in terms of environment cumulants. We show that such completely positive (CP) and trace preserving (TP) maps can be unraveled by linear stochastic Schr\"odinger equations (SSEs) characterized by sets of colored stochastic processes (with n-th order cumulants). We obtain both the conditions such that the SSEs provide CPTP dynamics, and those for unraveling an open system dynamics. We then focus on Gaussian non-Markovian unravellings, whose known structure displays a functional derivative. We provide a novel description that replaces the functional derivative with a recursive operatorial structure. Moreover, for the family of quadratic bosonic Hamiltonians, we are able to provide an explicit operatorial dependence for the unravelling.