Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Nonlinear light absorption in many-electron systems excited by an instantaneous electric field: a non-perturbative approach

Applications of low-cost non-perturbative approaches in real time, such as time-dependent density functional theory, for the study of nonlinear optical properties of large and complex systems are gaining increasing popularity. However, their assessment still requires the analysis and understanding of elementary dynamical processes in simple model systems. Motivated by the aim of simulating optical nonlinearities in molecules, here exemplified by the case of the quaterthiophene oligomer, we investigate light absorption in many-electron interacting systems beyond the linear regime by using a single broadband impulse of an electric field; i.e. an electrical impulse in the instantaneous limit. We determine non-pertubatively the absorption cross section from the Fourier transform of the time-dependent induced dipole moment, which can be obtained from the time evolution of the wavefunction. We discuss the dependence of the resulting cross section on the magnitude of the impulse and we highlight the advantages of this method in comparison with perturbation theory by working on a one-dimensional model system for which numerically exact solutions are accessible. Thus, we demonstrate that the considered non-pertubative approach provides us with an effective tool for investigating fluence-dependent nonlinear optical excitations.

Ultrafast dynamical Lifshitz transition

Fermi surface is at the heart of our understanding of metals and strongly correlated many-body systems. An abrupt change in the Fermi surface topology, also called Lifshitz transition, can lead to the emergence of fascinating phenomena like colossal magnetoresistance and superconductivity. While Lifshitz transitions have been demonstrated for a broad range of materials by equilibrium tuning of macroscopic parameters such as strain, doping, pressure, and temperature, a nonequilibrium dynamical route toward ultrafast modification of the Fermi surface topology has not been experimentally demonstrated. Combining time-resolved multidimensional photoemission spectroscopy with state-of-the-art TDDFT+U simulations, we introduce a scheme for driving an ultrafast Lifshitz transition in the correlated type-II Weyl semimetal T _d-MoTe₂ We demonstrate that this nonequilibrium topological electronic transition finds its microscopic origin in the dynamical modification of the effective electronic correlations. These results shed light on a previously unexplored ultrafast scheme for controlling the Fermi surface topology in correlated quantum materials.

Simulating Vibronic Spectra without Born-Oppenheimer Surfaces

We show how linear vibronic spectra in molecular systems can be simulated efficiently using first-principles approaches without relying on the explicit use of multiple Born-Oppenheimer potential energy surfaces. We demonstrate and analyze the performance of mean-field and beyond-mean-field dynamics techniques for the H₂ molecule in one dimension, in the later case capturing the vibronic structure quite accurately, including quantum Franck-Condon effects. In a practical application of this methodology we simulate the absorption spectrum of benzene in full dimensionality using time-dependent density functional theory at the multitrajectory Ehrenfest level, finding good qualitative agreement with experiment and significant spectral reweighting compared to commonly used single-trajectory Ehrenfest dynamics. These results form the foundation for nonlinear spectral calculations and show promise for future application in capturing phenomena associated with vibronic coupling in more complex molecular and potentially condensed phase systems.

From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape

Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.

Hybridized Defects in Solid-State Materials as Artificial Molecules

Two-dimensional materials can be crafted with structural precision approaching the atomic scale, enabling quantum defects-by-design. These defects are frequently described as "artificial atoms" and are emerging optically addressable spin qubits. However, interactions and coupling of such artificial atoms with each other, in the presence of the lattice, warrants further investigation. Here we present the formation of "artificial molecules" in solids, introducing a chemical degree of freedom in control of quantum optoelectronic materials. Specifically, in monolayer hexagonal boron nitride as our model system, we observe configuration- and distance-dependent dissociation curves and hybridization of defect orbitals within the bandgap into bonding and antibonding orbitals, with splitting energies ranging from ∼10 meV to nearly 1 eV. We calculate the energetics of cis and trans out-of-plane defect pairs CH_B-CH_B against an in-plane defect pair C_B-C_B and find that in-plane defect pair interacts more strongly than out-of-plane pairs. We demonstrate an application of this chemical degree of freedom by varying the distance between C_B and V_N of C_BV_N and observe changes in the predicted peak absorption wavelength from the visible to the near-infrared spectral band. We envision leveraging this chemical degree of freedom of defect complexes to precisely control and tune defect properties toward engineering robust quantum memories and quantum emitters for quantum information science.

