LayerPCM: An implicit scheme for dielectric screening from layered substrates

Ab Initio Modeling of Mixed-Dimensional Heterostructures: A Path Forward

Understanding the electronic structure of mixed-dimensional heterostructures is essential for maximizing their application potential. However, accurately modeling such interfaces is challenging due to the complex interplay between the subsystems. We employ a computational framework integrating first-principles methods, including GW, density functional theory (DFT), and the polarizable continuum model, to elucidate the electronic structure of mixed-dimensional heterojunctions formed by free-base phthalocyanines and monolayer molybdenum disulfide. We assess the impact of dielectric screening across various scenarios, from isolated molecules to organic films on a substrate-supported monolayer. Our findings show that while polarization effects cause significant renormalization of molecular energy levels, band energies and alignments in the most relevant setup can be accurately predicted through DFT simulations of the individual subsystems. Additionally, we analyze orbital hybridization, revealing potential pathways for interfacial charge transfer. This study offers new insights into hybrid inorganic/organic interfaces and provides a practical computational protocol suitable for scaled-up studies.

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Charge Localization in Acene Crystals from Ab Initio Electronic Structure

The performance of Koopmans-compliant hybrid functionals in reproducing the electronic structure of organic crystals is tested for a series of acene crystals. The calculated band gaps are found to be consistent with those achieved with the GW method at a fraction of the computational cost and in excellent accord with the experimental results at room temperature, when including the thermal renormalization. The energetics of excess holes and electrons reveals a struggle between polaronic localization and band-like delocalization. The consequences of these results on the transport properties of acene crystals are discussed.

Electronic Dynamics of a Molecular System Coupled to a Plasmonic Nanoparticle Combining the Polarizable Continuum Model and Many-Body Perturbation Theory

The efficiency of plasmonic metallic nanoparticles in harvesting and concentrating light energy in their proximity triggers a wealth of important and intriguing phenomena. For example, spectroscopies are able to reach single-molecule and intramolecule sensitivities, and important chemical reactions can be effectively photocatalyzed. For the real-time description of the coupled dynamics of a molecule's electronic system and of a plasmonic nanoparticle, a methodology has been recently proposed (J. Phys. Chem. C. 120, 2016, 28774-28781) which combines the classical description of the nanoparticle as a polarizable continuum medium with a quantum-mechanical description of the molecule treated at the time-dependent configuration interaction (TDCI) level. In this work, we extend this methodology by describing the molecule using many-body perturbation theory: the molecule's excitation energies, transition dipoles, and potentials computed at the GW/Bethe-Salpeter equation (BSE) level. This allows us to overcome current limitations of TDCI in terms of achievable accuracy without compromising on the accessible molecular sizes. We illustrate the developed scheme by characterizing the coupled nanoparticle/molecule dynamics of two prototype molecules, LiCN and p-nitroaniline.

Solvated Nuclear-Electronic Orbital Structure and Dynamics

Nonadiabatic dynamical processes such as proton-coupled electron transfer and excited state intramolecular proton transfer have been the subject of much research. One of the promising theoretical methods to describe these processes is the nuclear-electronic orbital (NEO) approach. This approach inherently accounts for nuclear quantum effects within quantum chemistry calculations, and it has recently been extended to directly simulate nonadiabatic processes with the development of real-time NEO methods. These processes can also be significantly dependent on the surrounding chemical environment, however, and capturing the effects of the environment is often necessary for analyzing experimentally relevant systems. This work couples the NEO density functional theory and real-time time-dependent density functional theory approaches with solvation through the polarizable continuum model. The effects of this coupling are investigated for ground state properties, solvent-dependent vibrational frequencies, and direct excited state intramolecular proton transfer dynamics.

Polarization Resolved Optical Excitation of Charge-Transfer Excitons in PEN:PFP Cocrystalline Films: Limits of Nonperiodic Modeling

Charge-transfer excitons (CTXs) at organic donor/acceptor interfaces are considered important intermediates for charge separation in photovoltaic devices. Crystalline model systems provide microscopic insights into the nature of such states as they enable microscopic structure-property investigations. Here, we use angular-resolved UV/vis absorption spectroscopy to characterize the CTXs of crystalline pentacene:perfluoro-pentacene (PEN:PFP) films allowing determination of the polarization of this state. This analysis is complemented by first-principles many-body calculations, performed on the three-dimensional PEN:PFP cocrystal, which confirm that the lowest-energy excitation is a CTX. Analogous simulations performed on bimolecular clusters are unable to reproduce this state. We ascribe this failure to the lack of long-range interactions and wave function periodicity in these cluster calculations, which appear to remain a valid tool for modeling properties of organic materials ruled by local intermolecular couplings.