Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERACTION RANGE”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Experimental realization of an extended Fermi-Hubbard model using a 2D lattice of dopant-based quantum dots

The Hubbard model is an essential tool for understanding many-body physics in condensed matter systems. Artificial lattices of dopants in silicon are a promising method for the analog quantum simulation of extended Fermi-Hubbard Hamiltonians in the strong interaction regime. However, complex atom-based device fabrication requirements have meant emulating a tunable two-dimensional Fermi-Hubbard Hamiltonian in silicon has not been achieved. Here, we fabricate 3 × 3 arrays of single/few-dopant quantum dots with finite disorder and demonstrate tuning of the electron ensemble using gates and probe the many-body states using quantum transport measurements. By controlling the lattice constants, we tune the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from metallic to Mott insulating behavior. We simulate thermally activated hopping and Hubbard band formation using increased temperatures. As atomically precise fabrication continues to improve, these results enable a new class of engineered artificial lattices to simulate interactive fermionic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Neural-network quantum states for ultra-cold Fermi gases

Abstract Ultra-cold Fermi gases exhibit a rich array of quantum mechanical properties, including the transition from a fermionic superfluid Bardeen-Cooper-Schrieffer (BCS) state to a bosonic superfluid Bose-Einstein condensate (BEC). While these properties can be precisely probed experimentally, accurately describing them poses significant theoretical challenges due to strong pairing correlations and the non-perturbative nature of particle interactions. In this work, we introduce a Pfaffian-Jastrow neural-network quantum state featuring a message-passing architecture to efficiently capture pairing and backflow correlations. We benchmark our approach on existing Slater-Jastrow frameworks and state-of-the-art diffusion Monte Carlo methods, demonstrating a performance advantage and the scalability of our scheme. We show that transfer learning stabilizes the training process in the presence of strong, short-ranged interactions, and allows for an effective exploration of the BCS-BEC crossover region. Our findings highlight the potential of neural-network quantum states as a promising strategy for investigating ultra-cold Fermi gases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

NNL.Fe.qSNAP-ZBL.2024.1: A Fe Spectral Neighbor Analysis Potential for Radiation Damage Simulations

The NNL.Fe.qSNAP-ZBL.2024.1 machine-learned potential (MLP) has been generated to support the development of an elemental body-centered cubic (BCC) Fe athermal recombination corrected neutron damage model and simulations of primary recoil atom (PRA) cascades in BCC Fe. This MLP is a quadratic spectral neighbor analysis potential (qSNAP) hybridized with the universal Ziegler-Beirsack-Littmark (ZBL) potential at short-range and is named according to Naval Nuclear Laboratory MLP naming conventions (NNL.material-system.MLP-type.year.version). Training set calculations for Fe are presented along with the subsequent MLP fitting procedure. A key criterion of the fitting procedure is that ZBL describes the short-range interaction with minimal impact on the MLP. The MLP is compared to density functional theory (DFT) predicted properties relevant to radiation damage simulation, including threshold displacement energies, for validation. The NNL.Fe.qSNAP-ZBL.2024.1 potential is considered suitable for molecular dynamics (MD) simulations of radiation defects up to 800 K and PRA cascades in BCC Fe up to around 10 keV. The potential can additionally be used on a limited basis for recoils of 10–20 keV, within which range the emergence of structures outside the training set in cascade simulations may cause system instabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1991

A statistical approach to the brittle fracture of a multi-phase solid

A stochastic damage model is proposed to quantify the inherent statistical distribution of the fracture toughness of a brittle, multi-phase solid. The model, based on the macrocrack-microcrack interaction, incorporates uncertainties in locations and orientations of microcracks. Due to the high concentration of microcracks near the macro-tip, a higher order analysis based on traction boundary integral equations is formulated first for an arbitrary array of cracks. The effects of uncertainties in locations and orientations of microcracks at a macro-tip are analyzed quantitatively by using the boundary integral equations method in conjunction with the computer simulation of the random microcrack array. The short range interactions resulting from surrounding microcracks closet to the main crack tip are investigated. The effects of microcrack density parameter are also explored in the present study. The validity of the present model is demonstrated by comparing its statistical output with the Neville distribution function, which gives correct fits to sets of experimental data from multi-phase solids.

