Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correspondence”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Stability of expanding accretion shocks for an arbitrary equation of state

We present a theoretical stability analysis for an expanding accretion shock that does not involve a rarefaction wave behind it. The dispersion equation that determines the eigenvalues of the problem and the explicit formulae for the corresponding eigenfunction profiles are presented for an arbitrary equation of state and finite-strength shocks. For spherically and cylindrically expanding steady shock waves, we demonstrate the possibility of instability in a literal sense, a power-law growth of shock-front perturbations with time, in the range of $h_c< h<1+2 {\mathcal {M}}_2$ , where $h$ is the D'yakov-Kontorovich parameter, $h_c$ is its critical value corresponding to the onset of the instability and ${\mathcal {M}}_2$ is the downstream Mach number. Shock divergence is a stabilizing factor and, therefore, instability is found for high angular mode numbers. As the parameter $h$ increases from $h_c$ to $1+2 {\mathcal {M}}_2$ , the instability power index grows from zero to infinity. This result contrasts with the classic theory applicable to planar isolated shocks, which predicts spontaneous acoustic emission associated with constant-amplitude oscillations of the perturbed shock in the range $h_c< h<1+2 {\mathcal {M}}_2$ . Examples are given for three different equations of state: ideal gas, van der Waals gas and three-terms constitutive equation for simple metals.

Huete, César (ORCID:0000000232278520)↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in a magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry–Horton system has been proposed by St-Onge ( J. Plasma Phys. , vol. 83, 2017, 905830504) as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$ , where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U'' \doteq {\rm d}^2 U/{\rm d} x^2 $ . Then, depending on $U''$ , the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $| {\rm d} U / {\rm d} x |$ . Our analytic predictions are in agreement with direct numerical simulations of the modified Terry–Horton system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
DOE OSTI2024

Nonlinear solution of classical three-wave interaction via finite-dimensional quantum model

The quantum three-wave interaction, the lowest-order nonlinear interaction in plasma physics, describes energy–momentum transfer between three resonant waves in the quantum regime. We describe how it may also act as a finite-degree-of-freedom approximation to the classical three-wave interaction in certain circumstances. By promoting the field variables to operators, we quantize the classical system, show that the quantum system has more free parameters than the classical system and explain how these parameters may be selected to optimize either initial or long-term correspondence. We then numerically compare the long-time quantum–classical correspondence far from the fixed point dynamics. We discuss the Poincaré recurrence of the system and the mitigation of quantum scrambling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Geometrical optics without singularities: using the ray time as the coordinate space

Geometrical optics (GO) is widely used for reduced modelling of waves in plasmas, but it fails near reflection points, where it predicts a spurious singularity of the wave amplitude. We show how to avoid this singularity by adopting a different representation of the wave equation. Instead of the physical coordinate 𝑥 and the wavevector 𝑘, we use the ray time 𝜏 as the new canonical coordinate and the ray energy ℎ as the associated canonical momentum. To derive the envelope equation in the 𝜏-representation, we construct the Weyl symbol calculus on the (𝜏,ℎ) space and show that the corresponding Weyl symbols are related to their (𝑥,𝑘) counterparts by the Airy transform. This allows us to express the coefficients in the envelope equation through the known properties of the original dispersion operator. When necessary, solutions of this equation can be mapped to the 𝑥-space using a generalised metaplectic transform. However, the field per se might not even be needed in practice. Instead, knowing the corresponding Wigner function usually suffices for linear and quasilinear calculations. As a Weyl symbol itself, the Wigner function can be mapped analytically, using the aforementioned Airy transform. We show that the standard Airy patterns that form in regions where conventional GO fails are successfully reproduced within metaplectic GO (MGO) simply by remapping the field from the 𝜏-space to the 𝑥-space. An extension to mode-converting waves is also presented. This formulation, which we call generalised MGO, can be particularly useful, for example, for reduced modelling of the O–X conversion in inhomogeneous plasma near the critical density, an effect that is important for fusion applications and also occurs in the ionosphere. Overall, MGO can replace GO for any practical purposes, because it better handles cutoffs and is similar otherwise.

