Search NASA⌕ Search

SEARCH · Search NASA

Results for “motion simulation”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗
DOE OSTI2025

The effect of African Easterly Wave suppression by periodicity on Atlantic tropical cyclones

Research has shown that suppressing African Easterly Waves (AEWs) does not reduce basin-wide North Atlantic tropical cyclone (TC) frequency but can enhance TC environmental favorability. We investigated the AEW-TC relationship further by examining the effects of suppressing the two AEW periodicities individually on TC activity. Using regional model simulations, AEWs were prescribed or suppressed in the 2–6 d or 6–10 d ranges through the lateral boundary conditions. Seasonal TC frequency increased significantly when either AEW periodicity was suppressed, with a larger increase when the 2–6 d waves were suppressed. We also found that suppressing the 2–6 d waves increased mid-tropospheric moisture by up to 8%, as well as overall atmospheric instability, near the western coast of northern Africa. Furthermore, the convective disturbances that developed into TCs exhibited stronger rotation, increased ascending motion, and higher rainfall. Our results suggest that reduced 2–6 d AEW activity may trigger a more active TC season.

Danso, Derrick K↗
DOE OSTI2024

Synergistic Integration of Multiple Wave Energy Converters with Adaptive Resonance and Offshore Floating Wind Turbines through Bayesian Optimization

We developed a synergistic ocean renewable system where an array of Wave Energy Converters (WEC) with adaptive resonance was collocated with a Floating Offshore Wind Turbine (FOWT) such that the WECs, capturing wave energy through the resonance adapting to varying irregular waves, consequently reduced FOWFT loads and turbine motions. Combining Surface-Riding WECs (SR-WEC) individually designed to feasibly relocate their natural frequency at the peak of the wave excitation spectrum for each sea state, and to obtain the highest capture width ratio at one of the frequent sea states for annual average power in a tens of kilowatts scale with a 15 MW FOWT based on a semi-submersible, Bayesian Optimization is implemented to determine the arrangement of WECs that minimize the annual representation of FOWT’s wave excitation spectra. The time-domain simulation of the system in the optimized arrangement is performed, including two sets of interactions: one set is the wind turbine dynamics, mooring lines, and floating body dynamics for FOWT, and the other set is the nonlinear power-take-off dynamics, linear mooring, and individual WECs’ floating body dynamics. Those two sets of interactions are further coupled through the hydrodynamics of diffraction and radiation. For sea states comprising Annual Energy Production, we investigate the capture width ratio of WECs, wave excitation on FOWT, and nacelle acceleration of the turbine compared to their single unit operations. We find that the optimally arranged SR-WECs reduce the wave excitation spectral area of FOWT by up to 60% and lower the turbine’s peak nacelle acceleration by nearly 44% in highly occurring sea states, while multiple WECs often produce more than the single operation, achieving adaptive resonance with a larger wave excitation spectra for those sea states. The synergistic system improves the total Annual Energy Production (AEP) by 1440 MWh, and we address which costs of Levelized Cost Of Energy (LCOE) can be reduced by the collocation.

Engineering↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2025

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

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

Beam Dynamics of the Muon $g\textrm{-}2$ Experiment

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Elucidating the mechanical and thermal response of nanotwinned Ni alloys

Transformative advances in nanoscale materials synthesis, characterization and modeling are enabling the synthesis of materials with unprecedent properties and greater understanding of the nanoscale mechanisms that underpin these properties. Nanotwinned Cu alloys have received considerable attention due to their impressive balance of strength and ductility, but they have limited microstructural stability. Recent instantiation of nanotwins in sputter deposited Ni-Mo-W, and several commercial Ni-based superalloys, point to the potential development of a new class of high temperature materials with a very beneficial suite of mechanical and physical properties. The experimental study described in this report was undertaken to elucidate the nanoscale origins of the thermal and mechanical behavior of nanotwinned Ni alloys. Micropillar compression experiments have shown nanotwinned Ni 85 Mo 15-x W x alloys to possess unusually high strengths above 3.5GPa and enhanced microstructural stability. The strength of these nanotwinned alloys is determined by the abrupt formation of shear bands. This study focused on identifying the nanoscale trigger, or triggers, for shear banding in nanotwinned materials using in situ experiments and atomic-scale postmortem analysis. Contrasting and comparing the mechanical response and postmortem nanostructures of “Mo-rich” and “W-rich” nanotwinned specimens elucidated the importance of the lateral motion of easy glide defects that results in erasure of twins and local coarsening of the nanotwinned microstructure. This twin coarsening was then associated with very localized plasticity, strain softening, and shear band formation. The difference between alloys was further studied by considering the role of various material factors: stacking fault energy, coherent twin boundary spacing and flatness, grain size, and the interactions of solute atoms with twin and grain boundaries. The effect of alloy composition and sputtering power and temperature were also investigated and used to control twin spacing and geometry. Combined with atomic-scale simulations, these experiments insights should allow us to identify strategies for achieving concomitant ultrahigh strength and ductility in nanotwinned Ni alloys. In parallel but synergistic studies, our work was buoyed by collaborative state-of-the-art nanoscale orientation and strain mapping at the DOE National Center for Electron Microscopy (NCEM). For example, novel 4D-STEM techniques were used to acquire nanoscale strain maps. The extremely fine twin spacing limited our measurements of atomic-scale thermal expansion within twins, but recent results from NCEM suggest that it will soon be possible to conduct such measurements and to unravel the origins of the novel thermal expansion characteristics of nanotwinned Ni alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mitigating Band Tailing in Kesterite Solar Absorbers: Ab Initio Quantum Dynamics

