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.

392 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Direct Observation of the ππ* to nπ* Transition in 2-Thiouracil via Time-Resolved NEXAFS Spectroscopy

The photophysics of nucleobases has been the subject of both theoretical and experimental studies over the past decades due to the challenges posed by resolving the steps of their radiationless relaxation dynamics, which cannot be described in the framework of the Born–Oppenheimer approximation (BOA). In this context, the ultrafast dynamics of 2-thiouracil has been investigated with a time-resolved NEXAFS study at the Free Electron Laser FLASH. Near Edge X-ray Absorption Fine Structure spectroscopy (NEXAFS) can be used to observe electronic transitions in ultrafast molecular relaxation. We performed time-resolved UV-pump/X-ray probe absorption measurements at the sulfur 2s (L1) and 2p (L2/3) edges. We are able to identify absorption features corresponding to the S2 (ππ*) and S1 (nπ*) electronic states. We observe a delay of 102 ± 11 fs in the population of the nπ* state with respect to the initial optical excitation and interpret the delay as the time scale for the S2 → S1 internal conversion. We furthermore identify oscillations in the absorption signal that match a similar observation in our previous X-ray photoelectron spectroscopy study on the same molecule.

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

Time-resolved momentum imaging of UV photodynamics in structural isomers of iodopropane probed by site-selective XUV ionization

The photodynamics of 1- and 2-iodopropane (1 and 2-IP) were studied in a time-resolved scheme incorporating ultraviolet (UV) excitation and extreme ultraviolet (XUV) probing, which initiates photoionization selectively from the I 4d core orbital. UV absorption in the A-band of both isomers leads to prompt C–I bond fission, with significant disposal of internal energy into the propyl radical product. Site-selective ionization enables a range of charge transfer (CT) processes between the nascent highly charged iodine ions and neutral propyl radicals, dependent on the interfragment distance at the instant of ionization. Subtle differences in the dynamics of these CT processes between the two isomers are observed. In 1-IP, the kinetic energies of iodine ions produced by UV photodissociation and subsequent XUV multiple ionization increased notably over the first few hundred femtoseconds, which could be understood in terms of differing gradients along the photodissociation coordinates of the neutral and polycationic states involved in the pump and probe steps, respectively. Led by a recent report of HI elimination in UV photoexcited 2-IP [Todt et al., Phys. Chem. Chem. Phys., 22(46), 27338 (2020)], we also model the most likely signatures of this process in the present experiment, and can identify signal in the 2-IP data (that is absent or significantly weaker in the data from the unbranched 1-IP isomer) that is consistent with such a process occurring on ultrafast timescales.

