Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerated molecular dynamics”

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 343 records · Page 19

Sulfur-enhanced dynamics of coinage metal(111) surfaces: Step edges versus terraces as locations for metal-sulfur complex formation

The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.

36 MATERIALS SCIENCE↗

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Structural Insights into the Dynamics of Water in SOD1 Catalysis and Drug Interactions

Superoxide dismutase 1 (SOD1) is a crucial enzyme that protects cells from oxidative damage by converting superoxide radicals into H 2 O 2 and O 2 . This detoxification process, essential for cellular homeostasis, relies on a precisely orchestrated catalytic mechanism involving the copper cation, while the zinc cation contributes to the structural integrity of the enzyme. This study presents the 2.3 Å crystal structure of human SOD1 (PDB ID: 9IYK), revealing an assembly of six homodimers and twelve distinct active sites. The water molecules form a complex hydrogen-bonding network that drives proton transfer and sustains active site dynamics. Our structure also uncovers subtle conformational changes that highlight the intrinsic flexibility of SOD1, which is essential for its function. Additionally, we observe how these dynamic structural features may be linked to pathological mutations associated with amyotrophic lateral sclerosis (ALS). By advancing our understanding of hSOD1’s mechanistic intricacies and the influence of water coordination, this study offers valuable insights for developing therapeutic strategies targeting ALS. Our structure’s unique conformations and active site interactions illuminate new facets of hSOD1 function, underscoring the critical role of structural dynamics in enzyme catalysis. Moreover, we conducted a molecular docking analysis using SOD1 for potential radical scavengers and Abelson non-receptor tyrosine kinase (c-Abl, Abl1) inhibitors targeting misfolded SOD1 aggregation along with oxidative stress and apoptosis, respectively. The results showed that CHEMBL1075867, a free radical scavenger derivative, showed the most promising docking results and interactions at the binding site of hSOD1, highlighting its promising role for further studies against SOD1-mediated ALS.

60 APPLIED LIFE SCIENCES↗

Orientation-dependent enhanced ionization in acetylene revealed by ultrafast cross-polarized pulse pairs

We investigate the orientation dependence of enhanced ionization (EI) during strong-field-driven nuclear motion in acetylene (C 2 ⁢H 2 ). Here, we both initiate and probe molecular dynamics in acetylene with intense 6-fs cross-polarized pulse pairs, separated by a variable delay. Following multiple ionization by the first pulse, acetylene undergoes simultaneous elongation of the carbon-carbon and carbon-hydrogen bonds, enabling further ionization by the second pulse and the formation of a very highly charged state, [C 2 ⁢H 2 ] 6+ . At small interpulse delays (< 20 fs), this enhancement occurs when the molecule is aligned to the probe pulse. Conversely, at large delays (> 40 fs), formation of [C 2 ⁢H 2 ] 6+ occurs when the molecule is aligned to the pump pulse. By analyzing the polarization and time dependence of sequentially ionized [C 2 ⁢H 2 ] 6+ , we resolve two transient alignments that both contribute to a large increase in the multiple ionization yield. In conclusion, this cross-polarized pulse pair scheme uniquely enables selective probing of deeply bound orbitals, providing new insights on orientation-dependent EI in highly charged hydrocarbons.

Atomic & molecular clusters↗

The Impact of Naphthenic Acids on Dynamic Fluid–Fluid Interactions: Implication for Enhanced Oil Recovery