Liu, W. K.↗
DOE OSTI2021

Variational design principles for nonequilibrium colloidal assembly

Using large deviation theory and principles of stochastic optimal control, we show that rare molecular dynamics trajectories conditioned on assembling a specific target structure encode a set of interactions and external forces that lead to enhanced stability of that structure. Such a relationship can be formulated into a variational principle, for which we have developed an associated optimization algorithm and have used it to determine optimal forces for targeted self-assembly within nonequilibrium steady-states. We illustrate this perspective on inverse design in a model of colloidal cluster assembly within linear shear flow. We find that colloidal clusters can be assembled with high yield using specific short-range interactions of tunable complexity. Shear decreases the yields of rigid clusters, while small values of shear increase the yields of nonrigid clusters. The enhancement or suppression of the yield due to shear is rationalized with a generalized linear response theory. Furthermore, by studying 21 unique clusters made of six, seven, or eight particles, we uncover basic design principles for targeted assembly out of equilibrium.

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

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Many-body chaos in the antiferromagnetic quantum critical metal

Here, we compute the scrambling rate at the antiferromagnetic (AFM) quantum critical point, using the fixed point theory found by Schlief, Lunts, and Lee. At this strongly coupled fixed point, there is an emergent control parameter $w \ll 1$ that is a ratio of natural parameters of the theory. The strong coupling is unequally felt by the two degrees of freedom: the bosonic AFM collective mode is heavily dressed by interactions with the electrons, while the electron is only marginally renormalized. We find that the scrambling rates act as a measure of the “degree of integrability” of each sector of the theory: the Lyapunov exponent for the boson $λ^{(B)}_L \sim \mathcal{O} (\sqrt{w})k_BT/ℏ$ is significantly larger than the fermion one $λ^{(F)}_L \sim \mathcal{O}(w^2)k_BT/ℏ$, where $\textit{T}$ is the temperature. Although the interaction strength in the theory is of order unity, the larger Lyapunov exponent is still parametrically smaller than the universal upper bound of $λ_L = 2πk_BT/ℏ$. We also compute the spatial spread of chaos by the boson operator, whose low-energy propagator is highly nonlocal. We find that this nonlocality leads to a scrambled region that grows exponentially fast at intermediate distances, giving an infinite “butterfly velocity” of the chaos front, a result that has also been found in lattice models with long-range interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. Here, we also present results for the energy as a function of temperature in the continuum limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

INTERACTING COSMIC RAYS WITH MOLECULAR CLOUDS: A BREMSSTRAHLUNG ORIGIN OF DIFFUSE HIGH-ENERGY EMISSION FROM THE INNER 2 Degree-Sign Multiplication-Sign 1 Degree-Sign OF THE GALACTIC CENTER

The high-energy activity in the inner few degrees of the Galactic center is traced by diffuse radio, X-ray, and {gamma}-ray emission. The physical relationship between different components of diffuse gas emitting at multiple wavelengths is a focus of this work. We first present radio continuum observations using the Green Bank Telescope and model the nonthermal spectrum in terms of a broken power-law distribution of {approx}GeV electrons emitting synchrotron radiation. We show that the emission detected by Fermi is primarily due to nonthermal bremsstrahlung produced by the population of synchrotron emitting electrons in the GeV energy range interacting with neutral gas. The extrapolation of the electron population measured from radio data to low and high energies can also explain the origin of Fe I 6.4 keV line and diffuse TeV emission, as observed with Suzaku, XMM-Newton, Chandra, and the H.E.S.S. observatories. The inferred physical quantities from modeling multiwavelength emission in the context of bremsstrahlung emission from the inner {approx}300 Multiplication-Sign 120 pc of the Galactic center are constrained to have the cosmic-ray ionization rate {approx}1-10 Multiplication-Sign 10{sup -15} s{sup -1}, molecular gas heating rate elevating the gas temperature to 75-200 K, fractional ionization of molecular gas 10{sup -6}-10{sup -5}, large-scale magnetic field 10-20 {mu}G, the density of diffuse and dense molecular gas {approx}100 and {approx}10{sup 3} cm{sup -3} over 300 pc and 50 pc path lengths, and the variability of Fe I K{alpha} 6.4 keV line emission on yearly timescales. Important implications of our study are that GeV electrons emitting in radio can explain the GeV {gamma}-rays detected by Fermi and that the cosmic-ray irradiation model, like the model of the X-ray irradiation triggered by past activity of Sgr A*, can also explain the origin of the variable 6.4 keV emission from Galactic center molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Role of weakly bound complexes in temperature-dependence and relative rates of M x O y – + H 2 O ( M = Mo, W) reactions