plasma waves↗
DOE OSTI2020

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

NASICON Na 3 V 2 (PO 4 ) 3 Enables Quasi-Two-Stage Na + and Zn 2+ Intercalation for Multivalent Zinc Batteries

Identifying positive electrode materials capable of reversible multivalent electrochemistry in electrolytes containing divalent ions such as Mg 2+ , Ca 2+ , and Zn 2+ at high operating potentials remains an ongoing challenge in “beyond lithium-ion” research. In this paper, we explore the Zn 2+ charge-storage mechanism of a vanadium-based Na + superionic conductor (NASICON), Na 3 V 2 (PO 4 ) 3 . By using X-ray synchrotron techniques to unravel potential-dependent structure–property relationships, we ascribe the reversible electrochemical behavior of Na 3 V 2 (PO 4 ) 3 to a quasi-two-stage intercalation process that involves both Na+ and Zn 2+ . Initial charging of Na 3 V 2 (PO 4 ) 3 leads to a Na + -extracted phase corresponding to Na 3 V 2 (PO 4 ) 3 , whereas subsequent discharge results predominantly in Na + intercalation followed by Zn 2+ intercalation. Operando X-ray diffraction of Na 3 V 2 (PO 4 ) 3 was used to study the phase changes associated with the first charge/discharge process, and ex situ measurements were used to precisely link the changes in the crystal structure to a quasi-two-stage intercalation of Na + and Zn 2+ . The corresponding changes in the V-oxidation state, V-O coordination, and the presence of Zn 2+ were confirmed by X-ray absorption spectroscopy. The results of this work present a comprehensive understanding of the charge-storage properties for a well-established NASICON structure that confers both the high capacity (~100 mA h g -1 ) and high potential (1.35 and 1.1 V vs Zn/Zn 2+ ).

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

Localized High Concentration Electrolytes for High Voltage Lithium-Metal Batteries: Correlation between the Electrolyte Composition and Its Reductive/Oxidative Stability

We demonstrate a first-principles screening methodology as an effective tool to explore electrolyte formulations for the new generation of high energy density rechargeable batteries. We study the liquid structure and electronic properties in dilute electrolytes, high concentration electrolytes (HCE), and localized high concentration electrolytes (LHCE), with focus on electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and bis(2,2,2-trifluoroethyl) ether (BTFE) as a diluent. We describe the solvation complexes in the dilute electrolyte and explore structural changes triggered by the increase in lithium salt concentration for HCEs and the diluent effects in LHCEs. In HCE formulations, there is a 4-fold coordination environment of lithium-ions as in the dilute electrolyte, but the number of lithium-ion interactions with O atoms from FSI- anions dominates. In these solutions, the ability of the FSI- anions to interact with multiple lithium-ions allows complex 3D network formation and influences the reductive/oxidative behavior of the electrolyte. Interestingly, in LHCEs, the BTFE diluent molecules do not change the 3D solution structure when diluting the HCE formulation from 5.49 to 3.83 M. However, there is a composition threshold where the structural and electronic behavior may change. We show that diluting the HCE electrolyte with BTFE down to 1.77 M breaks the three-dimensional solution structure into an island-like solvation complex. We relate these structural changes to the electronic properties of the electrolytes finding a causal relationship between the reductive/oxidative behavior and the lithium-oxygen interaction mechanisms in the solvated complexes. The coordination with lithium-ions lowers the electrolyte LUMO and HOMO levels: the higher is the number of interactions with lithium-ions, the more likely the solvent molecule, FSI- anion, or diluent molecule is to be reduced and the less likely it is to become oxidized. The evolution of the solvated ion structure in HCE and LHCE suggests a close connection to a corresponding change in the lithium-ion transport mechanisms for these electrolytes.