Previous coreflooding results and wettability analyses in our group show that injection of naphthenic-acid-enriched water can improve oil recovery over traditional waterflooding. This observation is still a subject of research efforts without a definitive explanation. Naphthenic acids (NA) have been reported to drive wettability alteration and increase the water–oil interface elasticity. These alterations depend on the NA carbon number and aqueous-phase salinity, among other conditions, as reported in the literature. Smart-water flooding (SWF) research often links recovery to the initial wettability condition, being higher for initially oil-wet rock. SWF refers to a technique in which the aqueous-phase ion composition or/and salinity are changed to maximize oil recovery. Given NAs’ complex solution behavior, selecting acid combinations that prompt oil recovery is a difficult objective. The aim of this research is to determine the effects of select naphthenic acids on the oil–water interfacial rheology and wettability alteration and how these interfacial effects are associated with oil recovery under spontaneous imbibition. NAs were selected based on their carbon number, molecular structure, and solubility in the saline solution used in this research. We aimed at exploring which NAs should be used to regulate interfacial properties so as to either increase oil recovery or accelerate production. Time-domain nuclear magnetic resonance, interfacial dilatational rheology, and liquid-bridge experiments, i.e., proxy of snap-off, were conducted. A baseline was established using results obtained with a previously tested sulfate-rich aqueous phase, shown to be effective in recovering oil. Results show that NA14 and N18 increase the water–oil interfacial viscoelasticity and induce interfacial healing but led to different recovery factors. N10, while effective at inducing water wetness in oil-wet rock, is ineffective at increasing the recovery factor. We concluded that wettability and oil–water interfacial rheology are not exclusive, and instead they can synergistically favor EOR benefits. Moreover, oil recovery benefits under spontaneous imbibition are shown to depend strongly on the initial wettability conditions.

