Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic structure factor”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Nanometer-Scale Correlations in Aqueous Salt Solutions

The intermediate-range (1–4 nm) correlation of cations, an- ions, and water in aqueous alkaline earth salt solutions is measured using synchrotron X-ray diffraction. We differentiate from the entire solution structure factor, S(Q), a separate region at low Q (<1 °A-1) containing local diffraction maxima (pre-peaks) that indicate nanometer-scale oscillatory behavior. These features are quantitatively reproduced with classical molecular dynamics simulations. At high concentrations, the pre-peaks emerge from correlations arising from existence of small quasi close-packed lattice- like structures comprised of cation hydration spheres. We also analyze the concentration dependence of the pre-peak and discuss the overall results in light of the rich literature dealing with intermediate-range correlations underlying universal phenomena in concentrated electrolytes. The authors were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated by Battelle for the U.S. Department of Energy under Contract No. DE-AC05-76RL01830. Computing resources were generously allocated by PNNL's Institutional Computing program. DOE Contract No. DEAC02- 06CH11357 supports operation of the Advanced Photon Source at Argonne National Laboratory.

Chemical physics, nucleation, Madelung potential↗
DOE OSTI2021

Paradigms of frustration in superionic solid electrolytes

Superionic solid electrolytes have widespread use in energy devices, but the fundamental motivations for fast ion conduction are often elusive. In this Perspective, we draw upon atomistic simulations of a wide range of superionic conductors to illustrate some ways frustration can lower diffusion cation barriers in solids. Based on our studies of halides, oxides, sulfides and hydroborates and a survey of published reports, we classify three types of frustration that create competition between different local atomic preferences, thereby flattening the diffusive energy landscape. These include chemical frustration, which derives from competing factors in the anion–cation interaction; structural frustration, which arises from lattice arrangements that induce site distortion or prevent cation ordering; and dynamical frustration, which is associated with temporary fluctuations in the energy landscape due to anion reorientation or cation reconfiguration. For each class of frustration, we provide detailed simulation analyses of various materials to show how ion mobility is facilitated, resulting in stabilizing factors that are both entropic and enthalpic in origin. We propose the use of these categories as a general construct for classifying frustration in superionic conductors and discuss implications for future development of suitable descriptors and improvement strategies.

25 ENERGY STORAGE↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomic-scale 3D structural dynamics and functional degradation of Pt alloy nanocatalysts during the oxygen reduction reaction

Pt-based electrocatalysts are the primary choice for fuel cells due to their superior oxygen reduction reaction (ORR) activity. To enhance ORR performance and durability, extensive studies have investigated transition metal alloying, doping, and shape control to optimize the three key governing factors for ORR: geometry, local chemistry, and strain of their surface and subsurface. However, systematic optimization remains incomplete, as it requires an atomic-scale understanding of these factors and their dynamics over potential cycling, as well as their relationship to ORR activity. Here, we implement neural network-assisted atomic electron tomography to measure the 3D atomic structural dynamics and their effects on the functional degradation of PtNi alloy catalysts. Our results reveal that PtNi catalysts undergo shape changes, surface alloying, and strain relaxation during cycling, which can be effectively mitigated by Ga doping. By combining geometry, local chemistry, and strain analysis, we calculated the changes in ORR activity over thousands of cycles and observed that Ga doping leads to higher initial activity and greater stability. These findings offer a pathway to understanding 3D atomic structural dynamics and their relation to ORR activity during cycling, paving the way for the systematic design of durable, high-efficiency nanocatalysts.

Jeong, Chaehwa↗
DOE OSTI2024

Metabolic interactions shape emergent biofilm structures in a conceptual model of gut mucosal bacterial communities

Abstract The gut microbiome plays a major role in human health; however, little is known about the structural arrangement of microbes and factors governing their distribution. In this work, we present an in silico agent-based model (ABM) to conceptually simulate the dynamics of gut mucosal bacterial communities. We explored how various types of metabolic interactions, including competition, neutralism, commensalism, and mutualism, affect community structure, through nutrient consumption and metabolite exchange. Results showed that, across scenarios with different initial species abundances, cross-feeding promotes species coexistence. Morphologically, competition and neutralism resulted in segregation, while mutualism and commensalism fostered high intermixing. In addition, cooperative relations resulted in community properties with little sensitivity to the selective uptake of metabolites produced by the host. Moreover, metabolic interactions strongly influenced colonization success following the invasion of newcomer species. These results provide important insights into the utility of ABM in deciphering complex microbiome patterns.

