Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular orbitals”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Using projection operators with maximum overlap methods to simplify challenging self-consistent field optimization

Maximum overlap methods are effective tools for optimizing challenging ground- and excited-state wave functions using self-consistent field models such as Hartree-Fock and Kohn-Sham density functional theory. Nevertheless, such models have shown significant sensitivity to the user-defined initial guess of the target wave function. Here in this work, a projection operator framework is defined and used to provide a metric for non-aufbau orbital selection in maximum-overlap-methods. The resulting algorithms, termed the Projection-based Maximum Overlap Method (PMOM) and Projection-based Initial Maximum Overlap Method (PIMOM), are shown to perform exceptionally well when using simple user-defined target solutions based on occupied/virtual molecular orbital permutations. This work also presents a new metric that provides a simple and conceptually convenient measure of agreement between the desired target and the current or final SCF results during a calculation employing a maximum-overlap method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Plutonium Hybrid Materials: A Platform to Explore Assembly and Metal–Ligand Bonding

In this work, we report the synthesis of five new hybrid materials containing the [PuCl 6 ] 2- anion and charge balancing, non-covalent interaction donating 4-X-pyridinium (X = H, Cl, Br, I) cations. Single crystals of the title compounds were grown and harvested from acidic, chloride-rich, aqueous media and their structures were determined via X-ray diffraction. Compounds 1-4, (4XPyH) 2 [PuCl 6 ] and 5, (4IPyH) 4 [PuCl 6 ] · 2Cl, exhibit two distinct sheet-like structure types. Structurally relevant non-covalent interactions were tabulated from crystallographic data and verified computationally using electrostatic surface potential maps and the quantum theory of atoms in molecules (QTAIM) approach. The strength of the hydrogen and halogen bonds was quantified using Kohn-Sham density functional theory and a hierarchy of acceptor-donor pairings established. In turn, the PuIV-Cl bonds were studied using the QTAIM and natural localized molecular orbital (NLMO) approaches to delineate the underlying bond mechanism and hybrid atomic orbital contributions therein. Energy decomposition (ED) and natural ED analyses were also explored to probe the bond mechanism and, more broadly, explore the efficacy of these techniques in studying these anionic systems. The results of the PuIV-Cl bond analyses were compared across composition via analogous treatments of previously reported [PuO 2 Cl 4 ] 2- and [PuCl 3 (H 2 O) 5 ] molecular units. In summary, our study indicates that the Pu-Cl bonds are predominately ionic, yet exhibit small varying degrees of covalent character that increase from [PuCl 3 (H 2 O) 5 ], [PuO 2 Cl 4 ] 2- , to [PuCl 6 ] 2- , while the participation of the Pu based s/d and f orbitals concurrently decrease and increase, respectively.

transuranic↗
DOE OSTI2020

Resonant x-ray absorption of strong-field-ionized CF 3 Br

We report on an experimental and theoretical study of strong-field laser ionization of CF 3 Br followed by resonant x-ray absorption at the Br K-edge. Distinct 1s -> 4p, 5p Rydberg transitions of Br q+ (q = 1-4) atomic ions are observed and identified with Hartree-Fock-Slater and relativistic configuration interaction calculations. Time-dependent density functional theory and ab initio molecular dynamics calculations were performed to simulate the dissociative ionization process and the molecular orbitals for the q = 1-4 charge states. Measurements were made with both parallel and perpendicular linear polarizations of the laser and x-rays, but dichroism was not observed, indicating negligible alignment by the laser ionization process. This result is explained by calculations on atomic Br and the molecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

High sensitivity pressure and temperature quantum sensing in pentacene-doped p-terphenyl single crystals

Quantum sensors’ responsiveness to their physical environment enables detection of variables such as temperature (T), pressure (P), and strain. We present a molecular platform for PT sensing using para-terphenyl crystals doped with pentacene (PDP), leveraging optically detected magnetic resonance (ODMR) of photoexcited triplet electron spins. We observe maximal frequency variations of df/dP=1.8 MHz/bar from 0-8 bar and df/dT=247 kHz/K from 79–330 K, over 1200 times and threefold greater, respectively, than those seen with nitrogen-vacancy centers in diamond and > 85-fold greater pressure sensitivity over the previous record. Density functional theory calculations indicate picometer-level PT-induced molecular orbital shifts are measurable via ODMR. PDP offers additional advantages including high sensor doping levels, narrow ODMR linewidths, high contrast, and low-cost single crystal growth. Overall, this work reports low-cost, optically-interrogated PT sensors and lays the foundation for increased versatility of quantum sensors through synthetic molecular design.

Magnetic properties and materials↗
DOE OSTI2021

Modeling dark- and light-induced crystal structures and single-crystal optical absorption spectra of ruthenium-based complexes that undergo SO 2 -linkage photoisomerization