Medina-Rodriguez, Bryan X. (ORCID:0000000262922831↗

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic’s evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (M pro ), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. Here, this work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic “warhead” to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of M pro with micromolar affinities (K I of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

60 APPLIED LIFE SCIENCES↗

Molecular-scale mechanisms of CO 2 mineralization in nanoscale interfacial water films

The calamitous impacts of untethered carbon emission from fossil-fuel-burning energy infrastructure calls for accelerated development of large-scale CO 2 capture, utilization, and storage technologies that are underpinned by fundamental understanding of molecular-level chemical processes. In the subsurface, rocks rich in divalent metals can react with CO 2 , permanently sequestering it in the form of stable metal carbonate minerals, with the CO 2 -H 2 O composition of the post-injection pore fluid acting as a primary control variable. Herein, we compare mechanistic reaction pathways for aqueous-mediated carbonation with carbon mineralization occurring in nanoscale adsorbed water films. In the extreme of pores filled with a CO 2 -dominant fluid, carbonation reactions are confined to nanometer-thick water films coating mineral surfaces, which enable metal cation release, transport, nucleation, and crystallization of metal carbonate minerals. Though seemingly counterintuitive, laboratory studies have demonstrated facile carbonation rates in these low water environments, which in recent years has begun to be better understood in mechanistic detail. The overarching objective of this review is to delineate the unique underlying molecular-scale reaction mechanisms that govern CO 2 mineralization in these reactive and dynamic quasi-2D interfaces. Here. we highlight the importance of understanding unique properties in thin water films, such as how water dielectric properties, and consequently ion solvation/hydration behavior, can change under nanoconfinement. We conclude by identifying important frontiers for future work and opportunities to exploit these fundamental chemical insights for decarbonization technologies in the 21st century.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In-vivo neuronal dysfunction by Aβ and tau overlaps with brain-wide inflammatory mechanisms in Alzheimer’s disease

The molecular mechanisms underlying neuronal dysfunction in Alzheimer’s disease (AD) remain uncharacterized. Here, we identify genes, molecular pathways and cellular components associated with whole-brain dysregulation caused by amyloid-beta (Aβ) and tau deposits in the living human brain. We obtained in-vivo resting-state functional MRI (rs-fMRI), Aβ- and tau-PET for 47 cognitively unimpaired and 16 AD participants from the Translational Biomarkers in Aging and Dementia cohort. Adverse neuronal activity impacts by Aβ and tau were quantified with personalized dynamical models by fitting pathology-mediated computational signals to the participant’s real rs-fMRIs. Then, we detected robust brain-wide associations between the spatial profiles of Aβ-tau impacts and gene expression in the neurotypical transcriptome (Allen Human Brain Atlas). Within the obtained distinctive signature of in-vivo neuronal dysfunction, several genes have prominent roles in microglial activation and in interactions with Aβ and tau. Moreover, cellular vulnerability estimations revealed strong association of microglial expression patterns with Aβ and tau’s synergistic impact on neuronal activity (q < 0.001). These results further support the central role of the immune system and neuroinflammatory pathways in AD pathogenesis. Neuronal dysregulation by AD pathologies also associated with neurotypical synaptic and developmental processes. In addition, we identified drug candidates from the vast LINCS library to halt or reduce the observed Aβ-tau effects on neuronal activity. Top-ranked pharmacological interventions target inflammatory, cancer and cardiovascular pathways, including specific medications undergoing clinical evaluation in AD. Our findings, based on the examination of molecular-pathological-functional interactions in humans, may accelerate the process of bringing effective therapies into clinical practice.

60 APPLIED LIFE SCIENCES↗

Structural and Mechanistic Advances in the Chemistry of Methyl-Coenzyme M Reductase (MCR)

Methane represents 34% of U.S. energy consumption and is a major greenhouse gas related to the global carbon cycle and energy production. However, current industrial practices significantly increase atmospheric methane levels, necessitating a deeper understanding of its biosynthesis and oxidation. Methyl-coenzyme M reductase (MCR) is central to biological methane metabolism, catalyzing the final step of methanogenesis and the first step in anaerobic methane oxidation. It is also a key target for strategies to capture and transform methane into value-added chemicals.The active site of MCR is a buried Ni-based cofactor only accessible by the substrates via a 50 Å long tunnel. Although the Ni(I) state is required to initiate catalysis, capturing this state remains a challenge for the current structural techniques. Recent advances in structural biology using X-ray Free-Electron Laser serial crystallography have provided insights into MCR's inactive Ni(II) state at room temperature and show promise for capturing its active Ni(I) form.Our team has established several critical aspects of the MCR mechanism using a combination of experimental and computational studies. MCR uses CH 3 -SCoM and CoBSH as substrates, producing methane and a disulfide product CoMSSCoB. Kinetic analysis showed that productive substrate binding requires CH 3 -SCoM to bind first, inducing conformational changes that optimize the active site for subsequent CoBSH binding. Following substrate binding, four proposed methane production/oxidation mechanisms were examined, establishing whether the reaction proceeds through an organometallic methyl-nickel(III), methyl anion ion, or methyl radical intermediate. Experimental measurements using CoBSH analogs successfully slowed the reaction, allowing for mechanistic insight that demonstrated the methyl radical pathway, where the initial interactions involve homolytic cleavage of the methyl-sulfur bond, generating a methyl radical that quickly abstracts the thiol hydrogen atom of CoBSH to form methane. Computational studies further confirmed that, compared to other mechanisms, the methyl radical mechanism is thermodynamically more favorable and accessible under physiological conditions.Spectroscopic and computational studies challenged the conventional understanding of substrate binding in MCR by proposing an alternative positioning of CH 3 -SCoM and CoMSSCoB in the active site pocket. The research suggested that CH 3 -SCoM (substrate) and CoMSSCoB (product) bind via their sulfonate groups to the Ni(I) center of cofactor F 430 . This binding allows for the reaction without substrate reorganization in the pocket but would require a long-range electron transfer. Altogether, the work summarized in this review reflects our current understanding of the enzyme's catalytic mechanism and structural dynamics. This is essential for developing efficient methane conversion technologies that could mitigate its environmental impact while harnessing energy-storage properties.

03 NATURAL GAS↗

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗

Evolution of carbapenemase activity in the class C β-lactamase ADC-1

Antibiotic resistance in bacteria poses a significant threat to public health. Among dozens of available antimicrobial agents, carbapenems are used as drugs of choice for the treatment of serious infections caused by pathogens resistant to other antibiotics. However, their usefulness has been severely compromised due to the emergence and wide spread of carbapenem-resistant clinical isolates worldwide. High-level resistance to carbapenems in bacteria is mediated by the production of β-lactamases from three molecular classes, A, B, and D, but not by class C enzymes. In this study, we selected a triple mutant of the intrinsic class C Acinetobacter-derived cephalosporinase ADC-1 (ADC-1 TM ) that confers high-level resistance to the carbapenems meropenem, ertapenem, and doripenem. Kinetic experiments demonstrated that the apparent binding affinity, along with the acylation and deacylation rates, were all improved for the mutant enzyme. X-ray crystallography, molecular docking, and molecular dynamics simulations revealed that the amino acid substitutions in ADC-1 TM produce significant changes in the enzyme active site architecture and binding mode of the carbapenem ertapenem. These changes allow for better positioning of a deacylating water for nucleophilic attack, thus explaining the significantly improved rate of ertapenem deacylation by ADC-1 TM . In this study, we showed for the first time that a class C β-lactamase can produce high-level resistance to carbapenem antibiotics, which underlines the potential for enzymes of this class to evolve such resistance and could further exacerbate the problem of antibiotic resistance in bacteria.

59 BASIC BIOLOGICAL SCIENCES↗

Distinguishing Electron and Hole Dynamics in Functionalized CdSe/CdS Core/Shell Quantum Dots Using Complementary Ultrafast Spectroscopies and Kinetic Modeling

The evolution of excitation energy and photogenerated charges in semiconductor quantum dots (QDs) functionalized with molecular acceptors can be probed on ultrafast time scales using techniques such as transient absorption (TA) spectroscopy. However, historical interpretations that the 1S(e)-1S 3/2 (h) transition in Cd-chalcogenide QDs is fully attributable to electrons may be misleading, and multiexponential models used to fit TA kinetics do not correspond directly to specific photophysical processes. Here, we present visible-wavelength and mid-IR TA and time-resolved photoluminescence measurements to inform a comprehensive kinetic model of the photoexcited CdSe/CdS core/shell QDs functionalized with passivating oleic acid (OA), hole-accepting ferrocene, or electron-accepting naphthalene bisimide (NBI). We show that ~30% of the 1S signal and 72% of the IR signal can originate from holes in well-passivated core/shell QDs. We also demonstrate evidence of electron trapping in OA-capped core/shell QDs, with additional electron transfer and hole trapping in the QDs functionalized with NBI. Electron (hole) trapping and detrapping occur in 450 ± 100 ps (430 ± 70 ps) and 340 ± 100 ps (1.1 ± 0.4 ns) respectively, while the time constant for electron transfer to NBI is ~1.8 ns. The comprehensive picture of photophysical processes provided by the complementary ultrafast techniques and kinetic modeling can accelerate both the fundamental science and application development of nanostructured and molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast structural dynamics of UV photoexcited cis , cis -1,3-cyclooctadiene observed with time-resolved electron diffraction

Conjugated diene molecules are highly reactive upon photoexcitation and can relax through multiple reaction channels that depend on the position of the double bonds and the degree of molecular rigidity. Understanding the photoinduced dynamics of these molecules is crucial for establishing general rules governing the relaxation and product formation. Here, in this study, we investigate the femtosecond time-resolved photoinduced excited-state structural dynamics of cis,cis-1,3-cyclooctadiene, a large-flexible cyclic conjugated diene molecule, upon excitation with 200 nm using mega-electron-volt ultrafast electron diffraction and trajectory surface hopping dynamics simulations. We tracked the photoinduced structural changes from the Franck–Condon region through the conical intersection seam to the ground state. Our findings revealed a novel primary reaction coordinate involving ring distortion, where the ring stretches along one axis and compresses along the perpendicular axis. The nuclear wavepacket remains compact along this reaction coordinate until it reaches the conical intersection seam, and it rapidly spreads as it approaches the ground state, where multiple products are formed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser Optogalvanic Spectroscopy pf Neon and Argon in a Discharge Plasma and its Significance for Microgravity Combustion

A detailed study of combustion mechanisms in flames, employing laser-based diagnostics, has provided good knowledge and understanding of the physical phenomena, and led to better characterization of the dynamical and chemical combustion processes, both under low-gravity (in space) and normal gravity (in ground based facilities, e.g. drop towers). Laser induced fluorescence (LIF), laser-induced incandescence (LII) and LIF thermometry have been widely used to perform nonintrusive measurements and to better understand combustion phenomena. Laser optogalvanic (LOG) spectroscopy has well-established applications in ion mobility measurements, atomic and molecular spectroscopy, ionization rates, recombination rates, velocity measurements and as a combustion probe for trace element detection. Absorption spectra of atomic and molecular species in flames can be obtained via LOG spectroscopy by measuring the voltage and current changes induced by laser irradiation. There are different kinds of processes that contribute to a discharge current, namely: (1) electron impact ionization, (2) collisions among the excited atoms of the discharge species and (3) Penning ionization. In general, at higher discharge currents, the mechanism of electron impact ionization dominates over Penning ionization, whereby the latter is hardly noticeable. In a plasma, whenever the wavelength of a laser coincides with the absorption of an atomic or molecular species, the rate of ionization of the species momentarily increases or decreases due to laser-assisted acceleration of collisional ionization. Such a rate of change in the ionization is monitored as a variation in the transient current by inserting a high voltage electrode into the plasma. Optogalvanic spectroscopy in discharges has been useful for characterizing laser line-widths and for providing convenient calibration lines for tunable dye lasers in the ultraviolet, visible and infrared wavelength regions. Different kinds of quantitative information, such as the electron collisional ionization rate, can be extracted from the complex processes occurring within the discharge. In the optogalvanic effect (OGE), there is no problem of overlap from background emissions, and hence even weak signals can be detected with a high signal-to-noise ratio, which makes the optogalvanic effect sensitive enough to resolve vibrational changes in molecular bonds and differences in energy levels brought about by different electron spins. For calibration purposes, neon and argon gaseous discharges have been employed most extensively, because these gases are commonly used as buffer gases within hollow-cathode lamps and provide an acceptable density of calibration lines. In the present work, our main aim has been to understand the dominant physical processes responsible for the production of the OGE signal, based on the extensive time resolved optogalvanic waveforms recorded, and also to extract quantitative information on the rates of excited state collisional processes.

Misra, Prabhakar↗

ALMA Spectral Line and Imaging Survey of a Low and a High Mass-Loss Rate AGB Star Between 335 and 362 GHz

Context. Low and intermediate mass stars are known to power strong stellar winds when evolving through the asymptotic giant branch (AGB) phase. Initial mass, luminosity, temperature, and composition determine the pulsation characteristics of the star and the dust species formed in the pulsating photospheric layers. Radiation pressure on these grains triggers the onset of a stellar wind. However, as of today, we still cannot predict the wind mass-loss rates and wind velocities from first principles neither do we know which species are the first to condense in the upper atmospheric regions. Aims. We aim to characterise the dominant physical, dynamical, and chemical processes in the inner wind region of two archetypical oxygen-rich (C/O < 1) AGB stars, that is, the low mass-loss rate AGB star R Dor (M’ ∼ 1 × 10(exp −7) M/yr) and the high mass-loss rate AGB star IK Tau (M’ ∼ 5 × 10(exp −6) M/yr). The purpose of this study is to observe the key molecular species contributing to the formation of dust grains and to cross-link the observed line brightnesses of several species to the global and local properties of the star and its wind. Methods. A spectral line and imaging survey of IK Tau and R Dor was made with ALMA between 335 and 362 GHz (band 7) at a spatial resolution of ∼150 mas, which corresponds to the locus of the main dust formation region of both targets. Results. Some two hundred spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states (up to v = 5 for SiO). Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v∞, previously derived for both objects (36 km/s versus v(∞) ∼ 17.7 km/s for IK Tau and 23 km/s versus v(∞) ∼ 5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ∼400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ∼3.7 × 10(exp −7) Mand ∼2 × 10(exp −8) M for IK Tau and R Dor, respectively. Conclusions. The observations presented here provide important constraints on the properties of these two oxygen-dominated AGB stellar winds. In particular, the ALMA data prove that both the dynamical and chemical properties are vastly different for this high mass-loss rate (IK Tau) and low mass-loss rate (R Dor) star.

L. Decin↗