Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “General Physics”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Rate Theory of Ion Pairing at the Water Liquid–Vapor Interface

There is overwhelming evidence that certain ions are present near the vapor–liquid interface of aqueous salt solutions. Despite their importance in many chemical reactive phenomena, how ion–ion interactions are affected by interfaces and their influence on kinetic processes is not well understood. Molecular simulations were carried out to exam the thermodynamics and kinetics of small alkali halide ions in the bulk and near the water vapor–liquid interface. We calculated dissociation rates using classical transition state theory, and corrected them with transmission coefficients determined by the reactive flux method and Grote-Hynes theory. Our results show that, in addition to affecting the free energy of ions in solution, the interfacial environments significantly influence the kinetics of ion pairing. The results obtained from the reactive flux method and Grote-Hynes theory on the relaxation time present an unequivocal picture of the interface suppressing ion dissociation. This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. The calculations were carried out using computer resources provided by the Office of Basic Energy Sciences.

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

A First-Principles Study of Hydrogen Diffusivity and Dissociation on δ-Pu (100) and (111) Surfaces

The diffusivity and chemical reactivity of hydrogen on the δ-Pu (100) and (111) surfaces has been studied with density functional theory using both spin-polarization and spin–orbit coupling calculations. Comparison of the total electronic density of states and atomic hydrogen diffusion energies indicates that spin-polarization yields accurate results, with spin–orbit coupling yielding slightly smaller barriers to hopping between adsorption sites. On the (100) surface, both sets of calculations indicate that the dissociation reaction for molecular hydrogen is highly active at ambient conditions and results in the hydrogen ions bonded within the Pu surface, similar to hydride formation. In contrast, calculations on the (111) surface indicate a lower barrier for dissociation without the formation of the hydride-like end product. Our results help quantify the effect of spin–orbit coupling on hydrogen surface chemistry on δ-Pu. Additionally, we observe that the degree of resulting corrosion could depend on the specific geometry of the plutonium surface, which could have ramifications for future engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Automatic Prediction of Surface Phase Diagrams Using Ab Initio Grand Canonical Monte Carlo

The properties of a material are often strongly influenced by its surfaces. Depending on the nature of the chemical bonding in a material, its surface can undergo a variety of stabilizing reconstructions that dramatically alter the chemical reactivity, light absorption, and electronic band offsets. For decades, ab initio thermodynamics has been the method of choice for computationally determining the surface phase diagram of a material under different conditions. The surfaces considered for these studies, however, are often human-selected and too few in number, leading both to insufficient exploration of all possible surfaces and to biases toward portions of the composition–structure phase space that often do not encompass the most stable surfaces. To overcome these limitations and automate the discovery of realistic surfaces, we combine density functional theory and grand canonical Monte Carlo (GCMC) into “ab initio GCMC.” This paper presents the successful application of ab initio GCMC to the study of oxide overlayers on Ag(111), which, for many years, mystified experts in surface science and catalysis. Specifically, we report that ab initio GCMC is able to reproduce the surface phase diagram of Ag(111) with no preconceived notions about the system. Using nonlinear, random forest regression, we discover that the Ag coordination number with O and the surface O–Ag–O bond angles are good descriptors of the surface energy. Additionally, using the composition–structure evolution histories produced by ab initio GCMC, we deduce a mechanism for the formation of oxide overlayers based on the Ag 3 O 4 pyramid motif that is common to many reconstructions of Ag(111). Lastly, ab initio GCMC is a promising tool for the discovery of realistic surfaces that can then be used to study phenomena on complex surfaces such as heterogeneous catalysis and material growth, enabling reliable and insightful interpretations of experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Towards Prediction of Non-Radiative Decay Pathways in Organic Compounds I: The Case of Naphthalene Quantum Yields

Many emerging technologies depend on human’s ability to control and manipulate the excited-state properties of molecular systems. These technologies include fluorescent labeling in biomedical imaging, light harvesting in photovoltaics, and electroluminescence in light-emitting devices. All of these systems suffer from non-radiative loss pathways that dissipate electronic energy as heat, which causes the overall system efficiency to be directly linked to quantum yield (Φ) of the molecular excited state. Unfortunately, Φ is very difficult to predict from first principles because the description of a slow non-radiative decay mechanism requires an accurate description of long-timescale excited-state quantum dynamics. In the present study, we introduce an efficient semi-empirical method of calculating the fluorescence quantum yield (Φ fl ) for molecular chromophores, which converts simple electronic energies computed using time-dependent density functional theory (TDDFT) into an estimate of Φ fl . As with all machine learning strategies, the algorithm needs to be trained on fluorescent dyes for which Φ fl ’s are known, so as to provide a black-box method which can later predict Φfl’s for chemically similar chromophores that have not been studied experimentally. Here, as a first illustration of how our proposed algorithm can be trained, we examine a family of 25 naphthalene derivatives. The simplest application of the energy gap law is found to be inadequate to explain the rates of internal conversion (IC) or intersystem crossing (ISC) – the electronic properties of at least one higher-lying electronic state (S n or T n ) or one far-from-equilibrium geometry are typically needed to obtain accurate results. Indeed, the key descriptors turn out to be the transition state between the Franck–Condon minimum and a distorted local minimum near an S 0 /S 1 conical intersection (which governs IC) and the magnitude of the spin-orbit coupling (which governs ISC). The resulting Φ fl ’s are predicted with reasonable accuracy (~ 22%), making our approach a promising ingredient for high-throughput screening and rational design of the molecular excited states with desired Φ fl ’s. We thus conclude that our model, while semi-empirical in nature, does in fact extract sound physical insight into the challenge of describing non-radiative relaxations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Entangled Electrons Drive a Non-superexchange Mechanism in a Cobalt Quinoid Dimer Complex

