Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean square displacement”

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

Dynamics of bottlebrush polymers

Bottlebrushes are an interesting class of polymers which shows intriguing material properties often associated with dynamics. While dynamical phenomena in linear polymers are well understood and existing theories can describe them in a good way, bottlebrush dynamics have only rarely been investigated. Therefore, we performed dielectric spectroscopy and quasi-elastic neutron scattering to study the dynamics of polydimethylsiloxane-based bottlebrush polymers, PDMS-g-PDMS focusing mostly on the segmental dynamics of the side chains. Comparing the relaxation times of the α – relaxation, tracked with dielectric spectroscopy, of bottlebrush polymers with those of their respective linear side chains show a slowing down once the side chains are attached to the backbone. This effect diminishes and finally vanishes with increasing side chain length. The time and length scale, offered by quasi-elastic neutron scattering, fits for the segmental dynamics together with faster processes. The Q -dependence of the segmental relaxation times allows to classify bottlebrush polymers as heterogenous including a non-Gaussian character. For such a dynamical system, the mean square displacement needs to be separated into single processes before an overall mean square displacement can be generated by applying the time temperature superposition principle.

36 MATERIALS SCIENCE↗

New techniques for diffusing-wave spectroscopy

We present two new types of measurements that can be made with diffusing-wave spectroscopy (DWS), a form of dynamic light scattering that applies in limit of strong multiple scattering. The first application is to measure the frequency-dependent linear viscoelastic moduli of complex fluids using light scattering. This is accomplished by measuring the mean square displacement of probe particles using DWS. Their response to thermal fluctuations is determined by the fluctuation-dissipation relation, and is controlled by the response of the surrounding complex fluid. This response can be described in terms of a memory function, which is directly related to the complex elastic modulus of the system. Thus by measuring the mean square displacement, we are able to determine the frequency dependent modulus. The second application is the measurement of shape fluctuations of scattering particles. This is achieved by generalizing the theory for DWS to incorporate the effects if amplitude fluctuations in the scattering intensity of the particles. We apply this new method to study the thermally induced fluctuations in the shape of spherical emulsion droplets whose geometry is controlled by surface tension.

Mason, T. G.↗

Atomistic Insights into the Reactive Diffusion of CO 2 in Guanidine-Based Facilitated Transport Membranes

The pressing need to address climate change has led to significant advancements in carbon dioxide (CO 2 ) capture technologies. Notably, facilitated transport membranes (FTMs) are distinguished by their exceptional selectivity and permeance, attributed to their reversible chemical reactions with CO 2 . This study, for the first time, sheds light on the reactive diffusion mechanism of CO 2 in FTMs, utilizing 1,1,3,3-tetramethylguanidine (TMG) as a mobile carrier. Specifically, state-ofthe- art molecular dynamics (MD) simulations, augmented by a reparameterized reactive force field (ReaxFF) capable of describing atomistic interactions and reaction pathways, are conducted to investigate the transport of CO 2 in TMG. The analysis of mean squared displacement (MSD) and diffusion coefficients reveals a clear hierarchy in the mobility of reaction components. Our findings highlight a unique hopping diffusion mechanism between bicarbonate ions and TMG molecules, increasing the diffusivity of reacted CO 2 by 1.4 times. The hopping events observed not only enhance our understanding of molecular mobility but also offer a means to boost the performance of FTMs in CO 2 capture applications. Overall, this research lays the groundwork for the future design of FTMs with optimal carrier properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Side Chain Dynamics of Poly(norbornene)- g -Poly(propylene oxide) Bottlebrush Polymers