Allum, Felix [Deutsches Elektronen-Synchrotron (DE↗
DOE OSTI2024

Pseudogap in electron-doped cuprates: Strong correlation leading to band splitting

The pseudogap phenomena have been a long-standing mystery of the cuprate high-temperature superconductors. The pseudogap in the electron-doped cuprates has been attributed to band folding due to antiferromagnetic (AFM) long-range order or short-range correlation. We performed an angle-resolved photoemission spectroscopy study of the electron-doped cuprates Pr 1.3−x La 0.7 Ce x CuO 4 showing spin-glass, disordered AFM behaviors, and superconductivity at low temperatures and, by measurements with fine momentum cuts, found that the gap opens on the unfolded Fermi surface rather than the AFM Brillouin zone boundary. The gap did not show a node, following the full symmetry of the Brillouin zone, and its magnitude decreased from the zone-diagonal to (π,0) directions, opposite to the hole-doped case. These observations were reproduced by cluster dynamical-mean-field-theory calculation, which took into account electron correlation precisely within a (CuO 2 ) 4 cluster. In conclusion, the present experimental and theoretical results are consistent with the mechanism that electron or hole doping into a Mott insulator creates an in-gap band that is separated from the upper or lower Hubbard band by the pseudogap.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Mechanistic investigation of SARS-CoV-2 main protease to accelerate design of covalent inhibitors

Targeted covalent inhibition represents one possible strategy to block the function of SARS-CoV-2 Main Protease (M PRO ), an enzyme that plays a critical role in the replication of the novel SARS-CoV-2. Toward the design of covalent inhibitors, we built a covalent inhibitor dataset using deep learning models followed by high throughput virtual screening of these candidates against M PRO . Two top-ranking inhibitors were selected for mechanistic investigations—one with an activated ester warhead that has a piperazine core and the other with an acrylamide warhead. Specifically, we performed a detailed analysis of the free energetics of covalent inhibition by hybrid quantum mechanics/molecular mechanics simulations. Cleavage of a fragment of the non-structured protein (NSP) from the SARS-CoV-2 genome was also simulated for reference. Simulations show that both candidates form more stable enzyme-inhibitor (E-I) complexes than the chosen NSP. It was found that both the NSP fragment and the activated ester inhibitor react with CYS145 of M PRO in a concerted manner, whereas the acrylamide inhibitor follows a stepwise mechanism. Most importantly, the reversible reaction and the subsequent hydrolysis reaction from E-I complexes are less probable when compared to the reactions with an NSP fragment, showing promise for these candidates to be the base for efficient M PRO inhibitors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Local halide heterogeneity drives surface wrinkling in mixed-halide wide-bandgap perovskites

Compositional heterogeneity in wide-bandgap (1.8 - 2.1 eV) mixed-halide perovskites is a key bottleneck in the processing of high-quality solution-processed thin films and prevents their application in efficient multijunction solar cells. Notably, mixed-cation (formamidinium-methylammonium) wide-bandgap perovskite films are prone to form micrometer-scale wrinkles which can interfere with the smooth surfaces ideal for multijunction devices. Here, we study the formation dynamics of wrinkled mixed-halide perovskite films and its impact on the local composition and optoelectronic properties. We use in situ X-ray scattering during perovskite film formation to show that crystallization of bromide-rich perovskites precedes that of mixed-halide phases in wrinkled films cast using an antisolvent-based process. Using nanoscopic X­-ray fluorescence and hyperspectral photoluminescence imaging, we also demonstrate the formation of iodide- and bromide-rich phases in the wrinkled domains. This intrinsic spatial halide segregation results in an increased local bandgap variation and Urbach energy. Morphological disorder and compositional heterogeneity also aggravate the formation of sub-bandgap electronic defects, reducing photostability and accelerating light-induced segregation of iodide and bromide ions in thin films and solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2025

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

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

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2025

Visualizing the strong field–induced molecular breakup of C 60 via x-ray diffraction

Laser-driven dynamics in polyatomic molecules poses a complex many-body problem. Understanding intense light-matter interaction is crucial for steering intramolecular quantum dynamical processes. Here, we record time-resolved x-ray diffraction images of C 60 molecules during and after their interaction with intense near-infrared fields, giving direct access to structural changes of the molecules and their fragmentation in real time. Tuning the intensity of the excitation pulses, we uncover a transition from a weak-field regime of excited but stable molecules to a high-field regime dominated by Coulomb explosion. In the transition region, the molecules expand by up to 50% of their initial size within just 140 fs, with major fragmentation only setting in afterward. This work demonstrates that x-ray diffractive imaging is capable of retrieving time-resolved structural information of large molecules reshaped by intense laser fields. Laser-driven fragmentation is a first step toward observing molecular processes modified by laser fields of increasing intensity.

Schnorr, Kirsten [Paul Scherrer Inst. (PSI), Villi↗
NASA NTRS2011

NASA's Lunar Atmosphere and Dust Environment Explorer (LADEE)

Nearly 40 years have passed since the last Apollo missions investigated the mysteries of the lunar atmosphere and the question of levitated lunar dust. The most important questions remain: what is the composition, structure and variability of the tenuous lunar exosphere? What are its origins, transport mechanisms, and loss processes? Is lofted lunar dust the cause of the horizon glow observed by the Surveyor missions and Apollo astronauts? How does such levitated dust arise and move, what is its density, and what is its ultimate fate? The US National Academy of Sciences/National Research Council decadal surveys and the recent "Scientific Context for Exploration of the Moon" (SCEM) reports have identified studies of the pristine state of the lunar atmosphere and dust environment as among the leading priorities for future lunar science missions. These measurements have become particularly important since recent observations by the Lunar Crater Observation and Sensing Satellite (LCROSS) mission point to significant amounts of water and other volatiles sequestered within polar lunar cold traps. Moreover Chandrayaan/M3, EPOXI and Cassini/VIMS have identified molecular water and hydroxyl on lunar surface regolith grains. Variability in concentration suggests these species are likely to be present in the exosphere, and thus constitute a source for the cold traps. NASA s Lunar Atmosphere and Dust Environment Explorer (LADEE) is currently under development to address these goals. LADEE will determine the composition of the lunar atmosphere and investigate the processes that control its distribution and variability, including sources, sinks, and surface interactions. LADEE will also determine whether dust is present in the lunar exosphere, and reveal its sources and variability. LADEE s results are relevant to surface boundary exospheres and dust processes throughout the solar system, will address questions regarding the origin and evolution of lunar volatiles, and will have implications for future exploration activities. LADEE will be the first mission based on the Ames Common Bus design. LADEE employs a high heritage instrument payload: a Neutral Mass Spectrometer (NMS), an Ultraviolet/Visible Spectrometer (UVS), and the Lunar Dust Experiment (LDEX). It will also carry a space terminal as part of the Lunar Laser Communication Demonstration (LLCD), which is a technology demonstration. LLCD will also supply a ground terminal. LLCD is funded by the Space Operations Mission Directorate (SOMD), managed by GSFC, and built by MIT Lincoln Lab. NMS was directed to the Goddard Space Flight Center (GSFC) and UVS to Ames Research Center (ARC). LDEX was selected through the Stand Alone Missions of Opportunity Notice (SALMON) Acquisition Process, and is provided by the University of Colorado at Boulder. The LADEE NMS covers a m/z range of 2-150 and draws its design from mass spectrometers developed at GSFC for the MSL/SAM, Cassini Orbiter, CONTOUR, and MAVEN missions. The UVS instrument is a next-generation, high-reliability version of the LCROSS UV-Vis spectrometer, spanning 250-800 nm wavelength, with high (<1 nm) spectral resolution. UVS will also perform dust occultation measurements via a solar viewer optic. LDEX senses dust impacts in situ, at LADEE orbital altitudes of 50 km and below, with a particle size range of between 100 nm and 5 micron. Dust particle impacts on a large hemispherical target create electron and ion pairs. The latter are focused and accelerated in an electric field and detected at a microchannel plate. LADEE is an important part of NASA s portfolio of near-term lunar missions; launch is planned for May, 2013. The lunar atmosphere is the most accessible example of a surface boundary exosphere, and may reveal the sources and cycling of volatiles. Dynamic dust activity must be accounted for in the design and operation of lunar surface operations.

Elphic, Richard↗
DOE OSTI2025

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

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

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2026

Structural Mechanism of an Efficacy Photoswitch Targeting the β 2 ‐adrenergic Receptor

The field of photopharmacology develops light-responsive drugs that can modulate protein activity, enabling precise and dynamic investigations of their roles in health and disease. Adrenergic receptors are prominent targets for this approach because they are prototypical G protein-coupled receptors with high clinical relevance in bronchial and cardiovascular diseases. Here, we employed the azobenzene-based compound photoazolol-1 in combination with time-resolved serial crystallography at X-ray free-electron lasers to resolve the molecular mechanisms by which photoswitchable β-blockers modulate activity of the β 2 -adrenoceptor (β 2 AR). Time-resolved structures of the receptor bound to trans-photoazolol-1 (pre-photoconversion), a strained intermediate in the nanosecond range, and the fully photoisomerized cis-photoazolol-1 reveal how isomerization of the azobenzene moiety induces distinct conformational changes within the orthosteric ligand binding pocket. Within seconds, light-excited photoazolol-1 adopts a new binding pose, altering interactions with extracellular loop 2 and shifting the positions of transmembrane helices 5, 6, and 7. Functional assays of β 2 AR in cellular membranes show that photoazolol-1 acts as an efficacy photoswitch, changing from an inverse agonist to a neutral antagonist upon isomerization without leaving the binding pocket. In combination, these findings suggest a molecular mechanism for activity modulation via efficacy photoswitches and provide a framework for designing ligands that exploit light-driven transitions within the binding pocket to achieve spatiotemporal control of receptor function.

G protein-coupled receptors↗
NASA NTRS2002

Motion of Drops on Surfaces with Wettability Gradients

A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Subramanian, R. Shankar↗
DOE OSTI2025

Wavelength-dependent photodissociation of iodomethylbutane

Ultrashort XUV pulses of the Free-Electron-LASer in Hamburg (FLASH) were used to investigate laser-induced fragmentation patterns of the prototypical chiral molecule 1-iodo-2-methyl-butane (C 5 H 11 I) in a pump-probe scheme. Ion velocity-map images and mass spectra of optical-laser-induced fragmentation were obtained for subsequent FEL exposure with photon energies of 63 eV and 75 eV. These energies specifically address the iodine 4d edge of neutral and singly charged iodine, respectively. The presented ion spectra for two optical pump-laser wavelengths, i.e., 800 nm and 267 nm, reveal substantially different cationic fragment yields in dependence on the wavelength and intensity. For the case of 800-nm-initiated fragmentation, the molecule dissociates notably slower than for the 267 nm pump. The results underscore the importance of considering optical-laser wavelength and intensity in the dissociation dynamics of this prototypical chiral molecule that is a promising candidate for future studies of its asymmetric nature.

74 ATOMIC AND MOLECULAR PHYSICS↗