A central theme in chemistry is the understanding of the mechanisms that drive chemical transformations. A well-known, highly cited mechanism in organometallic chemistry is the superexchange mechanism in which unpaired electrons on two or more metal centers interact through an electron pair of the bridging ligand. We use a combination of novel synthesis and computation to show that such interactions may in fact occur by a more direct mechanism than superexchange that is based on direct quantum entanglement of the two metal centers. Specifically, we synthesize and experimentally characterize a novel cobalt dimer complex with benzoquinoid bridging ligands and investigate its electronic structure with the variational two-electron reduced density matrix method using large active spaces. The result draws novel connections between inorganic mechanisms and quantum entanglement, thereby opening new possibilities for the design of strongly correlated organometallic compounds whose magnetic and spin properties have applications in superconductors, energy storage, thermoelectrics, and spintronics.

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

Rotational Modulation of à 2 A State Photodissociation of HCO via Renner-Teller Nonadiabatic Transitions

By examining the product-state distribution of a prototypical non- adiabatic predissociation system, HCO($Ã^2A''–\tilde X^2A'$), we demonstrate here that the dissociation dynamics is strongly modulated by parent rotational quantum numbers. The predissociation of the nominal ($ν_{\text{C–H}} = 0, ν_{\text{bend}}, ν_{\text{C–O}} = 1$) vibronic levels of the $AÃ^2A''$ state surprisingly gives rise to both vibrational ground and excited states of the CO product, despite the assumed spectator nature of the CO moiety. This anomaly is attributed to the dependence of the lifetime of the vibronic resonance facilitated by the Renner–Teller interaction on the parent rotational angular momentum quantum numbers coupled with transient intensity borrowing from nearby vibronic resonances with $ν_{\text{C–O}} = 0$. This unique phenomenon is a purely quantum mechanical behavior that has no classical analogue.

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

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Beyond Thermodynamics: Assessing the Dynamical Softness of Hydrated Ions from First Principles

Ion (de)hydration is a key rate-determining step in interfacial processes from corrosion to electrochemical energy storage. However, predicting the kinetics of ion (de)hydration remains challenging, prompting the use of static proxies such as hydration energy and valence. While useful for assessing thermodynamic preferences, such descriptors cannot fully capture the dynamical softness of the hydration shell that dictates kinetics. Accordingly, we use first-principles molecular dynamics to analyze hydration shell softness for a diverse set of metal cations. Additionally, three dynamic metrics are introduced to intuitively describe the bond rigidity, shape deformability, and exchange fluidity of the solvation shell. Together, these metrics capture the relevant physics in the static descriptors, while offering a far more complete and efficient representation for the overall propensity for (de)hydration. Application to the hydrated ion set demonstrates a weak connection between dynamical softness and hydration energy, confirming that dynamical descriptors of hydration are key for correctly describing ion transfer processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantification of Exciton Fine Structure Splitting in a Two-Dimensional Perovskite Compound

Applications of two-dimensional (2D) perovskites have significantly outpaced the understanding of many fundamental aspects of their photophysics. The optical response of 2D lead halide perovskites is dominated by strongly bound excitonic states. However, a comprehensive experimental verification of the exciton fine structure splitting and associated transition symmetries remains elusive. Here we employ low temperature magneto-optical spectroscopy to reveal the exciton fine structure of (PEA) 2 PbI 4 (here PEA is phenylethylammonium) single crystals. We observe two orthogonally polarized bright in-plane free exciton (FX) states, both accompanied by a manifold of phonon-dressed states that preserve the polarization of the corresponding FX state. Introducing a magnetic field perpendicular to the 2D plane, we resolve the lowest energy dark exciton state, which although theoretically predicted, has systematically escaped experimental observation (in Faraday configuration) until now. These results corroborate standard multiband, effective-mass theories for the exciton fine structure in 2D perovskites and provide valuable quantification of the fine structure splitting in (PEA) 2 PbI 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solvent Induced Proton Polarization within the Nuclear–Electronic Orbital Framework

To explicitly account for nuclear quantum effects and solvent environments in simulations of chemical processes, the nuclear-electronic orbital approach is coupled with a polarizable continuum model (PCM). This NEO–PCM approach is used to explore the influence of solvation on nuclear polarization through applications to a water dimer and a set of protonated water tetramers. Nuclear polarization in these species is analyzed in terms of changes in proton density and oxygen–hydrogen bond length. Solvation is shown to enhance nuclear polarization with increasing dielectric constant. For the water dimer, the internal, hydrogen-bonded proton is shown to polarize more than the external, free proton. Moreover, proton quantization leads to greater solvent polarization through their mutual polarization. Furthermore, these calculations highlight the complex interplay among electronic, nuclear, and solvent polarization in chemical systems.

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

Spectral Signatures of a Negative Polaron in a Doped Polymer Semiconductor: Energy Levels and Hubbard U Interactions

The modern picture of negative charge carriers on conjugated polymers invokes the formation of a singly occupied (spin-up/spin-down) level within the polymer gap and a corresponding unoccupied level above the polymer conduction band edge. The energy splitting between these sublevels is related to on-site Coulomb interactions between electrons, commonly termed Hubbard U. However, spectral evidence for both sublevels and experimental access to the U value is still missing. Here, in this work, we provide evidence by n-doping the polymer P(NDI 2 OD-T 2 ) with [RhCp*Cp] 2 , [N-DMBI] 2 , and cesium. Changes in the electronic structure after doping are studied with ultraviolet photoelectron and low-energy inverse photoemission spectroscopies (UPS, LEIPES). UPS data show an additional density of states (DOS) in the former empty polymer gap while LEIPES data show an additional DOS above the conduction band edge. These DOS are assigned to the singly occupied and unoccupied sublevels, allowing determination of a U value of ~1 eV.

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

Nonadiabatic Derivative Couplings through Multiple Franck–Condon Modes Dictate the Energy Gap Law for Near and Short-Wave Infrared Dye Molecules

Near infrared (NIR, 700–1000 nm) and short-wave infrared (SWIR, 1000–2000 nm) dye molecules exhibit significant nonradiative decay rates from the first singlet excited state to the ground state. While these trends can be empirically explained by a simple energy gap law, detailed mechanisms of nearly universal behavior have remained unsettled for many cases. Theoretical and experimental results for two representative NIR/SWIR dye molecules reported here clarify the key mechanism for the observed energy gap law behavior. It is shown that the first derivative nonadiabatic coupling terms serve as major coupling pathways for nonadiabatic decay processes from the first excited singlet state to the ground state for these NIR and SWIR dye molecules and that vibrational modes other than the highest frequency modes also make significant contributions to the rate. Finally, this assessment is corroborated by further theoretical comparison with possible alternative mechanisms of intersystem crossing to triplet states and also by comparison with experimental data for deuterated molecules.

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

Exact-Factorization-Based Surface Hopping without Velocity Adjustment

While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. Furthermore, the resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.

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

A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters

Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. Here, by performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.

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

Understanding Strain and Failure of a Knot in Polyethylene Using Molecular Dynamics with Machine-Learned Potentials

A neural network potential (NNP) has been developed by fitting to ab initio electronic structure data on hydrocarbons and is used to study failure of linear and knotted polyethylene (PE) chains. A linear PE chain must be highly strained before breaking as the stress is equally distributed across the chain. In contrast, the stress in a PE chain with a 31 or overhand knot, accumulates at the knot’s entrance/exit. We find the strain energy is greatest when the bond length and angle are strained simultaneously, and that the knot weakens the chain by increasing the variance of the C–C–C angle, thereby allowing rupture at lower bond strains. Here, we extend our analysis to both 51 and 52 knots and find that both break at the entrance/exit of a loop. Notably, molecular scale PE knots exhibit many of the same characteristics as knots in a macroscopic rope, with stick–slip phenomena upon tightening and similar points of failure.

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

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Mechanistic Regimes of Vibronic Transport in a Heterodimer and the Design Principle of Incoherent Vibronic Transport in Phycobiliproteins

Following the observation of coherent oscillations in nonlinear spectra of photosynthetic pigment protein complexes, in particular, phycobilliproteins such as PC645, coherent vibronic transport has been suggested as a design principle for novel light-harvesting materials. Vibronic transport between energetically remote pigments is coherent when the presence of a vibration resonant with the electronic energy gap supports transient delocalization between the electronic excited states. We establish the mechanism of vibronic transport for a model heterodimer across a wide range of molecular parameter values. The resulting mechanistic map demonstrates that the molecular parameters of phycobiliproteins in fact support incoherent vibronic transport. This result points to an important design principle: Incoherent vibronic transport is more efficient than a coherent mechanism when energetic disorder exceeds the coupling between the donor and vibrationally excited acceptor states. Finally, our results suggest that the role of coherent vibronic transport in pigment protein complexes should be reevaluated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