The segmental dynamics of the side chains of poly(norbornene)- g -poly(propylene oxide) (PNB- g -PPO) bottlebrush polymer in comparison to PPO is studied by quasi-elastic neutron scattering. Having experimental time and length scale information simultaneously allows to extract spatial information in addition to relaxation time. Tethering one end of the PPO side chain onto a stiff PNB backbone slows down the segmental relaxation, over the length and time scales investigated. The power law dependence of the relaxation time on the momentum transfer, Q , indicates a more heterogeneous relaxation pattern for the bottlebrush polymer, whereas the linear PPO has less deviations from a homogenous relaxation. Similar conclusions can be drawn from the time dependent mean square displacement, $\langle$ r 2 ( t )$\rangle$, and the non-Gaussian parameter, α 2 ( t ). Herein, the bottlebrush polymer shows a more restricted dynamics, whereas the linear PPO reaches $\langle$ r 2 ( t )$\rangle$∝ t 0.5 at the highest temperature. The deviations from Gaussian behavior are evident at the α 2 ( t ). Both samples show a decaying α 2 ( t ). Finally, the non-Gaussian parameter supports the results from the power law dependence of the relaxation times, with lower α 2 ( t ) values for PPO compared to those for PNB-g-PPO, pointing to less deviations from Gaussian behavior.

36 MATERIALS SCIENCE↗

Atomically Precise Nanoclusters as Co‐Catalysts for Light‐Activated Microswimmer Motility

Microswimmers are self-propelled particles that navigate fluid environments, offering significant potential for applications in environmental pollutant decomposition, biosensing, and targeted drug delivery. Their performance relies on engineered catalytic surfaces. Gold nanoclusters (AuNCs), with atomically precise structures, tunable optical properties, and high surface area-to-volume ratio, provide a new optimal catalyst for enhancing microswimmer propulsion. Unlike bulk gold or nanoparticles, AuNCs may deliver tunable photocatalytic activity and increased catalytic specificity, making them ideal co-catalysts for hybrid microswimmers. For the first time, this study combines AuNCs with TiO 2 /Cr 2 O 3 Janus microswimmers, combining the unique properties of both materials. This hybrid system capitalizes on the tuned optical properties of AuNCs and their role as co-catalysts with TiO 2 , driving enhanced photocatalytic performance under ultraviolet (UV) excitation. Using motion analysis, it is shown that the AuNC-microswimmers exhibit significantly greater propulsion and mean squared displacement (MSD) as compared to controls. These findings suggest that the integration of nanoclusters with semiconductor materials enables state of the art, light-switchable microswimmers. These AuNC-microswimmer systems may thus offer new opportunities for environmental catalysis and other applications, providing precise control over catalytic and motile behaviors at the microscale.

36 MATERIALS SCIENCE↗

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗

Large negative thermal expansion in GdFe(CN) 6 driven by unusual low-frequency modes

Understanding the negative thermal expansion (NTE) mechanism is of great importance. In this work, we consider the new NTE compound GdFe(CN) 6 (αv = -34.2×10 -6 K -1 ) as a case study to investigate the NTE mechanism from the perspective of the lattice vibrational dynamics. The atomic mean-square displacements suggest that the NTE of GdFe(CN) 6 comes from the strong tension effect induced by the transverse vibrations of the atomic –Fe–Ctriple bondN–Gd– linkages, with the largest contribution given by N atoms. Lattice dynamics calculations show that three low-frequency optical modes at about 50 cm -1 show the largest negative Grüneisen parameters thus providing the largest contribution to the NTE. In conclusion, the existence of these unusual low-frequency vibrational modes can be ascribed to the presence of GdN 6 trigonal prisms in the framework structure of GdFe(CN) 6 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transferable predictions of energetic and structural properties for refractory solid solution alloys across chemical compositions

We present a data-efficient approach to train graph neural networks (GNNs) on density functional theory (DFT) data for accurate and transferable predictions of energetic and structural properties of refractory solid solution alloys in the niobium-tantalum-vanadium (Nb-Ta-V) chemical space. We start by training the GNN model only on DFT data that describes refractory binary alloys niobium-tantalum (Nb-Ta), niobium-vanadium (Nb-V), and tantalum-vanadium (Ta-V) to predict formation enthalpy and root mean squared displacement. Once trained, the GNN predictions are tested on DFT data describing refractory ternary alloys Nb-Ta-V. While, unsurprisingly, direct transferability from binary to ternary is not sufficiently accurate, augmenting the training with only 1% of the available ternary data (uniformly distributed across the entire range of chemical compositions) improves significantly the quality of the GNN predictions. For comparison, we assess the transferability in the opposite direction by training GNN models on ternary Nb-Ta-V data and making predictions on binaries Nb-Ta, Nb-V, and Ta-V, which exhibits notably higher predictive errors. The proposed methodology, which favors transferability from lower-component to higher-component alloys, offers an efficient path towards avoiding the curse of dimensionality incurred when collecting DFT data for discovery and design of multi-component disordered alloys.

