Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex systems theory”

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 937 records · Page 52

Strongly correlated itinerant magnetism near superconductivity in NiTa 4 Se 8

Metallic ferromagnets with strongly interacting electrons often exhibit remarkable electronic phases such as ferromagnetic superconductivity, complex spin textures, and nontrivial topology. Here, in this report, we discuss the synthesis of a layered magnetic metal NiTa 4 Se 8 (or Ni 1/4 TaSe 2 ) with a Curie temperature of 58 Kelvin. Magnetization data and density functional theory calculations indicate that the nickel atoms host uniaxial ferromagnetic order of about 0.7µ B per atom, while an even smaller moment is generated in the itinerant tantalum conduction electrons. Strong correlations are evident in flat bands near the Fermi level, a high heat capacity coefficient, and a high Kadowaki-Woods ratio. Density functional theory calculations suggest that electron and hole Fermi surfaces in the ferromagnetic phase are associated with opposite spin polarization. When the system is diluted of magnetic ions, the samples become superconducting below about 2 Kelvin. We discuss possible mechanisms for superconductivity in this family.

30 DIRECT ENERGY CONVERSION↗

Forwardscattering corrections for optical extinction measurements in aerosol media. II - Polydispersions

The paper presents a parametric study of the forwardscattering corrections for experimentally measured optical extinction coefficients in polydisperse particulate media, since some forward scattered light invariably enters, along with the direct beam, into the finite aperture of the detector. Forwardscattering corrections are computed by two methods: (1) using the exact Mie theory, and (2) the approximate Rayleigh diffraction formula for spherical particles. A parametric study of the dependence of the corrections on mode radii, real and imaginary parts of the complex refractive index, and half-angle of the detector's view cone has been carried out for three different size distribution functions of the modified gamma type. In addition, a study has been carried out to investigate the range of these parameters in which the approximate formulation is valid. The agreement is especially good for small-view cone angles and large particles, which improves significantly for slightly absorbing aerosol particles. Also discussed is the dependence of these corrections on the experimental design of the transmissometer systems.

Deepak, A.↗

A combined experimental and computational study on the reaction dynamics of the 1-propynyl (CH 3 CC, X 2 A 1 ) – propylene (CH 3 CHCH 2 , X 1 A') system: formation of 1,3-dimethylvinylacetylene (CH 3 CCCHCHCH 3 , X 1 A') under single collision conditions