anions, molecules, oxygen, solvents, electrolytes↗
DOE OSTI2021

Anionic Oxygen Redox in the High-Lithium Material Li 8 SnO 6

Lithiated transition metal oxide cathode materials that combine oxygen-anion redox [O 2 – /(O 2 ) n– ] with cation redox offer substantial capacities and higher voltage than cathodes without oxygen redox, but oxygen-anion redox is usually accompanied by the formation of peroxo or superoxo species that may cause oxygen release. To clarify the relationship between the oxygen release and peroxide and superoxide dimers without the complication of transition-metal redox, we performed density functional theory calculations to study the layered Li 8 SnO 6 cathode because it has only oxygen-anion redox during delithiation. Features with O–O distances of ~1.5 Å (corresponding to the formation of a peroxide) were observed at 75% lithium concentration, and further delithiation to 62.5% lithium concentration shortened the O–O bond distance to ~1.3 Å (corresponding to the formation of superoxide). The solely anionic oxygen redox stabilizes the delithiated structure of Li 8–x SnO 6 , and it does not cause structural disorder or release oxygen gas. Ab initio molecular dynamics calculations in the supercell at 300 K indicate that oxygen dimerization may occur at 75% lithium concentration. On the (0001) surface with a high lithium concentration (98%), oxygen dimerization is unfavorable. The calculations by HSE06 show that oxygen redox in Li 8 SnO 6 offers an initial voltage (vs Li/Li + ) greater than 4.2 V. The GGA+U calculations indicate that the Li + -polaron in Li 95 Sn 12 O 72 migrates with a migration barrier of only 0.43 eV. Here, this work also quantifies a solely anionic oxygen redox mechanism in Li 8 SnO 6 , and it provides a deeper understanding of anionic oxygen redox in Li-excess oxide cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2021

Nanoscale Phase Separation and Large Refrigerant Capacity in Magnetocaloric Material LaFe 11.5 Si 1.5

Here, ferromagnetic transitions, the magnetocaloric effect, and the atomic-scale microstructure in an annealed bulk material of LaFe 11.5 Si 1.5 were investigated using magnetic measurements, Mössbauer spectroscopy, high-resolution X-ray diffraction analysis, and high-resolution transmission electron microscopy. The results provide evidence of the coexistence of two ferromagnetic phases with different Curie temperatures. The phase with a Curie temperature of 216 K corresponds to a Si-rich phase of a NaZn 13 -type structure with a small lattice volume, whereas the phase with a Curie temperature of 185 K corresponds to a Si-poor phase of a NaZn 13 -type structure with a large lattice volume. This phase coexistence is observed on a nanometer scale and it can account for inverse thermal hysteresis of the ferromagnetic transition in the Si-rich phase during warming. Furthermore, the phase coexistence helps the bulk material achieve a refrigeration capacity of 170 J kg -1 due to a combination of the magnetocaloric effect in each phase. It is determined that the magnetocaloric effect in the Si-rich and the Si-poor phase brings about a maximum of the total entropy change of -3.55 and -6.80 and J kg -1 K -1 at respective Curie temperatures, respectively. Fundamentally, the phase coexistence indicates coloring of Si atoms in the lattice of a NaZn 13 -type structure and can be attributed to nanoscale phase separation of a disordered precursor during annealing of the bulk material. The nanoscale phase separation may occur in other magnetocaloric materials, thus allowing for improvement of their refrigeration capacity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reducing Surface Recombination Velocity of Methylammonium-Free Mixed-Cation Mixed-Halide Perovskites via Surface Passivation

