Search NASA⌕ Search

SEARCH · Search NASA

Results for “agreement dynamics”

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

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

At least 73 records · Page 4

Detailed study of the eikonal reaction theory for the breakup of one-neutron halo nuclei

Background: One-neutron removal reactions are used to study the single-particle structure of unstable nuclei, and in particular the exotic halo nuclei. The eikonal reaction theory (ERT) has been developed by Yahiro, Ogata, and Minomo [Prog. Theor. Phys. 126, 167 (2011)] to include dynamical effects, which are missing in the usual eikonal description of these reactions. Encouraging results have been obtained for total breakup cross sections in comparison to more elaborate reaction models. Purpose: Here, we extend these comparisons to more differential breakup cross sections expressed as functions of the relative energy or parallel momentum between the core and halo neutron. Method: ERT predictions of these cross sections are compared to state-of-the-art calculations. Results: The hypotheses upon which the ERT is based are confirmed and their range of validity is made clearer. The actual ordering of the evolution operators affects ERT differential cross sections and a specific choice leads to excellent agreement with the reference calculation. Dynamical effects in the treatment of the neutron-target interaction can be significant in the parallel-momentum observable. Conclusions: The role of the different interactions in the dynamics of breakup reactions of one-neutron halo nuclei are better understood and improvements to the ERT are suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nature of the nonequilibrium phase transition in the non-Markovian driven Dicke model

The Dicke model famously exhibits a phase transition to a superradiant phase with a macroscopic population of photons and is realized in multiple settings in open quantum systems. In this paper, we study a variant of the Dicke model where the cavity mode is lossy due to the coupling to a Markovian environment while the atomic mode is coupled to a colored bath. We analytically investigate this model by inspecting its low-frequency behavior via the Schwinger-Keldysh field theory and carefully examine the nature of the corresponding superradiant phase transition. Integrating out the fast modes, we can identify a simple effective theory allowing us to derive analytical expressions for various critical exponents including the dynamical exponent. We find excellent agreement with previous numerical results when the non-Markovian bath is at zero temperature; however, contrary to these studies, our low-frequency approach reveals that the same exponents govern the critical behavior when the colored bath is at finite temperature unless the chemical potential is zero. Furthermore, we show that the superradiant phase transition is classical in nature, while it is genuinely nonequilibrium. In this work, we derive a fractional Langevin equation and conjecture the associated fractional Fokker-Planck equation that captures the system's long-time memory as well as its nonequilibrium behavior. Finally, we consider finite-size effects at the phase transition and identify the finite-size scaling exponents, unlocking a rich behavior in both statics and dynamics of the photonic and atomic observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗

Extension of Argonne’s Plant Dynamics Code (PDC) Capabilities to TerraPower Direct-Cycle Supercritical CO2 (“Pascal”) Nuclear Reactor Concept. Final CRADA report

Argonne National Laboratory (the Contractor) and TerraPower, LLC, (the Participant), entered into a Cooperative Research and Development Agreement (CRADA) to extend the Plant Dynamics Code (PDC) to model the Pascal reactor plant. PDC is the system-level computer code developed at Argonne, primarily for design and transient analysis of supercritical carbon dioxide (sCO2) energy conversion systems for various power plants. Pascal is a supercritical carbon dioxide (sCO2)-cooled direct-cycle nuclear power reactor under development by TerraPower. New features were added to the code to model components and systems specific to the Pascal plant and to make the code readily available as a tool for design, control, and safety analysis for this and other direct-cycle power plants. TerraPower provided design information of the Pascal reactor and the requirements for the development of the new PDC features. Argonne modified the PDC capabilities to introduce the new features.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Atomic fluctuations in electronic materials revealed by dephasing

The microscopic origin and timescale of the fluctuations of the energies of electronic states has a significant impact on the properties of interest of electronic materials, with implication in fields ranging from photovoltaic devices to quantum information processing. Spectroscopic investigations of coherent dynamics provide a direct measurement of electronic fluctuations. Modern multidimensional spectroscopy techniques allow the mapping of coherent processes along multiple time or frequency axes and thus allow unprecedented discrimination between different sources of electronic dephasing. Exploiting modern abilities in coherence mapping in both amplitude and phase, we unravel dissipative processes of electronic coherences in the model system of CdSe quantum dots (QDs). The method allows the assignment of the nature of the observed coherence as vibrational or electronic. The expected coherence maps are obtained for the coherent longitudinal optical (LO) phonon, which serves as an internal standard and confirms the sensitivity of the technique. Fast dephasing is observed between the first two exciton states, despite their shared electron state and common environment. This result is contrary to predictions of the standard effective mass model for these materials, in which the exciton levels are strongly correlated through a common size dependence. In contrast, the experiment is in agreement with ab initio molecular dynamics of a single QD. Electronic dephasing in these materials is thus dominated by the realistic electronic structure arising from fluctuations at the atomic level rather than static size distribution. The analysis of electronic dephasing thereby uniquely enables the study of electronic fluctuations in complex materials.

36 MATERIALS SCIENCE↗

