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

Cloud adjustments from large-scale smoke–circulation interactions strongly modulate the southeastern Atlantic stratocumulus-to-cumulus transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (~1–2 K) and greatly enhanced moisture (>2 g kg -1 ) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

Cloud Adjustments From Large-Scale Smoke–Circulation Interactions Strongly Modulate the Southeastern Atlantic Stratocumulus-to-Cumulus Transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (∼1–2 K) and greatly enhanced moisture (>2 g kg−1) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

southeast Atlantic Ocean↗
DOE OSTI2020

Disorder effects in nitride semiconductors: impact on fundamental and device properties

Abstract Semiconductor structures used for fundamental or device applications most often incorporate alloy materials. In “usual” or “common” III–V alloys, based on the InGaAsP or InGaAlAs material systems, the effects of compositional disorder on the electronic properties can be treated in a perturbative approach. This is not the case in the more recent nitride-based GaInAlN alloys, where the potential changes associated with the various atoms induce strong localization effects, which cannot be described perturbatively. Since the early studies of these materials and devices, disorder effects have indeed been identified to play a major role in their properties. Although many studies have been performed on the structural characterization of materials, on intrinsic electronic localization properties, and on the impact of disorder on device operation, there are still many open questions on all these topics. Taking disorder into account also leads to unmanageable problems in simulations. As a prerequisite to address material and device simulations, a critical examination of experiments must be considered to ensure that one measures intrinsic parameters as these materials are difficult to grow with low defect densities. A specific property of nitride semiconductors that can obscure intrinsic properties is the strong spontaneous and piezoelectric fields. We outline in this review the remaining challenges faced when attempting to fully describe nitride-based material systems, taking the examples of LEDs. The objectives of a better understanding of disorder phenomena are to explain the hidden phenomena often forcing one to use ad hoc parameters, or additional poorly defined concepts, to make simulations agree with experiments. Finally, we describe a novel simulation tool based on a mathematical breakthrough to solve the Schrödinger equation in disordered potentials that facilitates 3D simulations that include alloy disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Observation of quantum effects on radiation reaction in strong fields

Radiation reaction, the force experienced by an accelerated charge due to radiation emission, has long been the subject of extensive theoretical and experimental research. Experimental verification of a quantum, strong-field description of radiation reaction is fundamentally important, and has wide-ranging implications for astrophysics, laser-driven particle acceleration, next-generation particle colliders and inverse-Compton photon sources for medical and industrial applications. However, the difficulty of accessing regimes where strong field and quantum effects dominate inhibited previous efforts to observe quantum radiation reaction in charged particle dynamics with high significance. We report the first high significance (> 5σ) observation of strong-field radiation reaction on electron spectra where quantum effects are substantial. We obtain the first, quantitative, strong evidence favouring the quantum-continuous and quantum-stochastic models over the classical model; the quantum models perform comparably. The lower electron energy losses predicted by the quantum models accounts for their improved performance. Model comparison was performed using a novel Bayesian framework which has widespread utility for laser-particle collision experiments, including those utilising conventional accelerators, where some collision parameters cannot be measured directly.

