Ground states of linear rotor chains via the density matrix renormalization group

Ground states of planar dipolar rotor chains with recurrent neural networks

In this contribution, we employ a recurrent neural network (RNN) architecture in a variational optimization to obtain the ground state of linear chains of planar, dipolar rotors. We test different local basis sets and discuss their impact on the sign structure of the many-body ground state wavefunction. It is demonstrated that the RNN ansatz we employ is able to treat systems with and without a sign problem in the ground state. For larger chains with up to 50 rotors, accurate properties, such as correlation functions and Binder parameters, are calculated. By employing quantum annealing, we show that precise entanglement properties can be obtained. All these properties allow one to identify a quantum phase transition between a paraelectric and a ferroelectric quantum phase.

Quantum criticality in chains of planar rotors with dipolar interactions

In this work, we perform a density matrix renormalization group study of chains of planar rotors interacting via dipolar interactions. By exploring the ground state from weakly to strongly interacting rotors, we find the occurrence of a quantum phase transition between a disordered and a dipole-ordered quantum state. We show that the nature of the ordered state changes from ferroelectric to antiferroelectric when the relative orientation of the rotor planes varies and that this change requires no modification of the overall symmetry. The observed quantum phase transitions are characterized by critical exponents and central charges, which reveal different universality classes ranging from that of the (1 + 1)D Ising model to the 2D classical XY model.

Quantum Criticality and Universal Behavior in Molecular Dipolar Lattices of Endofullerenes

Fullerene cages allow the confinement of single molecules and the construction of molecular assemblies whose properties strongly differ from those of free species. In this work, we employ the density-matrix renormalization group method to show that chains of fullerenes filled with polar molecules (LiF, HF, and H₂O) can form dipole-ordered quantum phases. In symmetry broken environments, these ordered phases are ferroelectric, a property that makes them promising candidates for quantum devices. We demonstrate that for a given guest molecule, the occurrence of these quantum phases can be enforced or influenced either by changing the effective electric dipole moment or by isotopic substitution. In the ordered phase, all systems under consideration are characterized by universal behavior that depends only on the ratio of the effective electric dipole and of the rotational constant. A phase diagram is derived, and further molecules are proposed as candidates for dipole-ordered endofullerene chains.

Excited-State DMRG Made Simple with FEAST

We introduce DMRG[FEAST], a new method for optimizing excited-state many-body wave functions with the density matrix renormalization group (DMRG) algorithm. Our approach applies the FEAST algorithm, originally designed for large-scale diagonalization problems, to matrix product state wave functions. We show that DMRG[FEAST] enables the stable optimization of both low- and high-energy eigenstates, therefore overcoming the limitations of state-of-the-art excited-state DMRG algorithms. We demonstrate the reliability of DMRG[FEAST] by calculating anharmonic vibrational excitation energies of molecules with up to 30 fully coupled degrees of freedom.

A path integral ground state approach for asymmetric top rotors with nuclear spin symmetry: Application to water chains

A path integral ground state (PIGS) approach for the simulation of asymmetric top rotors is presented. The method is based on Monte Carlo sampling of angular degrees of freedom. A symmetry-adapted rotational density matrix is used to account for nuclear spin statistics. To illustrate the method, ground-state properties of collections of para-water molecules confined to a one-dimensional lattice are computed. Those include energetic and structural observables. An advantage of the PIGS method is that expectation values can be obtained directly since the square of the wavefunction is sampled during a simulation. To benchmark the method, ground state energies and orientational distributions are computed using exact diagonalization for a single para-water molecule in an external field using a finite basis of symmetric top eigenfunctions. Benchmark results are also provided for N = 2 para-water molecules pinned to lattice sites at various distances to sample the crossover from hydrogen bonding to the dipole-dipole interaction regime. Excellent agreement between the PIGS results and the finite basis set calculations is observed. A thorough analysis of the convergence in terms of the imaginary time propagation length and systematic Trotter error is performed. The PIGS approach is then applied to a chain of N = 11 water molecules, and an equation of state is constructed in terms of the intermolecular separation. Ordering effects are also studied, and a transition between hydrogen bonding to dipole-dipole alignment is observed. The method is scalable and can also be applied in higher dimensions.

Comparison of the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalization group (DMRG) for ground state properties of linear rotor chains

We demonstrate the applicability of the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method to the problem of computing ground states of one-dimensional chains of linear rotors with dipolar interactions. Specifically, we successfully obtain energies, entanglement entropies, and orientational correlations that are in agreement with the Density Matrix Renormalization Group (DMRG), which has been previously used for this system. We find that the entropies calculated by ML-MCTDH for larger system sizes contain nonmonotonicity, as expected in the vicinity of a second-order quantum phase transition between ordered and disordered rotor states. We observe that this effect remains when all couplings besides nearest-neighbor are omitted from the Hamiltonian, which suggests that it is not sensitive to the rate of decay of the interactions. In contrast to DMRG, which is tailored to the one-dimensional case, ML-MCTDH (as implemented in the Heidelberg MCTDH package) requires more computational time and memory, although the requirements are still within reach of commodity hardware. The numerical convergence and computational demand of two practical implementations of ML-MCTDH and DMRG are presented in detail for various combinations of system parameters.

Vibrational Raman Shifts of Spin Isomer Combinations of Hydrogen Dimers and Isotopologues

Binding energies for para-para, ortho-para, and ortho-ortho hydrogen dimers (H₂)₂ are calculated using the six-dimensional (6D) interaction potential developed by Hinde [ J. Chem. Phys. 2008, 128, 154308]. The eigenstates of the dimers are computed by diagonalization using, as a basis, products of the rovibrational states of the monomers, a radial grid for the distance between the monomers, and spherical harmonics for the end-over-end rotation of the dimer. We describe the overall nuclear spin symmetry and use these properties to determine the relative population of various states, making use of a Boltzmann factor for each spin isomer to assess the effect of temperature. A predicted Raman spectrum in the Q(0) and Q(1) region of the hydrogen dimer is produced. To assess the accuracy of our model, we verify our produced shifts with experimental results obtained previously by Montero et al. [ Eur. Phys. J. D 2009, 52, 31-34] and find good agreement. These results are extended to other cases involving the deuterium (D₂)₂ and tritium dimer (T₂)₂ isotopologues, to predict Raman shifts.

Perspective: Accurate treatment of the quantum dynamics of light molecules inside fullerene cages: Translation-rotation states, spectroscopy, and symmetry breaking

In this perspective, I review the current status of the theoretical investigations of the quantum translation-rotation (TR) dynamics and spectroscopy of light molecules encapsulated inside fullerenes, mostly C₆₀ and C₇₀. The methodologies developed in the past decade allow accurate quantum calculations of the TR eigenstates of one and two nanoconfined molecules and have led to deep insights into the nature of the underlying dynamics. Combining these bound-state methodologies with the formalism of inelastic neutron scattering (INS) has resulted in the novel and powerful approach for the quantum calculation of the INS spectra of a diatomic molecule in a nanocavity with an arbitrary geometry. These simulations have not only become indispensable for the interpretation and assignment of the experimental spectra but are also behind the surprising discovery of the INS selection rule for diatomics in near-spherical nanocavities. Promising directions for future research are discussed.