Results of a systematic comparison of the Mo x O y – + H 2 O and W x O y – + H 2 O reaction rate coefficients are reported and compared to previous experimental and computational studies on these reactions. W x O y – clusters undergo more direct oxidation by water to yield W x O y+1 – + H 2 , while for Mo x O y – clusters, production of Mo x O y H 2 – (trapped intermediates in the oxidation reaction) is comparatively more prevalent. However, Mo x O y – clusters generally have higher rate coefficients than analogous W x O y – clusters if Mo x O y+1 H 2 – formation is included. Results of calculations on the M 2 O y – + H 2 O (M = Mo, W; y = 4, 5) reaction entrance channel are reported. They include charge-dipole complexes formed from long-range interactions, and the requisite conversion to a Lewis acid-base complex that leads to M x O y+1 H 2 – formation. The results predict that the Lewis acid-base complex is more strongly bound for Mo x O y – clusters than for W x O y – clusters. The calculated free energies along this portion of the reaction path are also consistent with the modest anti-Arrhenius temperature dependence measured for most Mo x O y – + H 2 O reactions, and the W x O y – + H 2 O reaction rate coefficients generally being constant over the temperature range sampled in this study. For clusters that exhibit evidence of both water addition and oxidation reactions, increasing the temperature increases the branching ratio toward oxidation for both species. In conclusion, a more direct reaction path to H 2 production may therefore become accessible at modest temperatures for certain cluster stoichiometries and structures.

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

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗
DOE OSTI2022

Chemical order-disorder nanodomains in Fe 3 Pt bulk alloy

Chemical ordering is a common phenomenon and highly correlated with the properties of solid materials. By means of the redistribution of atoms and chemical bonds, it invokes an effective lattice adjustment and tailors corresponding physical properties. To date, however, directly probing the 3D interfacial interactions of chemical ordering remains a big challenge. In this work, we deciphered the interlaced distribution of nanosized domains with chemical order/disorder in Fe 3 Pt bulk alloy. HAADF-STEM images evidence the existence of such nanodomains. The reverse Monte Carlo method with the X-ray pair distribution function data reveal the 3D distribution of local structures and the tensile effect in the disordered domains at the single-atomic level. The chemical bonding around the domain boundary changes the bonding feature in the disordered side and reduces the local magnetic moment of Fe atoms. This results in a suppressed negative thermal expansion and extended temperature range in Fe 3 Pt bulk alloy with nanodomains. Our study demonstrates a local revelation for the chemical order/disorder nanodomains in bulk alloy. The understanding gained from atomic short-range interactions within the domain boundaries provides useful insights with regard to designing new functional compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Tan relations in one dimension

We derive exact relations that connect the universal C/k{sup 4}-decay of the momentum distribution at large k with both thermodynamic properties and correlation functions of two-component Fermi gases in one dimension with contact interactions. The relations are analogous to those obtained by Tan in the three-dimensional case and are derived from an operator product expansion of the one- and two-particle density matrix. They extend earlier results by Olshanii and Dunjko (2003) for the bosonic Lieb-Liniger gas. As an application, we calculate the pair distribution function at short distances and the dimensionless contact in the limit of infinite repulsion. The ground state energy approaches a universal constant in this limit, a behavior that also holds in the three-dimensional case. In both one and three dimensions, a Stoner instability to a saturated ferromagnet for repulsive fermions with zero range interactions is ruled out at any finite coupling. - Highlights: > We derive universal relations for the two-component, contact-interacting 1D Fermi gas. > These relations connect the tail of the momentum distribution to thermodynamic properties. > There is no saturated ferromagnetism at finite, repulsive couplings for the 3D model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Realistic Ion Dynamics through Charge Renormalization in Nonaqueous Electrolytes

