Search NASA⌕ Search

SEARCH · Search NASA

Results for “sharp features”

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

Understanding the Oxygen Evolution Reaction Mechanism on CoO x using Operando Ambient-Pressure X-ray Photoelectron Spectroscopy

Photoelectrochemical water splitting is a promising approach for renewable production of hydrogen from solar energy and requires interfacing advanced water-splitting catalysts with semiconductors. Understanding the mechanism of function of such electrocatalysts at the atomic scale and under realistic working conditions is a challenging, yet important, task for advancing efficient and stable function. This is particularly true for the case of oxygen evolution catalysts and, here, we study a highly active Co 3 O 4 /Co(OH) 2 biphasic electrocatalyst on Si by means of operando ambient-pressure X-ray photoelectron spectroscopy performed at the solid/liquid electrified interface. Spectral simulation and multiplet fitting reveal that the catalyst undergoes chemical-structural transformations as a function of the applied anodic potential, with complete conversion of the Co(OH) 2 and partial conversion of the spinel Co 3 O 4 phases to CoO(OH) under precatalytic electrochemical conditions. Furthermore, we observe new spectral features in both Co 2p and O 1s core-level regions to emerge under oxygen evolution reaction conditions on CoO(OH). The operando photoelectron spectra support assignment of these newly observed features to highly active Co 4+ centers under catalytic conditions. Comparison of these results to those from a pure phase spinel Co 3 O 4 catalyst supports this interpretation and reveals that the presence of Co(OH) 2 enhances catalytic activity by promoting transformations to CoO(OH). The direct investigation of electrified interfaces presented in this work can be extended to different materials under realistic catalytic conditions, thereby providing a powerful tool for mechanism discovery and an enabling capability for catalyst design.

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

Investigation of hardening mechanisms and size effects in proton-irradiated HT-9 steels

Ferritic/martensitic steels, such as HT-9, are known for their complex microstructural features and mechanical properties. Here in this paper, in-situ micro-tensile tests and traditional fractography methods were utilized to study the fracture behavior of proton-irradiated HT-9 steels. First, to evaluate the viability of micro-tensile tests for nuclear material qualification process, meso-tensile tests on as-received HT-9 steels were performed. Fracture mechanisms of unirradiated HT-9 steels at both length scales were compared and underlying mechanisms discussed. The direct comparison of micro- and meso-scale data shows a distinctive size effect demonstrated by the increase in yield stress (YS). Upon completion of initial assessment, specimens were irradiated with 4 MeV + protons to three fluences, all of which were lower than 0.01 displacements per atom (dpa). As expected, the YS increases with irradiation. However, at 7 × 10 -3 dpa, the reversal of the trend was observed, and the YS exhibited sharp decline. We demonstrate that at lower length scales, grain structure has a more profound impact on the mechanical properties of irradiated materials, which provides information needed to fill in the gap in current understanding of the HT-9 fracture at different length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Injection of oxygen vacancies in the bulk lattice of layered cathodes

Surfaces, interfaces, and grain boundaries are, classically, known to be sinks at which defects generated within bulk lattice will preferentially migrate to and get annihilated, therefore rejuvenating a crystalline lattice to restore the intended functionalities.1, 2, 3 Here we report an inverse case by which the defects generated at the particle surface are continuously pumped into the bulk lattice and consequentially compromising the intended functionalities. We discovered that during the operation of a rechargeable battery, oxygen vacancies generated at cathode particle surfaces migrate toward inside of a perfect lattice in lithium-rich layered cathodes, featuring electrochemical driven surface vacancy creation, injection and condensation phenomenon and revealing an unrecognized surface process induced bulk lattice degradation featured with void formation and lattice reorganization. This process is associated with a high-cutoff voltage at which the anionic redox is activated. Coupling with calculations, we substantiate that triggering of anionic redox leads to sharply lowering of both formation energy of oxygen vacancy and migration barrier of oxidized oxide ions, therefore enabling migration of oxygen vacancy via oxidized oxide ion and leading to the injection of oxygen vacancy generated at the particle surface into the bulk lattice. This work unveils unpredicted behavior of oxidized oxide ion, leading to essential insight for designing of high capacity cathode for better batteries through coupled cationic and anionic redox dynamics in layered cathodes. In a broad term, this work provides insight on the intimate coupling between interfacial process and bulk lattice behavior for governing electrode functionality.

Yan, Pengfei↗
DOE OSTI2024

On the computation of moments in the Super-Transition-Arrays model for radiative opacity calculations

In the Super-Transition-Array statistical method for the computation of radiative opacity of hot dense matter, the moments of the absorption or emission features involve partition functions with reduced degeneracies, occurring through the calculation of averages of products of subshell populations. Here, in the present work, we discuss several aspects of the computation of such peculiar partition functions, insisting on the precautions that must be taken in order to avoid numerical difficulties. In a previous work, we derived a formula for supershell partition functions, which takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients for which we obtained a recursion relation and an explicit formula. We show that such an expansion can be combined with the recurrence relation for shifted partition functions. We also propose, neglecting the effect of fine structure as a first step, a positive-definite formula for the Super-Transition-Array moments of any order, providing an insight into the asymmetry and sharpness of the latter. The corresponding formulas are free of alternating sums. Several ways to speed up the calculations are also presented.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

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

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2018