Density functional theory calculations↗

Microdynamics of active particles in defect-rich colloidal crystals

Hypothesis. Because they are self-propulsive, active colloidal particles can interact with their environment in ways that differ from passive, Brownian particles. Here, we explore how interactions in different microstructural regions may contribute to colloidal crystal annealing. Experiments. We investigate active particles propagating in a quasi-2D colloidal crystal monolayer produced by alternating current electric fields (active-to-passive particle ratio ~ 1:720). The active particle is a platinum Janus sphere propelled by asymmetric decomposition of hydrogen peroxide. Crystals are characterized for changes in void properties. The mean-squared-displacement of Janus particles are measured to determine how active microdynamics depend on the local microstructure, which is comprised of void regions, void-adjacent regions (defined as within three particle diameters of a void), and interstitial regions. Findings. At active particle energy E A = 2.55 k B T, the average void size increases as much as three times and the average void anisotropy increases about 40% relative to the passive case. Further, the average microdynamical enhancement, <δ(t)>, of Janus particles in the crystal relative to an equivalent passive Janus particle is reduced compared to that of a free, active particle (<δ(t) > is 1.88 ± 0.04 and 2.66 ± 0.08, respectively). The concentration of active particles is enriched in void and void-adjacent regions. Active particles exhibit the greatest change in dynamics relative to the passive control in void-adjacent regions (<δ(t)> = 2.58 ± 0.06). The results support the conjecture that active particle microdynamical enhancement in crystal lattices is affected by local defect structure.

36 MATERIALS SCIENCE↗

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Molecular dynamics study of interstitial He clusters in nickel

This study presents a molecular dynamics analysis focusing on the behavior of interstitial helium (He) clusters in nickel (Ni), examining their formation, stability, and migration energetics. Consistent with previous research, we found that the binding energies of interstitial helium within a helium cluster are positive and increase with the cluster size, indicating a preference for helium atoms to cluster together. However, our findings also reveal that while the formation energy increases monotonically with cluster size, the increase in binding energy is non-monotonic. Importantly, small He clusters were observed to be thermally unstable at reactor operational temperatures (approximately 600 K), with the He 2 cluster exhibiting instability even at room temperature. With a binding energy of 0.49 eV for a He 4 cluster, we hypothesize that for helium bubbles to form via homogeneous nucleation (i.e., through trap mutation) at reactor operating temperatures, the helium concentration must be high enough to facilitate the formation of helium clusters of at least size 4 or larger. As expected, interstitial helium and small helium clusters are highly mobile. This mobility was observed not only at room temperature but also at temperatures as low as approximately 200 K. Furthermore, the mean squared displacement method has been utilized to determine the migration barriers and the corresponding prefactors for clusters ranging from He 1 to He 6

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Atomistic Simulations of Polydisperse Lignin Melts Using Simple Polydisperse Residue Input Generator

Understanding the physics of lignin will help rationalize its function in plant cell walls as well as aiding practical applications such as deriving biofuels and bioproducts. Here, in this work, we present SPRIG (Simple Polydisperse Residue Input Generator), a program for generating atomic-detail models of random polydisperse lignin copolymer melts i.e., the state most commonly found in nature. Using these models, we use all-atom molecular dynamics (MD) simulations to investigate the conformational and dynamic properties of polydisperse melts representative of switchgrass (Panicum virgatum L.) lignin. Polydispersity, branching and monolignol sequence are found to not affect the calculated glass transition temperature, T g . The Flory–Huggins scaling parameter for the segmental radius of gyration is 0.42 ± 0.02, indicating that the chains exhibit statistics that lie between a globular chain and an ideal Gaussian chain. Below T g the atomic mean squared displacements are independent of molecular weight. In contrast, above T g , they decrease with increasing molecular weight. Therefore, a monodisperse lignin melt is a good approximation to this polydisperse lignin when only static properties are probed, whereas the molecular weight distribution needs to be considered while analyzing lignin dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion of Atoms in Glassy Mixtures of Deuterium and Lithium

