Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical force fields”

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

Perspective on Some Recent and Future Developments in Casimir Interactions

Here, we present a critical review of recent developments in Casimir physics motivated by discoveries of novel materials. Specifically, topologically nontrivial properties of the graphene family, Chern and topological insulators, and Weyl semimetals have diverse manifestations in the distance dependence, presence of fundamental constants, magnitude, and sign of the Casimir interaction. Limited studies of the role of nonlinear optical properties in the interaction are also reviewed. We show that, since many new materials have greatly enhanced the nonlinear optical response, new efficient pathways for investigation of the characteristic regimes of the Casimir force need to be explored, which are expected to lead to new discoveries. Recent progress in the dynamical Casimir effect is also reviewed and we argue that nonlinear media can open up new directions in this field as well.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Large ordered moment with strong easy-plane anisotropy and vortex-domain pattern in the kagome ferromagnet Fe 3 Sn

We report the magnetic anisotropy of kagome bilayer ferromagnet Fe 3 Sn probed by the bulk magnetometry and magnetic force microscopy (MFM) on high-quality single crystals. The dependence of magnetization on the orientation of the external magnetic field reveals strong easyplane magnetocrystalline anisotropy and anisotropy of the saturation magnetization. The leading magnetocrystalline anisotropy constant shows a monotonous increase from $\mathrm{K_1$≈ –1.0 \times 10^6 J/m^3}$ at 300 K to $\mathrm{–1.3 \times 10^6 J/m^3}$ at 2 K. Our ab initio electronic structure calculations yield the value of total magnetic moment of 7.1$μ_Β$/$f.u.$ and a magnetocrystalline anisotropy energy density of –0.57 $\mathrm{meV/}$$f.u.$ $\mathrm{(–1.62 \times 10^6J/m^3)}$ both being in reasonable agreement with the experimental values. The MFM imaging reveals micrometer-scale magnetic vortices with weakly pinned cores that vanish at the saturation field of ~3T applied perpendicular to the kagome plane. The observed vortex-domain structure is well reproduced by the micromagnetic simulations, using the experimentally determined value of the anisotropy and exchange stiffness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamic, Adaptive, Systems and Materials: Complex, Simple and Emergent Behaviors