Biotechnology & Applied Microbiology↗
DOE OSTI2026

Delving into the depths of NGC 3783 with XRISM V. Broad-band modelling of ionized outflows

The Seyfert 1 galaxy NGC 3783 hosts a multi-phase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000−2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We performed joint spectral modelling of the XMM-Newton /RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derived the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigated their thermal stability and absorption measure distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compared these results with the 2000−2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log ξ = 1.08 − 3.38 and outflow velocities of 480−1230 km s −1 . The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000−2001, but the earlier data contained more low-ionization, high-velocity components. The total column density of all the outflows in 2024 is 1.5 times larger than in 2000−2001, which means that it has been replenished by fresh material. The dominant unresolved transition array (UTA) absorber (component B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Quantum mereology: Factorizing Hilbert space into subsystems with quasiclassical dynamics

We study the question of how to decompose Hilbert space into a preferred tensor-product factorization without any preexisting structure other than a Hamiltonian operator, in particular the case of a bipartite decomposition into “system” and “environment.” Such a decomposition can be defined by looking for subsystems that exhibit quasiclassical behavior. The correct decomposition is one in which pointer states of the system are relatively robust against environmental monitoring (their entanglement with the environment does not continually and dramatically increase) and remain localized around approximately classical trajectories. We present an in-principle algorithm for finding such a decomposition by minimizing a combination of entanglement growth and internal spreading of the system. Both of these properties are related to locality in different ways. Furthermore, this formalism is relevant to questions in the foundations of quantum mechanics and the emergence of spacetime from quantum entanglement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Targeted modulation of pennycress lipid droplet proteins impacts droplet morphology and seed oil content

SUMMARY Lipid droplets (LDs) are unusual organelles that have a phospholipid monolayer surface and a hydrophobic matrix. In oilseeds, this matrix is nearly always composed of triacylglycerols (TGs) for efficient storage of carbon and energy. Various proteins play a role in their assembly, stability and turnover, and even though the major structural oleosin proteins in seed LDs have been known for decades, the factors influencing LD formation and dynamics are still being uncovered mostly in the “model oilseed” Arabidopsis. Here we identified several key LD biogenesis proteins in the seeds of pennycress, a potential biofuel crop, that were correlated previously with seed oil content and characterized here for their participation in LD formation in transient expression assays and stable transgenics. One pennycress protein, the lipid droplet associated protein‐interacting protein (LDIP), was able to functionally complement the Arabidopsis ldip mutant, emphasizing the close conservation of lipid storage among these two Brassicas. Moreover, loss‐of‐function ldip mutants in pennycress exhibited increased seed oil content without compromising plant growth, raising the possibility that LDIP or other LD biogenesis factors may be suitable targets for improving yields in oilseed crops more broadly.

Guzha, Athanas↗
DOE OSTI2023

A Moist Potential Vorticity Model for Midlatitude Long-Lived Mesoscale Convective Systems over Land

Mesoscale convective systems (MCSs) bring large amounts of rainfall and strong wind gusts to the midlatitude land regions, with significant impacts on local weather and hydrologic cycle. However, weather and climate models face a huge challenge in accurately modeling the MCS life cycle and the associated precipitation, highlighting an urgent need for a better understanding of the underlying mechanisms of MCS initiation and propagation. From a theoretical perspective, a suitable model to capture the realistic properties of MCSs and isolate the bare-bones mechanisms for their initiation, intensification, and eastward propagation is still lacking. To simulate midlatitude MCSs over land, we develop a simple moist potential vorticity (PV) model that readily describes the interactions among PV perturbations, air moisture, and soil moisture. Multiple experiments with or without various environmental factors and external forcing are used to investigate their impacts on MCS dynamics and mesoscale circulation vertical structures. The result shows that mechanical forcing can induce lower-level updraft and cooling, providing favorable conditions for MCS initiation. A positive feedback among surface winds, evaporation rate, and air moisture similar to the wind-induced surface heat exchange over tropical ocean is found to support MCS intensification. Both background surface westerlies and vertical westerly wind shear are shown to provide favorable conditions for the eastward propagation of MCSs. Last, our result highlights the crucial role of stratiform heating in shaping mesoscale circulation response. In conclusion, the model should serve as a useful tool for understanding the fundamental mechanisms of MCS dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

SISGR: The regulation of carbon fixation in plant and green algae: Rubisco activase and the origin of heat inactivation of CO 2 assimilation

Rubisco activase (Rca) is a critical AAA+ ATPase protein complex that remodels and promotes the Rubisco enzyme, a key player in photosynthetic performance and carbon fixation. The assembly and function of the Rca protein complex are regulated by a range of factors, including subunit concentration, nucleotide-binding states, thermal conditions, metal-ion coordination, and post-translational modifications, such as phosphorylation. Despite its importance in photosynthesis, the detailed molecular mechanisms underlying the regulation of plant Rca and how it activates Rubisco remain elusive. This project aims to bridge this knowledge gap by integrating sophisticated enzymology tools with single-molecule methods and high-resolution electron microscopy to elucidate the structure and function of plant Rca. Through these multiple approaches, we have systematically investigated how the activity of plant Rca is impacted by various factors, such as phosphorylation and metal-ion coordination. The Rca complex assembly/disassembly dynamics were captured using anti-Brownian electrokinetic (ABEL) trap-based measurements, providing unprecedented insight into its structural flexibility and diverse assembly states. Furthermore, the structural analysis of Rca through electron crystallography and single-particle cryogenic electron microscopy (cryo-EM) reveals novel assembly states of the spinach Rca, providing insight into the mechanistic action for Rubisco remodeling. By combining cutting-edge tools and approaches, this work uncovers critical aspects of Rca’s regulation and assembly, paving the way for a deeper understanding of its role in photosynthetic efficiency and the potential for enhancing carbon fixation in crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Maximizing ion dynamics and electrochemical performance of ionic liquid-acetonitrile electrolyte in Ti 3 C 2 T x MXene

Modification of the structure and morphology of MXene electrodes and the formulation of the electrolytes used in their supercapacitor configurations are significant factors affecting the performance of electrochemical devices. In this study, we investigated the electrochemical performance and ion dynamics of 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [EmimTFSI], ionic liquid in the presence of acetonitrile (ACN) at different concentrations in Ti 3 C 2 T x MXene supercapacitor. We found an optimum concentration of ACN, at which more cations from the ionic liquid attach to the MXene electrode surface, providing higher electrochemical performance. This higher capacitance is also associated with increased microscopic dynamics of the cation away from the pore wall. These findings give a guideline to optimize the performance of MXene-based supercapacitors using organic solvents-ionic liquid-based electrolyte systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structure of disordered TiO 2 phases from ab initio based deep neural network simulations

Amorphous TiO 2 (a-TiO 2 ) is widely used in many fields, ranging from photo-electrochemistry to bio-engineering, hence detailed knowledge of its atomic structure is of scientific and technological interest. Here we use an ab initio-based deep neural network potential (DP) to simulate large scale atomic models of crystalline and disordered TiO 2 with molecular dynamics. Our DP reproduces the structural properties of all (11) TiO 2 crystalline phases, predicts the densities and structure factors of molten and amorphous TiO 2 with only a few percent deviation from experiments, and describes the pressure dependence of the amorphous structure in agreement with recent observations. Furthermore, it can be extended to model additional structures and compositions and can be thus of great value in the study of TiO 2 -based (nano-)materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dissolving salt is not equivalent to applying a pressure on water

Abstract Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

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

A MEMS dynamic mechanical analyzer for in situ viscoelastic characterization of 3D printed nanostructures

Cellular metamaterial structures with sub-micron features have shown the ability to become excellent energy absorbing materials for impact mitigation due to the enhanced mechanical properties of materials at the nanoscale. However, in order to optimize the design of these energy absorbing metamaterial structures we need to be able to measure the dynamic properties of the sub-micron features such as storage and loss moduli and the loss factor. Therefore, at scale testing is required to capture the scale, temperature, and strain rate dependent material properties of these nanoscale materials. This study presents the design, fabrication, and calibration of a MEMS-based dynamic mechanical analyzer (DMA) that can be directly integrated with the two photon lithography (TPL) process commonly used to fabricate metamaterial structures with nanoscale features. The MEMS-based DMA consists of a chevron style thermal actuator used to generate a tensile load on the structure and two differential capacitive sensors on each side of the structure used to measure load and displacement. This design demonstrated 1.5 ± 0.75 nm displacement resolution and 104 ± 52 nN load resolution, respectively. Dynamic mechanical analysis was successfully conducted on a single nanowire feature printed between the load and displacement stages of the MEMS device with testing frequencies ranging between 0.01–10 Hz and testing temperatures ranging between 22 °C–47 °C. These initial tests on an exemplar TPL part demonstrate that the printed nanowire behaves as a viscoelastic material wherein the transition from glassy to viscous behavior has already set in at the room temperature.

42 ENGINEERING↗
DOE OSTI2021

Collective Nanoparticle Dynamics Associated with Bridging Network Formation in Model Polymer Nanocomposites

The addition of nanoparticles (NPs) to polymers is a powerful method to improve the mechanical and other properties of macromolecular materials. Such hybrid polymer–particle systems are also rich in fundamental soft matter physics. Among several factors contributing to mechanical reinforcement, a polymer-mediated NP network is considered to be the most important in polymer nanocomposites (PNCs). Here, we present an integrated experimental–theoretical study of the collective NP dynamics in model PNCs using X-ray photon correlation spectroscopy and microscopic statistical mechanics theory. Silica NPs dispersed in unentangled or entangled poly(2-vinylpyridine) matrices over a range of NP loadings are used. Static collective structure factors of the NP subsystems at temperatures above the bulk glass transition temperature reveal the formation of a network-like microstructure via polymer-mediated bridges at high NP loadings above the percolation threshold. The NP collective relaxation times are up to 3 orders of magnitude longer than the self-diffusion limit of isolated NPs and display a rich dependence with observation wavevector and NP loading. A mode-coupling theory dynamical analysis that incorporates the static polymer-mediated bridging structure and collective motions of NPs is performed. It captures well both the observed scattering wavevector and NP loading dependences of the collective NP dynamics in the unentangled polymer matrix, with modest quantitative deviations emerging for the entangled PNC samples. Here, we identify an unusual and weak temperature dependence of collective NP dynamics, in qualitative contrast with the mechanical response. Hence, the present study has revealed key aspects of the collective motions of NPs connected by polymer bridges in contact with a viscous adsorbing polymer medium and identifies some outstanding remaining challenges for the theoretical understanding of these complex soft materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

FLAMES─Fast, Low-Storage, Accurate, and Memory-Efficient Adaptive Sampling─Approach to Resolve Spatially Dependent Dynamics of Molecular Liquids

Many critical phenomena in soft matter occur at large length scales, necessitating the resolution of their structure and dynamics at low wavenumbers. However, resolving wavenumber-dependent dynamics computationally via molecular dynamics simulations presents significant challenges, as these phenomena span several orders of magnitude in both time and length scales, resulting in high computational costs and memory demands. Here, this work highlights the computational and memory challenges associated with analyzing molecular trajectories in reciprocal space and demonstrates a method to address them. We introduce FLAMESFast, Low-storage, Accurate, and Memory-Efficient adaptive Sampling, which is a direct method for calculation of structure factors, allowing us to select only the required number of wavevectors for binning. We also use wavenumber-dependent time steps to extract dynamics. Our FLAMES approach effectively mitigates computational and memory/storage bottlenecks. We demonstrate the method using simulations of a model system, liquid octane, at various temperatures. Comparisons with experimental data and real space computation show that the FLAMES technique achieves high accuracy in resolving temperature- and spatially dependent dynamics while being significantly more computationally efficient and requiring less memory and storage than methods based on a uniform wavevector grid and fixed temporal spacing.

Chen, Guang [Argonne National Laboratory (ANL), Ar↗
DOE OSTI2022

Effects of Salt Aggregation in Perfluoroether Electrolytes

Electrolytes comprised of polymers mixed with salts have great potential for enabling the use of Li metal anodes in batteries for increased safety. Ionic conductivity is one of the key performance metrics of these polymer electrolytes and achieving high room-temperature conductivity remains a challenge to date. For a bottom-up design of the polymer electrolytes, we must first understand how the structure of polyelectrolytes on a molecular level determines their properties. Here, we use classical molecular dynamics to study the solvation structure and ion diffusion in electrolytes composed of a short-chain perfluoroether with LiFSI or LiTFSI salts. Density functional theory is also used to provide some insights into the structures and energies of the salt interactions with the perfluoroether. We observe the formation of aggregates of salts in the fluorinated systems even at low salt concentrations. The fluorine-fluorine attraction in the solvent is the governing factor for creating the salt aggregates. The aggregates’ size and lifetime change with concentration and anion. These simulations provide an insight into the structure and dynamics of perfluoroether based electrolytes that can be used to improve Li-ion batteries.

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

Probing remdesivir nucleotide analogue insertion to SARS-CoV-2 RNA dependent RNA polymerase in viral replication

Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. Here in this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (~1.5 kcal mol -1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (~2.6 kcal mol -1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.

36 MATERIALS SCIENCE↗