Simulation of photo-electron spectrum and electron scattering by dual time propagation

A dual time propagation approach is introduced to describe electron scattering and ionization. The space is divided into two regions, a central region with a full time-dependent Hamiltonian and an outer region where the kinetic operator and the laser field dominate. The two regions are connected by a source term. Time-dependent density functional theory calculations of wave packet scattering on molecules and photoelectron spectrum due to circularly polarized laser are presented to illustrate the efficiency and applicability of the approach.

Light-matter interaction of a molecule in a dissipative cavity from first principles

Cavity-mediated light-matter coupling can dramatically alter opto-electronic and physico-chemical properties of a molecule. Ab initio theoretical predictions of these systems need to combine non-perturbative, many-body electronic structure theory-based methods with cavity quantum electrodynamics and theories of open-quantum systems. Here, we generalize quantum-electrodynamical density functional theory to account for dissipative dynamics of the cavity and describe coupled cavity-single molecule interactions in the weak-to-strong-coupling regimes. Specifically, to establish this generalized technique, we study excited-state dynamics and spectral responses of benzene and toluene under weak-to-strong light-matter coupling. By tuning the coupling, we achieve cavity-mediated energy transfer between electronically excited states. This generalized ab initio quantum-electrodynamical density functional theory treatment can be naturally extended to describe cavity-mediated interactions in arbitrary electromagnetic environments, accessing correlated light-matter observables and thereby closing the gap between electronic structure theory, quantum optics, and nanophotonics.

Intermolecular interactions in optical cavities: An ab initio QED study

Intermolecular bonds are weak compared to covalent bonds, but they are strong enough to influence the properties of large molecular systems. In this work, we investigate how strong light-matter coupling inside an optical cavity can modify intermolecular forces and illustrate the varying necessity of correlation in their description. The electromagnetic field inside the cavity can modulate the ground state properties of weakly bound complexes. Tuning the field polarization and cavity frequency, the interactions can be stabilized or destabilized, and electron densities, dipole moments, and polarizabilities can be altered. We demonstrate that electron-photon correlation is fundamental to describe intermolecular interactions in strong light-matter coupling. This work proposes optical cavities as a novel tool to manipulate and control ground state properties, solvent effects, and intermolecular interactions for molecules and materials.

Role of exchange and correlation in high-harmonic generation spectra of H2, N2, and CO2: Real-time time-dependent electronic-structure approaches

This study arises from the attempt to answer the following question: how different descriptions of electronic exchange and correlation affect the high-harmonic generation (HHG) spectroscopy of H₂, N₂, and CO₂ molecules? We compare HHG spectra for H₂, N₂, and CO₂ with different ab initio electronic structure methods: real-time time-dependent configuration interaction and real-time time-dependent density functional theory (RT-TDDFT) using truncated basis sets composed of correlated wave functions expanded on Gaussian basis sets. In the framework of RT-TDDFT, we employ Perdew-Burke-Ernzerhof (PBE) and long-range corrected Perdew-Burke-Ernzerhof (LC-ωPBE) functionals. We study HHG spectroscopy by disentangling the effect of electronic exchange and correlation. We first analyze the electronic exchange alone, and in the case of RT-TDDFT with LC-ωPBE, we use ω = 0.3 and ω = 0.4 to tune the percentage of long-range Hartree-Fock exchange and short-range exchange PBE. Then, we added the correlation as described by the PBE functional. All the methods give very similar HHG spectra, and they seem not to be particularly sensitive to the different description of exchange and correlation or to the correct asymptotic behavior of the Coulomb potential. Despite this general trend, some differences are found in the region connecting the cutoff and the background. Here, the harmonics can be resolved with different accuracy depending on the theoretical schemes used. We believe that the investigation of the molecular continuum and its coupling with strong fields merits further theoretical investigations in the near future.

