Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR FLOW”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Adaptive simulations enable computational design of electron beam processing of nanomaterials with supersonic micro-jet precursor

Focused Electron Beam Induced Processing (FEBIP) is a powerful tool for the “direct-write” of nanomaterials with the possibility of atomistic control on suspended 2D material substrates. FEBIP capabilities have been significantly expanded by using a localized jet-based delivery of precursors, and especially when a thermally energized supersonic micro-jet enhances the delivery of mass flux along with controlling the far-from-equilibrium thermodynamic state of adsorbed adatoms. The possibilities of growing nanomaterials with “dialed-in” composition with ultra-high growth rates and aspect ratios up to 100:1 have been demonstrated using the supersonic micro-jet-FEBIP. Bringing this scientific discovery to the level of maturity required for practical applications in additive nanomanufacturing requires simulation tools that are capable of capturing the complex flow physics of micro-jet-substrate interactions that bridge a wide range of flow regimes from the high-density gas micro-jet expanding into the vacuum environment of FEBIP. To address this significant computational challenge, a new approach has been developed and described in this work for a multiscale adaptive DSMC (Direct Simulation Monte Carlo) algorithm to predict the micro-jet gas dynamics in an FEBIP environment, spanning the full range of flow regimes from low Knudsen (Kn) number O(0.01) continuum flow to high Kn of O(10) for the molecular flow within a unified computational framework. Here, the fundamental principles of the adaptive DSMC algorithm are described, its viability as a DSMC technique is demonstrated, and computational improvements for the benchmark cases are discussed in the context of advancing 3D nanofabrication with FEBIP. Ultimately, combining the first principle simulations via adaptive DSMC with the complementary experimental data will enable the creation of the powerful CAD tools for in silico design and optimal operation of micro-jet-FEBIP.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carbon nanotube nanofluidics

Fluid flow under extreme spatial confinement exhibits unusual physical behaviors. This nanofluidic transport regime is relevant to a variety of mass transport, separation, and energy production processes in biological and industrial systems. Carbon nanotubes (CNTs) offer a nearly ideal platform for exploring nanofluidic transport because of their extremely narrow, smooth, hydrophobic inner pores, which enable very fast molecular flow while providing strong selectivity. In this review, we aim to provide a comprehensive understanding of nanofluidics in CNTs, focusing on the basic physics of mass transport in CNTs, various experimental platforms developed to investigate these phenomena, and key results on the permeation of water, protons, and ions. We focus on the critical factors that influence transport efficiency and selectivity, such as slip flow and charge regulation in CNTs, and the roles of entrance effects, dehydration processes and ion–charge interactions at the CNT entrances. We also explore the confinement effects, highlighting how the unique one-dimensional structure of CNTs imposes distinct constraints on fluid behavior and leads to novel single-file transport phenomena. Finally, we address current challenges and future directions of CNT nanofluidics.

