Search NASA⌕ Search

SEARCH · Search NASA

Results for “Langevin diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

B-DeepONet: An enhanced Bayesian DeepONet for solving noisy parametric PDEs using accelerated replica exchange SGLD

Here, the Deep Operator Network (DeepONet) is a neural network architecture used to approximate operators, including the solution operator of parametric PDEs. DeepONets have shown remarkable approximation ability. However, the performance of DeepONets deteriorates when the training data is polluted with noise, a scenario that occurs in practice. To handle noisy data, we propose a Bayesian DeepONet based on replica exchange Langevin diffusion (reLD). Replica exchange uses two particles. The first particle trains a DeepONet to exploit the loss landscape and make predictions. The other particle trains a different DeepONet to explore the loss landscape and escape local minima via swapping. Compared to DeepONets trained with state-of-the-art gradient-based algorithms (e.g., Adam), the proposed Bayesian DeepONet greatly improves the training convergence for noisy scenarios and accurately estimates the uncertainty. To further reduce the high computational cost of the reLD training of DeepONets, we propose (1) an accelerated training framework that exploits the DeepONet's architecture to reduce its computational cost up to 25% without compromising performance and (2) a transfer learning strategy that accelerates training DeepONets for PDEs with different parameter values. Finally, we illustrate the effectiveness of the proposed Bayesian DeepONet using four parametric PDE problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Transport coefficients of the quark-gluon plasma at the critical point and across the first-order line

A bottom-up Einstein-Maxwell-dilaton holographic model is used to compute, for the first time, the behavior of several transport coefficients of the hot and baryon-rich strongly coupled quark-gluon plasma at the critical point and also across the first-order phase transition line in the phase diagram. The observables under study are the shear and bulk viscosities, the baryon and thermal conductivities, the baryon diffusion, the jet quenching parameter $\hat{q}$, as well as the heavy-quark drag force and the Langevin diffusion coefficients. These calculations provide a phenomenologically promising estimate for these coefficients, given that our model quantitatively reproduces lattice QCD thermodynamics results, both at zero and finite baryon density, besides naturally incorporating the nearly perfect fluidity of the quark-gluon plasma. Here we find that the diffusion of baryon charge, and also the shear and bulk viscosities, are suppressed with increasing baryon density, indicating that the medium becomes even closer to perfect fluidity at large densities. On the other hand, the jet quenching parameter and the heavy-quark momentum diffusion are enhanced with increasing density. The observables display a discontinuity gap when crossing the first-order phase transition line, while developing an infinite slope at the critical point. The transition temperatures associated with different transport coefficients differ in the crossover region but are found to converge at the critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Equation of State and Energy Loss of Hot and Dense Quark-Gluon matter from Holographic Black Holes

By using gravity/gauge correspondence, we construct a holographic model, constrained to mimic the lattice QCD equation of state at zero density, to investigate the temperature and baryon chemical potential dependence of the equation of state. We also obtained the energy loss of light and heavy partons within the hot and dense plasma represented by the heavy quark drag force, Langevin diffusion coefficients and jet quenching parameter at the critical point and across the first-order transition line predicted by the model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Topology induced crossover between Langevin, subdiffusion, and Brownian diffusion regimes in supercooled water

Despite extensive studies of supercooled water, it remains challenging to understand its peculiar dynamic anomalous properties. In this work, we integrated full atomistic simulations of supercooled water over the temperature range of room temperature to 200 K using a quantum-mechanics-based polarizable force field with the dressed dynamics method that couples fast collision events and slow reorganization dynamics of hydrogen-bond networks. Our analysis unveils the salient multiscale features in the transient relaxation dynamics of supercooled water. A classical Langevin behavior dominates at fast timescales, while long-time relaxations unveil two different activation barriers in two temperature regions: below and above 230 K. The modulation of the entropy spectrum by temperature is elucidated in terms of a three-state model underlined by the complexity of the water dynamics associated with a topological transition of a strong hydrogen-bond network. This state-dependent network topology is quantitatively characterized by power-law exponents of inverse network connectivity from 200 to 298 K. Finally, this work provides valuable guidance for further studies on the transient relaxation dynamics of supercooled water.

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

Coarse-Graining Conformational Dynamics with Multidimensional Generalized Langevin Equation: How, When, and Why

A data-driven ab initio generalized Langevin equation (AIGLE) approach is developed to learn and simulate high-dimensional, heterogeneous, coarse-grained (CG) conformational dynamics. Constrained by the fluctuation–dissipation theorem, the approach can build CG models in dynamical consistency (DC) with all-atom molecular dynamics. Here, we also propose practical criteria for AIGLE to enforce long-term DC. Case studies of a toy polymer, with 20 CG sites, and the alanine dipeptide, with two dihedral angles, elucidate why one should adopt AIGLE or its Markovian limit for modeling CG conformational dynamics in practice.

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

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2024

Excitonic Shockley-Read-Hall recombination in organic semiconductors