This program has been funded by DoE/BES for twenty years. It has moved into and out of various subjects as it has developed, but it has retained its focus on complexity and complex systems. The project has evolved in the following way: Self-Assembly and Biomimetic Self-Assembly: All self-assembling systems depend upon a minimum of two types of interaction: a repulsion and an attraction. For the familiar molecular systems, attractive interactions are typically hydrogen bonds and electrostatic interactions. Repulsive interactions include steric effects, hydrophobic effects (in biological systems), and charge-charge repulsion. We have expanded this repertoire to include surface interactions, magnetic interactions, and others. I list these systems in the order in which we have explored them: i) A key emphasis in current work is in understanding how the movement of ions in a magnetic field (the Lorentz effect) interacts with catalytic systems. We have demonstrated that an acceleration in rate of reduction of CO 2 to CO can be accomplished by applying an external magnetic field. This acceleration is largely due to the application of the Lorentz effect on mass transport at the catalyst’s surface. ii) We have also extensively explored the influence of electrostatics, as exhibited in self-assembling systems, by tribocharging. iii) Another key system involves surface tension effects; examples include interactions between heavy particles floating at a liquid-air interface, and interacting by changes in surface area; interactions of bubbles and bubble rafts, behaviors of bubble trains in microfluidic networks, and behaviors of microorganisms in constraining environments. iv) This work has intentionally de-emphasized biological systems; but it does include some work on protein-ligand interactions and interactions among microorganisms. v) We have also explored applications of some of these effects, these explorations include bubble rafts as diffraction gratings, exploration of the structures that can be obtained by tribocharging and uses of these structures in exploring nucleation and melting of crystals. vi) Although not a major focus of this work, several other topics have emerged and offer opportunities for future work. These include the behavior of bubble trains and bubble rafts in microfluidic systems. A particularly interesting example is the formation of bubble trains that repeat in the alteration of large and small bubbles according to rules we do not presently understand, but are uniquely large-period oscillating systems. These systems offer a new route into understanding the instabilities of the type represented by oscillations. vii) We have also begun exploratory projects on magnetic levitation (especially to determine molecular density), and information storage (in molecules). Magnetic Levitation: Self-assembling and biomimetic systems require both attraction and repulsion. We have used electrostatics (tribocharging), interfacial free-energies (surface tension and related forces) and others. Potential uses include reconfigurable diffraction gratings and liquid lenses; exploration of mechanisms in tribocharging; tunneling in EGaIn junctions; and bubble trains (especially in micro-fluidic systems). Examples of systems representing these topics is included in the following papers: Complexity: Disks rotating at a water-air interface; Benard-Marangoni effects; Vortex-Crystals from spinning magnetic disks (Marangoni effects); EGaIn Electrode to study quantum tunneling; Self-Assembly of electrostatically-charged metallic spheres (electrets); Dynamically reconfigurable lens; Using computational designs of ligands for enzymes; Electrostatic self-assembly by tribocharging; Monodisperse bubble trains in microchannel systems; Inverted dripping faucet; Flames; Printing of micro-organisms to regenerate the “ink” of printing device; Using micro-organisms to move loads (“microoxen”); Motion of bacterial swarms near surfaces; Making monodisperse particles in microfluidic systems; Coding/decoding of information stored in droplet trains in microfluidic networks; Magnetic levitation; and Information storage. i) Tribocharging. The change in focus of this work on electrets from the fundamentals of charging to applications of these materials in studying self-assembly using electrostatic interactions. ii) Bubbles in Microchannels. The realization that systems of bubbles in microchannels represented a major opportunity to study complexity in a very tractable system, and the development of a semi-quantitative theory of this subject. iii) Flames. The growth of “flames” remains an exploratory subject for the research, although their currently relatively little active work involving it ongoing. iv) Systems with Microorganisms. The removal of work in biological systems from this project. Based on work supported in this program, we now have a significant project on the development of microfluidic tools for studying C. elegans (a nematode), but this work was not appropriate for a program focused on complexity, and we developed separate support for it. (It is, however, an example of successful seeding of a new area by BES.) The work on electrets has gone through a period in which a part of the program was the subject of a MURI; the focus of this work was to develop materials that did not charge electrostatically on friction or contact. The MURI is now over, and the work on dynamic self-assembly (supported by BES) is the major focus. “Flames” has also enjoyed synergistic support, in terms of a project supported by DARPA on flame suppression (in the absence of extinguishing agents, using acoustic and electrostatic interactions). This work was helpful in understanding some of the basics of flames, but is entirely distinct from the BES focus in complexity. A growing interest is in the Lorentz effect. The Lorentz effect is the force exerted on charged particles (electrons, ions, charged molecules) when they move through a perpendicular magnetic field. The Lorentz effect is almost ubiquitous in modern technology: examples of applications include electric motors, dynamos, cathode ray tubes, many batteries, and most systems that control electrical currents with magnetic forces. We have begun to explore the Lorentz effect in electrochemical systems and heterogeneous catalytic systems involving charged organic species and inorganic ions. This work is still at an early stage, but initial studies that Lorentz effects can be large when ions move through magnetic fields, or magnetic fields move in the presence of ions.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Gauged Q-balls

Classical non-topological soliton configurations are considered within the theory of a complex scalar field with a gauged U symmetry. Their existence and stability against dispersion are demonstrated and some of their properties are investigated analytically and numerically. The soliton configuration is such that inside the soliton the local U symmetry is broken, the gauge field becomes massive and for a range of values of the coupling constants the soliton becomes a superconductor pushing the charge to the surface. Furthermore, because of the repulsive Coulomb force, there is a maximum size for these objects, making impossible the existence of Q-matter in bulk form. Also briefly discussed are solitons with fermions in a U gauge theory.

