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

PyUNIxMD: A Python-based excited state molecular dynamics package

Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.

Many-electron dynamics in laser-driven molecules: wavefunction theory vs. density functional theory

With recent experimental advances in laser-driven electron dynamics in polyatomic molecules, the need arises for their reliable theoretical modelling. Among efficient, yet fairly accurate methods for many-electron dynamics are Time-Dependent Configuration Interaction Singles (TD-CIS) (a Wave Function Theory (WFT) method), and Real-Time Time-Dependent Density Functional Theory (RT-TD-DFT), respectively. Here we compare TD-CIS combined with extended Atomic Orbital (AO) bases, TD-CIS/AO, with RT-TD-DFT in a grid representation of the Kohn-Sham orbitals, RT-TD-DFT/Grid. Possible ionization losses are treated by complex absorbing potentials in energy space (for TD-CIS/AO) or real space (for RT-TD-DFT), respectively. The comparison is made for two test cases: (i) state-to-state transitions using resonant lasers (π-pulses), i.e., bound electron motion, and (ii) large-amplitude electron motion leading to High Harmonic Generation (HHG). Test systems are a H2 molecule and cis- and trans-1,2-dichlorethene, C2H2Cl2, (DCE). From time-dependent electronic energies, dipole moments and from HHG spectra, the following observations are made: first, for bound state-to-state transitions enforced by π-pulses, TD-CIS nicely accounts for the expected population inversion in contrast to RT-TD-DFT, in agreement with earlier findings. Secondly, when using laser pulses under non-resonant conditions, dipole moments and lower harmonics in HHG spectra are obtained by TD-CIS/AO which are in good agreement with those obtained with RT-TD-DFT/Grid. Deviations become larger for higher harmonics and at low laser intensities, i.e., for low-intensity HHG signals. We also carefully test effects of basis sets for TD-CIS/AO and grid size for RT-TD-DFT/Grid, different exchange-correlation functionals in RT-TD-DFT, and absorbing boundaries. Finally, for the present examples, TD-CIS/AO is observed to be at least an order of magnitude more computationally efficient than RT-TD-DFT/Grid.

LayerPCM: An implicit scheme for dielectric screening from layered substrates

We present LayerPCM, an extension of the polarizable-continuum model coupled to real-time time-dependent density-functional theory, for an efficient and accurate description of the electrostatic interactions between molecules and multilayered dielectric substrates on which they are physisorbed. The former are modeled quantum-mechanically, while the latter are treated as polarizable continua characterized by their dielectric constants. The proposed approach is purposely designed to simulate complex hybrid heterostructures with nano-engineered substrates including a stack of anisotropic layers. LayerPCM is suitable for describing the polarization-induced renormalization of frontier energy levels of the adsorbates in the static regime. Moreover, it can be reliably applied to simulating laser-induced ultrafast dynamics of molecules through the inclusion of electric fields generated by Fresnel-reflection at the substrate. Depending on the complexity of the underlying layer structure, such reflected fields can assume non-trivial shapes and profoundly affect the dynamics of the photo-excited charge carriers in the molecule. In particular, the interaction with the substrate can give rise to strong delayed fields, which lead to interference effects resembling those of multi-pulse-based spectroscopy. The robustness of the implementation and the above-mentioned features are demonstrated with a number of examples, ranging from intuitive models to realistic systems.

Phonoritons as Hybridized Exciton-Photon-Phonon Excitations in a Monolayer h-BN Optical Cavity

A phonoriton is an elementary excitation that is predicted to emerge from hybridization between exciton, phonon, and photon. Besides the intriguing many-particle structure, phonoritons are of interest as they could serve as functional nodes in devices that utilize electronic, phononic, and photonic elements for energy conversion and thermal transport applications. Although phonoritons are predicted to emerge in an excitonic medium under intense electromagnetic wave irradiation, the stringent condition for their existence has eluded direct observation in solids. In particular, on-resonance, intense pumping schemes have been proposed, but excessive photoexcitation of carriers prevents optical detection. Here, we theoretically predict the appearance of phonoritonic features in monolayer hexagonal boron nitride (h-BN) embedded in an optical cavity. The coherent superposition nature of phonoriton states is evidenced by the hybridization of exciton-polariton branches with phonon replicas that is tunable by the cavity-matter coupling strength. This finding simultaneously provides an experimental pathway for observing the predicted phonoritons and opens a new avenue for tuning materials properties.

Nuclear dynamics of singlet exciton fission in pentacene single crystals

Singlet exciton fission (SEF) is a key process for developing efficient optoelectronic devices. An aspect rarely probed directly, yet with tremendous impact on SEF properties, is the nuclear structure and dynamics involved in this process. Here, we directly observe the nuclear dynamics accompanying the SEF process in single crystal pentacene using femtosecond electron diffraction. The data reveal coherent atomic motions at 1 THz, incoherent motions, and an anisotropic lattice distortion representing the polaronic character of the triplet excitons. Combining molecular dynamics simulations, time-dependent density-functional theory, and experimental structure factor analysis, the coherent motions are identified as collective sliding motions of the pentacene molecules along their long axis. Such motions modify the excitonic coupling between adjacent molecules. Our findings reveal that long-range motions play a decisive part in the electronic decoupling of the electronically correlated triplet pairs and shed light on why SEF occurs on ultrafast time scales.

Ab initio dynamics simulation of laser-induced photodissociation of phenol

We theoretically investigated the photodissociation dynamics of phenol molecules steered by a sequence of temporally shaped femtosecond laser pulses with high intensity and ultrashort duration, via the real-time Time-Dependent Density Functional Theory (rt-TDDFT) combined with a Molecular Dynamics (MD) simulation. The principal findings of this research are that the phenol photodissociation can take place in 50 fs; the bonds broke sequentially; the degree of phenol molecular dissociation has a strong linear correlation with the intensity. For an incident laser being 800 nm-40 fs (wavelength-pulse duration), the threshold intensity is 7 × 10¹⁴ W cm^-2 and the products are hydrogen from OH1 (phenolic hydroxyl group) and C₆H₅O-fragments. More fragments will be found at stronger intensity, shorter wavelengths, and longer pulse duration. More accurately, we estimated the critical values of bond cleavage of an isolated phenol molecule are 1.779 Å for O-H1 and 2.184 Å for C-Hs via Electron localization function (ELF) analysis. The photodissociation of the phenol molecule was triggered via the excitation of electrons and the dissociation process of phenol here is in good agreement with the characteristics of field-assisted dissociation (FAD) theory. Orthogonal tests with an L₉ (3⁴) matrix and threshold intensity decrease tests were conducted to confirm the mechanism. Our research gives an insight into the photodissociation experiment of phenol and provides a simple yet effective way to understand the photochemical experiments of more complex organic pollutants with toxicity.

Time-dependent density functional theory applied to average atom opacity

We focus on studying the opacity of iron, chromium, and nickel plasmas at conditions relevant to experiments carried out at Sandia National Laboratories [J. E. Bailey et al., Nature (London) 517, 56 (2015)NATUAS0028-083610.1038/nature14048]. We calculate the photoabsorption cross sections and subsequent opacity for plasmas using linear-response time-dependent density functional theory (TD-DFT). Our results indicate that the physics of channel mixing accounted for in linear-response TD-DFT leads to an increase in the opacity in the bound-free quasicontinuum, where the Sandia experiments indicate that models underpredict iron opacity. However, the increase seen in our calculations is only in the range of 5%-10%. Further, we do not see any change in this trend for chromium and nickel. This behavior indicates that channel mixing effects do not explain the trends in opacity observed in the Sandia experiments.

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.