Open-circuit voltage deficits are limiting factors in kesterite solar cells. Addressing this issue by suppressing band tailing and nonradiative charge recombination is essential for enhancing the performance. We employ ab initio nonadiabatic molecular dynamics to elucidate the origin of band tailing and charge losses and propose a mitigation strategy. The simulations show that Cu–Zn disorder, associated with antisite defect clusters [Cu Zn +Zn Cu ], is a significant source of band tailing in kesterites, as evidenced by the much larger Urbach energy in disordered than ordered kesterites. Cu–Zn disorder gives rise to new sulfur-centered coordination polyhedra, increases structural inhomogeneity, changes electrostatic potential at sulfur centers, and shifts the S(3p) orbital energy. Differences in the S(3p)/Cu(3d) and S(3p)/Sn(5s) hybridization strengths and the S(3p) orbital energy shift reduce the band gap by 0.37 eV. Furthermore, Cu–Zn disorder enhances vibrational motion of sulfur anions and surrounding cations, increasing band gap fluctuations by 15 meV. The stronger electron–phonon interactions reduce charge carrier lifetimes and limit the kesterite solar cell efficiency. Partial substitution of Zn with Cd facilitates structural ordering and significantly suppresses band tailing, particularly in disordered systems. The improvement can be attributed to the larger atomic radius and mass of Cd, which weakens bonding around the anion, suppresses S-related vibrations within the covalent tetrahedra, and reduces nonadiabatic coupling, thereby increasing charge carrier lifetimes. The reported results establish the key influence of cation disorder on band tailing and reduced charge carrier lifetimes in kesterites and highlight cation disorder engineering as a strategy to achieve high-efficiency kesterite solar cells.

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

First Sagittarius A* Event Horizon Telescope Results. VIII. Physical Interpretation of the Polarized Ring

Abstract In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗
DOE OSTI2025

Minimal cyclic behavior in sheared amorphous solids

Although jammed packings of soft spheres exist in potential energy landscapes with a vast number of minima, when subjected to cyclic shear they may revisit the same configurations repeatedly. Simple hysteretic spin models, in which particle rearrangements are represented by interacting spin flips called hysterons, capture many features of this periodic behavior. Yet it has been unclear to what extent individual rearrangements can be described by such binary objects and how such objects interact with one another. Using a particularly sensitive algorithm, we identify rearrangements in simulated jammed packings and select pairs of rearrangements that undo one another to create periodic cyclic behavior. We find that the rearrangement pairs surprisingly persist down to the smallest increments in strain, even in the smallest systems we can study. We explore the statistics of these rearrangement pairs and find that there is a relation between the amount of hysteresis and the energy drop and mean-square displacement of the particles; these results are inconsistent with the scaling found in models that treat rearrangements as localized buckling events. Finally, our analysis shows that there is no clean distinction between the particle motions that represent the identity of a single, individual rearrangement and the particle motions that lead to interactions between separated rearrangements or hysterons. These results offer insight into how complex systems such as amorphous solids can reach a limit cycle.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Additive friction stir deposition induced stress ripples in aluminum alloy

In additive manufacturing, modification of the viscosity of softened or molten metal plays a significant role in the development of mechanical properties distribution. Similar to soft polymers, molten or softened metal usually exhibits vastly different viscoelastic properties. During scratching tests, viscoelastic polymers undergo spatiotemporal stress development and relaxation, causing stress ripples to form on their surfaces. Here, it is observed that friction stir-based solid-state manufacturing can yield similar stress ripples on softened aluminum alloys due to the linear motion of the deposition tool. Softened aluminum alloys exhibit viscoelastic and viscoplastic responses similar to polymeric material during the solid-state deposition process. In the final product, residual stress ripples appear on the surface of the metal alloy due to periodic residual stress relaxation, leading to oscillatory behavior from the bulk mechanical properties of the material. The orientation of the in-built stress pattern formed in the bulk of the additively manufactured metal alloy underneath the visually evident onion ring-like surface pattern formation is diametrically in the opposite direction. The experimentally observed results are consistent with analytical calculation and numerical simulations in terms of behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A molecular conveyor belt-associated protein controls the rotational direction of the bacterial type 9 secretion system