Lee, Kimyeong↗
NASA NTRS2013

Lessons Learned in the Selection and Development of Test Cases for the Aeroelastic Prediction Workshop: Rectangular Supercritical Wing

The Aeroelastic Prediction Workshop brought together an international community of computational fluid dynamicists as a step in defining the state of the art in computational aeroelasticity. The Rectangular Supercritical Wing (RSW) was chosen as the first configuration to study due to its geometric simplicity, perceived simple flow field at transonic conditions and availability of an experimental data set containing forced oscillation response data. Six teams performed analyses of the RSW; they used Reynolds-Averaged Navier-Stokes flow solvers exercised assuming that the wing had a rigid structure. Both steady-state and forced oscillation computations were performed by each team. The results of these calculations were compared with each other and with the experimental data. The steady-state results from the computations capture many of the flow features of a classical supercritical airfoil pressure distribution. The most dominant feature of the oscillatory results is the upper surface shock dynamics. Substantial variations were observed among the computational solutions as well as differences relative to the experimental data. Contributing issues to these differences include substantial wind tunnel wall effects and diverse choices in the analysis parameters.

Heeg, Jennifer↗
DOE OSTI2020

On the proper derivation of the Floquet-based quantum classical Liouville equation and surface hopping describing a molecule or material subject to an external field

We investigate different approaches to derive the proper Floquet-based quantum–classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis and then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form that is similar to the usual QCLE with two modifications: (1) The electronic degrees of freedom are expanded to infinite dimension and (2) the nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a fictitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., [J. Chem. Phys. 140, 084104 (2014)] who found that transforming to the adiabatic representation must always be the last operation applied, although now we have extended this result to a time-dependent Hamiltonian. Finally, this paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.

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

Tuning interfacial Dzyaloshinskii-Moriya interactions in thin amorphous ferrimagnetic alloys

Skyrmions can be stabilized in magnetic systems with broken inversion symmetry and chiral interactions, such as Dzyaloshinskii-Moriya interactions (DMI). Further, compensation of magnetic moments in ferrimagnetic materials can significantly reduce magnetic dipolar interactions, which tend to favor large skyrmions. Tuning DMI is essential to control skyrmion properties, with symmetry breaking at interfaces offering the greatest flexibility. However, in contrast to the ferromagnet case, few studies have investigated interfacial DMI in ferrimagnets. Here we present a systematic study of DMI in ferrimagnetic CoGd films by Brillouin light scattering. We demonstrate the ability to control DMI by the CoGd cap layer composition, the stack symmetry and the ferrimagnetic layer thickness. The DMI thickness dependence confirms its interfacial nature. In addition, magnetic force microscopy reveals the ability to tune DMI in a range that stabilizes sub-100 nm skyrmions at room temperature in zero field. Our work opens new paths for controlling interfacial DMI in ferrimagnets to nucleate and manipulate skyrmions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air

Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.

42 ENGINEERING↗
NASA NTRS2009

Lightweight Bulldozer Attachment for Construction and Excavation on the Lunar Surface

A lightweight bulldozer blade prototype has been designed and built to be used as an excavation implement in conjunction with the NASA Chariot lunar mobility platform prototype. The combined system was then used in a variety of field tests in order to characterize structural loads, excavation performance and learn about the operational behavior of lunar excavation in geotechnical lunar simulants. The purpose of this effort was to evaluate the feasibility of lunar excavation for site preparation at a planned NASA lunar outpost. Once the feasibility has been determined then the technology will become available as a candidate element in the NASA Lunar Surface Systems Architecture. In addition to NASA experimental testing of the LANCE blade, NASA engineers completed analytical work on the expected draft forces using classical soil mechanics methods. The Colorado School of Mines (CSM) team utilized finite element analysis (FEA) to study the interaction between the cutting edge of the LANCE blade and the surface of soil. FEA was also used to examine various load cases and their effect on the lightweight structure of the LANCE blade. Overall it has been determined that a lunar bulldozer blade is a viable technology for lunar outpost site preparation, but further work is required to characterize the behavior in 1/6th G and actual lunar regolith in a vacuum lunar environment.

