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.
The HDIS/COUPi discrete element method modeling system was used to simulate the interaction between various debris and the Research Debris Diversion Platform (RDDP)
Viscously driven magnetospheric models, as well as a model involving interconnection between the geomagnetic field and the magnetic field in the solar wind, have been proposed to describe the driving force for magnetospheric convection. Lack of a satisfactory theory for the interconnection in the latter model and, in the case of the viscous interaction models, inadequacies in predicting the quantity of the driving force, make these two classes of models less than successful. Accordingly, a mechanically driven magnetospheric model is proposed: solar wind plasma enters the magnetosphere around the neutral points, covers the inner surface of the magnetopause and subsequently expands, driving convection as it escapes from the open tail.
The recently obtained spin-crossover nanoparticles are possible candidates for applications in the recording media industry as materials for data storage, or as pressure and temperature sensors. For these applications, the intermolecular interactions and interactions between spin-crossover nanoparticles are extremely important, as they may be essential factors in triggering the transition between the two stable phases: the high-spin and low-spin ones. In order to find correlations between the distributions in size and interactions and the transition temperatures distribution, we apply the FORC (First Order Reversal Curves) method, using simulations based on a mechanoelastic model applied to 2D triangular lattices composed of molecules linked by springs and embedded in a surfactant. We consider two Gaussian distributions: one is the size of the nanoparticles and another is the elastic interactions between edge spin-crossover molecules and the surfactant molecules. In order to disentangle the kinetic and non-kinetic parts of the FORC distributions, we compare the results obtained for different temperature sweeping rates. We also show that the presence of few larger particles in a distribution centered around much smaller particles dramatically increases the hysteresis width.
Another concern for the vehicle during its design trajectory was the separation of the first stage solid rocket booster from the upper stage component after it had depleted its solid fuel propellant. There has been some concern about the interstage of the first stage from clearing the nozzle of the J2-X engine. A detailed separation aerodynamic wind tunnel investigation was conducted in the AEDC VKF Tunnel A to help to investigate the interaction aerodynamic effects5. A comparison of the separation plane details between the Ares I architecture and the Ares I-X demonstration flight architecture is shown in figure 12. The Ares I design requires a more complex separation sequence and requires better control in order to avoid contact with the nozzle of the upper stage engine. The interstage, which houses the J2-X engine for the Ares I vehicle, must be able to separate cleanly to avoid contact of the J2-X engine. There is only about approximately 18 inches of buffer inside the interstage on each size of the nozzle so this is a challenging controlled separation event. This complex experimental investigation required two separate Ares I models (upper stage and first stage with interstage attached) with independent strain gauge balances installed in each model. It also required the Captive Trajectory System (CTS) that was needed to precisely locate the components in space relative to each other to fill out the planned test matrix. The model setup in the AEDC VKF Tunnel A is shown in figure 13. The CTS remotely positioned the first stage at the required x, y, and z positions and was able to provide interactions within 0.2" of the upper stage. A sample of the axial force on the first stage booster is shown in figure 14. These results, as a function of separation distance between the two stages, are compared to pre-test CFD results. Since this is a very challenging, highly unsteady flow field for CFD to correctly model, the experimental results have been utilized by GN&C discipline to more accurately represent the interaction aerodynamics. In addition to the integrated forces and moments obtained from the test, flow visualization data was obtained from this test in the form of Schlieren photographs, as shown in figure 15, which show the shock structure and interaction effects after the two stages separate during flight. This separation test was crucial in the successful flight test of the Ares I-X vehicle and provided the GN&C discipline with the unpowered proximity aerodynamic effect for a separation of the Ares I vehicle.
Fe 5-x GeTe 2 is a centrosymmetric, layered van der Waals (vdW) ferromagnet that displays Curie temperatures T c (270–330 K) that are within the useful range for spintronic applications. However, little is known about the interplay between its topological spin textures (e.g., merons, skyrmions) with technologically relevant transport properties such as the topological Hall effect (THE) or topological thermal transport. Here, via high-resolution Lorentz transmission electron microscopy, it is shown that merons and anti-meron pairs coexist with Néel skyrmions in Fe 5-x GeTe 2 over a wide range of temperatures and probe their effects on thermal and electrical transport. A THE is detected, even at room T, that senses merons at higher T's, as well as their coexistence with skyrmions as T is lowered, indicating an on-demand thermally driven formation of either type of spin texture. Remarkably, an unconventional THE is also observed in absence of Lorentz force, and it is attributed to the interaction between charge carriers and magnetic field-induced chiral spin textures. Here these results expose Fe 5-x GeTe 2 as a promising candidate for the development of applications in skyrmionics/meronics due to the interplay between distinct but coexisting topological magnetic textures and unconventional transport of charge/heat carriers.
Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.
Nonlocal modeling has come a long way. Researchers in the continuum mechanics and computational mechanics communities increasingly recognize that nonlocality is critical in realistic mathematical models of many aspects of the physical world. Physical interaction over a finite distance is fundamental at the atomic and nanoscale level, in which atoms and molecules interact through multibody potentials. Long-range forces partially determine the mechanics of surfaces and the behavior of dissolved molecules and suspended particles in a fluid. Nonlocality is therefore a vital feature of any continuum model that represents these physical systems at small length scales.
Amyloid β precursor protein binding family B member 1(APBB1) was first identified as a binding partner of amyloid precursor protein during brain development, but its function in the context of cancer remain unclear. Here we show for the first time that APBB1 is partly associated with intensifying cancer stem cell(CSC) and epithelial-to-mesenchymal transition (EMT) and enhancing radiation-resistant properties of lung cancer cells. We found that APBB1 was highly expressed in ALDH1{sup high} CSC-like cells sorted from A549 lung cancer cells. In APBB1-deficient H460 cells with forced overexpression of APBB1, the protein directly interacted with IGF1Rβ, enhanced phosphorylation of IGF1Rβ/PI3K/AKT pathway(activation) and subsequently induced the phosphorylation of GSK3β(inactivation). This phosphorylation stabilized Snail1, a negative regulator of E-cadherin expression, and regulated β-catenin–mediated ALDH1 expression, which are representative markers for EMT and CSCs, respectively. In contrast, suppression of APBB1 expression with siRNA yielded the opposite effects in APBB1-rich A549 cells. We concluded that APBB1 partly regulates the expression of ALDH1. We also found that APBB1 regulates activation of nuclear factor-κB, which is involved in reducing various stresses including oxidative stress, which suggests that APBB1 is associated with γ-radiation sensitivity. Our findings imply that APBB1 plays an important role in the maintenance of EMT-associated CSC-like properties and γ-radiation resistance via activation of IGF1Rβ/AKT/GSK3β pathway in lung cancer cells, highlighting APBB1 as a potential target for therapeutic cancer treatment. - Highlights: • APBB1 was highly expressed in ALDH1{sup high} cells sorted from A549 cells. • APBB1 is associated with intensifying CSC and EMT and enhancing radiation-resistant properties of lung cancer cells. • APBB1 regulates EMT-associated CSC-like properties and γ-radiation resistance via activation of IGF1Rβ/AKT/GSK3β pathway. • APBB1 partly regulates the expression of ALDH1 via AKT/GSK3β-associated β-catenin pathway.
The pressing need to address climate change has led to significant advancements in carbon dioxide (CO 2 ) capture technologies. Notably, facilitated transport membranes (FTMs) are distinguished by their exceptional selectivity and permeance, attributed to their reversible chemical reactions with CO 2 . This study, for the first time, sheds light on the reactive diffusion mechanism of CO 2 in FTMs, utilizing 1,1,3,3-tetramethylguanidine (TMG) as a mobile carrier. Specifically, state-ofthe- art molecular dynamics (MD) simulations, augmented by a reparameterized reactive force field (ReaxFF) capable of describing atomistic interactions and reaction pathways, are conducted to investigate the transport of CO 2 in TMG. The analysis of mean squared displacement (MSD) and diffusion coefficients reveals a clear hierarchy in the mobility of reaction components. Our findings highlight a unique hopping diffusion mechanism between bicarbonate ions and TMG molecules, increasing the diffusivity of reacted CO 2 by 1.4 times. The hopping events observed not only enhance our understanding of molecular mobility but also offer a means to boost the performance of FTMs in CO 2 capture applications. Overall, this research lays the groundwork for the future design of FTMs with optimal carrier properties.
We report that nanoporous materials such as metal–organic frameworks (MOFs) have attractive properties for selective capture of chemical warfare agents (CWAs). For obvious reasons, most research on adsorption of CWAs is performed with simulant molecules rather than real agents. This paper examines how effectively common CWA simulants mimic the adsorption properties of sarin and soman. To this end, we perform molecular simulations in the dilute adsorption limit for four simulants [dimethyl methylphosphonate (DMMP), diethyl chlorophosphate (DCP), diisopropyl fluorophosphate, and dimethyl p-nitrophenyl phosphate (DMNP)] and sarin and soman in a set of 2969 MOFs with experimentally known crystal structures. To establish the robustness of the conclusions with respect to the force field used in these simulations, each system was examined with two independent force fields, a “generic” force field and a density functional theory (DFT)-derived force field we established based on extensive dispersion-corrected DFT calculations of adsorption in the well-known MOF UiO-66. Our results show that when judging the performance of adsorbents using the heat of adsorption, DCP and DMMP are the best simulants for the adsorption of sarin, while DMNP is the best simulant for soman. The adsorption properties of DCP or DMMP show a strong correlation with sarin over a range of MOFs, but the correlation between DMNP and soman is considerably weaker. Finally, comparisons of results with both force fields indicate that our main conclusions are robust with respect to the force field used to define adsorbate–MOF interactions.
Atomistic-scale insights into the growth of a continuous, atomically thin hexagonal boron nitride (hBN) lattice from elemental boron and nitrogen on Ni substrates were obtained from multiscale modeling combining density functional theory (DFT) and reactive molecular dynamics. The quantum mechanical calculations focused on the adsorption and reaction energetics for the hBN building-block species, i.e., atomic B, N, BxNy (x, y = 1, 2), on Ni(111) and Ni(211), and the diffusion pathways of elemental B and N on these slab model surfaces and in the sublayer. B can diffuse competitively on both the surface and in the sublayer, while N diffuses strictly on the substrate surface. The DFT data were then used to generate a classical description of the Ni–B and Ni–N pair interactions within the formulation of the reactive force field, ReaxFF. Using the potential developed from this work, the elementary nucleation and growth process of an hBN monolayer structure from elemental B and N is shown at the atomistic scale. The nucleation initiates from the growth of linear BN chains, which evolve into branched and then hexagonal lattices. Subsequent DFT calculations confirmed the structure evolution energetically and validate the self-consistency of this multiscale modeling framework. On the basis of this framework, the fundamental aspects regarding crystal quality and the role of temperature and substrates used during hBN growth can also be understood.
Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.
Abstract The Coulomb-free 1 S 0 proton-proton ( p - p ) scattering length relies heavily on numerous and distinct theoretical techniques to remove the Coulomb contribution. Here, it has been determined from the half-off-the-energy-shell p - p scattering cross section measured at center-of-mass energies below 1 MeV using the quasi-free p + d → p + p + n reaction. A Bayesian data-fitting approach using the expression of the s-wave nucleon-nucleon scattering cross section returned a p - p scattering length $${a}_{pp}=-18.1{7}_{-0.58}^{+0.52}{| }_{stat}\pm 0.0{1}_{syst}$$ a p p = − 18.1 7 − 0.58 + 0.52 ∣ s t a t ± 0.0 1 s y s t fm and effective range r 0 = 2.80 ± 0.05 s t a t ± 0.001 s y s t fm. A model based on universality concepts has been developed to interpret this result. It accounts for the short-range interaction as a whole, nuclear and residual electromagnetic, according to what the s-wave phase-shift δ does in the description of low-energy nucleon-nucleon scattering data. We conclude that our parameters are representative of the short-range physics and propose to assess the charge symmetry breaking of the short-range interaction instead of just the nuclear interaction. This is consistent with the current understanding that the charge dependence of nuclear forces is due to different masses of up-down quarks and their electromagnetic interactions. This achievement suggests that these properties have a lesser than expected impact in the context of the charge symmetry breaking.
The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of Δ G ele , Δ G vDW , Δ G polar , Δ G non-polar , T Δ S total , Δ G deform , Δ G con , and Δ G ion . To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (Δ G deform ), the energy cost from translational and rotational entropic contributions ( T Δ S tran+rot ), the total electrostatic interactions (Δ G polar-PB/GB + Δ G ele ) of incorporation, the intramolecular electrostatic interactions (Δ G ele ) and electrostatic polar solvation interactions (Δ G polar-PB/GB ) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (Δ G vDW ), the non-polar solvation interaction (Δ G non-polar ), the vibrational entropic contribution ( T Δ S vib ), and the standard concentration dependent free energies of DOX (Δ G con ) and the ionic solution (Δ G ion ). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX–DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. Here, we have calculated the total binding free energy (BFE) (Δ G t-sim ) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX–DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX–DNA1 and DOX–DNA2 in the 1 : 1 complex is –9.1 and –5.1 kcal mol –1 respectively. The same quantities from the modified MM-PBSA method are –12.74 and –8.35 kcal mol –1 respectively. The value of the total BFE Δ G t-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of –7.7 ± 0.3 kcal mol –1 .
Silicates are the most abundant materials in the earth's crust. In recent years, two-dimensional (2D) versions of them grown on metal supports (known as bilayer silicates) have allowed their study in detail down to the atomic scale. These structures are self-containing. Furthermore, they are not covalently bound to the metal support but interact with it through van der Waals forces. Like their three-dimensional counterparts, the 2D-silicates can form both crystalline and vitreous structures. Furthermore, the interconversion between vitreous to crystalline structures has been experimentally observed at the nanoscale. While theoretical work has been carried out to try to understand these transformations, a limitation for ab initio methods, and even molecular dynamics methods, is the computational cost of studying large systems and long timescales. In this work, we present a simple and computationally inexpensive approach, that can be used to represent the evolution of bilayer silicates using a bond-switching algorithm. This approach allows reaching equilibrium ring size distributions as a function of a parameter that can be related to the ratio between temperature and the energy required for the bond-switching event. The ring size distributions are compared to experimental data available in the literature.
Here, we numerically examine a system of monodisperse self-propelled particles interacting with each other via simple steric forces and aligning torques moving through a periodic array of obstacles. Without obstacles, this system shows a transition to a polarized or aligned state for critical alignment parameters. In the presence of obstacles, there is still a polarization transition, but for dense enough arrays, the polarization is locked to the symmetry directions of the substrate. When the obstacle array is made anisotropic, at low densities the particles can form a quasi-isotropic state where the system can be polarized in any of the dominant symmetry directions. For intermediate anisotropy, the particles self-organize into a coherent lane state with one-dimensional polarization. In this phase, a small number of highly packed lanes are adjacent to less dense lanes that have the same polarization, but lanes further away can have the opposite polarization, so that global polarization is lost. For the highest anisotropy, hopping between lanes is suppressed, and the system forms uniformly dense uncoupled but polarized lanes.
The mechanism and quality of the welding of single-crystal (SC) and amorphous gold nanowires (NWs) with head-to-head contact are studied using molecular dynamics simulations based on the second-moment approximation of the many-body tight-binding potential. The results are discussed in terms of atomic trajectories, slip vectors, stress, and radial distribution function. Simulation results show that the alignment for the amorphous NWs during welding is easier than that for the SC NWs due to the former's relatively stable geometry. A few dislocations nucleate and propagate on the (111) close-packed plane (slip plane) inside the SC NWs during the welding and stretching processes. During welding, an incomplete jointing area first forms through the interactions of the van der Waals attractive force, and the jointing area increases with increasing extent of contact between the two NWs. A crystallization transition region forms in the jointing area for the welding of SC-amorphous or amorphous-SC NWs. With increasing interference, an amorphous gold NW shortens more than does a SC gold NW due to the former's relatively poor strength. The pressure required for welding decreases with increasing temperature.
The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. Finally, the patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.