Features of Short-Period Oscillations of Microwave Emission from the Solar Active Region before the Flare

The radiation of microwave sources above sunspots at a frequency of 17 GHz gives information about the parameters of solar plasma in the regions where the magneti-field strength is B ~ 2000 G in the transition region between the chromosphere and corona. Short-period oscillations (with a period of several minutes) of microwave emission from solar active regions (ARs) reflect wave processes in magnetic flux tubes of sunspots. Short-period oscillations of microwave emission from AR NOAA 12242 before two flares on December 17, 2014 are analyzed. This analysis is based on solar radio images obtained by means of the Nobeyama Radio Heliograph with a 10″−15″ two-dimensional spatial resolution. The radio maps of the whole solar disk were synthesized in a nonstandard mode with a cadence of 10 s and an averaging time of 10 s. An increase in the power of about ten-minute oscillations of microwave radiation approximately 40 to 50 min before the M1.5 flare (01: 00 UT) is found. On the same day, an increase in the power of ten-minute oscillations is observed about 60 min before the M8.7 flare (04: 42 UT). This effect is similar to the effect found earlier by two groups of authors independently for three-minute oscillations—namely, they observed a sharp increase in three-minute oscillations 15 to 20 min before the radio burst accompanying the flare. The effect in question may be interpreted as an relationship ofMHD waves propagating along the magnetic flux tube of a sunspot and the onset of the solar flare.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Overcharge investigation of large format lithium-ion pouch cells with Li(Ni 0.6 Co 0.2 Mn 0.2 )O 2 cathode for electric vehicles: Thermal runaway features and safety management method

The overcharge-induced thermal runaway features of large format commercial lithium-ion batteries with Li(Ni 0.6 Co 0.2 Mn 0.2 )O 2 (NCM622) cathode for electric vehicles under different current rates (C-rates) have been systematically studied at ambient temperature. The overcharge process is characterized as four stages. The temperature rise and the maximum temperature of the battery surface don't increase in proportion to the applied C-rates. However, with the increase of C-rates, the crest voltage of voltage curve rises linearly. When the voltage reaches approximately 5.1 V, a new voltage plateau appears in the cases below 2C. It is not sufficient that the temperature sensor is placed only near the terminal tab for most battery packs of EVs. In addition, the accumulated heat analysis demonstrates that side reactions dominate the temperature rise and contribute to most of the accumulated heat before thermal runaway. To mitigate the impact of overcharge and avoid the thermal runaway risk, a safety management method is proposed. Furthermore, the sharp drop in voltage before thermal runaway also provides a feasible approach to forewarn the users of the impending risk. These results are important for building safer batteries and providing information for the safety monitoring function of the battery management system (BMS).

25 ENERGY STORAGE↗
DOE OSTI2023

Unusual magnetotransport in twisted bilayer graphene from strain-induced open Fermi surfaces

Anisotropic hopping in a toy Hofstadter model was recently invoked to explain a rich and surprising Landau spectrum measured in twisted bilayer graphene away from the magic angle. Suspecting that such anisotropy could arise from unintended uniaxial strain, we extend the Bistritzer–MacDonald model to include uniaxial heterostrain and present a detailed analysis of its impact on band structure and magnetotransport. We find that such strain strongly influences band structure, shifting the three otherwise-degenerate van Hove points to different energies. Coupled to a Boltzmann magnetotransport calculation, this reproduces previously unexplained nonsaturating B 2 magnetoresistance over broad ranges of density near filling ν = ± 2 and predicts subtler features that had not been noticed in the experimental data. In contrast to these distinctive signatures in longitudinal resistivity, the Hall coefficient is barely influenced by strain, to the extent that it still shows a single sign change on each side of the charge neutrality point—surprisingly, this sign change no longer occurs at a van Hove point. The theory also predicts a marked rotation of the electrical transport principal axes as a function of filling even for fixed strain and for rigid bands. More careful examination of interaction-induced nematic order versus strain effects in twisted bilayer graphene could thus be in order.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Mechanism and Kinetics for Reaction of the Chemical Warfare Agent Simulant, DMMP( g ), with Zirconium(IV) MOFs: An Ultrahigh-Vacuum and DFT Study