ABSTRACT Many bacteria utilize the type 9 secretion system (T9SS) for gliding motility, surface colonization, and pathogenesis. This dual-function motor supports both gliding motility and protein secretion, where rotation of the T9SS plays a central role. Fueled by the energy of the stored proton motive force and transmitted through the torque of membrane-anchored stator units, the rotary T9SS propels an adhesin-coated conveyor belt along the bacterial outer membrane like a molecular snowmobile, thereby enabling gliding motion. However, the mechanisms controlling the rotational direction and gliding motility of T9SS remain elusive. Shedding light on this mechanism, we find that in the gliding bacteriumFlavobacterium johnsoniae, deletion of the C-terminus of the conveyor belt-associated protein GldJ controls and, in fact, reverses the rotational direction of T9SS from counterclockwise (CCW) to clockwise (CW). This suggests that the interface between the conveyor belt-associated protein GldJ and the T9SS ring protein GldK plays an important role in controlling the directionality of T9SS, potentially by modulating its interaction with the stator complex GldLM, which drives motor rotation. Combined with MD simulation of the T9SS stator units GldLM, we suggest a “tri-component gearset” model where GldJ controls the rotational direction of its driver, the T9SS, thus providing adaptive sensory feedback to influence the motility of the gliding bacterium. IMPORTANCE The type 9 secretion system (T9SS) is fundamental to bacterial gliding motility, pathogenesis, and surface colonization. Our findings reveal that the C-terminal region of the conveyor belt-associated protein GldJ functions as a molecular switch which is capable of reversing the rotational direction of T9SS. Through the coordinated actions of the T9SS stator units (akin to a driving motor), the GldK ring (the gear that converts rotational energy into linear movement), and GldJ, this machinery forms a smart conveyor belt system reminiscent of flexible or cognitive mechanical conveyors. Such advanced conveyors can alter their direction to adapt to shifting demands. Here, we show that the bacterial T9SS similarly adjusts its rotational bias based on feedback from the conveyor belt-associated protein GldJ. This dual-role feedback mechanism underscores an evolved, controllable biological snowmobile, offering new avenues for studying how bacteria fine-tune motility in dynamic environments.

Microbiology↗
DOE OSTI2024

The reaction mechanism of the Ideonella sakaiensis PETase enzyme

Abstract Polyethylene terephthalate (PET), the most abundantly produced polyester plastic, can be depolymerized by the Ideonella sakaiensis PETase enzyme. Based on multiple PETase crystal structures, the reaction has been proposed to proceed via a two-step serine hydrolase mechanism mediated by a serine-histidine-aspartate catalytic triad. To elucidate the multi-step PETase catalytic mechanism, we use transition path sampling and likelihood maximization to identify optimal reaction coordinates for the PETase enzyme. We predict that deacylation is likely rate-limiting, and the reaction coordinates for both steps include elements describing nucleophilic attack, ester bond cleavage, and the “moving-histidine” mechanism. We find that the flexibility of Trp185 promotes the reaction, providing an explanation for decreased activity observed in mutations that restrict Trp185 motion. Overall, this study uses unbiased computational approaches to reveal the detailed reaction mechanism necessary for further engineering of an important class of enzymes for plastics bioconversion.

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

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗
DOE OSTI2024

Elucidating the Mechanism of Helium Evaporation from Liquid Water

Here, we investigate the evaporation of trace amounts of helium solvated in liquid water using molecular dynamics simulations and theory. Consistent with experimental observations, we find a super-Maxwellian distribution of kinetic energies of evaporated helium. This excess of kinetic energy over typical thermal expectations is explained by an effective continuum theory of evaporation based on a Fokker–Planck equation, parametrized molecularly by a potential of mean force and position-dependent friction. Using this description, we find that helium evaporation is strongly influenced by the friction near the interface, which is anomalously small near the Gibbs dividing surface due to the ability of the liquid–vapor interface to deform around the gas particle. Our reduced description provides a mechanistic interpretation of trace gas evaporation as the motion of an underdamped particle in a potential subject to a viscous environment that varies rapidly across the air–water interface. From it we predict the temperature dependence of the excess kinetic energy of evaporation, which is yet to be measured.

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

Ab Initio Polariton Spectra of ZnTPP Molecules Collectively Coupled Inside an Optical Cavity

Exciton-polaritons are quasi-particles formed by the quantum mechanical hybridization of electronic and photonic excitations. Despite extensive investigations, a fundamental understanding of molecular polariton spectra and the polariton delocalization from an ab initio theoretical perspective remains elusive. We simulate experimentally measured linear transmission spectroscopy of many Zinc(II) tetraphenylporphyrin (ZnTPP) molecules collectively coupled to a cavity from first principles. Our theoretical approach incorporates many low-lying electronic excitations in ZnTPP molecules, as well as collective light-matter couplings between ZnTPP and the quantized radiation modes, both of which are shown to be the key to accurately recovering the experimental spectra. We further analyzed to what extent the polariton and dark states are delocalized over many molecules, for the first time, using fully ab initio descriptions of the molecules. We finally investigate the line width as a function of detuning, providing new theoretical insights into the experimentally observed motional narrowing behavior. Our work presents first-ofits- kind theoretical studies on molecular polariton spectra, offering a new perspective on molecular polariton formation in realistic ab initio molecular systems whose rich, many-state nature provides spectral features enabled by the high density of electronic states beyond simple quantum optics models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