Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Force on Force”

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 451 records · Page 25

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Improving the precision of forces in real-space pseudopotential density functional theory

The high-order finite difference real-space pseudopotential density functional theory (DFT) approach is a valuable method for large-scale, massively parallel DFT calculations. A significant challenge in the approach is the oscillating “egg-box” error introduced by aliasing associated with a coarse grid spacing. To address this issue while minimizing computational cost, we developed a finite difference interpolation (FDI) scheme [Roller et al., J. Chem. Theory Comput. 19, 3889 (2023)] as a means of exploiting the high resolution of the pseudopotential to reduce egg-box effects systematically. Here, we show an implementation of this method in the PARSEC code and examine the practical utility of the combination of FDI with additional methods for improving force precision and/or reducing its computational cost, including orbital-based forces, compensating charges (namely, adding and subtracting a judiciously chosen charge density such that the total density is unaltered), and a modified spatial domain in which the real-space grid is defined. Using selected small molecules, as well as metallic Li, as test cases, we show that a combination of all four aspects leads to a significant reduction in computational cost while retaining a high level of precision that supports accurate structures and vibrational spectra, as well as stable and accurate molecular dynamics runs.

Chemistry↗

Exciting DeePMD: Learning excited-state energies, forces, and non-adiabatic couplings

We extend the DeePMD neural network architecture to predict electronic structure properties necessary to perform non-adiabatic dynamics simulations. While learning the excited state energies and forces follows a straightforward extension of the DeePMD approach for ground-state energies and forces, how to learn the map between the non-adiabatic coupling vectors (NACV) and the local chemical environment descriptors of DeePMD is less trivial. Most implementations of machine-learning-based non-adiabatic dynamics inherently approximate the NACVs, with an underlying assumption that the energy-difference-scaled NACVs are conservative fields. We overcome this approximation, implementing the method recently introduced by Richardson [J. Chem. Phys. 158, 011102 (2023)], which learns the symmetric dyad of the energy-difference-scaled NACV. Furthermore, the efficiency and accuracy of our neural network architecture are demonstrated through the example of the methaniminium cation CH 2 NH 2 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Astrophysical observations of a dark matter-Baryon fifth force

We consider the effects of an attractive, long-range Yukawa interaction between baryons and dark matter (DM), focusing in particular on temperature and pulsar timing observations of neutron stars (NSs). We show that such a fifth force, with strength modestly stronger than gravity at ranges greater than tens of kilometers (corresponding to mediator masses less than 10 -11 eV), can dramatically enhance dark matter kinetic heating, capture, and pulsar timing Doppler shifts relative to gravity plus short range interactions alone. Using the coldest observed NS and pulsar timing array (PTA) data, we derive limits on fifth force strength over a DM mass range spanning light dark matter up to order solar mass composite DM objects. We also consider an indirect limit by combining bullet cluster limits on the DM self-interaction with weak equivalence principle test limits on baryonic self-interactions. Furthermore, we find the combined indirect limits are moderately stronger than kinetic heating and PTA limits, except when considering a DM subcomponent.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomic force microscopy phase imaging of epitaxial graphene films

Abstract Dynamic mode atomic force microscopy phase imaging is known to produce distinct contrast between graphene areas of different atomic thickness. But the intrinsic complexity of the processes controlling the tip motion and the phase angle shift excludes its use as an independent technique for a quantitative type of analysis. By investigating the relationship between the phase shift, the tip-surface interaction, and the thickness of the epitaxial graphene areas grown on silicon carbide, we shed light on the origin of such phase contrast, and on the complex energy dissipation processes underlying phase imaging. In particular, we study the behavior of phase shift and energy dissipation when imaging the interfacial buffer layer, single-layer, and bilayer graphene regions as a function of the tip-surface separation and the interaction forces. Finally, we compare these results with those obtained on differently-grown quasi free standing single- and bilayer graphene samples.

Lavini, Francesco (ORCID:0000000257223040)↗

Decoherence from long-range forces in atom interferometry