A family of coordination complexes of the type [Ru(SO 2 )(NH 3 ) 4 X] m+ Y n – (m, n = 1 or 2) exhibit optical switching capabilities in their single-crystal states. This striking effect is caused by the light-induced formation of SO 2 -linkage photoisomers, which are metastable if kept at suitably cool temperatures. We modeled the dark- and light-induced states of these large crystalline complexes via plane-wave (PW)- and molecular-orbital (MO)-based density functional theory (DFT) and time-dependent DFT in order to calculate their structural and optical properties; the calculated results are compared with experimental data. We show that the PW-DFT-based periodic models replicate the structural properties of these complexes more effectively than the MO-DFT-based molecular-fragment models, observing only small deviations in key bond lengths relative to the experimentally derived crystal structures. The periodic models were also found to more effectively simulate trends seen in experimental optical absorption spectra, with optical absorbance and coverage of the visible region increasing with the formation of the photoinduced geometries. The contribution of the metastable photoisomeric species more heavily focuses on the lower-energy end of the spectra. Spectra generated from the molecular-fragment models are limited by the geometry of the fragment used and the number of excited-state roots considered in those calculations. In general, periodic models outperform the molecular-fragment models owing to their ability to better appreciate the periodic phenomena that are present in these crystalline materials as opposed to MO approaches, which are finite methods. We thus demonstrate that PW-DFT-based periodic models should be considered as a more than viable method for simulating the optical and electronic properties of these single-crystal optical switches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Trimerized small-molecule acceptors enable high-performance organic solar cells with high open-circuit voltage and prolonged life-time

Although the recent development of Y-series small-molecule acceptors (SMAs) has led to a dramatic increase in the power conversion efficiency (PCE) of organic solar cells (OSCs), the operational stability (device lifetime) of OSCs is inadequate for commercialization. In this study, we develop a new trimer acceptor (TYT), consisting of three Y-based molecules linked by electron-donating spacers, to realize an OSC with high-performance (PCE > 18%) and -stability (t 80% lifetime > 8000 h under 1 sun illumination, in which t 80% lifetime is the time required for the PCE of the OSC to reach 80% of its initial value). We demonstrate that the trimerization approach affords an acceptor, TYT, with an upshifted lowest unoccupied molecular orbital energy level, which, in turn, affords an efficient OSC with a high open-circuit voltage (0.964 V). In addition, the glass-transition temperature (T g ) of TYT (217 °C) is significantly higher than those of monomer (MYT, T g = 80 °C) and dimer (DYT, T g = 127 °C) acceptors, which effectively suppresses the molecular diffusion of TYT in a blend film with a polymer donor. Accordingly, a TYT-based OSC demonstrates a high PCE (18.2%) and long t 80% lifetime under 1 sun illumination (8454 h), outperforming MYT- and DYT-based OSCs that exhibit PCEs and t 80% lifetimes of 16.4% and 35 h, and 17.3% and 2551 h, respectively. Therefore, this study provides important guidelines for the design of electron acceptors in achieving OSCs with high performance and stability close to a commercial level.

14 SOLAR ENERGY↗
DOE OSTI2021

Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction

While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory (κ-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of κ-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The aug-cc-pwCVTZ basis set is used for all calculations. The use of κ-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing two multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a κ-OOMP2 orbital reference (CCSD(T):κ-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm -1 , and 8.50 cm -1 , and 8.75 cm -1 respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets shows these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. CCSD(T) with RHF orbitals is seen to improve upon UHF for the closed-shell species, while spatial symmetry breaking in a number of restricted open-shell HF (ROHF) references leads CCSD(T) with ROHF reference orbitals to exhibit the poorest statistical performance of all methods surveyed for open-shell species. The use of κ-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Orbital Engineering Band Degeneracy in a Dual-Square Carbon-Oxide Framework

Electron band degeneracies in momentum space give rise to exotic quantum phenomena that have sparked intense interest in condensed matter physics and materials science. Nodal-lines─isolines in k-space formed by the incidental touching of two bands that share the same energy but belong to discrete eigenstates─arise in the presence of symmetries that preclude effective hybridization. Despite recent advances in the design, bottom-up assembly, and engineering of exotic electronic states in graphene nanomaterials, the extension of this approach to access synthetic two-dimensional (2D) quantum materials derived from metal- or covalent-organic frameworks (COFs) has lagged behind. Here we present a molecular orbital engineering approach for designing and fabricating an edge-centered dual square lattice within a π-conjugated 2D-tetraoxa[8]circulene (2D-TOC) COF. First-principles calculations and scanning tunnelling spectroscopy reveal the emergence of Frontier states at the center of a 3 × 3 lattice that give rise to Dirac nodal-lines in 2D-TOC. Our findings not only provide a general guide for the design of conjugated COFs with custom tailored electronic properties from molecular fragments but enable the exploration of emergent topological phenomena in synthetic 2D materials with potential application for high-speed, low-power data processing, transmission, and storage.

Dirac nodal line semimetal↗
DOE OSTI2024

Manipulating symmetry-breaking charge separation employing molecular recognition

The exploration of symmetry-breaking charge separation (SB-CS) is imperative when designing functional light-harvesting materials. Past explorations, however, have been confined to covalent systems, more often than not requiring complicated/demanding syntheses and facing inconvenient regulation of charge transfer processes. Here, in this work, we present a concept that regulates the efficiency of SB-CS through molecular recognition utilizing a pyridinium-based cyclophane as a host. This host undergoes photo-driven excited-state SB-CS. By employing different guests with distinct frontier molecular orbital energy levels, we have achieved comprehensive control of electron transfer pathways in the cyclophane, modulating between accelerated (>10-fold) intramolecular SB-CS involving superexchange and direct intermolecular electron transfer between the host and guest. The improvement in SB-CS efficiency results in catalytic activity for the photo-oxidation of a sulfur-mustard simulant. This research offers an opportunity for tuning SB-CS by utilizing molecular recognition, which holds the potential for achieving precise regulation without complicated organic syntheses.

charge transfer↗
DOE OSTI2021

InteraChem: Virtual Reality Visualizer for Reactive Interactive Molecular Dynamics

Interactive molecular dynamics in virtual reality (IMD-VR) simulations provide a digital molecular playground for students as an alternative or complement to traditional molecular modeling kits or 2D illustrations. Previous IMD-VR studies have used molecular mechanics to enable simulations of macromolecules such as proteins and nanostructures for the classroom setting with considerable success. Here, we present the InteraChem molecular visualizer, intended for reactive IMD-VR simulation using semiempirical and ab initio methods. InteraChem visualizes not only the molecular geometry but also (1) isosurfaces such as molecular orbitals and electrostatic potentials and (2) two-dimensional graphs of time-varying simulation quantities such as kinetic/potential energy, internal coordinates, and user-applied force. Additionally, InteraChem employs speech recognition to facilitate user interaction and introduces a novel “atom happiness” visualization using emojis to indicate the energetic feasibility of a particular bonding arrangement. Furthermore, we include a set of accompanying exercises that we have used to teach chemical reactivity in small molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Excited-State Dynamics during Primary C–I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Conjugated-linker dependence of the photophysical properties and electronic structure of chlorin dyads

The synthesis, photophysical properties and electronic structure of seven new chlorin dyads and associated benchmark monomers are described. Each dyad contains two identical chlorins linked at the macrocycle β-pyrrole 13-position. The extent of electronic communication between chlorin constituents depends on the nature of the conjugated linker. The communication is assessed by modification of prominent ground-state absorption and redox properties, rate constants and yields of excited-state decay processes, and molecular-orbital characteristics. Relative to the benchmark monomers, the chlorin dyads in toluene exhibit a substantial bathochromic shift of the long-wavelength absorption band (30 nm average), two-fold increased radiative rate constant [average (10 ns) -1 vs. (22 ns) -1 ], reduced singlet excited-state lifetimes (average 5.0 ns vs. 8.2 ns), and increased fluorescence quantum yields (average 0.56 vs. 0.42). The excited-state lifetime and fluorescence yield for the chlorin dyad with a benzothiadiazole linker are reduced substantially in benzonitrile vs. toluene due largely to ~25-fold accelerated internal conversion. The results aid design strategies for molecular architectures that may find utility in solar-energy conversion and photomedicine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Molecular-Scale Characterization of Photoinduced Charge Separation in Mixed-Dimensional InSe–Organic van der Waals Heterostructures

Layered indium selenide (InSe) is an emerging two-dimensional semiconductor that has shown significant promise for high-performance transistors and photodetectors. The range of optoelectronic applications for InSe can potentially be broadened by forming mixed-dimensional van der Waals heterostructures with zero-dimensional molecular systems that are widely employed in organic electronics and photovoltaics. Here, we report the spatially resolved investigation of photoinduced charge separation between InSe and two molecules (C 70 and C 8 -BTBT) using scanning tunneling microscopy combined with laser illumination. We experimentally and computationally show that InSe forms type-II and type-I heterojunctions with C 70 and C 8 -BTBT, respectively, due to an interplay of charge transfer and dielectric screening at the interface. Laser-excited scanning tunneling spectroscopy reveals a ~0.25 eV decrease in the energy of the lowest unoccupied molecular orbital of C 70 with optical illumination. Furthermore, photoluminescence spectroscopy and Kelvin probe force microscopy indicate that electron transfer from InSe to C 70 in the type-II heterojunction induces a photovoltage that quantitatively matches the observed downshift in the tunneling spectra. In contrast, no significant changes are observed upon optical illumination in the type-I heterojunction formed between InSe and C 8 -BTBT. Density functional theory calculations further show that, despite the weak coupling between the molecular species and InSe, the band alignment of these mixed-dimensional heterostructures strongly differs from the one suggested by the ionization potential and electronic affinities of the isolated components. Self-energy-corrected density functional theory indicates that these effects are the result of the combination of charge redistribution at the interface and heterogeneous dielectric screening of the electron–electron interactions in the heterostructure. In addition to providing specific insight for mixed-dimensional InSe–organic van der Waals heterostructures, this work establishes a general experimental methodology for studying localized charge transfer at the molecular scale that is applicable to other photoactive nanoscale systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

High-performance n-Type Polymers Based on Multiple Electron-withdrawing Groups Decorated ( E )-1,2-Di(thiophen-2-yl)ethene Building Blocks

The performance of n-type conjugated polymers lags far behind that of p-type polymers, which significantly restricts the development of organic electronics. The (E)-1,2-di(thiophen-2-yl)ethene (TVT) unit, owing to its unique advantages, has been widely applied in the design of p-type polymer semiconductors. Previous studies have demonstrated that introducing electron-withdrawing groups can lower the frontier orbital energy levels of polymers and enhance electron injection/transporting capabilities. Based on this, we proposed incorporating multiple electronwithdrawing groups, such as amide groups, fluorine atoms, and cyano groups, into the polymer backbones of TVT-based polymer to facilitate the electron transport. Here, we successfully designed and synthesized the polymers TVTDA-4FTVT and TVTDA-2F2CNTVT. Both polymers exhibited low frontier orbital energy levels. Due to its significantly higher crystallization tendency and favorable intermolecular packing structure, the organic field-effect transistor (OFET) device based on TVTDA-4FTVT demonstrated an electron mobility one order of magnitude higher than that of TVTDA-2F2CNTVT. TVTDA-4FTVT showed the highest electron mobility of 0.87 cm 2 ·V −1 ·s −1 , while TVTDA-2F2CNTVT exhibited the highest electron mobility of 0.049 cm 2 ·V −1 ·s −1 . Owing to its deeper lowest unoccupied molecular orbital (LUMO) level, the OFET devices based on TVTDA-2F2CNTVT showed good air stability after being placed in a natural environment for 15 d.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular engineering of ethereal electrolyte for ultrastable Si-based high voltage full cells

The successful application of Si-based high-energy Li-ion batteries (LIBs) depends on our ability to tailor electrolyte properties to achieve long-term stability and reliable performance. In this work, we demonstrate our rationale for the molecular design of ethereal solvents to address low anodic stability issues and produce a highly electrochemically stable electrolyte for Si‖LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) high-energy full cells. Unlike the trimethylsilyl group, the trifluoromethyl (–CF 3 ) group exerts a very strong electron-withdrawing effect on the glycol ether backbone, reducing the highest occupied molecular orbital (HOMO) energy level of the fluorinated glycol ether (FGE) and significantly enhancing its oxidation potential. The FGE-based electrolyte enables stable cycling of Si‖NMC811 full cells, delivering high specific capacity (900 mA h g −1 ) and coulombic efficiency (>99.78%) over extended (500) cycles. The improved electrochemical performance originates from the terminal fluorination of the diglyme backbone, which strengthens anion coordination in the solvation structure, leading to the preferential reduction of the FSI anion and the formation of robust solid electrolyte interphases (SEIs) on the Si surface. Through molecular engineering of ethereal solvents, we have discovered a promising candidate for a next-generation stable electrolyte, paving the way for the design of practical and commercially viable Si batteries.

Silicon anode↗
DOE OSTI2024

Synthesis and characterization of low-dimensional N-heterocyclic carbene lattices

The covalent interaction of N-heterocyclic carbenes (NHCs) with transition metal atoms gives rise to distinctive frontier molecular orbitals (FMOs). These emergent electronic states have spurred the widespread adoption of NHC ligands in chemical catalysis and functional materials. Although formation of carbene-metal complexes in self-assembled monolayers on surfaces has been explored, design and electronic structure characterization of extended low-dimensional NHC-metal lattices remains elusive. Here, in this work, we demonstrate a modular approach to engineering one-dimensional (1D) metal-organic chains and two-dimensional (2D) Kagome lattices using the FMOs of NHC–Au–NHC junctions to create low-dimensional molecular networks exhibiting intrinsic metallicity. Scanning tunneling spectroscopy and first-principles density functional theory reveal the contribution of C–Au–C π-bonding states to dispersive bands that imbue 1D- and 2D-NHC lattices with exceptionally small work functions.

Qie, Boyu↗
DOE OSTI2021

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