Quantum Information and Algorithms for Correlated Quantum Matter

Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.

Direct Measurement of Electron-Phonon Coupling with Time-Resolved ARPES

Time- and angular- resolved photoelectron spectroscopy is a powerful technique to measure electron dynamics in solids. Recent advances in this technique have facilitated band and energy resolved observations of the effect that excited phonons, have on the electronic structure. Here, we show with the help of ab initio simulations that the Fourier analysis of the time-resolved measurements of solids with excited phonon modes enables the determination of the band- and mode-resolved electron-phonon coupling directly from the experimental data without any additional input from theory. Such an observation is not restricted to regions of strong electron-phonon coupling and does not require strongly excited or hot phonons, but can be employed to monitor the dynamical renormalization of phonons in driven phases of matter.

Light-Matter Hybrid-Orbital-Based First-Principles Methods: The Influence of Polariton Statistics

A detailed understanding of strong matter–photon interactions requires first-principle methods that can solve the fundamental Pauli–Fierz Hamiltonian of nonrelativistic quantum electrodynamics efficiently. A possible way to extend well-established electronic-structure methods to this situation is to embed the Pauli–Fierz Hamiltonian in a higher-dimensional light–matter hybrid auxiliary configuration space. In this work we show the importance of the resulting hybrid Fermi–Bose statistics of the polaritons, which are the new fundamental particles of the “photon-dressed” Pauli–Fierz Hamiltonian for systems in cavities. We show that violations of these statistics can lead to unphysical results. We present an efficient way to ensure the correct statistics by enforcing representability conditions on the dressed one-body reduced density matrix. We further present a general prescription how to extend a given first-principles approach to polaritons and as an example introduce polaritonic Hartree–Fock theory. While being a single-reference method in polariton space, polaritonic Hartree–Fock is a multireference method in the electronic space, i.e., it describes electronic correlations. We also discuss possible applications to polaritonic QEDFT. We apply this theory to a lattice model and find that, the more delocalized the bound-state wave function of the particles is, the stronger it reacts to photons. The main reason is that within a small energy range, many states with different electronic configurations are available as opposed to a strongly bound (and hence energetically separated) ground-state wave function. This indicates that under certain conditions coupling to the quantum vacuum of a cavity can indeed modify ground state properties.

Ab initio polaritonic potential-energy surfaces for excited-state nanophotonics and polaritonic chemistry

Advances in nanophotonics, quantum optics, and low-dimensional materials have enabled precise control of light-matter interactions down to the nanoscale. Combining concepts from each of these fields, there is now an opportunity to create and manipulate photonic matter via strong coupling of molecules to the electromagnetic field. Toward this goal, here we demonstrate a first principles framework to calculate polaritonic excited-state potential-energy surfaces, transition dipole moments, and transition densities for strongly coupled light-matter systems. In particular, we demonstrate the applicability of our methodology by calculating the polaritonic excited-state manifold of a formaldehyde molecule strongly coupled to an optical cavity. This proof-of-concept calculation shows how strong coupling can be exploited to alter photochemical reaction pathways by influencing avoided crossings with tuning of the cavity frequency and coupling strength. Therefore, by introducing an ab initio method to calculate excited-state potential-energy surfaces, our work opens a new avenue for the field of polaritonic chemistry.

Astrochemical relevance of VUV ionization of large PAH cations. ASTRONOMY AND ASTROPHYSICS 2020;641:A98. [PMID: 33154599 PMCID: PMC7116310 DOI: 10.1051/0004-6361/202038139]

CONTEXT

As a part of interstellar dust, polycyclic aromatic hydrocarbons (PAHs) are processed by the interaction with vacuum ultra-violet (VUV) photons that are emitted by hot young stars. This interaction leads to the emission of the well-known aromatic infrared bands but also of electrons, which can significantly contribute to the heating of the interstellar gas.

AIMS

Our aim is to investigate the impact of molecular size on the photoionization properties of cationic PAHs.

METHODS