Los, Eva E. [Imperial College London (United Kingd↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗
DOE OSTI2021

Effect of chemical substituents attached to the zwitterion cation on dielectric constant

Materials with high dielectric constant, ε s , are desirable in a wide range of applications including energy storage and actuators. Recently, zwitterionic liquids have been reported to have the largest ε s of any liquid and, thus, have the potential to replace inorganic fillers to modulate the material ε s . Although the large ε s for zwitterionic liquids is attributed to their large molecular dipole, the role of chemical substituents attached to the zwitterion cation on ε s is not fully understood, which is necessary to enhance the performance of soft energy materials. Here, we report the impact of zwitterionic liquid cation chemical substituents on ε s (50 < ε s < 300 at room temperature). Dielectric relaxation spectroscopy reveals that molecular reorientation is the main contributor to the high ε s . The low Kirkwood factor g calculated for zwitterionic liquids (e.g., 0.1–0.2) suggests the tendency for the antiparallel zwitterion dipole alignment expected from the strong electrostatic intermolecular interactions. With octyl cation substituents, the g is decreased due to the formation of hydrophobic-rich domains that restrict molecular reorientation under applied electric fields. In contrast, when zwitterion cations are functionalized with ethylene oxide (EO) segments, g increases due to the EO segments interacting with the cations, allowing more zwitterion rotation in response to the applied field. The reported results suggest that high ε s zwitterionic liquids require a large molecular dipole, compositionally homogeneous liquids (e.g., no aggregation), a maximized zwitterion number density, and a high g, which is achievable by incorporating polar chemical substituents onto the zwitterion cations.

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

Confinement‐Driven Inverse Domain Scaling in Polycrystalline ErMnO 3

Abstract The research on topological phenomena in ferroelectric materials has revolutionized the way people understand polar order. Intriguing examples are polar skyrmions, vortex/anti‐vortex structures, and ferroelectric incommensurabilties, which promote emergent physical properties ranging from electric‐field‐controllable chirality to negative capacitance effects. Here, the impact of topologically protected vortices on the domain formation in improper ferroelectric ErMnO 3 polycrystals is studied, demonstrating inverted domain scaling behavior compared to classical ferroelectrics. It is observed that as the grain size increases, smaller domains are formed. Phase field simulations reveal that elastic strain fields drive the annihilation of vortex/anti‐vortex pairs within the grains and individual vortices at the grain boundaries. The inversion of the domain scaling behavior has far‐reaching implications, providing fundamentally new opportunities for topology‐based domain engineering and the tuning of the electromechanical and dielectric performance of ferroelectrics in general.

36 MATERIALS SCIENCE↗
NASA NTRS2007

SIMBOL-X: A Formation Flying Mission on HEO for Exploring the Universe

SIMBOL-X is a high energy new generation telescope covering by a single instrument a continuous energy range starting at classical X-rays and extending to hard X-rays, i.e. from 0.5 to 80 keV. It is using in this field a focalizing payload which until now was used for energy below 10 keV only, via the construction of a telescope distributed on two satellites flying in formation. SIMBOL-X permits a gain of two orders of magnitude in sensibility and spatial resolution in comparison to state of the art hard X-rays instruments. The mirror satellite will be in free flight on a high elliptical orbit and will target the object to observe very precisely, thus focusing the hard X-ray emission thanks to this mirror module. At the focal point area which is situated 20 meters behind the mirror satellite, the detector satellite maintains its position on a forced orbit thanks to a radio link with the mirror satellite and a lateral displacement sensor using a beam emitted onboard the mirror satellite. This configuration is said "formation flying". The location of the detector satellite shall be very finely tuned as it carries the focal plane of this distributed telescope. To provide science measurements, the Simbol-X orbit has been chosen High elliptic (HEO), which means elliptical orbit with a high perigee altitude. Preliminary studies where made with an orbit with an altitude of the perigee of 44000km and altitude of the apogee of 253000km. The orbit was seven days ground track repeated in order to maintain a perigee pass over the Malindi ground station to download scientific telemetry. But as studies went on, difficulties in mass budget, link budget, perigee maintenance and formation flying maintenance were raised. This was mainly due to the vicinity of the Moon and its disturbing effect on the satellites orbits. Alternative orbits have been proposed in order to demonstrate the feasibility of the mission. The problematic of bringing the two satellites from their injection orbit to their operational orbit 20 m apart from each other and then maintain this configuration is very challenging. It requires theoretical development of the relative motion between two satellites in high eccentric orbit with large differential disturbance on the two bodies. This paper will present the mission analysis for the Simbol-X satellites with the complex problematic of doing formation flying in high elliptic orbit.

Gamet, Philippe↗
DOE OSTI2023

Hydrogen irradiation-driven computational surface chemistry of lithium oxide and hydroxide

We have investigated, using molecular dynamics, the surface chemistry of hydrogen incident on the amorphous and crystalline lithium oxide and lithium hydroxide surfaces upon being slowed down by a collision cascade and retained in the amorphous surface of either Li 2 O or LiOH. We looked for the bonding of H to the resident structures in the surface to understand a possible chain of chemical reactions that can lead to surface transformation upon H atom impact. Our findings, using Density-Functional Theory (DFT) trained ReaxFF force field/electronegativity equalization method potentials, stress the importance of inclusion of polarization in the dynamics of a Li–O–H system, which is also illustrated by DFT energy minimization and quantum–classical molecular dynamics using tight binding DFT. The resulting polar-covalent chemistry of the studied systems is complex and very sensitive to the instantaneous positions of all atoms as well as the ratio of concentrations of various resident atoms in the surface.

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

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2025

Quantum confining excitons with an electrostatic moiré superlattice

Quantum confining excitons has been a persistent challenge in the pursuit of strong exciton interactions and quantum light generation. Unlike electrons, which can be readily controlled via electric fields, imposing strong nanoscale potentials on excitons to enable quantum confinement has proven challenging. In this study, we utilize piezoelectric force microscopy to image the domain structures of twisted hexagonal boron nitride (h – BN), revealing evidence of strong in-plane electric fields at the domain boundaries. By placing a monolayer MoSe 2 only one to two nanometers away from the twisted h – BN interface, we observe energy splitting of neutral excitons and Fermi polarons by several millielectronvolts at the moiré domain boundaries. By directly correlating local structural and optical properties, we attribute such observations to excitons confined in a nanoscale one-dimensional electrostatic potential created by the strong in-plane electric fields at the moiré domain boundaries. Intriguingly, this 1D quantum confinement results in pronounced polarization anisotropy in the excitons’ reflection and emission, persistent to temperatures as high as ~80 Kelvins. Furthermore, these findings open new avenues for exploring and controlling strongly interacting excitons for classical and quantum optoelectronics.

2-dimensional systems↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2021

Reconstructing phase-resolved hysteresis loops from first-order reversal curves

The first order reversal curve (FORC) method is a magnetometry based technique used to capture nanoscale magnetic phase separation and interactions with macroscopic measurements using minor hysteresis loop analysis. This makes the FORC technique a powerful tool in the analysis of complex systems which cannot be effectively probed using localized techniques. However, recovering quantitative details about the identified phases which can be compared to traditionally measured metrics remains an enigmatic challenge. We demonstrate a technique to reconstruct phase-resolved magnetic hysteresis loops by selectively integrating the measured FORC distribution. From these minor loops, the traditional metrics—including the coercivity and saturation field, and the remanent and saturation magnetization—can be determined. In order to perform this analysis, special consideration must be paid to the accurate quantitative management of the so-called reversible features. This technique is demonstrated on three representative materials systems, high anisotropy FeCuPt thin-films, Fe nanodots, and SmCo/Fe exchange spring magnet films, and shows excellent agreement with the direct measured major loop, as well as the phase separated loops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electro-optic transduction in silicon via gigahertz-frequency nanomechanics

Interfacing electronics with optical fiber networks is key to the long-distance transfer of classical and quantum information. Piezo-optomechanical transducers enable such interfaces by using gigahertz-frequency acoustic vibrations as mediators for converting microwave photons to optical photons via the combination of optomechanical and piezoelectric interactions. However, despite successful demonstrations, efficient quantum transduction remains out of reach due to the challenges associated with hybrid material integration and increased loss from piezoelectric materials when operating in the quantum regime. Here, we demonstrate an alternative approach in which we actuate 5-GHz phonons in a conventional silicon-on-insulator platform. In our experiment, microwave photons resonantly drive a phononic crystal oscillator via the electrostatic force realized in a charge-biased narrow-gap capacitor. The mechanical vibrations are subsequently transferred via a phonon waveguide to an optomechanical cavity, where they transform into optical photons in the sideband of a pump laser field. Operating at room temperature and atmospheric pressure, we measure a microwave-to-optical photon conversion efficiency of 1.72±0.14×10 −7 in a 3.3 MHz bandwidth. Our results mark a stepping stone towards quantum transduction with integrated devices made from crystalline silicon, which promise efficient high-bandwidth operation and integration with superconducting qubits. Additionally, the lack of need for piezoelectricity or other intrinsic nonlinearities makes our approach applicable to a wide range of materials for potential applications beyond quantum technologies.

42 ENGINEERING↗
DOE OSTI2022

The influence of particle size on the fluid dynamics of a laser-induced plasma

The interaction of a laser-induced shock wave with nanoparticles and microparticles of aluminum oxide is investigated through experiments and modeling. The chemistry and physics of the interaction between the particles and plasma generated from laser ablation shows similarities and discrete differences for the two particle sizes. For both particle sizes, early stage (< 10 μs) ionization was dominant and evidenced by higher concentrations of Al II. While both sizes exhibit ionization over the same duration, the intensity of emission was greater for nanoparticles indicating greater concentrations of ionized species. Moreover, the dispersion of species was notably more elongated for microparticles while radial dispersion was more pronounced for nanoparticles with elevated drag forces. At later stages (i.e., > 10 μs), oxidation reactions were dominant for both particle sizes, but the same distinctions in flow field were observed and attributed to particle drag. In all stages of interaction, microparticles expand axially with less drag that suppresses their radial expansion. As a result, the dispersion of reactive species was mapped over an up to 80% larger area for nanoparticles relative to microparticles. Results shown here can be applied toward advancing experimental diagnostics and particle-shock wave modeling and simulation efforts for energetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Fast Quantum Algorithm for Predicting Descriptive Statistics of Stochastic Processes

Stochastic processes are used as a modeling tool in several sub-fields of physics, biology, and finance. Analytic understanding of the long term behavior of such processes is only tractable for very simple types of stochastic processes such as Markovian processes. However, in real world applications more complex stochastic processes often arise. In physics, the complicating factor might be nonlinearities; in biology it might be memory effects; and in finance is might be the non-random intentional behavior of participants in a market. In the absence of analytic insight, one is forced to understand these more complex stochastic processes via numerical simulation techniques. In this paper we present a quantum algorithm for performing such simulations. In particular, we show how a quantum algorithm can predict arbitrary descriptive statistics (moments) of N-step stochastic processes in just O(square root of N) time. That is, the quantum complexity is the square root of the classical complexity for performing such simulations. This is a significant speedup in comparison to the current state of the art.

Williams Colin P.↗
NASA NTRS2014

Data-Intensive Science meets Inquiry-Driven Pedagogy: Interactive Big Data Exploration, Threshold Concepts, and Liminality

Threshold concepts in any discipline are the core concepts an individual must understand in order to master a discipline. By their very nature, these concepts are troublesome, irreversible, integrative, bounded, discursive, and reconstitutive. Although grasping threshold concepts can be extremely challenging for each learner as s/he moves through stages of cognitive development relative to a given discipline, the learner's grasp of these concepts determines the extent to which s/he is prepared to work competently and creatively within the field itself. The movement of individuals from a state of ignorance of these core concepts to one of mastery occurs not along a linear path but in iterative cycles of knowledge creation and adjustment in liminal spaces - conceptual spaces through which learners move from the vaguest awareness of concepts to mastery, accompanied by understanding of their relevance, connectivity, and usefulness relative to questions and constructs in a given discipline. For example, challenges in the teaching and learning of atmospheric science can be traced to threshold concepts in fluid dynamics. In particular, Dynamic Meteorology is one of the most challenging courses for graduate students and undergraduates majoring in Atmospheric Science. Dynamic Meteorology introduces threshold concepts - those that prove troublesome for the majority of students but that are essential, associated with fundamental relationships between forces and motion in the atmosphere and requiring the application of basic classical statics, dynamics, and thermodynamic principles to the three dimensionally varying atmospheric structure. With the explosive growth of data available in atmospheric science, driven largely by satellite Earth observations and high-resolution numerical simulations, paradigms such as that of dataintensive science have emerged. These paradigm shifts are based on the growing realization that current infrastructure, tools and processes will not allow us to analyze and fully utilize the complex and voluminous data that is being gathered. In this emerging paradigm, the scientific discovery process is driven by knowledge extracted from large volumes of data. In this presentation, we contend that this paradigm naturally lends to inquiry-driven pedagogy where knowledge is discovered through inductive engagement with large volumes of data rather than reached through traditional, deductive, hypothesis-driven analyses. In particular, data-intensive techniques married with an inductive methodology allow for exploration on a scale that is not possible in the traditional classroom with its typical problem sets and static, limited data samples. In addition, we identify existing gaps and possible solutions for addressing the infrastructure and tools as well as a pedagogical framework through which to implement this inductive approach.

Ramachandran, Rahul↗
DOE OSTI2020

Seeing is believing: visualization of antiferromagnetic domains

Abstract Understanding and utilizing novel antiferromagnetic (AFM) materials has been recently one of the central issues in condensed matter physics, as well as in materials science and engineering. The relevant contemporary topics include multiferroicity, topological magnetism and AFM spintronics. The ability to image magnetic domains in AFM materials is of key importance for the success of these exciting fields. While imaging techniques of magnetic domains on the surfaces of ferro-(ferri)magnetic materials with, for example, magneto-optical Kerr microscopy and magnetic force microscopy have been available for a number of decades, AFM domain imaging is a relatively new development. We review various experimental techniques utilizing scanning, optical, and synchrotron X-ray probes to visualize AFM domains and domain walls, and to unveil their physical properties. We also discuss the existing challenges and opportunities in these techniques, especially with further increase of spatial and temporal resolution.

36 MATERIALS SCIENCE↗