The diffusion coefficients of D and Li were calculated in amorphous, glassy Li+D mixtures for various concentrations of D in lithium over the temperature range 100–1000 K. The densities of the mixtures as a function of temperature were also determined. The diffusion coefficients were obtained by the analysis of the mean-squared displacement using molecular dynamics with ReaxFF, a reactive force field. In conclusion, the diffusion and density data obtained for Li+D mixtures were compared with the available experimental and calculated data reported in the literature.

Diffusion↗

Challenges in Pulsed-Field Gradient Nuclear Magnetic Resonance on Magnetically Heterogeneous Interfaces: Sequence and Field-Dependent Apparent Diffusion Coefficients

It is well known that the internal gradient (gi) that exists within pores haunts the diffusion coefficient (D) as measured by the pulsed-field gradient (PFG) nuclear magnetic resonance (NMR). Several PFG-NMR methods developed to determine an accurate D were not successful. Then, the steady-state diffusion coefficient (Dapp,8) for the cation [C4mim]+ of [C4mim][Tf2N]; [1-butyl-3-methylimidazolium][bis(trifluoromethylsulfonyl)imde] ionic liquid confined in ordered mesoporous carbon (OMC) were determined by comparing Dapp,8 obtained from 1H PFG-NMR performed with three different stimulated echo sequences: STE, APFG, and MPFG under the two external magnetic field strength, B0 = 9.4 and 14.1 Tesla. The measured Dapp,8 which is an order of magnitude smaller than D of bulk [C4mim][Tf2N], is in good agreement between APFG and MPFG both in B0 = 9.4 and 14.1 Tesla. However, the strong gi artifact, which caused apparent diffusion coefficient (Dapp) depending strongly and weakly on B0 and temperature, respectively, in diffusion-time dependent Dapp, Dapp(?) obtained from a sequence with monopolar gradients (STE) was suppressed by using sequences employing bipolar gradients (APFG and MPFG) in the region of steady-state diffusion. But incompletely suppressed gi artifact resulting in the different behaviors of the early part of Dapp(?) between the sequences leads a ˜ 0.6 and 0.9 in MPFG and APFG, respectively, in the relationship between mean squared displacement and diffusion time: = 2Dta, where a = 0.5 and 1 for 1-dimensional single file diffusion and 3-dimensional bulk diffusion, respectively. The above observations clearly show that the diffusion behavior of ions/molecules within the pores and pore structure, such as the surface-to-volume ratio? (D?_app (?)=D_0 [1-4/(9vp) S/V v(D_0 ?)]) and tortuosity (T = D0/Dapp,8), are possible to be misunderstood, especially in the systems with a non-negligible gi. This work demonstrates that it may be necessary to test several PFG sequences under multiple external magnetic fields for the correct determination of the diffusion behavior of ions/molecules in the pores with a larger internal gradient, gi.

Han, Kee Sung↗

Ab Initio Polariton Transport Dynamics with the Classical Path Approximation

We present an ab initio framework for simulating polariton transport dynamics based on the classical path approximation (CPA). The quantum dynamics of polariton transport involves simulating many electronic degrees of freedom, making a fully ab initio dynamics simulation computationally expensive. We demonstrate that the CPA, which removes the need for excited-state nuclear gradients, is well-suited for polaritonic systems because collective light–matter coupling leads to vanishing excited-state forces. Benchmark comparisons between CPA and full evaluation of the excited-state forces show excellent agreement for polariton transport results in model light–matter systems such as polariton group velocities and mean-squared displacements. Ab initio simulations of polariton transport using CPA reproduce key physical trends that are observed in experiments with BODIPY molecules. Our work establishes the CPA as a highly efficient tool for ab initio investigations of transport and energy flow in hybrid light–matter systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