Extracting Inelastic Scattering Cross Sections for Finite and Aperiodic Materials from Electronic Dynamics Simulations

Explicit time-dependent electronic structure theory methods are increasingly prevalent in the areas of condensed matter physics and quantum chemistry, with the broad-band optical absorptivity of molecular and small condensed-phase systems nowadays routinely studied with such approaches. Here, in this paper, it is demonstrated that electronic dynamics simulations can similarly be employed to study cross sections for the scattering-induced electronic excitations probed in nonresonant inelastic X-ray scattering and momentum-resolved electron energy loss spectroscopies. A method is put forth for evaluating the electronic dynamic structure factor, which involves the application of a momentum boost-type perturbation and transformation of the resulting reciprocal space density fluctuations into the frequency domain. Good agreement is first demonstrated between the dynamic structure factor extracted from these electronic dynamics simulations and the corresponding transition matrix elements from linear response theory. The method is then applied to some extended (quasi)one-dimensional systems, for which the wave vector becomes a good quantum number in the thermodynamic limit. Finally, the dispersion of many-body excitations in a series of hydrogen-terminated graphene flakes (and twisted bilayers thereof) is investigated to highlight the utility of the presented approach for capturing morphology-dependent effects in the inelastic scattering cross sections of nanostructured and/or noncrystalline materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-temperature active oxidation of nanocrystalline silicon-carbide: A reactive force-field molecular dynamics study

Flexible woven SiC ceramics are prone to accelerated fiber embrittlement under high temperature oxidation in dynamic oxygen environments. The nanocrystalline structure of the constituent fibers impacts the reaction kinetics and phase transformations during active oxidation. However, fundamental understanding and quantification of grain boundary effects on oxidation behavior in nanocrystalline SiC remain elusive when temperatures exceed 1500 K. This study deploys large-scale molecular dynamics simulations with a reactive force-field to elucidate the complex roles of atomic oxygen reservoir conditions and grain size on oxidation kinetics and the nature of oxides produced in both monocrystalline and nanocrystalline 3C-SiC between 1100 K and 2000 K. The simulations with dynamically replenished oxygen provide good agreement with oxidation kinetics and activation energies for the monocrystalline Si(100) and C(100) orientations published in the available literature. This study reveals that, by contrast, nanocrystalline SiC samples exhibit two distinct oxidation kinetics with a transition point at 1500 K due to surface melting, which is supported by experimental evidence. The introduction of a grain-boundary network produces a two-fold decrease in oxidation activation energies compared to monocrystalline SiC below 1500 K. Above 1500 K, however, the activation energies rise substantially due to the formation of a liquid Si phase at the SiC/Si oxide interface. Further, it is shown that the stability of the interfacial liquid phase is promoted by incoherent grain boundaries in the crystalline SiC. These findings are important for the deployment of nanocrystalline SiC fibers in advanced thermal protection systems for high-temperature applications.

36 MATERIALS SCIENCE↗

QUIC-fire: A fast-running simulation tool for prescribed fire planning

Coupled fire-atmospheric modeling tools are increasingly used to understand the complex and dynamic behavior of wildland fires. Multiple research tools linking combustion to fluid flow use Navier-Stokes numerical solutions coupled to a thermodynamic model to understand fire-atmospheric feedbacks, but these computational fluid dynamics approaches require high-performance computing resources. We present a new simulation tool called QUIC-Fire to rapidly solve these feedbacks by coupling the mature 3-D rapid wind solver QUIC-URB to a physics-based cellular automata fire spread model Fire-CA. QUIC-Fire uses 3-D fuels inputs similar to the CFD model FIRETEC, allowing this tool to simulate effects of fuel structure on local winds and fire behavior. Results comparing fire behavior metrics to the computational fluid dynamic model FIRETEC show strong agreement. As a result, QUIC-Fire is the first tool intended to provide an opportunity for prescribed fire planners to compare, evaluate and design burn plans, including complex ignition patterns and coupled fire atmospheric feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Load-Displacement Relation and Gap Distribution Between Rough Surfaces: Partial Differential Equations Approach

We develop a theoretical model to predict the load-displacement relation and probability density function for gaps between contacting rough surfaces. We derive partial differential equations from the previous model (Joe et al., 2018), and extend them to the non-adhesive contact problem. The predictions of the present theory are compared with numerical results using the Green's function molecular dynamics algorithm, and a good agreement is obtained.

adhesion↗

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

Spinodal enhancement of fluctuations in nucleus-nucleus collisions

Subensemble Acceptance Method (SAM) [1, 2] is an essential link between measured event-by-event fluctuations and their grand canonical theoretical predictions such as lattice QCD. The method allows quantifying the global conservation law effects in fluctuations. In its basic formulation, SAM requires a sufficiently large system such as created in central nucleus-nucleus collisions and sufficient space-momentum correlations. Directly in the spinodal region of the First Order Phase Transition (FOPT) different approximations should be used that account for finite size effects. Thus, we present the generalization of SAM applicable in both the pure phases, metastable and unstable regions of the phase diagram [3]. Obtained analytic formulas indicate the enhancement of fluctuations due to crossing the spinodal region of FOPT and are tested using molecular dynamics simulations. A rather good agreement is observed. Using transport model calculations with interaction potential we show that the spinodal enhancement of fluctuations survives till the later stages of collision via the memory effect [4]. However, at low collision energies the space-momentum correlation is not strong enough for this signal to be transferred to second and third order cumulants measured in momentum subspace. This result agrees well with recent HADES data on proton number fluctuations at $\sqrt{S_{NN}}$ = 2.4 GeV which are found to be consistent with the binomial momentum space acceptance [5].