We control surface recombination in the mixed-cation, mixed-halide perovskite, FA0.83Cs0.17Pb(I0.85Br0.15)3, by passivating nonradiative defects with the polymerizable Lewis base (3-aminopropyl)trimethoxysilane (APTMS). We demonstrate average minority carrier lifetimes >4 µs, nearly single exponential monomolecular photoluminescence decays, and high external photoluminescence quantum efficiencies (>20%, corresponding to ~97% of the maximum theoretical quasi-Fermi-level splitting) at low excitation fluence. We confirm both the composition and valence band edge position of the FA0.83Cs0.17Pb(I0.85Br0.15)3 perovskite using multi-institutional, cross-validated, X-ray photoelectron spectroscopy and UV photoelectron spectroscopy measurements. We extend the APTMS surface passivation to higher bandgap double-cation (FA and Cs) compositions (1.7, 1.75, and 1.8 eV) as well as the widely used triple-cation (FA, MA, and Cs) composition. Finally, we demonstrate that the average surface recombination velocity decreases from ~1000 to ~10 cm/s post APTMS passivation for FA0.83Cs0.17Pb(I0.85Br0.15)3. Our results demonstrate that surface-mediated recombination is the primary nonradiative loss pathway in many methylammonium (MA)-free mixed-cation mixed-halide films with a range of different bandgaps, which is a problem observed for a wide range of perovskite active layers and reactive electrical contacts. Our study also provides insights to develop passivating molecules that help reduce surface recombination in MA-free mixed-cation mixed-halide films and indicates that surface passivation and contact engineering will enable near-theoretical device efficiencies with these materials.

14 SOLAR ENERGY↗
DOE OSTI2023

Revealing the Correlation between the Solvation Structures and the Transport Properties of Water-in-Salt Electrolytes

Water-in-salt (WIS) electrolytes containing 21 m lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) have been considered as a safe and environment-friendly alternative to common organic electrolytes used in lithium-ion batteries. However, the relation between the solvation structures and transport properties of these materials remains elusive. Here, for this paper, we performed small-angle X-ray scattering (SAXS), small-angle neutron scattering (SANS), and X-ray pair distribution function (PDF) measurements of LiTFSI aqueous solutions at a wide range of concentrations. Combined with molecular dynamics simulations, the detailed solvation structures from long to short length scale were resolved. We found that the TFSI– solvation structures consist of TFSI– solvated structures and TFSI– networks; the former corresponds to solvent separated ion pairs, while the latter corresponds to contact ion pairs and cation–anion aggregates. In addition, we found that the relaxation time in the q range associated with the anion network structure exhibits the same concentration dependence as the viscosity. By combining the results from the experiments and simulations, this study revealed a correlation between the solvation structures of LiTFSI and the transport properties of the solutions, which is critical to understand the relation between the transport properties and the dynamics of the ions for imide-based lithium-ion salt aqueous electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spectroscopic Investigation into P(NDI2OD-T2) Charge Localization

P(NDI2OD-T2), commonly referred to as N2200, stands out as a promising electron-transporting (n-type) polymer for low-cost, flexible (photo)electrochemical applications due to its reversible two-electron reduction and high electron mobility. UV–vis spectroelectrochemistry in the tetrabutyl ammonium hexafluorophosphate/acetonitrile electrolyte shows two sets of chemically reversible redox signals in the cyclic voltammetry corresponding to the reduction of the neutral polymer film to polaronic and bipolaronic species. These electrochemical signatures suggest a distinct electronic reorganization upon reduction (polaron/bipolaron formation), highlighting the need for molecular-level insights into how charges are accommodated within the polymer backbone. While it has been previously hypothesized that charge predominantly localizes on the naphthalene diimide (NDI) unit during reductive charging, specific changes in atomic environments that confirm this localization have not been characterized in n-type polymers. Herein, we use near-edge X-ray absorption fine structure (NEXAFS) spectroscopy to probe electronic transitions in an electrochemically charged polymer to deduce charge localization. The O K-edge (1s) spectra exhibit two distinct π* peaks; the intensity of the lower-energy π* a peak that corresponds to an excitation to a largely localized carbonyl state decreases with reductive potentials relative to the higher-energy π* b peak. We corroborate this with Raman spectroscopy at different potentials, which shows a decrease in intensity on the C–C/C═C and C═O stretching bands of NDI as well as a red shift of the carbonyl band due to the formation of a polaron on the NDI. Additionally, new Raman active NDI signals associated with elongated C═O and C═C bonds are observed at lower energy during the formation of charged states. Together with theoretical calculations, these findings show that the injected charge spatially localizes on the NDI units and is dominantly distributed on the carbonyl groups. In conclusion, the combination of NEXAFS, optical and vibrational spectroscopies, and theoretical calculations is generalizable to other pi-conjugated polymers and can identify charge localization for the further development of organic semiconductors.

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