Atom interferometers provide a powerful means of realizing quantum coherent systems with increasingly macroscopic extent in space and time. These systems provide an opportunity for a variety of novel tests of fundamental physics, including ultralight dark matter searches and tests of modifications of gravity, using long drop times and microgravity environments. However, as experiments operate with longer periods of free fall and become sensitive to smaller background effects, key questions start to emerge about the fundamental limits to future atom interferometry experiments. Here we study the effects on atomic coherence from hard-to-screen backgrounds due to baths of ambient particles with long-range forces, such as gravitating baths and charged cosmic rays. Our approach—working in the Heisenberg picture for the atomic motion—makes proper inclusion of the experimental apparatus feasible and clearly shows how to handle long-range forces. We find that these potential backgrounds are likely negligible for the next generation of interferometers, as aggressive estimates for the gravitational decoherence from a background bath of dark matter particles gives a decoherence timescale on the order of years.

74 ATOMIC AND MOLECULAR PHYSICS↗

Deep learning of accurate force field of ferroelectric HfO 2

The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo

The accurate computation of forces and other energy derivatives has been a long-standing challenge for quantum Monte Carlo methods. A number of technical obstacles contribute to this challenge. We discuss how these obstacles can be removed with the auxiliary-field quantum Monte Carlo (AFQMC) approach. AFQMC is a general, high-accuracy, many-body total-energy method for molecules and solids. The implementation of back-propagation for pure estimators allows direct calculation of gradients of the energy via the Hellmann-Feynman theorem. A planewave basis with norm-conserving pseudopotentials is used for the study of periodic bulk materials. Completeness of the planewave basis minimizes the effect of so-called Pulay terms. The ionic pseudopotentials, which can be incorporated in AFQMC in exactly the same manner as in standard independent-electron methods, regulate the force and stress estimators and eliminate any potential divergence of the Monte Carlo variances. The resulting approach allows applications of full geometry optimizations in bulk materials. As a result, it also paves the way for many-body computations of the phonon spectrum in solids.

36 MATERIALS SCIENCE↗

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Finite-size effects on the self-force

Electromagnetic and linear gravitational radiation do not solely propagate on the null cone in 3 + 1 dimensions in curved spacetimes, contrary to their well-known behavior in flat spacetime. Their additional propagation inside the null cone is known as the tail effect. A compact body will produce a signal whose tail will interact with its future worldline, thus producing a tail-induced self-force. In this work, we present new results for the tail-induced scalar, electromagnetic, and gravitational self-force for a test mass in orbit around a central mass, including effects from the internal structure of that body.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dark dimension and the grand unification of forces

The dark dimension scenario, predicting one extra mesoscopic dimension in the micron range, has emerged by applying various swampland principles to the dark energy. In this note we find that realizing the grand unification of gauge forces is highly constraining in this context. Without actually constructing any grand unified theory (GUT) models, we argue that the mere assumption of grand unification of forces in this scenario, together with the experimental bounds on massive replicas of the Standard Model gauge bosons, predicts an upper bound for the GUT scale, 𝑀 GUT ≲ 10 16 GeV. Combined with the experimental bound on the proton lifetime, this predicts that the 𝑋 gauge boson mediating proton decay is a 5D solitonic string of Planckian tension stretched across a length scale 𝐿 ∼ (1–10 TeV) −1 ending on gauge branes of the same diameter ∼𝐿. This leads to a mass of 𝑀 𝑋 ∼ 10 15 –10 16 GeV. In particular assuming grand unification in the dark dimension scenario results in a tower of Kaluza-Klein excitations of Standard Model gauge bosons on the gauge branes in the 1–10 TeV range. This suggests that the diameter/separation 𝐿 of the gauge branes correlates with both the weak scale ∼1/𝐿 near a TeV and the GUT scale ∼𝑀$^2_5$⁢𝐿 at 10 16 GeV.

branes↗

Friction force in strongly magnetized plasmas

A charged particle moving through a plasma experiences a friction force that commonly acts antiparallel to its velocity. It was recently predicted that in strongly magnetized plasmas, in which the plasma particle gyrofrequency exceeds the plasma frequency, the friction also includes a transverse component that is perpendicular to both the velocity and Lorentz force. Here, this prediction is confirmed using molecular-dynamics simulations, and it is shown that the relative magnitude of the transverse component increases with plasma coupling strength. Furthermore, this result influences single-particle motion and macroscopic transport in strongly magnetized plasmas found in a broad range of applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synchronization of dust acoustic waves in a forced Korteweg–de Vries–Burgers model

