Density matrix analysis and simulation of electronic excitations in conjugated and aggregated molecules

Cavity Manipulation of Attosecond Charge Migration in Conjugated Dendrimers

Dendrimers are branched polymers with wide applications to photosensitization, photocatalysis, photodynamic therapy, photovoltaic conversion, and light sensor amplification. The primary step of numerous photophysical and photochemical processes in many molecules involves ultrafast coherent electronic dynamics and charge oscillations triggered by photoexcitation. This electronic wavepacket motion at short times where the nuclei are frozen is known as attosecond charge migration. We show how charge migration in a dendrimer can be manipulated by placing it in an optical cavity and monitored by time-resolved X-ray diffraction. Our simulations demonstrate that the dendrimer charge migration modes and the character of photoexcited wave function can be significantly influenced by the strong light-matter interaction in the cavity. This presents a new avenue for modulating initial ultrafast charge dynamics and subsequently controlling coherent energy transfer in dendritic nanostructures.

Spectroscopic study of wemple-didomenico empirical formula and taucs model to determine the optical band gap of dye-doped polymer based on chitosan: Common poppy dye as a novel approach to reduce the optical band gap of biopolymer

This study investigates the effect of a natural dye extracted from common poppy (Papaver rhoeas) waste flowers on the optical properties of chitosan (CS) films. The extraction of natural dyes from waste flowers can be considered a new field for research in green chemistry. CS films are flexible and biodegradable but have low optical activity and band gap, limiting their applications in optical devices. The doped CS polymer with different concentrations of Papaver rhoeas dye exhibited enhanced optical properties. Also, 30 % glycerol was added as a plasticizer to omit film brittleness. The FTIR examinations is helpful to propose a mechanism that explains the interaction of the dye with the host polymer. The UV-vis spectroscopic examination establish that the optical characteristics of the films can be modified by adjusting the dye concentration. Furthermore, optical absorption properties are described using the Tauc non-direct transition model, revealing an approximate optical band gap of 1.64 eV. This band gap defines the energy required for electron transitions, elucidating the material's electronic characteristics. The extinction coefficient (k) and refractive index (n) of the CS-doped films' shows a dispersion behavior at visible regions of EM radiation. The Wemple-DiDomenico single oscillator model was used to investigate the n dispersion and determine the oscillator energy equivalent to the optical band gap. Additionally, calculations have been performed on optical dielectric properties and optical conductivity. The Urbach energy was measured and used to detect the structure of the films. The findings underscore the potential applications of these natural dye-doped CS films in eco-friendly materials and optical devices.

A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn

Analysis of electron wavefunction is a key component of quantum chemistry investigations and is indispensable for the practical research of many chemical problems. After more than ten years of active development, the wavefunction analysis program Multiwfn has accumulated very rich functions, and its application scope has covered numerous aspects of theoretical chemical research, including charge distribution, chemical bond, electron localization and delocalization, aromaticity, intramolecular and intermolecular interactions, electronic excitation, and response property. This article systematically introduces the features and functions of the latest version of Multiwfn and provides many representative examples. Through this article, readers will be able to fully understand the characteristics and recognize the unique value of Multiwfn. The source code and precompiled executable files of Multiwfn, as well as the manual containing a detailed introduction to theoretical backgrounds and very rich tutorials, can all be downloaded for free from the Multiwfn website (http://sobereva.com/multiwfn).

Time-Resolved Circular Dichroism in Molecules: Experimental and Theoretical Advances

Following changes in chirality can give access to relevant information on the function or reactivity of molecular systems. Time-resolved circular dichroism (TRCD) spectroscopy proves to be a valid tool to achieve this goal. Depending on the class of molecules, different temporal ranges, spanning from seconds to femtoseconds, need to be investigated to observe such chiroptical changes. Therefore, over the years, several approaches have been adopted to cover the timescale of interest, especially based on pump-probe schemes. Moreover, various theoretical approaches have been proposed to simulate and explain TRCD spectra, including linear and non-linear response methods as well as non-adiabatic molecular dynamics. In this review, an overview on both experimental and theoretical advances in the TRCD field is provided, together with selected applications. A discussion on future theoretical developments for TRCD is also given.

Transient-absorption spectroscopy of dendrimers via nonadiabatic excited-state dynamics simulations

The efficiency of light-harvesting and energy transfer in multi-chromophore ensembles underpins natural photosynthesis. Dendrimers are highly branched synthetic multi-chromophoric conjugated supra-molecules that mimic these natural processes. After photoexcitation, their repeated units participate in a number of intramolecular electronic energy relaxation and redistribution pathways that ultimately funnel to a sink. Here, a model four-branched dendrimer with a pyrene core is theoretically studied using nonadiabatic molecular dynamics simulations. We evaluate excited-state photoinduced dynamics of the dendrimer, and demonstrate on-the-fly simulations of its transient absorption pump-probe (TA-PP) spectra. We show how the evolutions of the simulated TA-PP spectra monitor in real time photoinduced energy relaxation and redistribution, and provide a detailed microscopic picture of the relevant energy-transfer pathways. To the best of our knowledge, this is the first of this kind of on-the-fly atomistic simulation of TA-PP signals reported for a large molecular system.

Ultrafast Excited-State Energy Transfer in Phenylene Ethynylene Dendrimer: Quantum Dynamics with the Tensor Network Method

Photoinduced excited-state energy transfer (EET) processes play an important role in solar energy conversions. Owing to their excellent photoharvesting and exciton-transport properties, phenylene ethynylene (PE) dendrimers display great potential for improving the efficiency of solar cells. In this work, we investigated the intramolecular EET dynamics in a dendrimer composed of two linear PE units (2-ring and 3-ring) using a fully quantum description based on the tensor network method. We first constructed a diabatic model Hamiltonian based on the electronic structure calculations. Using this diabatic vibronic coupling model, we tried to obtain the main features of the EET dynamics in terms of the several diabatic models with different numbers of vibrational modes (from 4 modes to 129 modes) and to explore the corresponding vibronic coupling interactions. The results show that the EET in this PE dendrimer is ultrafast. Four modes of A' symmetry play dominant roles in the dynamics; the remaining 86 modes of A' symmetry can dampen the electronic coherence; and the modes of A″ symmetry do not exhibit significant influence on the EET process. Overall, the first-order intrastate vibronic coupling terms show the dominant role in the EET dynamics, while the second-order intrastate vibronic coupling terms cause damping of the electronic coherence and slow down the overall EET process. This work provides a microscopic understanding of the EET dynamics in PE dendrimers.

Charge transfer excitons and directional fluorescence of in-plane lateral MoSe2-WSe2 heterostructures

In this article, the linear and nonlinear optical properties of in-plane lateral MoSe2-WSe2 heterostructures are theoretically investigated. The polarization-dependent strongest optical absorption in one-photon absorption occurs in charge transfer excited states, where electrons transfer from WSe2 to MoSe2. This phenomenon is supported by the LUMO (lowest unoccupied molecular orbital) and HUMO (highest occupied molecular orbital) imaging obtained through scanning tunneling microscopy. The charge difference density and transition density matrix are used to interpret the electronic transitions, and these interpretations rely on the concept of transition density. The optical properties of two-photon absorption in its nonlinear optical process are significantly different from the excitation in one-photon absorption, where the strongest optical absorption is contributed from direct transition from the ground state to the final state without going through an intermediate excited state, due to the very large difference of permanent dipole moments between the excited and ground states. Our results also reveal directional fluorescence and physical mechanism of in-plane lateral MoSe2-WSe2 heterostructures. Our work can provide insights into the physical mechanism of the optical properties of in-plane lateral MoSe2-WSe2 heterostructures.

Coherent Transient Localization Mechanism of Interchain Exciton Transport in Regioregular P3HT: A Quantum-Dynamical Study

Transient localization has been proposed as a transport mechanism in organic materials, for both charge carriers and excitons. Here, we characterize a quantum coherent transient localization mechanism using full quantum simulations of an H-aggregated model system representative of regioregular poly(3-hexylthiophene) (rrP3HT). A Frenkel-Holstein Hamiltonian parametrized from first principles is considered, including local high-frequency modes and anharmonic, site-correlated interchain modes. Quantum-dynamical calculations are carried out using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a 13-site system with 195 vibrational modes, under periodic boundary conditions. It is shown that temporary localization of exciton polarons alternates with resonant transfer driven by interchain modes. While the transport process is mainly determined by exciton-polarons at the low-energy band edge, persistent coupling with the excitonic manifold is observed, giving rise to a nonadiabatic excitonic flux. This elementary transport mechanism remains preserved for limited static disorder and gives way to Anderson localization when the static disorder becomes dominant.

Visible Light Induced Exciton Dynamics and trans-to-cis Isomerization in Azobenzene Aggregates: Insights from Surface Hopping/Semiempirical Configuration Interaction Molecular Dynamics Simulations

Assemblies of photochromic molecules feature exciton states, which govern photochemical and photophysical processes in multichromophoric systems. Understanding the photoinduced dynamics of the assemblies requires nonadiabatic treatment involving multiple exciton states and numerous nuclear degrees of freedom, thus posing a challenge for simulations. In this work, we address this challenge for aggregates of azobenzene, a prototypical molecular switch, performing on-the-fly surface hopping calculations combined with semiempirical configuration interaction electronic structure and augmented with transition density matrix analysis to characterize exciton evolution. Specifically, we consider excitation of azobenzene tetramers in the nπ* absorption band located in the visible (blue) part of the electromagnetic spectrum, thus extending our recent work on dynamics after ππ* excitation corresponding to the ultraviolet region [Titov, J. Phys. Chem. C2023, 127, 13678-13688]. We find that the nπ* excitons, which are initially strongly localized by ground-state conformational disorder, undergo further (very strong) localization during short-time photodynamics. This excited-state localization process is extremely ultrafast, occurring within the first 10 fs of photodynamics. We observe virtually no exciton transfer of the localized excitons in the nπ* manifold. However, the transfer may occur via secondary pathways involving ππ* states or the ground state. Moreover, we find that the nπ* quantum yields of the trans-to-cis isomerization are reduced in the aggregated state.

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.

Efficient Quasiparticle Determination beyond the Diagonal Approximation via Random Compression

Calculations of excited states in the Green's function formalism often invokes the diagonal approximation, in which the quasiparticle states are taken from a mean-field calculation. In this paper, we extend the stochastic approaches applied in the many-body perturbation theory and overcome this limitation for large systems in which we are interested in a small subset of states. We separate the problem into a core subspace whose coupling to the remainder of the system environment is stochastically sampled. This method is exemplified on computing hole injection energies into CO2 on an extended gold surface with nearly 3000 electrons. We find that in the extended system the size of the problem can be compressed up to 95% using stochastic sampling. This result provides a way forward for self-consistent stochastic methods and determination of Dyson orbitals in large systems.

Design and Exploration of Benzene Like Azobis Triazoles for Long-range Push-Pull Photo-Switching Attributes

This research paper presents a comprehensive study on the design and photovoltaic parameters of azobenzene type 24 photo switches (PSs) of triazole by density functional theory (DFT). The focus was on investigating how to create a long-range push-pull effect of different substituents on the PS properties for their application in photovoltaics by further substituent decoration. Their range of values for the maximum wavelength (λmax) ranged 315-556 nm while their HOMO-LUMO energies (Egaps) were 0.57-6.35eV. The stability of the PS was evaluated by measuring hardness (η) and softness (σ) values. Additionally, photovoltaic parameters such as open-circuit voltage (Voc), short-circuit current density (Jsc), fill factor (FF), and maximum power (Pmax) were calculated to assess the performance of the PS as photovoltaic materials. The results revealed that PSs 6 exhibited promising photovoltaic parameters to include Voc values ranging from 0.4-1.63eV, FF values ranging from 0.5438-0.929, Jsc values ranging from 19.27-50.75 mA/cm2, and Pmax values ranging from 14.72-75.91W. This indicates its potential as an efficient light-harvesting material for photovoltaic applications. Moreover, this study presents a pioneering investigation on the correlation between rotational velocity (R) and Mayer bond index (MBI) for the first time. The findings revealed a significant correlation between R and MBI, providing valuable insights into the structural dynamics of the PS. This novel finding opens up new avenues for understanding the structural dynamics of PS and their potential applications in various fields, including photovoltaics. The study provides valuable insights into the structure-property relationships of azobenzene-based PS and their suitability for photovoltaic devices. Further investigations are warranted to optimize the design of the PS, enhance their photovoltaic performance, and explore the underlying mechanisms of the correlation between R and MBIs.

Theory and modeling of light-matter interactions in chemistry: current and future

Light-matter interaction not only plays an instrumental role in characterizing materials' properties via various spectroscopic techniques but also provides a general strategy to manipulate material properties via the design of novel nanostructures. This perspective summarizes recent theoretical advances in modeling light-matter interactions in chemistry, mainly focusing on plasmon and polariton chemistry. The former utilizes the highly localized photon, plasmonic hot electrons, and local heat to drive chemical reactions. In contrast, polariton chemistry modifies the potential energy curvatures of bare electronic systems, and hence their chemistry, via forming light-matter hybrid states, so-called polaritons. The perspective starts with the basic background of light-matter interactions, molecular quantum electrodynamics theory, and the challenges of modeling light-matter interactions in chemistry. Then, the recent advances in modeling plasmon and polariton chemistry are described, and future directions toward multiscale simulations of light-matter interaction-mediated chemistry are discussed.

Promising impact of push-pull configuration into designed octacyclic naphthalene-based organic scaffolds for nonlinear optical amplitudes: a quantum chemical approach

In opto-electronics, non-fullerene (NF) derivatives are regarded as efficient non-linear optical (NLO) materials. The present investigation was based on designing NF naphthalene-based derivatives (PCMD1-D9) with D-π-A configuration from PCMR. DFT analysis at M06/6-311G (d,p) level was accomplished to explore the photonic behavior of PCMD1-D9 compounds. Various kind of analysis like; UV-Vis, density of state (DOS), natural bond orbitals (NBOs), transition density matrix (TDM) and frontier molecular orbitals (FMOs) analyses were accomplished to understand the NLO properties of said chromophores. The configuration change led to considerable charge distribution over highest occupied and lowest unoccupied molecular orbitals with minimum band difference. The energy gap trend for all the entitled compounds was observed as; PCMD8 < PCMD5 = PCMD9 < PCMD6 < PCMD7 < PCMD4 < PCMD3 < PCMD2 < PCMD1 with the least band gap of 2.048 eV in PCMD8 among all the compounds. The UV-Visible spectrum of the entitled chromophores manifested high values of λmax in derivatives contrary to PCMR. Additionally, NBO findings explored effective intramolecular charge transfer and maximum energy of stabilization (34.31 kcal/mol) for PCMD8 chromophore. The highest linear polarizability (<α>) and dipole moment (µtot) values were exhibited by PCMD5 at 2.712 × 10-22. and 1.995 × 10-17 esu, respectively. PCMD8 push-pull configured molecular entity exhibited highest first hyper-polarizability (βtot) at 4.747 × 10-27 esu and second hyper-polarizability at 6.867 × 10-32 esu. Overall, all the formulated chromophores exhibited significant NLO results contrary to PCMR. Hence, through this structural tailoring via various acceptors, effective NLO materials were obtained for optoelectronic applications.

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort-design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

Triplet Generation Through Singlet Fission in Metal-Organic Framework: An Alternative Route to Inefficient Singlet-Triplet Intersystem Crossing

High quantum yield triplets, populated by initially prepared excited singlets, are desired for various energy conversion schemes in solid working compositions like porous MOFs. However, a large disparity in the distribution of the excitonic center of mass, singlet-triplet intersystem crossing (ISC) in such assemblies is inhibited, so much so that a carboxy-coordinated zirconium heavy metal ion cannot effectively facilitate the ISC through spin-orbit coupling. Circumventing this sluggish ISC, singlet fission (SF) is explored as a viable route to generating triplets in solution-stable MOFs. Efficient SF is achieved through a high degree of interchromophoric coupling that facilitates electron super-exchange to generate triplet pairs. Here we show that a predesigned chromophoric linker with extremely poor ISC efficiency (kISC ) butE S 1 ≥ 2 E T 1 ${{E}_{{S}_{1}}\ge {2E}_{{T}_{1}}}$ form triplets in MOF in contrast to the frameworks that are built from linkers with sizable kISC butE S 1 ≤ 2 E T 1 ${{E}_{{S}_{1}}\le {2E}_{{T}_{1}}}$ . This work opens a new photophysical and photochemical avenue in MOF chemistry and utility in energy conversion schemes.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Investigating Molecular Exciton Polaritons Using Ab Initio Cavity Quantum Electrodynamics

Coupling molecules to the quantized radiation field inside an optical cavity creates a set of new photon-matter hybrid states called polariton states. We combine electronic structure theory with quantum electrodynamics (QED) to investigate molecular polaritons using ab initio simulations. This framework joins unperturbed electronic adiabatic states with the Fock state basis to compute the eigenstates of the QED Hamiltonian. The key feature of this "parametrized QED" approach is that it provides the exact molecule-cavity interactions, limited by only approximations made in the electronic structure. Using time-dependent density functional theory, we demonstrated comparable accuracy with QED coupled cluster benchmark results for predicting potential energy surfaces in the ground and excited states and showed selected applications to light-harvesting and light-emitting materials. We anticipate that this framework will provide a set of general and powerful tools that enable direct ab initio simulation of exciton polaritons in molecule-cavity hybrid systems.

On-the-Fly Nonadiabatic Dynamics Simulations of Single-Walled Carbon Nanotubes with Covalent Defects

Single-walled carbon nanotubes (SWCNTs) with covalent surface defects have been explored recently due to their promise for use in single-photon telecommunication emission and in spintronic applications. The all-atom dynamic evolution of electrostatically bound excitons (the primary electronic excitations) in these systems has only been loosely explored from a theoretical perspective due to the size limitations of these large systems (>500 atoms). In this work, we present computational modeling of nonradiative relaxation in a variety of SWCNT chiralities with single-defect functionalizations. Our excited-state dynamics modeling uses a trajectory surface hopping algorithm accounting for excitonic effects with a configuration interaction approach. We find a strong chirality and defect-composition dependence on the population relaxation (varying over 50-500 fs) between the primary nanotube band gap excitation E₁₁ and the defect-associated, single-photon-emitting E₁₁* state. These simulations give direct insight into the relaxation between the band-edge states and the localized excitonic state, in competition with dynamic trapping/detrapping processes observed in experiment. Engineering fast population decay into the quasi-two-level subsystem with weak coupling to higher-energy states increases the effectiveness and controllability of these quantum light emitters.

Rational design of non-fullerene acceptors via side-chain and terminal group engineering: a computational study

We investigated the optoelectronic and photovoltaic properties of three types of acceptor-donor-acceptor-based non-fullerene acceptor (NFA) molecules for organic solar cell (OSC) applications. Density functional theory and its time-dependent variant were employed to compute the quadrupole moment perpendicular to the π-system (Q₂₀), open circuit voltage (V_OC), and other relevant solar cell parameters. The role of functionalization in the acceptor unit on the overall device performance was explored by incorporating halogen and methoxy-based electron-withdrawing groups. The electronegativity differences between the halogen atoms and the methoxy group demonstrated contrasting effects on the energy levels, molecular orbitals, and absorption maximum. We observed a trade-off between short-circuit current (J_SC) and V_OC, which was further substantiated by an inverse correlation between Q₂₀ and V_OC. We found an optimum value of Q₂₀ in the range of 80 to 130 ea₀² to achieve an optimized solar cell performance. Among the designed systems, Se-derived NFAs with a small band gap, red-shifted absorption maximum, high-oscillator strength, small exciton binding energy, and optimum Q₂₀ turned out to be potential candidates for future applications. These criteria can be generalized to design and screen next-generation non-fullerene acceptors to achieve improved OSC performance.

Alteration of the central core of a DF-PCIC chromophore to boost the photovoltaic applications of non-fullerene acceptor based organic solar cells

Modifying the central core is a very efficient strategy to boost the performance of non-fullerene acceptors. Herein five non-fullerene acceptors (M1-M5) of A-D-D'-D-A type were designed by substituting the central acceptor core of the reference (A-D-A'-D-A type) with different strongly conjugated and electron donating cores (D') to enhance the photovoltaic attributes of OSCs. All the newly designed molecules were analyzed through quantum mechanical simulations to compute their optoelectronic, geometrical, and photovoltaic parameters and compare them to the reference. Theoretical simulations of all the structures were carried out through different functionals with a carefully selected 6-31G(d,p) basis set. Absorption spectra, charge mobility, dynamics of excitons, distribution pattern of electron density, reorganization energies, transition density matrices, natural transition orbitals and frontier molecular orbitals, respectively of the studied molecules were evaluated at this functional. Among the designed structures in various functionals, M5 showed the most improved optoelectronic properties, such as the lowest band gap (2.18 e V), highest maximum absorption (720 nm), and lowest binding energy (0.46 eV) in chloroform solvent. Although the highest photovoltaic aptitude as acceptors at the interface was perceived to be of M1, its highest band gap and lowest absorption maxima lowered its candidature as the best molecule. Thus, M5 with its lowest electron reorganization energy, highest light harvesting efficiency, and promising open-circuit voltage (better than the reference), amongst other favorable features, outperformed the others. Conclusively, each evaluated property commends the aptness of designed structures to augment the power conversion efficiency (PCE) in the field of optoelectronics in one way or another, which reveals that a central un-fused core having an electron-donating capability with terminal groups being significantly electron withdrawing, is an effective configuration for the attainment of promising optoelectronic parameters, and thus the proposed molecules could be utilized in future NFAs.

Optical Physical Mechanisms of Helicene Carbon Nanohoop with Möbius Topology

Optical and spectral properties of carbon nanohoop with Möbius topology is of great interest in nano-science and nano-technology. And it can be imagined that it has a lot of unexpected potential application prospects. However, theoretical calculations based on some figure-of-eight helicene carbon nanohoop with Möbius topology are still insufficient. Therefore, in this paper, we theoretically study the optical and spectral properties of figure-of-eight helicene carbon nanohoop with Möbius topology. Optical and spectral properties are analyzed with visualization method of transition density matrix and charge density difference, which reveal the unique characterization of carbon nanohoop with Möbius topology. Our results can not only deepen the understanding of the optical physical mechanisms of the nanorings with mobius carbons, but also provide deeper insight on optical properties and potential design on optical nanodevices.

Molecular Chirality and Its Monitoring by Ultrafast X-ray Pulses

Major advances in X-ray sources including the development of circularly polarized and orbital angular momentum pulses make it possible to probe matter chirality at unprecedented energy regimes and with Ångström and femtosecond spatiotemporal resolutions. We survey the theory of stationary and time-resolved nonlinear chiral measurements that can be carried out in the X-ray regime using tabletop X-ray sources or large scale (XFEL, synchrotron) facilities. A variety of possible signals and their information content are discussed.

Modeling Photolytic Decomposition of Energetically Functionalized Dodecanes

The photolytic stability of explosives and energetic functional groups is of importance for those who regularly handle or are exposed to explosives in typical environmental conditions. This study models the photolytic degradation of dodecane substituted with various energetic functional groups: azide, nitro, nitrate ester, and nitramine. For the studied molecules, it was found that excitons localize on the energetic functional group, no matter where they were initially formed, and thus, the predominant degradation pathway involves the degradation of the energetic functional group. The relative trends for both 4 and 8 eV excitation energies followed with what is expected from the relative stability of the energetic functional groups to thermal and sub-shock degradation. The one notable exception was the azide functional group; more work should be done to further understand the photolytic effects on the azide functional group.

Designing dibenzosilole core based, A2-π-A1-π-D-π-A1-π-A2 type donor molecules for promising photovoltaic parameters in organic photovoltaic cells

In this research work, four new molecules from the π-A-π-D-π-A-π type reference molecule "DBS-2PP", were designed for their potential application in organic solar cells by adding peripheral A2 acceptors to the reference. Under density functional theory, a comprehensive theoretical investigation was conducted to examine the structural geometries, along with the optical and photovoltaic parameters; comprising frontier molecular orbitals, density of states, light-harvesting effectiveness, excitation, binding, and reorganizational energies, molar absorption coefficient, dipole moment, as well as transition density matrix of all the molecules under study. In addition, some photo-voltaic characteristics (open circuit photo-voltage and fill factor) were also studied for these molecules. Although all the developed compounds (D1-D4) surpassed the reference molecule in the attributes mentioned above, D4 proved to be the best. D4 possessed the narrowest band-gap, as well as the highest absorption maxima and dipole moment of all the molecules in both the evaluated phases. Moreover, with PC61BM as the acceptor, D4 showed the maximum V OC and FF values. Furthermore, while D3 had the greatest hole mobility owing to its lowest value of hole reorganization energy, D4 exhibited the maximum electron mobility due to its lowermost value of electron reorganization energy. Overall, all the chromophores proposed in this study showed outstanding structural, optical, and photovoltaic features. Considering this, organic solar cell fabrication can be improved by using these newly derived donors at the donor-acceptor interfaces.

Effect of conformational disorder on exciton states of an azobenzene aggregate

Azobenzene is a prototypical molecular photoswitch, widely used to trigger a variety of transformations at different length scales. In systems like self-assembled monolayers or micelles, azobenzene chromophores may interact with each other, which gives rise to the emergence of exciton states. Here, using first-principles calculations, we investigate how conformational disorder (induced, e.g., by thermal fluctuations) affects localization of these states, on an example of an H-type azobenzene tetramer. We find that conformational disorder leads to (partial) exciton localization on a single-geometry level, whereas ensemble-averaging results in a delocalized picture. The ππ* and nπ* excitons at high and low temperatures are discussed.

Semi-Empirical Calculation of Bodipy Aggregate Spectroscopic Properties through Direct Sampling of Configurational Ensembles

Efficient prediction of the aggregation-induced callback of organic chromophores for utilization in molecular sensorics is a desirable development goal in modern computational chemistry. Dye aggregates are complicated to study when utilizing conventional quantum chemistry approaches, since they are usually composed of too many atoms to be effectively analyzed, even with high-throughput parallel systems. Here, we present a successful attempt to develop a protocol to assess the spectroscopic changes happening in BODIPY dyes upon aggregation from the first principles utilizing extended tight-binding (XTB) and Zerner's intermediate neglect of differential overlap (ZINDO) Hamiltonians. The developed sampling technique for aggregate configurational space scanning was found to be sufficient to both reproduce peculiarities and justify experimental data on the spectroscopic behavior of chromophore aggregates. The sTDA, sTD-DFT (GFN2-XTB) and CIS (ZINDO) approaches were assessed, and then sources of errors and benefits were outlined. Importantly, our goal was to keep any of the mentioned calculations within a computational cost feasible for a single workstation, whereas scaling was possible at any point in time. Finally, several aggregate structures were investigated in the external field to try to achieve distributions similar to the ones observed in the electrostatic potential of the air-water interface to assess the borderlines of practical applicability of the suggested scheme.

Infinitene: Computational Insights from Nonadiabatic Excited State Dynamics

Progress in organic synthesis opens exploration of a rich diversity of molecules with interesting new structural topologies. This is the case of a recently synthesized helically twisted figure-eight molecule coined infinitene. The molecule belongs to a numerous family of looped polyarenes, where the degree of π-conjugation is controlled by high strain energies and steric hindrances. A particular balance of these ingredients leads to unusual optoelectronic properties potentially suitable for a range of applications in nanoelectronics and photonics. Due to its recent discovery, the photophysical properties of infinitene remain unexplored. In this Letter, atomistic nonadiabatic excited state molecular dynamics modeling unveils unique features of intramolecular electronic and vibrational energy relaxation and redistribution that take place after molecular photoexcitation. Our results detail relationships between optical and electronic properties providing useful knowledge for future molecular designs related to infinitene.

Adiabatic Excited-State Molecular Dynamics with an Explicit Solvent: NEXMD-SANDER Implementation

We present a method to link the Nonadiabatic EXcited-state Molecular Dynamics (NEXMD) package to the SANDER package supplied by AMBERTOOLS to provide excited-state adiabatic quantum mechanics/molecular mechanics (QM/MM) simulations. NEXMD is a computational package particularly developed to perform simulations of the photoexcitation and subsequent nonadiabatic electronic and vibrational energy relaxation in large multichromophoric conjugated molecules involving several coupled electronic excited states. The NEXMD-SANDER exchange has been optimized in order to achieve excited-state adiabatic dynamics simulations of large conjugated materials in a QM/MM environment, such as an explicit solvent. Dynamics of a substituted polyphenylene vinylene oligomer (PPV3-NO₂) in vacuum and different explicit solvents has been used as a test case by performing comparative analysis of changes in its optical spectrum, state-dependent conformational changes, and quantum bond orderings. The method has been tested and compared with respect to previous implicit solvent implementations. Also, the impact on the expansion of the QM region by including a variable number of solvent molecules has been analyzed. Altogether, these results encourage future implementations of NEXMD simulations using the same combination of methods.

Extending machine learning beyond interatomic potentials for predicting molecular properties

Machine learning (ML) is becoming a method of choice for modelling complex chemical processes and materials. ML provides a surrogate model trained on a reference dataset that can be used to establish a relationship between a molecular structure and its chemical properties. This Review highlights developments in the use of ML to evaluate chemical properties such as partial atomic charges, dipole moments, spin and electron densities, and chemical bonding, as well as to obtain a reduced quantum-mechanical description. We overview several modern neural network architectures, their predictive capabilities, generality and transferability, and illustrate their applicability to various chemical properties. We emphasize that learned molecular representations resemble quantum-mechanical analogues, demonstrating the ability of the models to capture the underlying physics. We also discuss how ML models can describe non-local quantum effects. Finally, we conclude by compiling a list of available ML toolboxes, summarizing the unresolved challenges and presenting an outlook for future development. The observed trends demonstrate that this field is evolving towards physics-based models augmented by ML, which is accompanied by the development of new methods and the rapid growth of user-friendly ML frameworks for chemistry.

Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs. H-aggregate models

We report on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature, using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used, along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport, by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 x 10-3 cm2/s, showing a non-monotonous increase with temperature.

A Theoretical Evaluation of the Efficiencies of Metal-Free 1,3,4-Oxadiazole Dye-Sensitized Solar Cells: Insights from Electron–Hole Separation Distance Analysis

Herein, some novel metal-free 1,3,4-oxadiazole compounds O1–O7 were evaluated for their photovoltaic properties using density functional theory (DFT) and time-dependent density functional theory (TD-DFT) calculations to determine if they can serve as metal-free organic dyes in the use of dye-sensitized solar cells (DSSCs). To understand the trends in the relative efficiencies of the investigated compounds as dyes in DSSCs, their electron contributions, hole contributions, and electron–hole overlaps for each respective atom and fragment within the molecule were analyzed with a particular focus on the electron densities on the anchoring segments. As transition density matrices (TDM) provide details about the departure of each electron from its corresponding hole during excitations, which results in charge transfer (CT), the charge separation distance (Δr) between the electron and its corresponding hole was studied, in addition to the degree of electron–hole overlap (Λ). The latter, single-point excitation energy of each electron, the percentage electron contribution to the anchoring segments of each compound, the incident-photon-conversion-efficiency (IPCE), charge recombination, light harvesting efficiency (LHE), electron injection (Φinj), and charge collection efficiency (ncollect) were then compared to Δr to determine whether the expected relationships hold. Moreover, parameters such as diffusion constant (Dπ) and electron lifetime (t), amongst others, were also used to describe electron excitation processes. Since IPCE is the key parameter in determining the efficiency, O3 was found to be the best dye due to its highest value.

Impact of π-Conjugation Length on the Excited-State Dynamics of Star-Shaped Carbazole-π-Triazine Organic Chromophores

Correlating star-shaped donor-bridge-acceptor (DBA) molecular structures with intramolecular charge transfer (ICT) and intersystem crossing (ISC) is essential to their application in photocatalysis, photovoltaics, and organic light-emitting diodes (OLEDs). In this work, we report a systematic photophysical study on a series of star-shaped triazine-phenylene-carbazole DBA molecules with 0, 1, and 2 bridging phenylene units (pTCT-0P, pTCT-1P, pTCT-2P). Using a combination of steady-state and time-resolved spectroscopy with time-dependent density functional theory (TDDFT), we find that the bridge length can strongly impact the structural conformation, ICT, and ISC. Global target analysis of the time-resolved spectroscopy reveals that pTCT-0P has the most favorable ISC rate of 1.96 × 10^-4 ps^-1, which is competitive with a singlet relaxation rate of 1.92 × 10^-4 ps^-1. TDDFT aligns with spectroscopic results within an order of magnitude, predicting an ISC rate of 2.1 × 10^-5 ps^-1 and revealing that the donor/acceptor orthogonalization concomitantly suppresses singlet exciton recombination and lowers the singlet-triplet energy gap. The new fundamental insights gained from this work will help design the next generation of star-shaped DBA-type molecules for photocatalytic and photoelectronic applications.

Ultrafast coherent photoexcited dynamics in a trimeric dendrimer probed by X-ray stimulated-Raman signals

The photoinduced ultrafast coherent inter-chromophore energy redistribution in a triarylamine trimer is explored using nonadiabatic excited state molecular dynamics followed by simulations of X-ray Raman signals. The nitrogencentered system ensures strong interchromophore interactions and, thus, the presence of coherences. Nevertheless, the multitude of non-deterministic photoinduced pathways during the ultrafast inter-branch migration of the excitation results in random confinement on some branches and, therefore, spatial exciton scrambling and loss of phase information at long times. We show that the vibronic coherence dynamics evolving into the incoherent scrambling mechanism on ultrafast 50 fs timescale, is accurately probed by the TRUECARS X-ray stimulated Raman signal. In combination with previous results, where the technique has revealed long-lived coherences in a rigid heterodimer, the signal is most valuable for detecting ultrafast molecular coherences or their absence. We demonstrate that X-ray Raman spectroscopy is a useful tool in the chemical design of functional molecular building blocks.

The photoinduced ultrafast coherent inter-chromophore energy redistribution in a triarylamine trimer is explored using nonadiabatic excited state molecular dynamics followed by simulations of X-ray Raman signals.

Impact of packing arrangement on the optical properties of C60 cluster aggregates

The packing arrangement of organic π-conjugated molecules in a nanoscale material can have a strong impact on their optical properties. Here, using real-time-propagation time dependent density functional theory (rt-TDDFT) calculations with the support of transition contribution maps, we study how modifications in the packing arrangement (cubic-like and chain-like aggregates composed of eight C60 molecules) and packing density (assembled at close distances with center-to-center inter-fullerene distances (d) varying from 9 Å to 11 Å) of C60 molecules affect the optical properties of cluster aggregates. The important conclusions drawn from this work are summarized as follows. For d = 9 Å, the charge transfer excitons produced by cubic and chain-like C60 cluster aggregates have highly different optical characteristics, as evidenced by the transition contribution maps. On the other hand, for d = 10 Å and 11 Å, both kinds of aggregates produce qualitatively similar optical features with the emergence of Wannier-like delocalized excitons having distinct degrees of localization and spatial distribution. The theoretical findings in this study elucidate the optical excitations in C60 cluster aggregates and could help in the design of more efficient organic devices.

Superradiance and Directional Exciton Migration in Metal-Organic Frameworks

Crystalline metal-organic frameworks (MOFs) are promising synthetic analogues of photosynthetic light-harvesting complexes (LHCs). The precise assembly of linkers (organic chromophores) around the topology-defined pores offers the evolution of unique photophysical behaviors that are reminiscence of LHCs. These include MOF excited states with photoabsorbed energy that is spatially dispersed over multiple linkers defining the molecular excitons. The multilinker molecular excitons display superradiance─a hallmark of coupled oscillators seen in LHCs─with radiative rate constant (k_rad) exceeding that of a single linker. Our theoretical model and experimental results on three zirconium MOFs, namely, PCN-222(Zn), NU-1000, and SIU-100, with similar topology but varying linkers suggest that the size of such molecular excitons depends on the electronic symmetry of the linker. This multilinker exciton model effectively predicts the energy transfer rate constant; corresponding single-step exciton hopping time, ranging from a few picoseconds in SIU-100 and NU-1000 to a few hundreds of picoseconds in PCN-222(Zn), matches well with the experimental data. The model also predicts the anisotropy of exciton displacement with preferential migration along the crystallographic c-axis. Overall, these findings establish various missing links defining the exciton size and dynamics in MOF-assembled linkers. The understandings will provide design principles, especially, positioning the catalysts or electrode relative to the linker orientation for low-density solar energy conversion systems.

Rational designs of structurally similar TADF and HLCT emitters with benzo- or naphtho-carbazole units as electron donors

Structurally similar D–A type molecules with the combination of benzo- or naphtho-carbazole units as electron donors and tunable electron acceptors with different electron-withdrawing ability are designed to realize HLCT and TADF emissions.

Physical mechanisms of photoinduced charge transfer in neutral and charged donor-acceptor systems

In this paper, we provide visualization methods to reveal the physical mechanisms of photoinduced charge transfer in neutral and charged donor-acceptor systems. These visualization methods use the charge density difference and transition density matrix, which can promote deeper understanding of photoinduced charge transfer in donor-acceptor systems.

An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics

We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

Excitation Energy Transfer between bodipy Dyes in a Symmetric Molecular Excitonic Seesaw

We examine the redistribution of energy between electronic and vibrational degrees of freedom that takes place between a π-conjugated oligomer, a phenylene-butadiynylene, and two identical boron-dipyrromethene (bodipy) end-caps using femtosecond transient absorption spectroscopy, single-molecule spectroscopy, and nonadiabatic excited-state molecular dynamics (NEXMD) modeling techniques. The molecular structure represents an excitonic seesaw in that the excitation energy on the oligomer backbone can migrate to either one end-cap or the other, but not to both. The NEXMD simulations closely reproduce the characteristic time scale for redistribution of electronic and vibrational energy of 2.2 ps and uncover the vibrational modes contributing to the intramolecular relaxation. The calculations indicate that the dihedral angle between the bodipy dye and the oligomer change upon excitation of the oligomer. Single-molecule experiments reveal a difference in photoluminescence lifetime of the bodipy dyes depending on whether they are excited by direct absorption or by redistribution of energy from the backbone. This difference in lifetime may be attributed to the difference in dihedral angle. The simulations also suggest that a strong coupling can occur between the two end-caps, giving rise to a reversible shuttling of excitation energy between them. Strong coupling should lead to a pronounced loss in polarization memory of the fluorescence since the oligomer backbone tends to be slightly distorted and the two bodipy transition dipoles have different orientations. A sensitive single-molecule technique is presented to test for such coupling. However, although redistribution of electronic and vibrational energy between the end-caps can occur, it appears to be unidirectional and irreversible, suggesting that an additional localization mechanism is at play which is, as yet, not fully accounted for in the simulations.

Nonadiabatic dynamics with spin-flip vs linear-response time-dependent density functional theory: A case study for the protonated Schiff base C5H6NH2. J Chem Phys 2021;155:124111. [PMID: 34598550 DOI: 10.1063/5.0062757]

Nonadiabatic trajectory surface hopping simulations are reported for trans-C₅H₆NH₂ ⁺, a model of the rhodopsin chromophore, using the augmented fewest-switches algorithm. Electronic structure calculations were performed using time-dependent density functional theory (TDDFT) in both its conventional linear-response (LR) and its spin-flip (SF) formulations. In the SF-TDDFT case, spin contamination in the low-lying singlet states is removed by projecting out the lowest triplet component during iterative solution of the TDDFT eigenvalue problem. The results show that SF-TDDFT qualitatively describes the photoisomerization of trans-C₅H₆NH₂ ⁺, with favorable comparison to previous studies using multireference electronic structure methods. In contrast, conventional LR-TDDFT affords qualitatively different photodynamics due to an incorrect excited-state potential surface near the Franck-Condon region. In addition, the photochemistry (involving pre-twisting of the central double bond) appears to be different for SF- and LR-TDDFT, which may be a consequence of different conical intersection topographies afforded by these two methods. The present results contrast with previous surface-hopping studies suggesting that the LR-TDDFT method's incorrect topology around S₁/S₀ conical intersections is immaterial to the photodynamics.

Coupling between Emissive Defects on Carbon Nanotubes: Modeling Insights

Covalent functionalization of single-walled carbon nanotubes (SWCNTs) with organic molecules results in red-shifted emissive states associated with sp³-defects in the tube lattice, which facilitate their improved optical functionality, including single-photon emission. The energy of the defect-based electronic excitations (excitons) depends on the molecular adducts, the configuration of the defect, and concentration of defects. Here we model the interactions between two sp³-defects placed at various distances in the (6,5) SWCNT using time-dependent density functional theory. Calculations reveal that these interactions conform to the effective model of J-aggregates for well-spaced defects (>2 nm), leading to a red-shifted and optically allowed (bright) lowest energy exciton. H-aggregate behavior is not observed for any defect orientations, which is beneficial for emission. The splitting between the lowest energy bright and optically forbidden (dark) excitons and the pristine excitonic band are controlled by the single-defect configurations and their axial separation. These findings enable a synthetic design strategy for SWCNTs with tunable near-infrared emission.

Predicting phosphorescence energies and inferring wavefunction localization with machine learning

Phosphorescence is commonly utilized for applications including light-emitting diodes and photovoltaics. Machine learning (ML) approaches trained on ab initio datasets of singlet-triplet energy gaps may expedite the discovery of phosphorescent compounds with the desired emission energies. However, we show that standard ML approaches for modeling potential energy surfaces inaccurately predict singlet-triplet energy gaps due to the failure to account for spatial localities of spin transitions. To solve this, we introduce localization layers in a neural network model that weight atomic contributions to the energy gap, thereby allowing the model to isolate the most determinative chemical environments. Trained on the singlet-triplet energy gaps of organic molecules, we apply our method to an out-of-sample test set of large phosphorescent compounds and demonstrate the substantial improvement that localization layers have on predicting their phosphorescence energies. Remarkably, the inferred localization weights have a strong relationship with the ab initio spin density of the singlet-triplet transition, and thus infer localities of the molecule that determine the spin transition, despite the fact that no direct electronic information was provided during training. The use of localization layers is expected to improve the modeling of many localized, non-extensive phenomena and could be implemented in any atom-centered neural network model.

On the Low-Lying Electronically Excited States of Azobenzene Dimers: Transition Density Matrix Analysis

Azobenzene-containing molecules may associate with each other in systems such as self-assembled monolayers or micelles. The interaction between azobenzene units leads to a formation of exciton states in these molecular assemblies. Apart from local excitations of monomers, the electronic transitions to the exciton states may involve charge transfer excitations. Here, we perform quantum chemical calculations and apply transition density matrix analysis to quantify local and charge transfer contributions to the lowest electronic transitions in azobenzene dimers of various arrangements. We find that the transitions to the lowest exciton states of the considered dimers are dominated by local excitations, but charge transfer contributions become sizable for some of the lowest ππ* electronic transitions in stacked and slip-stacked dimers at short intermolecular distances. In addition, we assess different ways to partition the transition density matrix between fragments. In particular, we find that the inclusion of the atomic orbital overlap has a pronounced effect on quantifying charge transfer contributions if a large basis set is used.

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.