The reaction of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with propylene (CH 3 CHCH 2 ; X 1 A') was studied in a crossed molecular beam machine at a collision energy of 37 ± 1 kJ mol –1 . Experimental data combined with high-level electronic structure (CCSD(T)-F12/cc-pVTZ-F12//ωB97X-D/6-311G(d,p)) and RRKM calculations reveal the reaction mechanism. Here, the overall barrierless and exoergic reaction involves indirect reaction dynamics and commences preferentially with addition of 1-propynyl with its radical centre to the carbon–carbon double bond at the terminal carbon atom of propylene. This work focuses on molecular mass growth process (hydrogen loss channels) although theory suggests methyl loss as a prevalent channel. In these processes, the C6H9 collision complexes either emit atomic hydrogen or undergo isomerisation followed by atomic hydrogen loss to preferentially yield the cis/trans isomers of 1,3-dimethylvinylacetylene (2-hexen-4-yne) as the primary product. Analysis of reaction dynamics of 1-propynyl and ethynyl radicals with propylene along with their fractional abundance in deep space suggests formation of methyl- and dimethyl derivatives of vinylacetylene in cold molecular clouds. Once formed they may engage in fundamental molecular mass growth processes via the barrierless Hydrogen Abstraction Vinylacetylene Addition mechanism that leads to the formation of methyl- and dimethylnaphthalenes thus providing a versatile route to methyl-substituted PAHs in interstellar medium.

1-propynyl↗

Assessment of DFT methods for the prediction of detachment energies and electronic structures of complex and multiply charged anions

Here, in this work, we assess the ability of different density functional theory methods to reproduce the vertical/adiabatic detachment energies (VDEs/ADEs) of a series of 40 anions, including monoanion versus multiply charged anions (MCAs) and covalent versus noncovalent interaction. The statistical analysis of errors is firstly performed by comparing the theoretical values with the experimental benchmark data obtained from the negative ion photoelectron spectroscopy (NIPES). The proposed optimally tuned range-separated (OTRS) functionals are proved to not only well reproduce the experimental VDEs/ADEs, but also simulate experimental NIPES. The result implies a more serious issue of delocalization error for MCAs and a balanced description of electron-rich/- deficient region of anions is necessary. The radius of the spherically symmetric average electron localization function (ELF) region is demonstrated as a useful descriptor for the characterization of electronic structure of anionic systems and a quasilinear fitting model is proposed to efficiently obtain their OTRS parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Neural Network Force Fields to Ascertain the Quality of Ab Initio Simulations of Liquid Water

Accurately simulating the properties of bulk water, despite the apparent simplicity of the molecule, is still a challenge. In order to fully understand and reproduce its complex phase diagram, it is necessary to perform simulations at the ab initio level, including quantum mechanical effects both for electrons and nuclei. This comes at a high computational cost, given that the structural and dynamical properties tend to require long timescales and large simulation cells. In this work, we evaluate the errors that density functional theory (DFT)-based simulations routinely incur into due time- and size-scale limitations. These errors are evaluated using neural-network-trained force fields that are accurate at the level of DFT methods. We compare different exchange and correlation potentials for properties of bulk water that require large timescales. Here we show that structural properties are less dependent on the system size and that dynamical properties such as the diffusion coefficient have a strong dependence on the simulation size and timescale. Our results facilitate comparisons of DFT-based simulation results with experiments and offer a path to discriminate between model and convergence errors in these simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling prebiotic chemistries with quantum accuracy at classical costs

Molecular Dynamics (MD) simulations using classical force-fields are commonly employed in numerous scientific investigations. However, many natural processes involve bond breaking and quantum forces. This complexity is compounded by the presence of multiple competing length and timescales. For example, accurately modeling the thermodynamics and dynamics of a chemical reaction requires accounting for the concerted movements of numerous solvent molecules and ions with their own fast or slow timescales. While widely used static Density Functional Theory (DFT) calculations at 0 temperature can be beneficial for such investigations, they do not account for dynamics, and lack precision in describing the molecular environments. They particularly fail at correct, rigorous treatments of finite-temperature fluctuations, and thus generalization to experimentally relevant conditions. In PNAS Benayad et al develop a scalable, generalizable approach for designing Neural Network Potentials (NNPs) that can handle chemical reactivity in solvated systems with quantum accuracy at classical costs. Specifically, they study phosphoester bond formation and rupture, which is fundamentally relevant to the Phosphorus-Oxygen bond formation central to life, and especially for the RNA world hypothesis. The framework developed here has the potential to generalize to different chemical reactions of energy and biological relevance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Modeling of Floating Wind Arrays Including Current Loads and Seabed Bathymetry

Capabilities for modeling the effects of seabed bathymetry and current drag loads on a floating wind farm are now available in an open-source model for quasi-static analysis. In this model, mooring lines and dynamic cables are represented by a quasi-static solver that can quickly represent complex mooring/cabling arrangements and arrays of floating bodies. To account for seabed bathymetry, we expand the model to include a surface mesh that captures changes in water depth over a rectangular grid. We formulate modifications to the catenary equations that capture a mooring line's profile and tensions when contacting a slope seabed. To account for current drag loads on mooring lines and dynamic cables, we formulate a novel technique that rotates the reference frame so that the vector sum of the the weight and the current force are used in the catenary equations, while accounting for the seabed orientation. To complete the system, current drag loads on floating substructures are handled by inclusion of strip-theory drag calculations. These new capabilities are verified by comparing with results from the established offshore dynamics models MoorDyn and OrcaFlex in equivalent steady-state scenarios. The results show very good agreement for both sloped seabeds and current loads. With computation times of the quasi-static model typically under one second, the model additions are a useful capability toward rapidly evaluating a floating wind array's response to environmental loads under realistic site conditions.

bathymetry↗

A Sheaf Theoretical Approach to Uncertainty Quantification of Heterogeneous Geolocation Information

Integration of multiple, heterogeneous sensors is a challenging problem across a range of applications. Prominent among these are multi-target tracking, where one must combine observations from different sensor types in a meaningful and efficient way to track multiple targets. Because different sensors have differing error models, we seek a theoretically justified quantification of the agreement among ensembles of sensors, both overall for a sensor collection, and also at a fine-grained level specifying pairwise and multi-way interactions among sensors. We demonstrate that the theory of mathematical sheaves provides a unified answer to this need, supporting both quantitative and qualitative data. Furthermore, the theory provides algorithms to globalize data across the network of deployed sensors, and to diagnose issues when the data do not globalize cleanly. We demonstrate and illustrate the utility of sheaf-based tracking models based on experimental data of a wild population of black bears in Asheville, North Carolina. A measurement model involving four sensors deployed among the bears and the team of scientists charged with tracking their location is deployed. This provides a sheaf-based integration model which is small enough to fully interpret, but of sufficient complexity to demonstrate the sheaf’s ability to recover a holistic picture of the locations and behaviors of both individual bears and the bear-human tracking system. A statistical approach was developed in parallel for comparison, a dynamic linear model which was estimated using a Kalman filter. This approach also recovered bear and human locations and sensor accuracies. When the observations are normalized into a common coordinate system, the structure of the dynamic linear observation model recapitulates the structure of the sheaf model, demonstrating the canonicity of the sheaf-based approach. However, when the observations are not so normalized, the sheaf model still remains valid.

97 MATHEMATICS AND COMPUTING↗

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass and momentum perturbations

An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 11-12 July 1981 CCOPE case study indicated two episodes of coherent waves. While geostrophic adjustment, shearing instability, and terrain were all implicated separately or in combination as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to define the genesis processes from observations alone. The first part of this paper, 3D Numerical Modeling Studies of Terrain-Induced Mass/Momentum Perturbations, employs a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed gravity wave episode. The meso-beta scale numerical model is used to study various simulations of the role of multiple geostrophic adjustment processes in focusing a region for gravity wave genesis. The second part of this paper, Linear Theory and Theoretical Modeling, investigates the response of non-resting rotating homogeneous and continuously stratified Boussinesq models of the terrestrial atmosphere to temporally impulsive and uniformly propagating three-dimensional localized zonal momentum sources representative of midlatitude jet streaks. The methods of linear perturbation theory applied to the potential vorticity (PV) and wave field equations are used to study the geostrophic adjustment dynamics. The total zonal and meridional wind perturbations are separated into geostrophic and ageostrophic components in order to define and follow the evolution of both the primary and secondary mesocirculations accompanying midlatitude jetogenesis forced by geostrophic adjustment processes. This problem is addressed to help fill the gap in understanding the dynamics and structure of mesoscale inertia-gravity waves forced by geostrophic adjustment processes in simple two-dimensional quiescent current systems and those produced by mesoscale numerical models simulating the orographic and diabatic perturbation of three-dimensional quasi-geostrophically balanced synoptic scale jet streaks associated with complex baroclinic severe storm producing environments.

Lin, Yuh-Lang↗

A new stochastic systems approach to structural integrity

This paper develops improved stochastic models for the description of a large variety of fatigue crack growth phenomena that occur in components of considerable importance to the functionality and reliability of complex engineering structures. In essence, the models are based on the McGill-Markov and Closure-Lognormal stochastic processes. Not only do these models have the capability of predicting the statistical dispersion of crack growth rates, they also, by incorporating the concept of crack closure, have the capability of transferring stochastic crack growth properties measured under ideal laboratory conditions to situations of industrial significance, such as those occurring under adverse loading and/or environmental conditions. The primary data required in order to be in a position to estimate the pertinent parameters of these stochastic models are obtained from a statistically significant number of replicate tests. In this paper, both the theory and the experimental technique are illustrated using a Ti-6Al-4V alloy. Finally, important structural integrity, reliability, availability and maintainability concepts are developed and illustrated.

Provan, James W.↗

Development of a Comprehensive Physics-Based Model for Study of NASA Gateway Lunar Dust Contamination

NASA is committed to landing the first woman and first person of color on the Moon by 2025 to begin a sustaining presence. A major component of this mission is NASA’s orbiting lunar outpost, Gateway, that will be subjected to the harsh environment of space during its 7-day orbit around the Moon. One aspect of this environment that is not well quantified is microscopic lunar regolith particles, or simply, lunar dust. As Apollo 17 Mission Commander Eugene Cernan said, “I think dust is probably one of our greatest inhibitors to a nominal operation on the moon. I think we can overcome other physiological or physical or mechanical problems except dust” [1]. These dust particles, most of which are smaller than the width of two human hairs and of highly irregular shape, are composed mainly of dielectric materials. Such characteristics can allow the dust particles to collect electric charge from the surrounding environment, resulting in electrostatic interactions between the dust, electromagnetic fields, and electrically charged surfaces [2-3]. Due to the nature of these fields, dust can collect on and contaminate charged surfaces. Lunar dust introduced into the Gateway environment by a lunar ascent vehicle element returning from the surface of the Moon presents risks to Gateway hardware such as radiators, solar arrays, antennae, and docking mechanisms. To quantify this risk and inform NASA, International Partner, and Commercial Provider stakeholders, Booz Allen, in partnership with NASA, is developing a toolset to model the interaction of charged lunar dust particles with the cis-lunar and deep-space environments. This includes a comprehensive physics model of the space plasma and solar radiation environment and electromagnetic interaction of spacecraft and lunar dust. Spacecraft charging and plasma environment are solved using the open-source code from The French Aerospace Lab (ONERA) and the European Space Agency (ESA), known as Spacecraft Plasma Interaction Software (SPIS) [4]. The physics of particle charging and motion is handled by a time-dependent implementation of Orbital-Motion-Limited (OML) theory which accounts for plasma flows and positive potential grains [5-7], using Siemens STAR-CCM+ as the framework. Model validation includes lab experiments by NASA experts and academia, as well as future on-orbit dust detecting payloads on the exterior of Gateway and Commercial Lunar Payload Services (CLPS) missions to the lunar surface [7-9]. The results and analyses will inform Gateway Program system owners at risk for lunar dust contamination.

Complex Plasma↗

Efficient Modeling of Floating Wind Arrays Including Current Loads and Seabed Bathymetry: Preprint

Capabilities for modeling the effects of seabed bathymetry and current drag loads on a floating wind farm are now available in an open-source model for quasi-static analysis. In this model, mooring lines and dynamic cables are represented by a quasistatic solver that can quickly represent complex mooring/cabling arrangements and arrays of floating bodies. To account for seabed bathymetry, we expand the model to include a surface mesh that captures changes in water depth over a rectangular grid. We formulate modifications to the catenary equations that capture a mooring line's profile and tensions when contacting a slope seabed. To account for current drag loads on mooring lines and dynamic cables, we formulate a novel technique that rotates the reference frame so that the vector sum of the the weight and the current force are used in the catenary equations, while accounting for the seabed orientation. To complete the system, current drag loads on floating substructures are handled by inclusion of strip-theory drag calculations. These new capabilities are verified by comparing with results from the established offshore dynamics models MoorDyn and OrcaFlex in equivalent steady-state scenarios. The results show very good agreement for both sloped seabeds and current loads. With computation times of the quasi-static model typically under one second, the model additions are a useful capability toward rapidly evaluating a floating wind array's response to environmental loads under realistic site conditions.

bathymetry↗

Controlling Vibronic Coupling in Chlorophyll Proteins: The Effects of Excitonic Delocalization and Vibrational Localization

Vibrational-electronic (vibronic) coupling plays a critical role in excitation energy transfer in molecular aggregates and pigment-protein complexes (PPCs). But the interplay between excitonic delocalization and vibronic interactions is complex, often leaving even qualitative questions as to what conceptual framework (e.g., Redfield versus Forster theory) should be used to interpret experimental results. To shed light on this issue, we report here on the interplay between excitonic delocalization and vibronic coupling in site-directed mutants of the water soluble chlorophyll protein (WSCP), as reflected in 77 K fluorescence spectra. Experimentally, we find that in PPCs where excitonic delocalization is disrupted (either by mutagenesis or heterodimer formation), the relative intensity of the vibrational sideband (VSB) in fluorescence spectra is suppressed by up to 37% compared to the native protein. Furthermore, numerical simulations reveal that this effect results from the localization of high-frequency vibrations in the coupled system; while excitonic delocalization suppresses the purely electronic transition due to H-aggregate-like dipole-dipole interference, high-frequency vibrations are unaffected, leading to a relative enhancement of the VSB. By comparing VSB intensities of PPCs both in the presence and absence of excitonic delocalization, we extract a set of “local” Huang-Rhys (HR) factors for Chl a in WSCP. More generally, our results suggest a significant role for geometric effects in controlling energy-transfer rates (which depend sensitively on absorption/fluorescence line shapes) in molecular aggregates and PPCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stereomolecular interactions and microsystems in experimental protobiogenesis

A theory on the origin of the first reproducing (in the sense of Fox, 1973) cells on earth is presented. A degree of self-ordering has been observed in the selective thermal reactions of certain ones of the amino acids which have been found in samples of meteorites, lunar samples, and terrestrial lava, and which have been synthesized from formaldehyde and ammonia under conditions modeling those assumed to have been present in the primitive earth. The genetic code is supposed to have arisen from a nonrandom matrix of preproteins. The ordering of preproteins specified limitations in the identity of the pre-enzymes and primordial microsystems which gave rise to contemporary living systems. The evolution of progresively more complex system was largely determined by stereochemical forces involving the shapes of the molecules themselves and the regulatory effects of the confines of the microsystems.

Fox, S. W.↗

Approximate effect of parameter pseudonoise intensity on rate of convergence for EKF parameter estimators

When using parameter estimation methods based on extended Kalman filter (EKF) theory, it is common practice to assume that the unknown parameter values behave like a random process, such as a random walk, in order to guarantee their identifiability by the filter. The present work is the result of an ongoing effort to quantitatively describe the effect that the assumption of a fictitious noise (called pseudonoise) driving the unknown parameter values has on the parameter estimate convergence rate in filter-based parameter estimators. The initial approach is to examine a first-order system described by one state variable with one parameter to be estimated. The intent is to derive analytical results for this simple system that might offer insight into the effect of the pseudonoise assumption for more complex systems. Such results would make it possible to predict the estimator error convergence behavior as a function of the assumed pseudonoise intensity, and this leads to the natural application of the results to the design of filter-based parameter estimators. The results obtained show that the analytical description of the convergence behavior is very difficult.

Hill, Bryon K.↗

Soliton formation and topology manipulation of coupled spins via ultrafast re-magnetization

The major goal of the project was to explore the properties of magnetic order driven far out of equilibrium by optical excitations. These include fundamental questions related to the interplay of magnetic, structural and electronic degrees of freedom in materials that are optically excited. These problems break down in short term quenching of the magnetization as a results of the energy of the optical pulse studying how the energy of the optical pulse and the angular momentum of the magnetization flow between the different degrees of freedom. This is followed by the longer time re-emergence of the magnetism as the system cools. Within this goal we explored the formation of solitons, both with and without non-trivial topology, via rapid re-magnetization processes after optical-driven ultrafast demagnetization. Guided by theory, we predict that turbulence and modulational instabilities will drive the formation of solitons, including dispersive shock waves, magnon droplets, and skyrmions. In the case of topological defect formation, i.e. skyrmion generation, related processes have long been predicted based on general principles of symmetry-breaking phase transitions; the density of topological defects in a long-range-ordered phase can be controlled by varying the quench rate through the second-order phase transition, i.e. the Kibble-Zurek (KZ) mechanism. We, more broadly, had the goal to study a range of other emergent magnetic behaviors after optical excitation. While great attention has been given to ultrafast demagnetization, less is known of the subsequent spin dynamics and coupling to the lattice. What is increasingly appreciated is that ultrafast demagnetization leads to spin currents that can carry angular momentum from the rapidly demagnetized sample. These spin currents appear in many ways but fundamentally controls demagnetization, drives interactions between regions of the material, and control the spin resulting structure. They can be probed many ways including THz emission which, in turn, gives insight into the demagnetization processes. Simultaneously if there is magneto-elastic coupling ultrafast demagnetization can drive structural excitations that are not expected from thermal energy added to the lattice. We explored a range of complex phenomena that arise from photo-excitation of magnetic systems with the coupling of electronic, magnetic and structural degrees of freedom. We have made a number of fundamental discoveries that are reflected in our publications list with additional work still being prepared for publication. Highlights of this work include (i) ultra-efficient, nonlinear THz surface acoustics (ii) spin-current-mediated rapid magnon localization and coalescence, (iii) spin-wave soliton formation in ferromagnetic FePt nanoparticles, (iv) dynamic phonon coupling in elemental antiferromagnetic Cr, (v) THz emission from Co/Pt bilayers and FeRh/Pt bilayers, (vi) theoretical investigation of spin hydrodynamics, solitons and shock waves, and (vii) ultrafast perturbation of magnetic domains by optical pumping.

36 MATERIALS SCIENCE↗

Spectral Properties and Coding Transitions of Haar-Random Quantum Codes

A quantum error-correcting code with a nonzero error threshold undergoes a mixed-state phase transition when the error rate reaches that threshold. We explore this phase transition for Haar-random quantum codes, in which the logical information is encoded in a random subspace of the physical Hilbert space. We focus on the spectrum of the encoded system density matrix as a function of the rate of uncorrelated, single-qudit errors. For low error rates, this spectrum consists of well-separated bands, representing errors of different weights. As the error rate increases, the bands for high-weight errors merge. The evolution of these bands with increasing error rate is well described by a simple analytic ansatz. Using this ansatz, as well as an explicit calculation, we show that the threshold for Haar-random quantum codes saturates the hashing bound, and thus coincides with that for random stabilizer codes. For error rates that exceed the hashing bound, typical errors are uncorrectable, but postselected error correction remains possible until a much higher detection threshold. Postselection can in principle be implemented by projecting onto subspaces corresponding to low-weight errors, which remain correctable past the hashing bound.

decoherence↗