Structural Phase Transitions in the van der Waals Ferromagnets Fe x Pd y Te 2

Two-dimensional van der Waals (vdW) magnets are attracting significant attention, both as platforms for studying fundamental magnetic interactions and for the exciting possibility of utilizing them as building blocks in devices and heterostructures, which may lead to new physical phenomena and functionalities. Here, we provide a detailed study of the crystal structure and physical properties of the recently discovered vdW ferromagnet FePd 2 Te 2 . We find this compound has a relatively wide width of formation, and grow single crystals with compositions Fe x Pd y Te 2 where x ranges from 0.9 to 1.1 and y from 1.8 to 2.5, respectively. Temperature-dependent X-ray diffraction and transport measurements reveal that a first-order structural transition occurs in the range of T = 360–420 K, where the critical temperature, modulation wave vector, and corresponding room-temperature crystal structures all depend on chemical composition. Above the transition, the compounds with Pd fraction y > 2 adopt a disordered derivative of the tetragonal FeTe structure, with the Fe layer showing mixed Fe/Pd occupancy and the extra Pd atoms partially occupying interstitial sites. Below 370 K, the structure is incommensurately modulated, likely associated with the complex ordering of Pd/Fe atoms in the metal layers or the interstitial Pd in the vdW gaps. For y < 2, the composition Fe 1.1 Pd 1.8 Te 2 has monoclinic symmetry at room temperature that is consistent with the reported structure of FePd 2 Te 2 . This phase undergoes a structural transition at 420 K for which the high temperature structure is yet to be determined; however, based on the similarities with the y > 2 compounds, we speculate that its T > 420 K structure is also tetragonal. Importantly, the high temperature, symmetry-breaking structural transition observed here provides a likely explanation for the origin of the structural domains previously observed in FePd 2 Te 2 . All compounds investigated in the Fe x Pd y Te 2 series show metallic behavior, with magnetic characterization indicating that they are easy-plane, hard, ferromagnets with T C spanning 98–180 K. Both the critical temperature for the structural transition and the Curie temperature are moderately suppressed with increasing Pd fraction y and corresponding decreasing Fe fraction x, indicating that synthetic control over x and y paves way for the further exploration of these compounds.

crystal structure↗
DOE OSTI2020

Determination of the N–H Bond Dissociation Free Energy in a Pyridine(diimine)molybdenum Complex Prepared by Proton-Coupled Electron Transfer