While many practically important electrolytes contain lithium ions, interactions of these ions are particularly difficult to probe experimentally because of their small X-ray and neutron scattering cross sections and large neutron absorption cross sections. Molecular dynamics (MD) is a powerful tool for understanding the properties of nonaqueous electrolyte solutions from the atomic level, but the accuracy of this computational method crucially depends on the physics built into the classical force field. In this study, we demonstrate that several force fields for lithium bistriflimide (LiTFSI) in acetonitrile yield a solution structure that is consistent with the neutron scattering experiments, yet these models produce dramatically different ion dynamics in solution. Such glaring discrepancies indicate that inadequate representation of long-range interactions leads to excessive ionic association and ion-pair clustering. We show that reasonable agreement with the experimental observations can be achieved by renormalization of the ion charges using a “titration” method suggested herewith. This simple modification produces realistic concentration dependencies for ionic diffusion and conductivity in <2 M solutions, without loss in quality for simulation of the structure.

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

Small RNAs at the Interface of the Plant-Fungal Interactions

Fungi and plants interact in a myriad of ways. These interactions range from mutualism to parasitism, and yet, many plant species appear to use similar tools to deal with both. One recent observation is the role that small RNAs play in mediating the conversation between plant and fungus. Increasingly, many studies demonstrate these RNAs are pivotal modulators, reprogramming gene expression and cellular processes essential to the biogenesis of these inter-kingdom relationships.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis

Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parametrize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but do not always demonstrate, how cancellation of errors occurs among the remaining molecular interactions accounted for such as exchange repulsion, electrostatics, and polarization. In this work we present the first generation of the (many-body) MB-UCB force field that explicitly accounts for the decomposed molecular interactions commensurate with a variational energy decomposition analysis, including charge transfer, with force field design choices that reduce the computational expense of the MB-UCB potential while remaining accurate. We optimize parameters using only a single water molecule and water cluster data up through pentamers, with no fitting to condensed phase data, and we demonstrate that high accuracy is maintained when the force field is subsequently validated against conformational energies of larger water cluster data sets, radial distribution functions of the liquid phase, and the temperature dependence of thermodynamic and transport water properties. Finally, we conclude that MB-UCB is comparable in performance to MB-Pol but is less expensive and more transferable by eliminating the need to represent short-ranged interactions through large parameter fits to high order polynomials.

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

Role of a Multivalent Ion–Solvent Interaction on Restricted Mg 2+ Diffusion in Dimethoxyethane Electrolytes

The diffusion behavior of Mg 2+ in electrolytes is not as readily accessible as that from Li + or Na + utilizing PFG NMR, due to the low sensitivity, poor resolution, and rapid relaxation encountered when attempting 25 Mg NMR. In MgTFSI 2 /DME solutions, “bound” DME (coordinating to Mg 2+ ) and “free” DME (bulk) are distinguishable from 1 H NMR. With the exchange rates between them obtained from 2D 1 H EXSY NMR, we can extract the self-diffusivities of free DME and bound DME (which are equal to that of Mg 2+ ) before the exchange occurs using PFG diffusion NMR measurements coupled with analytical formulas describing diffusion under two-site exchange. Further, the high activation enthalpy for exhange (65–70 kJ/mol) can be explained by the structural change of bound DME as evidenced by its reduced C–H bond length. Comparison of the diffusion behaviors of Mg 2+ , TFSI – , DME, and Li + reveals a relative restriction to Mg 2+ diffusion that is caused by the long-range interaction between Mg 2+ and solvent molecules, especially those with suppressed motions at high concentrations and low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