The mechanism and kinetics of interactions between dimethyl methylphosphonate (DMMP), a key chemical warfare agent (CWA) simulant, and Zr 6 -based metal organic frameworks (MOFs) have been investigated with in situ infrared spectroscopy (IR), X-ray photoelectron spectroscopy (XPS), powder X-ray diffraction (PXRD), and DFT calculations. DMMP was found to adsorb molecularly to UiO-66 through the formation of hydrogen bonds between the phosphoryl oxygen and the free hydroxyl groups associated with Zr 6 nodes on the surface of crystallites and not within the bulk MOF structure. Unlike UiO-66, the infrared spectra for UiO-67 and MOF-808, recorded during DMMP exposure, suggest that uptake occurs through both physisorption and chemisorption. The XPS spectra of MOF-808 zirconium 3d electrons reveal a charge redistribution following exposure to DMMP. In addition, analysis of the phosphorus 2p electrons following exposure and thermal annealing to 600 K indicates that two types of stable phosphorus-containing species exist within the MOF. DFT calculations, used to guide the IR band assignments and to help interpret the XPS features, suggest that uptake is driven by nucleophilic addition of an OH group to DMMP with subsequent elimination of a methoxy substituent to form strongly bound methyl methylphosphonic acid (MMPA). The rates of product formation indicate that there are likely two distinct uptake processes, requiring rate constants that differ by approximately an order of magnitude. However, the rates of molecular uptake were found to be nearly identical to the rates of reaction, which strongly suggests that the reaction rates are diffusion-limited. Here, the final products were found to inhibit further reactions within the MOFs, and these products could not be thermally driven from the MOFs prior to decomposition of the MOFs themselves.

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

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗
DOE OSTI2024

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems—like fixation of cells coupled with fluorescence microscopy—provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. In this study, these snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

OzDES Reverberation Mapping Programme: the first Mg ii lags from 5 yr of monitoring

ABSTRACT Reverberation mapping is a robust method to measure the masses of supermassive black holes outside of the local Universe. Measurements of the radius–luminosity (R−L) relation using the Mg ii emission line are critical for determining these masses near the peak of quasar activity at z ≈ 1−2, and for calibrating secondary mass estimators based on Mg ii that can be applied to large samples with only single-epoch spectroscopy. We present the first nine Mg ii lags from our 5-yr Australian Dark Energy Survey reverberation mapping programme, which substantially improves the number and quality of Mg ii lag measurements. As the Mg ii feature is somewhat blended with iron emission, we model and subtract both the continuum and iron contamination from the multiepoch spectra before analysing the Mg ii line. We also develop a new method of quantifying correlated spectroscopic calibration errors based on our numerous, contemporaneous observations of F-stars. The lag measurements for seven of our nine sources are consistent with both the H β and Mg ii R−L relations reported by previous studies. Our simulations verify the lag reliability of our nine measurements, and we estimate that the median false positive rate of the lag measurements is $4{{\ \rm per\ cent}}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Prediction of DIII-D Pedestal Structure from Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. An experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (ne) and electron temperature (Te) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (Ip), toroidal magnetic field (Bφ), neutral beam heating power (PNBI) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. Here, an experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (n e ) and electron temperature (T e ) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (I p ), toroidal magnetic field (B Φ ), neutral beam heating power (P NBI ) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Coherent Vibrational Wavepacket Dynamics in Platinum(II) Dimers and Their Implications

Vibrational coherence in the metal-metal-to-ligand-charge transfer (MMLCT) excited state of cyclometalated platinum dimers with a pseudo C 2 symmetry was investigated where two nearly degenerate transitions from HOMO (metal-metal σ* orbital) to higher energy ligand π* orbitals could be simultaneously induced. We observed oscillatory features in femtosecond degenerate transient absorption (TA) signals from complexes [(ppy)Pt(μ- t Bu 2 pz)] 2 (1) and anti-[(ppy)Pt(μ-pyt)] 2 (2), which are attributed to coherent nuclear motions that modulate the HOMO (antibonding σ*) energy level, and hence the energy for the MMLCT transition. The characteristics of such coherent nuclear motions, such as the oscillatory frequency and the dephasing time, differ between 1 and 2 and are explained by mainly two structural factors that could influence the vibrational coherence, the Pt-Pt distance (2.97 Å for 1 vs. 2.85 Å for 2) and molecular shape (1 in an “A” frame vs. 2 in an “H” frame). Because the electronic coupling between the two Pt atoms determines the energy splitting of the bonding σ and antibonding σ* orbital, the Pt-Pt stretching mode coupled to the MMLCT transition changes the inter Pt distance from that of the ground state. Interestingly, while 1 shows a single PtPt stretching frequency of 120 cm-1 in the MMLCT state, 2 exhibits multiple downshifted frequencies (80 and 105 cm-1) in the MMLCT state along with a shorter vibrational dephasing time than 1. Based on the ground state optimized structures and Raman calculations, the changes evident in the vibrational wavepaket dynamics in 2 are closely correlated with the “H” framed geometry in 2 compared to the “A” frame in 1, leading to the sharp increase in - interaction between ppy ligands. Although the TA experiments do not directly reveal the ultrafast intersystem crossing (ISC) because of a strong coherent spike at early timescales, the dependence of the vibrational wavepacket dynamics on molecular geometry can be understood based on previously proposed potential energy surfaces (PES) as function of PtPt distance, suggesting that the interaction between the cyclometalating ligands can be a key factor in determining the PtPt shortening and the related energy relaxation dynamics in the Pt(II) dimers. Further experiments using femtosecond broadband TA and X-ray scattering spectroscopy are planned to investigate directly the ISC and PtPt contraction, to support the relationship between ground state molecular geometry and photoinduced structural changes in the Pt(II) dimers.

Kim, Pyosang↗