The pyridine(diimine) molybdenum bis(imido) complex, ( iPr PDI)Mo(=NTol) 2 (Tol = 4-methylphenyl) was synthesized by addition of two equivalents of 4-methylphenylazide to a corresponding molybdenum benzene derivative, ( iPr PDI)Mo(η 6 -C 6 H 6 ) ( iPr PDI = 2,6-(2,6-iPr 2 C 6 H 3 N=CMe) 2 C 5 H 3 N). Protonation of ( iPr PDI)Mo(=NTol) 2 with 2,6-lutinidum triflate yielded a cationic molybdenum amido complex, [( iPr PDI)Mo(NHTol)(=NTol)][OTf], which was further transformed into the neutral molybdenum amido complex, ( iPr PDI)Mo(NHTol)(=NTol) by reduction with zinc powder. A series of spectroscopic, synthetic and pK a determination studies along with electrochemical measurements by the protonation-reduction pathway were used to establish an N–H bond dissociation free energy (BDFE) between 65-69 kcal/mol for the molybdenum imido-amido compound, ( iPr PDI)Mo(NHTol)(=NTol). Full molecule DFT studies provided a computed value of 61 kcal/mol. By contrast, reduction of ( iPr PDI)Mo(=NTol) 2 with KC 8 afforded the corresponding anionic molybdenum complex, K[( iPr PDI)Mo(=NTol) 2 ] that has a potassium cation intercalated with the pyridine and the tolyl groups. Protonation of K[( iPr PDI)Mo(=NTol) 2 ] with the weak amidinium acids [TBD(H)][BArF 24 ] (TBD = triazabicyclodecene, BArF 24 = B[3,5-(CF 3 ) 2 C 6 H 3 ] 4 ) also produced the neutral molybdenum amido complex, ( iPr PDI)Mo(NHTol)(=NTol). Measurement of the pKa and oxidation potential of K[( iPr PDI)Mo(=NTol) 2 ] provided a range of 69-73 kcal/mol for the N–H BDFE of ( iPr PDI)Mo(NHTol)(=NTol), in good agreement with the protonation-reduction route and completing the square scheme. The similar pK a and redox potentials obtained from each pathway demonstrate that both sequences are energetically feasible for PCET events. This study on the determination of N–H BDFE of the molybdenum amido complex renders fundamental insight into the N 2 reduction cycle by proton-coupled electron transfer.

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

Aqueous TiO 2 Nanoparticles React by Proton-Coupled Electron Transfer

Redox reactions of aqueous colloidal TiO 2 4 nm nanoparticles (NPs) have been examined, including both citratecapped and uncapped NPs (c-TiO 2 and uc-TiO 2 ). Photoreduction gave stable blue colloidal c-TiO 2 R NPs with 10–60 electrons per particle. Equilibration of these reduced NPs with soluble redox reagents such as methylviologen (MV 2+ ) provided measurements of the colloid reduction potential as a function of pH. The potentials of c-TiO 2 from pH 2–9 varied linearly with pH, with a slope of –60 ± 5 mV/pH. Estimates of the potential at pH 12 were consistent with extrapolating that line to high pH. The reduction potentials did not correlate with the zeta potentials (ζ) or the surface charge of the NPs across this pH range. Similar reduction potentials were observed for c- and uc-TiO 2 at low pH even though they have quite different ζ potentials. These results show that the common surface-charging explanation of the pH dependence is not tenable in these systems. Oxidation of reduced c-TiO 2 R with the electron-transfer oxidant potassium triiodide (KI 3 ) occurred with a significant drop in pH, showing that protons were released when the electrons were removed from the NPs. Smaller pH drops were observed for the proton-coupled electron transfer (PCET) reagents O 2 (air) and 4-MeO-TEMPO (4-methoxy-2,2,6,6-tetramethylpiperine-1-oxy radical). The difference in the number of protons released with KI 3 vs O 2 and 4-MeO-TEMPO was roughly one proton per electron removed. Thus, the thermodynamically preferred reactivity of these colloidal TiO 2 NPs is PCET over the pH 2–13 range studied. The measured redox potentials refer to the chemical process TiO 2 + H + + e – → TiO 2 ·e – ,H + ; and therefore they do not correspond with an electronic energy such as a conduction band edge or flat band potential. The 1e – /1H + stoichiometry means that the TiO 2 reduction potentials correspond to a TiO 2 –H bond dissociation free energy (BDFE), determined to be 49 ± 2 kcal mol –1 . The PCET description is consistent with the pH dependence of E(TiO 2 /TiO 2 ·e – ,H + ), the release of protons upon oxidation, the lack of correlation with ζ potentials, the similarity of capped and uncapped NPs, and the small change in the potential and BDFE from the first to the last electron/proton pair (H atom) removed. Furthermore, this behavior is suggested to be the norm for redox-active oxide/water interfaces.

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

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

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

Localized Active Space Pair-Density Functional Theory

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