Absence of localization in interacting spin chains with a discrete symmetry.

Novel paradigms of strong ergodicity breaking have recently attracted significant attention in condensed matter physics. Understanding the exact conditions required for their emergence or breakdown not only sheds more light on thermalization and its absence in closed quantum many-body systems, but it also has potential benefits for applications in quantum information technology. A case of particular interest is many-body localization whose conditions are not yet fully settled. Here, we prove that spin chains symmetric under a combination of mirror and spin-flip symmetries and with a non-degenerate spectrum show finite spin transport at zero total magnetization and infinite temperature. We demonstrate this numerically using two prominent examples: the Stark many-body localization system (Stark-MBL) and the symmetrized many-body localization system (symmetrized-MBL). We provide evidence of delocalization at all energy densities and show that delocalization persists when the symmetry is broken. We use our results to construct two localized systems which, when coupled, delocalize each other. Our work demonstrates the dramatic effect symmetries can have on disordered systems, proves that the existence of exact resonances is not a sufficient condition for delocalization, and opens the door to generalization to higher spatial dimensions and different conservation laws.

Time evolution of ML-MCTDH wavefunctions. I. Gauge conditions, basis functions, and singularities.

We derive a family of equations-of-motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than those used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOM approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.

Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator.

The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3-4 orders of magnitude fewer Hamiltonian evaluations and 2-3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2-3/1-2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ∼1.3×109 parameters to obtain accurate dynamics.

Phonon-induced disorder in dynamics of optically pumped metals from nonlinear electron-phonon coupling.

The non-equilibrium dynamics of matter excited by light may produce electronic phases, such as laser-induced high-transition-temperature superconductivity, that do not exist in equilibrium. Here we simulate the dynamics of a metal driven at initial time by a spatially uniform pump that excites dipole-active vibrational modes which couple nonlinearly to electrons. We provide evidence for rapid loss of spatial coherence, leading to emergent effective disorder in the dynamics, which arises in a system unitarily evolving under a translation-invariant Hamiltonian, and dominates the electronic behavior as the system evolves towards a correlated electron-phonon long-time state, possibly explaining why transient superconductivity is not observed. Our framework provides a basis within which to understand correlation dynamics in current pump-probe experiments of vibrationally coupled electrons, highlight the importance of the evolution of phase coherence, and demonstrate that pumped electron-phonon systems provide a means of realizing dynamically induced disorder in translation-invariant systems.

Multiset Matrix Product State Calculations Reveal Mobile Franck-Condon Excitations Under Strong Holstein-Type Coupling.

We show that the dynamics of (vertical) Franck-Condon excitations in the regime where Holstein-coupled vibrational modes mix strongly with electronic degrees of freedom sharply contrasts with the known self-localized behavior of vibrationally relaxed excitations. Instead, the strongly coupled modes are found to periodically induce resonances between interacting electronic sites, during which effective excitation transfer occurs, allowing Franck-Condon excitations to attain substantial mean square displacements under conditions where relaxed excitations are essentially trapped to a single site. In demonstrating this behavior, we employ a multiset matrix product state formalism. We find this tensor network state method to be a remarkably efficient and accurate approach for the notoriously difficult problem posed by the Holstein model in the regime where the electronic coupling, the vibrational quantum, and the vibrational reorganization energy are comparable in magnitude.

Triplet Separation Drives Singlet Fission after Femtosecond Correlated Triplet Pair Production in Rubrene.

Singlet fission, a multistep molecular process in which one photon generates two triplet excitons, holds great technological promise. Here, by applying a combination of transient transmittance and two-dimensional electronic spectroscopy with 5 fs laser pulses, we resolve the full set of fission steps before the onset of spin dephasing. In addition to its role as a viable singlet fission material, single-crystalline rubrene is selected because its energetics and transition dipole alignment uniquely allow for the unambiguous identification of the various fission steps through their contributions to distinct spectroscopic features. The measurements reveal that the neighboring correlated triplet pair achieves its maximum population within 20 fs. Subsequent growth of the triplet signal on picosecond time scales is attributable to spatial separation of the triplets, proceeding nonadiabatically through weakly coupled but near-resonant states. As such, we provide evidence in crystalline rubrene for a singlet fission step that, until now, has not been convincingly observed.

Implementation of a novel projector-splitting integrator for the multi-configurational time-dependent Hartree approach.

The variational equations of motion of the multi-configuration time-dependent Hartree (MCTDH) approach contain the inverse of reduced density matrices which are typically ill-conditioned and therefore lead to small stepsizes for numerical time integration. This problem is usually dealt with via regularization of the density matrices, which works well in many cases but still calls for systematic improvement schemes. Recently this problem, its implications and possible solutions have become the subject of increased interest. Notably, a projector splitting integrator for the MCTDH approach that does not require the inversion of reduced density matrices has been proposed [C. Lubich, Appl. Math. Res. Express 2015, 311]. Here, we present the first implementation of this integration scheme. Results for low-dimensional benchmark systems are presented, and the case of initially unoccupied single-particle functions is discussed.

Density functional theory study of the conformation and optical properties of hybrid Au(n)-dithienylethene systems (n = 3, 19, 25).

We present a theoretical study of Aun-dithienylethene hybrid systems (n = 3, 19, 25), where the organic molecule is covalently linked to a nanometer-scaled gold nanoparticle (NP). We aim at gaining insights on the optical properties of such photochromic devices and proposing a size-limited gold aggregate model able to recover the optical properties of the experimental system. We thus present a DFT-based calculation scheme to model the ground-state (conformation, energetic parameters) and excited-state properties (UV-visible absorption spectra) of this type of hybrid systems. Within this framework, the structural parameters (adsorption site, orientation, and internal structure of the photochrome) are found to be slightly dependent on the size/shape of the gold aggregate. The influence of the gold fragment on the optical properties of the resulting hybrid system is then discussed with the help of TD-DFT combined with an analysis of the virtual orbitals involved in the photochromic transitions. We show that, for the open hybrid isomer, the number of gold atoms is the key parameter to recover the photoactive properties that are experimentally observed. On the contrary, for hybrid closed systems, the three-dimensional structure of the metallic aggregate is of high impact. We thus conclude that Au25 corresponds to the most appropriate fragment to model nanometer-sized NP-DTE hybrid device.