Trapped PAH cations of sizes between 30 and 48 carbon atoms were submitted to VUV photons in the range of 9 to 20 eV from the DESIRS beamline at the synchrotron SOLEIL. All resulting photoproducts including dications and fragment cations were mass-analyzed and recorded as a function of photon energy.

RESULTS

Photoionization is found to be predominant over dissociation at all energies, which differs from an earlier study on smaller PAHs. The photoionization branching ratio reaches 0.98 at 20 eV for the largest studied PAH. The photoionization threshold is observed to be between 9.1 and 10.2 eV, in agreement with the evolution of the ionization potential with size. Ionization cross sections were indirectly obtained and photoionization yields extracted from their ratio with theoretical photoabsorption cross sections, which were calculated using time-dependent density functional theory. An analytical function was derived to calculate this yield for a given molecular size.

CONCLUSIONS

Large PAH cations could be efficiently ionized in H i regions and provide a contribution to the heating of the gas by photoelectric effect. Also, at the border of or in H ii regions, PAHs could be exposed to photons of energy higher than 13.6 eV. Our work provides recipes to be used in astronomical models to quantify these points.

Understanding real-time time-dependent density-functional theory simulations of ultrafast laser-induced dynamics in organic molecules

Real-time time-dependent density functional theory, in conjunction with the Ehrenfest molecular dynamics scheme, is becoming a popular methodology to investigate ultrafast phenomena on the nanoscale. Thanks to recent developments, it is also possible to explicitly include in the simulations a time-dependent laser pulse, thereby accessing the transient excitation regime. However, the complexity entailed in these calculations calls for in-depth analysis of the accessible and yet approximate (either "dressed" or "bare") quantities in order to evaluate their ability to provide us with a realistic picture of the simulated processes. In this work, we analyze the ultrafast dynamics of three small molecules (ethylene, benzene, and thiophene) excited by a resonant laser pulse in the framework of the adiabatic local-density approximation. The electronic response to the laser perturbation in terms of induced dipole moment and excited-state population is compared to the results given by an exactly solvable two-level model. In this way, we can interpret the charge-carrier dynamics in terms of simple estimators, such as the number of excited electrons. From the computed transient absorption spectra, we unravel the appearance of nonlinear effects such as excited-state absorption and vibronic coupling. In this way, we observe that the laser excitation affects the vibrational spectrum by enhancing the anharmonicities therein, while the coherent vibrational motion contributes to stabilizing the electronic excitation already within a few tens of femtoseconds.

The CECAM electronic structure library and the modular software development paradigm

First-principles electronic structure calculations are now accessible to a very large community of users across many disciplines, thanks to many successful software packages, some of which are described in this special issue. The traditional coding paradigm for such packages is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, essentially from the compiler up, possibly with the exception of linear-algebra and message-passing libraries. This model has endured and been quite successful for decades. The successful evolution of the electronic structure methodology itself, however, has resulted in an increasing complexity and an ever longer list of features expected within all software packages, which implies a growing amount of replication between different packages, not only in the initial coding but, more importantly, every time a code needs to be re-engineered to adapt to the evolution of computer hardware architecture. The Electronic Structure Library (ESL) was initiated by CECAM (the European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure codes and redesign them as open-source libraries available to everybody. Such libraries include "heavy-duty" ones that have the potential for a high degree of parallelization and adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by, e.g., physicists and chemists when implementing new ideas. We envisage that this modular paradigm will improve overall coding efficiency and enable specialists (whether they be computer scientists or computational scientists) to use their skills more effectively and will lead to a more dynamic evolution of software in the community as well as lower barriers to entry for new developers. The model comes with new challenges, though. The building and compilation of a code based on many interdependent libraries (and their versions) is a much more complex task than that of a code delivered in a single self-contained package. Here, we describe the state of the ESL, the different libraries it now contains, the short- and mid-term plans for further libraries, and the way the new challenges are faced. The ESL is a community initiative into which several pre-existing codes and their developers have contributed with their software and efforts, from which several codes are already benefiting, and which remains open to the community.