Poberezhnyuk, Roman↗

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗

Phase Separation, Capillarity, and Odd-Surface Flows in Chiral Active Matter

Active phase separations evade canonical thermodynamic descriptions and have thus challenged our understanding of coexistence and interfacial phenomena. Considerable progress has been made towards a nonequilibrium theoretical description of these traditionally thermodynamic concepts. Spatial parity symmetry is conspicuously assumed in much of this progress, despite the ubiquity of chirality in experimentally realized systems. Here, in this Letter, we derive a theory for the phase coexistence and interfacial fluctuations of a system that microscopically violates spatial parity. We find suppression of the phase separation as chirality is increased as well as the development of steady-state currents tangential to the interface dividing the phases. These odd flows are irrelevant to stationary interfacial properties, with stability, capillary fluctuations, and surface area minimization determined entirely by the capillary surface tension. Using large-scale Brownian dynamics simulations, we find excellent agreement with our theoretical scaling predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uranium Carbide and Uranium-Metal Thermal Scattering Law Evaluations and Cross Sections [Slides]

New Uranium Carbide and Uranium Metal evaluations have been developed using the FLASSH platform. Ab initio LD/MD models converged electronic, magnetic, and dynamical properties and demonstrate good agreement with experimental data. Thermal scattering cross sections generate inelastic and elastic thermal neutron scattering libraries in FLASSH . New TSL evaluations occurred in the production of ENDF File 7 data for Uranium Carbide and Uranium Metal submitted to NNDC for ENDF/B-VIII.1 release.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonadiabatically Driven Quantum Interference Effects in the Ultracold K + KRb → Rb + K 2 Chemical Reaction

The K + KRb → Rb + K 2 chemical reaction is the first ultracold atom–diatom chemical reaction for which experimental results have been reported for temperatures below 1 μK more than a decade ago. The reaction occurs through coupling with an excited electronic state that is accessible even in the ultracold limit. A previous quantum dynamics study, excluding nonadiabatic effects, has reported a rate coefficient that is about 35% below the experimental value. Here, we report the first nonadiabatic quantum dynamics study of this reaction and obtain rate coefficients in better agreement with experiments. Our results show that short-range dynamics mediated by coupling with the excited electronic state introduces quantum interference effects that influence both the state-to-state rate coefficients and the overall reaction rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From ionic clusters dynamics to network constraints in ionic polymer solutions

Physical networks formed by ionizable polymers with ionic clusters as crosslinks are controlled by coupled dynamics that transcend from ionic clusters through chain motion to macroscopic response. Here, the coupled dynamics, across length scales, from the ionic clusters to the networks in toluene swollen polystyrene sulfonate networks, were directly correlated, as the electrostatic environment of the physical crosslinks was altered. The multiscale insight is attained by coupling neutron spin echo measurements with molecular dynamics simulations, carried out to times typical of relaxation of polymers in solutions. The experimental dynamic structure factor is in outstanding agreement with the one calculated from computer simulations, as the networks are perturbed by elevating the temperature and changing the electrostatic environment. In toluene, the long-lived clusters remain stable over hundreds of ns across a broad temperature range, while the polymer network remains dynamic. In conclusion, though the size of the clusters changes as the dielectric constant of the solvent is modified through the addition of ethanol, they remain stable but morph, enhancing the polymer chain dynamics.

36 MATERIALS SCIENCE↗

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗

Quasi-equilibrium predictions of water desorption kinetics from rapidly-heated metal oxide surfaces

Controlling sub-microsecond desorption of water and other impurities from electrode surfaces at high heating rates is crucial for pulsed power applications. Despite the short time scales involved, quasi-equilibrium ideas based on transition state theory (TST) and Arrhenius temperature dependence have been widely applied to fit desorption activation free energies. Here, we apply molecular dynamics (MD) simulations in conjunction with equilibrium potential-of-mean-force (PMF) techniques to directly compute the activation free energies (ΔG*) associated with desorption of intact water molecules from Fe 2 O 3 and Cr 2 O 3 (0001) surfaces. The desorption free energy profiles are diffuse, without maxima, and have substantial dependences on temperature and surface water coverage. Incorporating the predicted ΔG* into an analytical form gives rate equations that are in reasonable agreement with non-equilibrium molecular dynamics desorption simulations. We also show that different ΔG* analytical functional forms which give similar predictions at a particular heating rate can yield desorption times that differ by up to a factor of four or more when the ramp rate is extrapolated by 8 orders of magnitude. This highlights the importance of constructing a physically-motivated ΔG* functional form to predict fast desorption kinetics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