The synchronization of dust acoustic waves to an external periodic source is studied in the framework of a driven Korteweg–de Vries–Burgers equation that takes into account the appropriate nonlinear and dispersive nature of low-frequency waves in a dusty plasma medium. For a spatiotemporally varying source term, the system is shown to demonstrate harmonic (1:1) and superharmonic (1:2) synchronized states. Here, the existence domains of these states are delineated in the form of Arnold tongue diagrams in the parametric space of the forcing amplitude and forcing frequency and their resemblance to some past experimental results is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Primordial Black Holes from Long-Range Scalar Forces and Scalar Radiative Cooling

We describe a new scenario for the formation of primordial black holes (PBHs). In the early Universe, the long-range forces mediated by the scalar fields can lead to formation of halos of heavy particles even during the radiation-dominated era. The same interactions result in the emission of scalar radiation from the motion and close encounters of particles in such halos. Radiative cooling due the scalar radiation allows the halos to collapse to black holes. We illustrate this scenario on a simple model with fermions interacting via the Yukawa forces. The abundance and the mass function of PBHs are suitable to account for all dark matter, or for some gravitational wave events detected by LIGO. The model relates the mass of the dark-sector particles to the masses and abundance of dark matter PBHs in a way that can explain why the dark matter and the ordinary matter have similar mass densities. The model also predicts a small contribution to the number of effective light degrees of freedom, which can help reconcile different measurements of the Hubble constant.

79 ASTRONOMY AND ASTROPHYSICS↗

Chemical and steric effects in simulating noncontact atomic force microscopy images of organic molecules on a Cu (111) substrate

Noncontact mode of atomic force microscopy (nc-AFM) employing a CO-functionalized tip is a very powerful tool for studying molecular structures. However, interpreting nc-AFM images for nonplanar molecules can sometimes be problematic. To illustrate and resolve the nature of such problematic systems, we employ real-space pseudopotentials constructed within density functional theory to simulate nc-AFM images. In this work, we focus on several representative nonplanar organic molecules (pentacene, naphthanthrone, olympicene, and 6-phenylhexa-1,3,5-triynybenzene (PHTB)) on a typical substrate: the Cu (111) surface. This substrate results in significant distortions in the molecular geometries of pentacene and naphthanthrone. Including these distortions in simulated nc-AFM imaging notably improves the agreement between the simulated and measured images. In naphthanthrone, the relatively large interaction between the O atom and the Cu substrate offers a straightforward explanation for the absence of the C = O bond in the measured image. Nonplanar features such as tilting or twisting are also apparent in olympicene and PHTB. A “triangular” bright feature associated with the –CH 2 group in olympicene appears in simulated and measured nc-AFM images. This feature is directly related to the tilting angle of the molecule with respect to the Cu substrate. The “defective” benzene ring feature and the faint ellipsoidal C ≡ C feature in PHTB can be ascribed to its twisted nature. The ability to simulate such subatomic images in nc-AFM reflects the accuracy and efficiency of calculating quantum forces in real space.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cyclic phase transition of substrate-modulated two-dimensional dusty plasma driven by gyroscopic forces

Langevin dynamical simulations are performed to investigate the formation of clusters and voids of a two-dimensional periodic substrate (2DPS)-modulated two-dimensional dusty plasma (2DDP) driven by a gyroscopic force. It is discovered that, as the frequency of the gyroscopic force decreases gradually, the substrate-modulated 2DDP undergoes the cyclic transition between the ordered cluster and void phases. Between the observed ordered cluster and void phases, the studied 2DDP exhibits a more uniform arrangement of particles. The discovered cyclic transition is attributed to the symmetry of the time-averaged potential landscape induced by the 2DPS in the reference frame of the moving particle, as confirmed by superimposing the particle locations on the effective potential landscape under various conditions.

36 MATERIALS SCIENCE↗

Quadrupole forces between quark and gluon subsystems inside higher-spin particles

We generalize the mechanical interpretation of the forces between quark and gluon subsystems, previously studied for the nucleon, to arbitrary higher-spin particles. For spin-0 and spin-1/2 particles, this force is characterized by the nonconserved $\overline{𝑐}$⁡(𝑡) form factor. However, such an interpretation has not yet been established for higher-spin particles due to the intricate structure of the nonconserved energy-momentum tensor (EMT) form factors. By performing a multipole expansion, we identify the physically meaningful combinations of the nonconserved covariant EMT form factors and provide them with a clear mechanical interpretation.

Baryons↗