Li, Zhongwu [Lawrence Livermore National Laborator↗
DOE OSTI2020

Gas Transport in the Insensitive High Explosive PBX 9502

Small-scale laboratory experiments were performed to analyze gas transport mechanics in PBX 9502, an insensitive high explosive (IHE) composition. Two independent investigations are reported here. First, gas permeametry and dynamic pycnometry techniques were used to measure the molecular flow (Knudsen) coefficient and the internally interconnected void fraction of pristine and thermally damaged samples in two uniaxial pressing orientations. The permeability of PBX 9502 was found to be unmeasurably low (of the same order of magnitude as PTFE) with gas transport being diffusion-dominated. Secondly, a pressure vessel experiment was developed to measure quasi-static and dynamic gas generation as the explosive was heated to self-ignition (cookoff). The gas generation results and the permeability/diffusivity findings provide evidence that PBX 9502 remains impermeable until seconds prior to self-ignition. At ignition, internal void-pressure drives macro-scale cracking and the sample becomes uniformly incorporated in the following deconsolidating deflagration. These results are discussed within the context of previous observations of pressure-dependent cookoff behavior of PBX 9502 and provide a more complete description of thermal damage evolution in this explosive composition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

First operation and validation of simulations for the divertor cryo-vacuum pump in Wendelstein 7-X

Ten cryo-vacuum pumps (CVPs) were installed in the subdivertor region of each island divertor in the stellarator Wendelstein 7-X (W7-X) and operated for the first time during the recently completed plasma campaign OP2.1. A pumping speed of 70 ± 1 $\frac{m^3}{s}$ was measured during dedicated tests with known hydrogen gas injection. Based on a conductance model, the estimated pumping speed ranges from 86-93 $\frac{m^3}{s}$ for different sticking coefficients between 0.6 and 0.8. After completion of the initial tests the CVPs were operated successfully throughout the campaign, with regeneration performed once a week. Neutral gas pressures in the subdivertor in the range of 10 −4 mbar are well within the molecular flow regime and limit the particle exhaust capabilities of the CVPs. Simulations of the neutral gas pressure in the three-dimensional complex geometry of the subdivertor were performed using the DIVGAS code based on the direct simulation Monte Carlo method and a model implemented in the steady-state thermal package in ANSYS, which are in agreement with the measured values during plasma operation.

Cryo-vacuum pumping↗
DOE OSTI2020

Effect of 60 Co γ-radiation on the tetrafluoroethylene/perfluoro(methyl vinyl ether) copolymer

The chemical and paramagnetic composition of the surface of the tetrafluoroethylene / perfluoro(methyl vinyl ether) copolymer (MFA) was investigated before and after γ-irradiation. The content of fluorine and carbon decreases at radiolysis of MFA, as well as the proportion of oxygen increases as a result of the oxidation reaction of radicals formed during γ-irradiation. The electron paramagnetic resonance (EPR) spectrum of the copolymer irradiated at 77 K contains the radicals ~CF 2 C∙FCF 2 ~ (1) and >CFOC∙F 2 (2). The absence of ~CF 2 C∙F 2 radicals, which can be formed at radiolytic rupture of the C—C bond of the copolymer main chain is associated with the effect of recombination of the resulting radical pair. The EPR spectrum of irradiated MFA at 300 K only recorded radical (1), which transform into peroxide radicals when air oxygen is introduced into the irradiated sample. Thermomechanical curves of native and radiolyzed MFA show three structural states of the amorphous block (glass, transition region, plateau of highly elastic deformation), as well as melting and molecular flow of the copolymer. The radio thermoluminescence curve shows five maxima associated with the emission of light during the recombination of stabilized ions and radicals during heating of the pre-irradiated MFA at 77 K. The temperature of the luminescence peak maximum at 336 K is close to the temperature of the onset of the devitrification of the chains of the amorphous block of MFA on the thermomechanical curve at 341 K. Irradiation leads to an increase in crystallinity of MFA. The heat of fusion of irradiated MFA on the Differential Scanning Calorimeter (DSC) curve is 50 % higher than the fusion heat of native copolymer. During the heating of γ-irradiated MFA, the relative content of COF molecules in the thermal decomposition products increases by an order of magnitude.

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

Time reversal symmetry breaking and odd viscosity in active fluids: Green–Kubo and NEMD results

Active fluids, which are driven at the microscale by non-conservative forces, are known to exhibit novel transport phenomena due to the breaking of time reversal symmetry. Recently, Epstein and Mandadapu [arXiv:1907.10041 (2019)] obtained Green-Kubo relations for the full set of viscous coefficients governing isotropic chiral active fluids, including the so-called odd viscosity, invoking Onsager's regression hypothesis for the decay of fluctuations in active non-equilibrium steady states. In this Communication, we test these Green-Kubo relations using molecular dynamics simulations of a canonical model system consisting of actively torqued dumbbells. We find the resulting odd and shear viscosity values from the Green-Kubo relations to be in good agreement with values measured independently through non-equilibrium molecular dynamics flow simulations. This provides a test of the Green-Kubo relations and lends support to the application of the Onsager regression hypothesis in relation to viscous behaviors of active matter systems.

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

Ultra-High Vacuum Outgassing Characterization of Thermally Processed Low-Carbon Steel for Advanced Particle Accelerator and Gravitational Wave Detector Applications

This dissertation investigated AISI 1020 low-carbon steel as an alternative vacuum chamber material to conventional stainless steel for ultra-high vacuum (UHV) and extreme-high vacuum (XHV) applications. After a 400 °C/48 h bake, AISI 1020 tube chambers achieved a hydrogen outgassing rate of 2.4 × 10¿¹6 Torr·L·s¿¹·cm¿², approximately 2,300 times lower than the prebaked 316L stainless-steel comparator, among the lowest hydrogen outgassing rates ever reported for an uncoated metallic vacuum chamber. Bare and magnetite-coated AISI 1020 chambers were then compared using throughput and rate-of-rise methods. The magnetite coating yielded 5× lower water outgassing at room temperature, but this advantage disappeared after 80 °C baking. After full thermal conditioning (400 °C/48 h prebake followed by 150 °C/96 h and 200 °C/110 h), bare steel achieved 25× lower hydrogen outgassing than the magnetite-coated chamber (9.6 × 10¿¹6 Torr·L·s¿¹·cm¿²) and >99% H2 purity with carbon species below RGA detection. Monte Carlo molecular flow simulations of a CEBAF photogun beamline (96 scenarios) showed that replacing 304L stainless steel with AISI 1020 reduces equilibrium H2 pressure by a factor of 833; a single 304L electrode contributes 98.8% of the gas load despite occupying only 9.1% of the internal surface area. A 500-m Einstein Telescope beampipe screening showed that corrugated bellows contribute 18% of the gas load from only 0.7% of the surface area. A five-model adsorption isotherm framework applied to 22 pumpdown datasets (164 fits with AR(1)-GLS correction) established that the experimental protocol, not the material, controls isotherm identifiability: Dubinin–Radushkevich wins isothermal pipe pumpdowns; Langmuir wins thermally dominated chamber bakes. Cross-dataset joint fitting of the AISI 1020 pipe pumpdowns yielded an H2 diffusion activation energy Ed = 7.24 ± 1.28 kcal·mol¿¹, consistent with trap dominated diffusion in commercial low-carbon steels. Two companion innovations were developed: a Variable Conductance Device (VCD, patent pending IDF-00723) for XHV outgassing measurement, and VacuumDesignerPro (VDP), a MATLAB-based design tool validated against LIGO benchmarks.

Al-Allaq, Aiman H [Old Dominion University]↗
DOE OSTI2020

Realization of an Asymmetric Non‐Aqueous Redox Flow Battery through Molecular Design to Minimize Active Species Crossover and Decomposition

Abstract This communication presents a mechanism‐based approach to identify organic electrolytes for non‐aqueous redox flow batteries (RFBs). Symmetrical flow cell cycling of a pyridinium anolyte and a cyclopropenium catholyte resulted in extensive capacity fade due to competing decomposition of the pyridinium species. Characterization of this decomposition pathway enabled the rational design of next‐generation anolyte/catholyte pairs with dramatically enhanced cycling performance. Three factors were identified as critical for slowing capacity fade: (1) separating the anolyte–catholyte in an asymmetric flow cell using an anion exchange membrane (AEM); (2) moving from monomeric to oligomeric electrolytes to limit crossover through the AEM; and (3) removing the basic carbonyl moiety from the anolyte to slow the protonation‐induced decomposition pathway. Ultimately, these modifications led to a novel anolyte–catholyte pair that can be cycled in an AEM‐separated asymmetric RFB for 96 h with >95 % capacity retention at an open circuit voltage of 1.57 V.

Shrestha, Anuska↗
DOE OSTI2024

Properties of molecular clumps and cores in colliding magnetized flows

ABSTRACT We simulate the formation of molecular clouds in colliding flows of warm neutral medium with the adaptive mesh refinement code flash in eight simulations with varying initial magnetic field strength, between 0.01–5 μG. We include a chemical network to treat heating and cooling and to follow the formation of molecular gas. The initial magnetic field strength influences the fragmentation of the forming cloud because it prohibits motions perpendicular to the field direction and hence impacts the formation of large-scale filamentary structures. Molecular clump and core formation occurs anyhow. We identify 3D clumps and 3D cores, which are defined as connected, CO-rich regions. Additionally, 3D cores are heavily shielded. While we do not claim those 3D objects to be directly comparable to observations, this enables us to analyse their full virial state. With increasing field strength, we find more fragments with a smaller average mass; yet the dynamics of the forming clumps and cores only weakly depends on the initial magnetic field strength. The molecular clumps are mostly unbound, probably transient objects, which are weakly confined by ram pressure or thermal pressure, indicating that they are swept up by the turbulent flow. They experience significant fluctuations in the mass flux through their surface, such that the Eulerian reference frame shows a dominant time-dependent term due to their indistinct nature. We define the cores to encompass highly shielded molecular gas. Most cores are in gravitational-kinetic equipartition and are well described by the common virial parameter $\alpha _\mathrm{vir}$, while some undergo minor dispersion by kinetic surface effects.

Weis, M. (ORCID:0000000256838860)↗
DOE OSTI2020

A combined SECM and electrochemical AFM approach to probe interfacial processes affecting molecular reactivity at redox flow battery electrodes

Redox flow batteries are attractive technologies for grid energy storage since they use solutions of redox-active molecules that enable a superior scalability and the decoupling of power and energy density. However, the reaction mechanisms of the redox active components at RFB electrodes are complex, and there is currently a pressing need to understand how interfacial processes impact the kinetics and operational reversibility of RFB systems. Here, we developed a combined electrochemical imaging methodology rooted in scanning electrochemical microscopy (SECM) and atomic force microscopy (AFM) for exploring the impact of electrode structure and conditioning on the electron transfer properties of model redox-active dialkoxybenzene derivatives, 2,5-di-tert-butyl-1,4-bis(2-methoxyethoxy)benzene (C1) and 2,3-dimethyl-1,4-dialkoxybenzene (C7). Using AFM and secondary-ion mass spectrometry (SIMS), we observed the formation of interfacial films with distinct mechanical properties compared to those of cleaved graphitic surfaces, and exclusively during reduction of electrogenerated radical cations. These films had an impact on the median rate and distribution of the electron transfer rate constant at the basal plane of multilayer and single layer graphene electrodes, displaying kinetically-limited values that did not yield the activation expected per the Butler–Volmer model with a transfer coefficient ~0.5. These changes were dependent on redoxmer structure: furthermore, SECM showed strong attenuation of C7 kinetics by a surface layer on MLG and SLG, while C1 kinetics were only affected by SLG. SECM and AFM results together show that these limiting films operate exclusively on the basal plane of graphite, with the edge plane showing a relative insensitivity to cycling and operation potential. Ultimately, this integrated electrochemical imaging methodology creates new opportunities to understand the unique role of interfacial processes on the heterogeneous reactivity of redoxmers at electrodes for RFBs, with a future role in elucidating phenomena at high active concentrations and spatiotemporal variations in electrode dynamics.

25 ENERGY STORAGE↗
DOE OSTI2022

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

Here the flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear-stress ratio μ are relatively insensitive to system size for a wide range of inertial numbers I. Second-order effects in strain rate, such as normal stress differences, require large system sizes to accurately extract meaningful results, notably a nonmonotonic dependence in the first normal stress difference with strain rate. The rheological response represented by the μ(I) scalar model works well at describing the lower-order aspects of the rheology, except near the quasistatic limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems N≤1000 near the quasistatic limit can lead to the arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasistatic and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number, and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasistatic flow.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗
DOE OSTI2025

A spline-based method to obtain spatially dependent viscosity in confined flows

Coupling chemical physics to continuum theories is a critical step to understanding multi-scale phenomena. This paper will connect non-equilibrium molecular dynamics simulations to a continuum-based Navier-Stokes equation that has relaxed the assumption of spatial uniformity in viscosity. Using a form for viscosity based on spline interpolation, viscosity as a function of position is obtained from the least squares fit of the velocity profile measured from molecular simulations of flow in a nanochannel. Viscosity can vary widely, particularly near the channel boundaries, indicating that uniform viscosity is no longer appropriate. Variations of the viscosity near the channel surfaces imply that considering solution and surface chemistry could be necessary to rigorously understand molecular-scale flows in nanochannels.

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

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