Mueller, Robert↗
NASA NTRS1996

Design of a Superconducting Magnetic Suspension System for a Liquid Helium Flow Experiment

We discuss a preliminary design for a superconducting magnetic suspension system for measurement of drag on rotationally symmetric bodies in liquid helium. Superconducting materials are a natural choice for liquid helium studies, since temperatures are well below most critical temperatures, so that the resulting heat load is negligible. Also, due to its diamagnetic properties, a superconducting model (for example made or coated with Nb) is inherently stable against disturbances. Issues which we consider include model placement during initial cool-down, maintaining placement during anticipated drag and lift forces, and force measurement. This later can be achieved by a passive technique, where the body is allowed to deflect under the influence of drag from its neutral position. The resulting shift in flux is detected via a superconducting pickup coil. The pickup coil may be connected either to a SQUID, or a secondary loop wound around a Hall probe. Both options are discussed. The objective of this work is to gain a better understanding of the nature of turbulent fields in normal and superfluid helium for potential application to problems in classical high Reynolds number turbulence.

Smith, Michael R.↗
NASA NTRS2015

Lunar Orbit Stability for Small Satellite Mission Design

The irregular nature of the lunar gravity field will severely affect the orbit lifetime and behavior of future lunar small satellite missions. These spacecraft need stable orbits that do not require large deltaV budgets for station-keeping maneuvers. The initial classical elements of any lunar orbit are critical to address its stability and to comply with mission requirements. This publication identifies stable regions according to different initial conditions at the time of lunar orbit insertion (LOI). High fidelity numerical simulations with two different gravity models were performed. We focus in low altitude orbits where the dominant force in orbit propagation is the existence of unevenly distributed lunar mass concentrations. These orbits follow a periodic oscillation in some of the classical elements that is particularly useful for mission design. A set of orbital maintenance strategies for various mission concepts is presented.

Lunar Science↗
DOE OSTI2024

Gravitational scattering and beyond from extreme mass ratio effective field theory

We explore a recently proposed effective field theory describing electromagnetically or gravitationally interacting massive particles in an expansion about their mass ratio, also known as the self-force (SF) expansion. By integrating out the deviation of the heavy particle about its inertial trajectory, we obtain an effective action whose only degrees of freedom are the lighter particle together with the photon or graviton, all propagating in a Coulomb or Schwarzschild background. The 0SF dynamics are described by the usual background field method, which at 1SF is supplemented by a “recoil operator” that encodes the wobble of the heavy particle, and similarly computable corrections appearing at 2SF and higher. Our formalism exploits the fact that the analytic expressions for classical backgrounds and particle trajectories encode dynamical information to all orders in the couplings, and from them we extract multiloop integrands for perturbative scattering. As a check, we study the two-loop classical scattering of scalar particles in electromagnetism and gravity, verifying known results. We then present new calculations for the two-loop classical scattering of dyons, and of particles interacting with an additional scalar or vector field coupling directly to the lighter particle but only gravitationally to the heavier particle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy

Portlandite [Ca(OH) 2 ] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomic-scale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH) 2 ] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO 2 (OH) n ] 2–n (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum Computing for High-Energy Physics: State of the Art and Challenges

Quantum computers offer an intriguing path for a paradigmatic change of computing in the natural sciences and beyond, with the potential for achieving a so-called quantum advantage—namely, a significant (in some cases exponential) speedup of numerical simulations. The rapid development of hardware devices with various realizations of qubits enables the execution of small-scale but representative applications on quantum computers. In particular, the high-energy physics community plays a pivotal role in accessing the power of quantum computing, since the field is a driving source for challenging computational problems. This concerns, on the theoretical side, the exploration of models that are very hard or even impossible to address with classical techniques and, on the experimental side, the enormous data challenge of newly emerging experiments, such as the upgrade of the Large Hadron Collider. In this Roadmap paper, led by CERN, DESY, and IBM, we provide the status of high-energy physics quantum computations and give examples of theoretical and experimental target benchmark applications, which can be addressed in the near future. Having in mind hardware with about 100 qubits capable of executing several thousand two-qubit gates, where possible, we also provide resource estimates for the examples given using error-mitigated quantum computing. The ultimate declared goal of this task force is therefore to trigger further research in the high-energy physics community to develop interesting use cases for demonstrations on near-term quantum computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

First-principles molten salt phase diagrams through thermodynamic integration

Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. Here, we show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid–liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange–correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid–liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials.

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

Simulating the non-monotonic strain response of nanoporous multiferroic composites under electric field control

In this work, we simulate and analyze the mechanical response of a class of multiferroic materials consisting of a templated porous nanostructure made out of cobalt ferrite (CFO) partially filled by atomic layer deposition (ALD) with a ferroelectric phase of lead zirconate titanate (PZT). The strain in the device is measured when an electric field is applied for varying ALD thicknesses, displaying a non-monotonic dependence with a maximum strain achieved for a coating thickness of 3 nm. To understand this behavior, we apply finite element modeling to the smallest repeatable unit of the nanoporous template and simulate the mechanical response as a function of PZT coating thickness. We find that this non-monotonic response is caused by the interplay between two driving forces opposing one another. First, increased porosity works toward increasing the strain due to a reduced system stiffness. Second, decreased porosity involves a larger mass fraction of PZT, which drives the electro-mechanical response of the structure, thus leading to a larger strain. The balance between these two driving forces is controlled by the shear coupling at the CFO/PZT interface and the effective PZT cross section along the direction of the applied electric field. Here, our numerical results show that considering a nonlinear piezoelectric response for PZT leads to an improved agreement with the experimental data, consistent with ex situ poling of the nanostructure prior to magnetic measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

AFM-IR and s-SNOM-IR measurements of chemically addressable monolayers on Au nanoparticles

The performance of catalysts depends on their nanoscale properties, and local variations in structure and composition can have a dramatic impact on the catalytic reactivity. Therefore, probing the localized reactivity of catalytic surfaces using high spatial resolution vibrational spectroscopy, such as infrared (IR) nanospectroscopy and tip-enhanced Raman spectroscopy, is essential for mapping their reactivity pattern. Two fundamentally different scanning probe IR nanospectroscopy techniques, namely, scattering-type scanning near-field optical microscopy (sSNOM) and atomic force microscopy-infrared spectroscopy (AFM-IR), provide the capabilities for mapping the reactivity pattern of catalytic surfaces with a spatial resolution of ~20 nm. Herein, we compare these two techniques with regard to their applicability for probing the vibrational signature of reactive molecules on catalytic nanoparticles. For this purpose, we use chemically addressable self-assembled molecules on Au nanoparticles as model systems. We identified significant spectral differences depending on the measurement technique, which originate from the fundamentally different working principles of the applied methods. While AFM-IR spectra provided information from all the molecules that were positioned underneath the tip, the s-SNOM spectra were more orientation-sensitive. Due to its field-enhancement factor, the s-SNOM spectra showed higher vibrational signals for dipoles that were perpendicularly oriented to the surface. The s-SNOM sensitivity to the molecular orientation influenced the amplitude, position, and signal-to-noise ratio of the collected spectra. Ensemble-based IR measurements verified that differences in the localized IR spectra stem from the enhanced sensitivity of s-SNOM measurements to the adsorption geometry of the probed molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