SEARCH · Engineering Papers
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.
A preliminary investigation into the possible role of glass transition and glassy state behavior of lubricants in EHD contacts is reported. Measurements of the glass transition of lubricants as a function of pressure by two methods are presented along with a discussion indicating possible implications of the results to EHD lubrication.
Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.
Zonation due to alpha-decay damage in a natural single crystal of zircon from Sri Lanka is discussed. The zones vary in thickness on a scale from one to hundreds of microns. The uranium and thorium concentrations vary from zone to zone such that the alpha decay dose is between 0.2 x 10 to the 16th and 0.8 x 10 to the 16th alpha-events per milligram. The transition from the crystalline to the aperiodic metamict state occurs over this dose range. At doses greater than 0.8 x 10 to the 16th alpha events/mg there is no evidence for long-range order. This type of damage will accumulate in actinide-bearing, ceramic nuclear waste forms. The systematic pattern of fractures would occur in crystalline phases that are zoned with respect to actinide radionuclides.
We report an experimental investigation on the effect of ground-state charge transfer and its induced electric field on photoexcited charge transfer in van der Waals heterostructures. Two heterostructure samples were fabricated by stacking an undoped WSe 2 monolayer with either a Nb-doped or undoped MoSe 2 monolayer. While no ground-state charge transfer is expected in the MoSe 2 /WSe 2 heterostructure, the doped holes in the MoSe 2 :Nb/WSe 2 heterostructure can transfer to WSe 2 , creating a space-charge electric field. By comparing the photoluminescence and time-resolved differential reflectance of the two heterostructures, we find that photoexcited hole transfer from MoSe 2 to WSe 2 is largely blocked by this field, whereas photoexcited electron transfer from WSe 2 to MoSe 2 is less affected. Furthermore, these results provide insight into the impact of doping on the charge-transfer performance of van der Waals heterostructures.
A systematic structural investigation of Fe 2 Mo 3 O 8 as a function of pressure, temperature, and magnetic field reveals that the 𝑃6 3 𝑚𝑐 space group remains robust over a wide range of conditions. No changes in the long-range crystal structure are observed for pressures up to 10 GPa, temperatures between 11 K and 300 K, and magnetic fields up to 9 T. The magnetostructural response, quantified by Δ𝑐/𝑐, is determined for magnetic fields applied transverse to the crystallographic 𝑐 axis, demonstrating strong magnetoelastic coupling. The well-known magnetic-field-induced transition is confirmed to be first order and isostructural, occurring between two distinct altermagnetic states. Importantly, this transition can also be accessed using magnetic fields applied within the 𝑎𝑏 plane. When expressed in terms of the 𝑐/𝑎 ratio, the structural evolution under pressure (0–10 GPa) maps continuously onto that observed across the full Zn substitution range in (Fe 1−𝑦 Zn 𝑦 )2 Mo 3 O 8 (0 ≤ 𝑦 ≤ 1), indicating a common underlying structural control parameter. Furthermore, these results demonstrate that the 𝑐-axis lattice parameter is the primary structural degree of freedom governing magnetic behavior under external tuning parameters. More broadly, they indicate that magnetic order in this family of complex oxides (𝐴 2 Mo 3 O 8 , 𝐴 = Co,Mn,Ni) can be systematically tuned via pressure. This tunability highlights these materials as promising platforms for controllable magnetic switching, particularly in thin films integrated with piezoelectric substrates.
We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.
The Soil Moisture Active Passive (SMAP) observatory was launched January 31, 2015, and its L-band radiometer and radar instruments became operational during April 2015. This paper provides a summary of the quality assessment of its baseline soil moisture and freeze/thaw products as well as an overview of new products. The first new product explores the Backus Gilbert optimum interpolation based on the oversampling characteristics of the SMAP radiometer. The second one investigates the disaggregation of the SMAP radiometer data using the European Space Agency's Sentinel-1 C-band synthetic aperture radar (SAR) data to obtain soil moisture products at about 1 to 3 km resolution. In addition, SMAPs L-band data have been found useful for many scientific applications, including depictions of water cycles, vegetation opacity, ocean surface salinity and hurricane ocean surface wind mapping. Highlights of these new applications will be provided.The SMAP soil moisture, freeze/taw state and SSSprovide a synergistic view of water cycle. For example, Fig.7 illustrates the transition of freeze/thaw state, change of soilmoisture near the pole and SSS in the Arctic Ocean fromApril to October in 2015 and 2016. In April, most parts ofAlaska, Canada, and Siberia remained frozen. Melt onsetstarted in May. Alaska, Canada, and a big part of Siberia havebecome thawed at the end of May; some freshwater dischargecould be found near the mouth of Mackenzie in 2016, but notin 2015. The soil moisture appeared to be higher in the Oband Yenisei river basins in Siberia in 2015. As a result,freshwater discharge was more widespread in the Kara Seanear the mouths of both rivers in June 2015 than in 2016. TheNorth America and Siberia have become completely thawedin July. After June, the freshwater discharge from other riversinto the Arctic, indicated by blue, also became visible. Thefreeze-up started in September and the high latitude regionsin North America and Eurasia became frozen. Comparing thespread of freshwater in August 2015 and 2016 suggests thatthere was more discharge from Ob and Yenisei in 2015,which appeared to correspond to a higher soil moisturecontent in the Ob and Yenisei basins. In contrast, Mackenzieappeared to have more discharge in September 2016.
When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Lastly, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.
A recent work [Balram, Jain, and Barkeshli, Phys. Rev. Res. 2, 013349 (2020)] has suggested that an unconventional state describing Zn superconductivity of composite bosons, which supports excitations with charge 1/(3n) of the electron charge, is energetically better than the Laughlin wave function at \nu=7/3 in GaAs systems. All experiments to date, however, are consistent with the latter. To address this discrepancy, we study the effect of finite width on the ground state and predict a phase transition from an unconventional Zn state at small widths to the Laughlin state for widths exceeding ~1.5 magnetic lengths. Here, we also determine the parameter region where an unconventional state is stabilized in the one-third filled zeroth Landau level in bilayer graphene. The roles of Landau level mixing and spin are also considered.
In α-RuCl 3 , an external magnetic field applied within the honeycomb plane can induce a transition from a magnetically ordered state to a disordered state that is potentially related to the Kitaev quantum spin liquid. In zero field, single crystals with minimal stacking faults display a low-temperature state with in-plane zigzag antiferromagnetic order and a three-layer periodicity in the direction perpendicular to the honeycomb planes. In this work, we present angle-dependent magnetization, ac susceptibility, and thermal transport data that demonstrate the presence of an additional intermediate-field ordered state at fields below the transition to the disordered phase. Neutron-diffraction results show that the magnetic structure in this phase is characterized by a six-layer periodicity in the direction perpendicular to the honeycomb planes. Theoretically, the intermediate ordered phase can be accounted for by including spin-anisotropic couplings between the layers in a three-dimensional spin model. Together, this demonstrates the importance of interlayer exchange interactions in α-RuCl 3 .
The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024
Observations of the 1(01)-0(00) rotational transitions of A and E state acetaldehyde are reported. The transitions were detected, for the first time in interstellar space, in the cold dust clouds TMC-1 and L134N, and in Sgr B2. This is also the first time acetaldehyde has been found in a dust cloud and is the most complex oxygen-bearing molecule yet known in this environment. A column density of 6 x 10 to the 12th/sq cm in TMC-1, comparable to many other species detected there, and an approximately equal column density in L134N are formed. In the direction of Sgr B2, the CH3CHO profile appears to consist of broad emission features from the hot molecular cloud core, together with absorption features resulting from intervening colder material. The possible detection of HC9N toward IRC + 10 deg 216 through its J = 33-32 transition is also reported. Implications for cold dust cloud chemistry and excitation are discussed.
Activating metamagnetic transitions between ordered states in van der Waals magnets and devices bring great opportunities in spintronics. Here we show that external pressure, which enhances the inter-layer hopping without introducing chemical disorders, triggers multiple metamagnetic transitions upon cooling in the topological van der Waals magnets Mn(Bi 1–x Sb x ) 4 Te 7 , where the antiferro magnetic interlayer superexchange coupling competes with the ferromagnetic interlayer coupling mediated by the antisite Mn spins. The temperature-pressure phase diagrams reveal that while the ordering temperature from the paramagnetic to ordered states is almost pressure-independent, the metamagnetic transitions show non-trivial pressure and temperature dependence, even re-entrance. For these highly anisotropic magnets, we attribute the former to the ordering temperature being only weakly dependent on the intralayer parameters, the latter to the parametrically different pressure and temperature dependence of the two interlayer couplings. Our independent probing of these disparate magnetic interactions paves an avenue for efficient magnetic manipulations in van der Waals magnets.
A proposed highly accurate clock or oscillator would be based on the concept of an opto-electronic oscillator (OEO) stabilized to an atomic transition. Opto-electronic oscillators, which have been described in a number of prior NASA Tech Briefs articles, generate signals at frequencies in the gigahertz range characterized by high spectral purity but not by longterm stability or accuracy. On the other hand, the signals generated by previously developed atomic clocks are characterized by long-term stability and accuracy but not by spectral purity. The proposed atomic clock would provide high spectral purity plus long-term stability and accuracy a combination of characteristics needed to realize advanced developments in communications and navigation. In addition, it should be possible to miniaturize the proposed atomic clock. When a laser beam is modulated by a microwave signal and applied to a photodetector, the electrical output of the photodetector includes a component at the microwave frequency. In atomic clocks of a type known as Raman clocks or coherent-population-trapping (CPT) clocks, microwave outputs are obtained from laser beams modulated, in each case, to create two sidebands that differ in frequency by the amount of a hyperfine transition in the ground state of atoms of an element in vapor form in a cell. The combination of these sidebands produces a transparency in the population of a higher electronic level that can be reached from either of the two ground-state hyperfine levels by absorption of a photon. The beam is transmitted through the vapor to a photodetector. The components of light scattered or transmitted by the atoms in the two hyperfine levels mix in the photodetector and thereby give rise to a signal at the hyperfine- transition frequency. The proposed atomic clock would include an OEO and a rubidium- or cesium- vapor cell operating in the CPT/Raman regime (see figure). In the OEO portion of this atomic clock, as in a typical prior OEO, a laser beam would pass through an electro-optical modulator, the modulated beam would be fed into a fiber-optic delay line, and the delayed beam would be fed to a photodetector. The electrical output of the photodetector would be detected, amplified, filtered, and fed back to the microwave input port of the modulator. The laser would be chosen to have the same wavelength as that of the pertinent ground-state/higher-state transition of the atoms in the vapor. The modulator/ filter combination would be designed to operate at the microwave frequency of the hyperfine transition. Part of the laser beam would be tapped from the fiberoptic loop of the OEO and introduced into the vapor cell. After passing through the cell, this portion of the beam would be detected differentially with a tapped portion of the fiber-optically-delayed beam. The electrical output of the photodetector would be amplified and filtered in a loop that would control a DC bias applied to the modulator. In this manner, the long-term stability and accuracy of the atomic transition would be transferred to the OEO.
HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.
A steady-state photocapacitance (SSPC) setup directly connected to the beamline of a MeV ion implanter is utilized to study primary intrinsic defects in β – Ga 2 O 3 generated by He implantation at cryogenic temperatures (120 K). At low temperatures, the migration of defects is suppressed, and hence the generation of primary intrinsic defects is expected to prevail. SSPC measurements reveal defect-related optical transitions in halide vapor-phase epitaxy (HVPE) -grown β – Ga 2 O 3 thin films with onset energies at 1.3 ( T 1 ) , 1.7 ( T 2 ) , 1.9 ( T 3 ) , 2.6 ( T 4 ) , 3.7 ( T 5 ) , and 4.2 eV ( T 6 ) . T 2 , T 4 , T 5 , and T 6 were observed in as-received HVPE-grown β – Ga 2 O 3 thin films, whereby T 2 is only sporadically observed. The introduction rates for T 3 , T 4 , as well as T 6 indicate an origin related to primary intrinsic defects. Notably, T 1 and T 3 are only observed after He implantation at cryogenic temperatures. Hybrid-functional calculations were performed to estimate the optical absorption cross-section spectra for the gallium ( Ga i ) and oxygen ( O i ) interstitials as well as the corresponding vacancies ( V Ga and V O , respectively), and compared with the measured onsets for optical absorption found by SSPC measurements. Indeed, we propose T 3 to be associated with Ga i ( + / + 2 ) and/or V GaI ( - 3 / - 2 ) , while T 4 is suggested to be related to V O K ( 0 / + ) ( K = I , II , III ) and/or V GaII ( - 3 / - 2 ) . Additionally, several further charge-state transition levels associated with V GaI and V GaII may contribute to T 4 and T 6 . We further studied the kinetics of the defects created with He implantation by exposing the sample to room temperature. The kinetics observed for T 3 and T 4 further support the proposed assignments of the corresponding defect signatures.
Understanding protein folding is crucial for protein sciences. The conformational spaces and energy landscapes of cold (unfolded) protein states, as well as the associated transitions, are hardly explored. Furthermore, it is not known how structure relates to the cooperativity of cold transitions, if cold and heat unfolded states are thermodynamically similar, and if cold states play important roles for protein function. We created the cold unfolding 4-helix bundle DCUB1 with a de novo designed bipartite hydrophilic/hydrophobic core featuring a hydrogen bond network which extends across the bundle in order to study the relative importance of hydrophobic versus hydrophilic protein–water interactions for cold unfolding. Structural and thermodynamic characterization resulted in the discovery of a complex energy landscape for cold transitions, while the heat unfolded state is a random coil. Below ~0 °C, the core of DCUB1 disintegrates in a largely cooperative manner, while a near-native helical content is retained. The resulting cold core-unfolded state is compact and features extensive internal dynamics. Below -5 °C, two additional cold transitions are seen, that is, (i) the formation of a water-mediated, compact, and highly dynamic dimer, and (ii) the onset of cold helix unfolding decoupled from cold core unfolding. Our results suggest that cold unfolding is initiated by the intrusion of water into the hydrophilic core network and that cooperativity can be tuned by varying the number of core hydrogen bond networks. Protein design has proven to be invaluable to explore the energy landscapes of cold states and to robustly test related theories.
Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact electronic excitation of molecular hydrogen suitable for implementation in collisional-radiative models. The adiabatic-nuclei molecular convergent close-coupling method is used to calculate cross sections for excitation of all bound vibrational levels and dissociative excitation of the B 1 Σ u + , C 1 Π u , E F 1 Σ g + , B ′ 1 Σ u + , G K 1 Σ g + , I 1 Π g , J 1 Δ g , D 1 Π u , H 1 Σ g + , b 3 Σ u + , c 3 Π u , a 3 Σ g + , e 3 Σ u + , d 3 Π u , h 3 Σ g + , g 3 Σ g + , i 3 Π g , and j 3 Δ g electronic states from all –14 bound vibrational levels of the ground electronic ($\text{X}$ $^1Σ^+_g$) state. The data set consists of cross sections from threshold to 500 eV for over 5000 transitions, representing all possible electronic and vibrational transitions between the state and the –3 singlet and triplet states (where refers to the united-atoms-limit principle quantum number). The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files.