Trap-mediated recombination influences the performance of a wide range of electronic devices. The well-known Shockley-Read-Hall (SRH) expression for inorganic semiconductors is often invoked to describe the recombination rate in organic materials, although without a clear understanding of how its parameters relate to the underlying material properties or how it should be modified to account for the finite lifetime of exciton intermediates in, for example, the doped emissive layer of an organic light-emitting diode (OLED). Here, we formalize SRH recombination for organic semiconductors based on diffusive trapping and Langevin recombination. We show that including the exciton state suppresses the recombination rate in host-guest systems with type II energy level alignment whenever the interfacial gap between the host and guest molecular orbitals is comparable to the exciton energy. Furthermore, these results quantify the balance between bimolecular and trap-mediated recombination in doped OLED emissive layers, and indicate that devices with type II host-guest pairings can, in principle, beat the thermodynamic limit of their neat guest counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Generative Thermodynamic Computing

Here, we introduce a generative modeling framework for thermodynamic computing, in which structured data are synthesized from noise by the natural time evolution of a physical system governed by Langevin dynamics. While conventional diffusion models use neural networks to perform denoising, here the information needed to generate structure from noise is encoded by the dynamics of a thermodynamic system. Training proceeds by maximizing the probability with which the computer generates the reverse of a noising trajectory, which ensures that the computer generates data with minimal heat emission. We demonstrate this framework within a digital simulation of a thermodynamic computer. If realized in analog hardware, such a system would function as a generative model that produces structured samples without the need for artificially injected noise or active control of denoising.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2022

Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit

Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Diffusion-mediated passing of molecular species in linear nanopores constrained by orientational alignment

For diffusion-mediated catalytic conversion reactions in materials with narrow linear nanopores, e.g., mesoporous silica MCM-41, a key parameter is the propensity for product species to be able to pass reactant species and thus to efficiently exit the pore. For elongated species, this can require orientational alignment with the pore axis. In this work, we perform benchmark analyses for such solution-phase systems where one of these species is elongated in order to quantify the dependence of this passing propensity, P, on pore diameter and on the rotational diffusion coefficient, Dr, of the elongated species. In particular, we consider the passing of a spherical and an elongated spherocylindrical shaped species in a cylindrical pore in an implicit solvent, where these species cannot overlap. Passing is mediated by diffusive Brownian motion of these species as described by strongly damped Langevin dynamics. We quantify scaling of P for pore width just above the threshold where passing is sterically blocked, and also reveal a significant decrease in P for lower Dr. We also consider the dependence of P on the aspect ratio of the elongated species and obtain an exact result in the limiting regime of large aspect ratio.

Brownian motion↗
DOE OSTI2025

Simulation and generalized Langevin equation study of lipid subdiffusion in biomembrane phases

In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.

Malik, Sheeba [Oak Ridge National Laboratory (ORNL↗
NASA NTRS1983

Termolecular ionic recombination at high ambient gas density

It is shown that the importance of diffusion relative to drift in ion-ion recombination is a decreasing function of the separation between the ions (contrary to a common supposition). Langevin's formula for the recombination coefficient is obtained very simply using a physical model which allows for both diffusion and drift.

Bates, D. R.↗
DOE OSTI2024

Shot noise and universal Fano factor as a characterization of strongly correlated metals

Shot noise measures out-of-equilibrium current fluctuations and is a powerful tool to probe the nature of current-carrying excitations in quantum systems. Recent shot-noise measurements in the heavy-fermion strange metal YbRh 2 Si 2 exhibit a strong suppression of the Fano factor ( F )—the ratio of the current noise to the average current in the dc limit. This system is representative of metals in which electron correlations are extremely strong. Here we carry out the first theoretical study on the shot noise of diffusive metals in the regime of strong correlations. A Boltzmann-Langevin equation formulation is constructed in a quasiparticle description in the presence of strong correlations. We find that F = 3 / 4 in such a correlation regime. Thus, we establish the aforementioned Fano factor as universal to Fermi liquids, and we show that the Fano factor suppression observed in experiments on YbRh 2 Si 2 necessitates a loss of the quasiparticles. Our work opens the door to systematic theoretical studies of shot noise as a means of characterizing strongly correlated metallic phases and materials. Published by the American Physical Society 2024

Wang, Yiming (ORCID:0000000277582921)↗
DOE OSTI2026

Substrate-Mediated Evaporation and Stochastic Evolution of Supported Au Nanoparticles

Here, we use in situ transmission electron microscopy with automated tracking to study supported gold nanoparticles (NPs) during high-temperature vacuum annealing. The average mass loss per NP is governed by a flat, nearly size-independent substrate-mediated evaporation profile. On top of this mean shrinkage, individual NPs show significant fluctuations in apparent growth or shrinkage, and NP volume follows a random-walk-like trajectory. To rationalize both the ensemble-mean behavior and the particle-resolved variability, we develop a self-consistent theory that couples substrate-mediated evaporation to collective 2D Ostwald-type mass exchange through a shared adatom field, described in terms of a renormalized screening length and background concentration. In the experimentally relevant regime, the theory predicts an approximately size-independent mean shrinkage rate and clarifies how net mass loss suppresses classical coarsening. Superimposed on this deterministic drift, we quantify stochastic volume trajectories and capture their fluctuation spectrum with a minimal Langevin description consistent with intermittent adatom attachment and detachment events. In addition, we characterize the lateral diffusive motion of NPs, which is responsible for their coalescence. Altogether, our results highlight that stochasticity is intrinsic at the nanoscale and that predicting the evolution of supported NPs at early and intermediate times requires a unified framework combining substrate-mediated evaporation, collective mass exchange, and stochastic fluctuations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