SEARCH · Search NASA
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.
Amyloids are insoluble and misfolded fibrous protein aggregates and associated with more than 20 serious human diseases. We perform all-atom molecular dynamics simulations of amyloid fibril assembly and disassembly.
Abstract not provided
Calculation of elastic constants requires initializing a lattice of atoms, an h matrix defining the computational cell boundaries, and the parameters defining the interatomic potentials, then numerically integrating the equations of motion for the atoms and the computational cell boundaries, calculating strains and stresses at each timestep, and updating at each timestep the average quantities required for calculation of fluctuations. Enforcements of the DYNAMO FORTRAN code included addition of documentation, improved memory management and data flow, and enabling use of several interatomic potentials. Documentation included headers for each subroutine, which document the parameters input to the routine, algorithms employed, and output.
Previous work on the cubic direct interaction approximation applied to the truncated (in Fourier space) cubically nonlinear Schroedinger equation model of Langmuir turbulence is extended to more modes. In the undriven, undamped case, excellent agreement between the statistical theory and a numerical ensemble of solutions of the dynamic equations is obtained. In the driven, damped case, satisfactory agreement is obtained provided the dynamic ensemble is limited to initial conditions in the basin of attraction.
Here, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method developed in our previous work [Smith et al., J. Chem. Phys. 142, 244112 (2015)] for the determination of the dynamic structure factor of liquid para-hydrogen and ortho-deuterium at state points of (T = 20.0 K, n = 21.24 nm -3 ) and (T = 23.0 K, n = 24.61 nm -3 ), respectively. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor reported by Smith et al. [J. Chem. Phys. 140, 034501 (2014)] for all momentum transfers considered. Moreover, this shows that FK-QCW provides a substantial improvement over the Feynman-Kleinert linearized path-integral method, in which purely classical dynamics are used. Furthermore, for small momentum transfers, it is shown that FK-QCW provides nearly the same results as ring-polymer molecular dynamics (RPMD), thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since it is not formally limited to correlation functions involving linear operators.
libEnsemble is a Python-based toolkit for running dynamic ensembles, developed as part of the DOE Exascale Computing Project. The toolkit utilizes a unique generator–simulator–allocator paradigm, where generators produce input for simulators, simulators evaluate those inputs, and allocators decide whether and when a simulator or generator should be called. The generator steers the ensemble based on simulation results. Generators may, for example, apply methods for numerical optimization, machine learning, or statistical calibration. libEnsemble communicates between a manager and workers. Flexibility is provided through multiple manager–worker communication substrates each of which has different benefits. These include Python’s multiprocessing, mpi4py, and TCP. Multisite ensembles are supported using Balsam or Globus Compute. We overview the unique characteristics of libEnsemble as well as current and potential interoperability with other packages in the workflow ecosystem. We highlight libEnsemble’s dynamic resource features: libEnsemble can detect system resources, such as available nodes, cores, and GPUs, and assign these in a portable way. These features allow users to specify the number of processors and GPUs required for each simulation; and resources will be automatically assigned on a wide range of systems, including Frontier, Aurora, and Perlmutter. Such ensembles can include multiple simulation types, some using GPUs and others using only CPUs, sharing nodes for maximum efficiency. We also describe the benefits of libEnsemble’s generator–simulator coupling, which easily exposes to the user the ability to cancel, and portably kill, running simulations based on models that are updated with intermediate simulation output. We demonstrate libEnsemble’s capabilities, scalability, and scientific impact via a Gaussian process surrogate training problem for the longitudinal density profile at the exit of a plasma accelerator stage. In conclusion, the study uses gpCAM for the surrogate model and employs either Wake-T or WarpX simulations, highlighting efficient use of resources that can easily extend to exascale.
The unfolding dynamics of ubiquitin were studied using a combination of x-ray solution scattering (XSS) and molecular dynamics (MD) simulations. The kinetic analysis of the XSS ubiquitin signals showed that the protein unfolds through a two-state process, independent of the presence of destabilizing salts. In order to characterize the ensemble of unfolded states in atomic detail, the experimental XSS results were used as a constraint in the MD simulations through the incorporation of x-ray scattering derived potential to drive the folded ubiquitin structure toward sampling unfolded states consistent with the XSS signals. We detail how biased MD simulations provide insight into unfolded states that are otherwise difficult to resolve and underscore how experimental XSS data can be combined with MD to efficiently sample structures away from the native state. Our results indicate that ubiquitin samples unfolded in states with a high degree of loss in secondary structure yet without a collapse to a molten globule or fully solvated extended chain. Finally, we propose how using biased-MD can significantly decrease the computational time and resources required to sample experimentally relevant nonequilibrium states.
A new feature-based forecasting framework of abnormalities is presented for early warning for cyber-physical control systems. The abnormalities may refer to intelligent cyber-attacks or naturally occurring faults and failures. Techniques presented here are aimed at protecting against unauthorized intrusions as well as fault prevention. Time series signals from system monitoring nodes are converted to features using feature discovery techniques. The feature behavior for each monitoring node is characterized in the form of decision boundaries, separating normal and abnormal space with operating data collected from the plant or by running virtual models of the plant. A set of ensemble state-space models are constructed for representing feature evolution in time-domain, where the ensembles are selected using Gaussian Mixture Model (GMM) clustering. The forecasted outputs are anticipated time-evolution of features, computed by applying an adaptive Kalman predictor to each ensemble model. The overall features forecast is then obtained through dynamic ensemble averaging. This is done by projecting evolution of feature vector to future times in a receding horizon fashion. The features forecast are compared to the decision boundary to estimate if/when the feature vectors will cross the boundary. Simulation results in a high fidelity GE gas turbine platform, show the effectiveness of our approach for forecasting abnormalities.
Abstract Time-resolved small-angle X-ray experiments are reported here that capture and quantify a previously unknown rapid collapse of the unfolded oligonucleotide as an early step in the folding of hybrid 1 and hybrid 2 telomeric G-quadruplex structures. The rapid collapse, initiated by a pH jump, is characterized by an exponential decrease in the radius of gyration from 24.3 to 12.6 Å. The collapse is monophasic and is complete in <600 ms. Additional hand-mixing pH-jump kinetic studies show that slower kinetic steps follow the collapse. The folded and unfolded states at equilibrium were further characterized by SAXS studies and other biophysical tools, showing that G4 unfolding was complete at alkaline pH, but not in LiCl solution as is often claimed. The SAXS Ensemble Optimization Method analysis reveals models of the unfolded state as a dynamic ensemble of flexible oligonucleotide chains with a variety of transient hairpin structures. These results suggest a G4 folding pathway in which a rapid collapse, analogous to molten globule formation seen in proteins, is followed by a confined conformational search within the collapsed particle to form the native contacts ultimately found in the stable folded form.
During the 42-year period (1979–2020) of satellite measurements, four major winter (December–March) polynyas have been observed north of Greenland: one in December 1986 and three in the last decade, i.e., February of 2011, 2017, and 2018. The 2018 polynya was unparalleled in its magnitude and duration compared to the three previous events. Given the apparent recent increase in the occurrence of these extreme events, this study aims to examine their evolution and causality, in terms of forced versus natural variability. The limited weather station and remotely sensed sea ice data are analyzed combining with output from the fully coupled Regional Arctic System Model (RASM), including one hindcast and two ensemble simulations. We found that neither the accompanying anomalous warm surface air intrusion nor the ocean below had an impact (i.e., no significant ice melting) on the evolution of the observed winter open-water episodes in the region. Instead, the extreme atmospheric wind forcing resulted in greater sea ice deformation and transport offshore, accounting for the majority of sea ice loss in all four polynyas. Our analysis suggests that strong southerly winds (i.e., northward wind with speeds greater than 10 m s –1 ) blowing persistently over the study region for at least 2 d or more were required over the study region to mechanically redistribute some of the thickest Arctic sea ice out of the region and thus to create open-water areas (i.e., a latent heat polynya). To assess the role of internal variability versus external forcing of such events, we carried out and examined results from the two RASM ensembles dynamically downscaled with output from the Community Earth System Model (CESM) Decadal Prediction Large Ensemble (DPLE) simulations. Out of 100 winters in each of the two ensembles (initialized 30 years apart: one in December 1985 and another in December 2015), 17 and 16 winter polynyas were produced north of Greenland, respectively. The frequency of polynya occurrence had no apparent sensitivity to the initial sea ice thickness in the study area pointing to internal variability of atmospheric forcing as a dominant cause of winter polynyas north of Greenland. We assert that dynamical downscaling using a high-resolution regional climate model offers a robust tool for process-level examination in space and time, synthesis with limited observations, and probabilistic forecasts of Arctic events, such as the ones being investigated here and elsewhere.
State and parameter estimation of power transmission networks is important for monitoring power grid operating conditions and analyzing transient stability. Wind power generation depends on fluctuating input power levels, which are correlated in time and contribute to uncertainty in turbine dynamical models. The ensemble Kalman filter (EnKF), a standard state estimation technique, uses a deterministic forecast and does not explicitly model time-correlated noise in parameters such as mechanical input power. However, this uncertainty affects the probability of fault-induced transient instability and increased prediction bias. Here a novel approach is to model input power noise with time-correlated stochastic fluctuations, and integrate them with the network dynamics during the forecast. While the EnKF has been used to calibrate constant parameters in turbine dynamical models, the calibration of a statistical model for a time-correlated parameter has not been investigated. In this study, twin experiments on a standard transmission network test case are used to validate our time-correlated noise model framework for state estimation of unsteady operating conditions and transient stability analysis, and a methodology is proposed for the inference of the mechanical input power time-correlation length parameter using time-series data from PMUs monitoring power dynamics at generator buses.
Ensembles of actuated colloids are excellent model systems to explore emergent out-of-equilibrium structures, complex collective dynamics, and design rules for the next generation materials. Here, we demonstrate that ferromagnetic microparticles suspended at an air-water interface and energized by an external rotating magnetic field spontaneously form dynamic ensembles of synchronized spinners in a certain range of the excitation field parameters. Each spinner generates strong hydrodynamic flows, and collective interactions of the multiple spinners promote a formation of dynamic lattices. On the basis of experiments and simulations, we reveal structural transitions from liquid to nearly crystalline states in this novel active spinner material and demonstrate that dynamic spinner lattices are reconfigurable, capable of self-healing behavior and that the transport of embedded inert cargo particles can be remotely tuned by the parameters of the external excitation field. Our findings provide insights into the behavior of active spinner materials with reconfigurable structural order and tunable functionalities.
Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.
Magnetic skyrmions are nanoscale spin textures touted as next-generation computing elements. When subjected to lateral currents, skyrmions move at considerable speeds. Their topological charge results in an additional transverse deflection known as the skyrmion Hall effect (SkHE). While promising, their dynamic phenomenology with current, skyrmion size, geometric effects and disorder remain to be established. Here we report on the ensemble dynamics of individual skyrmions forming dense arrays in Pt/Co/MgO wires by examining over 20,000 instances of motion across currents and fields. The skyrmion speed reaches 24 m/s in the plastic flow regime and is surprisingly robust to positional and size variations. Meanwhile, the SkHE saturates at ~22 ° , is substantially reshaped by the wire edge, and crucially increases weakly with skyrmion size. Particle model simulations suggest that the SkHE size dependence - contrary to analytical predictions - arises from the interplay of intrinsic and pinning-driven effects. These results establish a robust framework to harness SkHE and achieve high-throughput skyrmion motion in wire devices.
The conformational ensemble of intrinsically disordered proteins, such as α-synuclein, are responsible for their function and malfunction. Misfolding of α-synuclein can lead to neurodegenerative diseases, and the ability to study their conformations and those of other intrinsically disordered proteins under varying physiological conditions can be crucial to understanding and preventing pathologies. In contrast to well-folded peptides, a consensus feature of IDPs is their low hydropathy and high charge, which makes their conformations sensitive to pH perturbation. We examine a prominent member of this subset of IDPs, α-synuclein, using a divide-and-conquer scheme that provides enhanced sampling of IDP structural ensembles. We constructed conformational ensembles of α-synuclein under neutral (pH ~ 7) and low (pH ~ 3) pH conditions and compared our results with available information obtained from smFRET, SAXS, and NMR studies. Specifically, α-synuclein has been found to in a more compact state at low pH conditions and the structural changes observed are consistent with those from experiments. We also characterize the conformational and dynamic differences between these ensembles and discussed the implication on promoting pathogenic fibril formation. We find that under low pH conditions, neutralization of negatively charged residues leads to compaction of the C-terminal portion of α-synuclein while internal reorganization allows α-synuclein to maintain its overall end-to-end distance. Here, we also observe different levels of intra-protein interaction between three regions of α-synuclein at varying pH and a shift towards more hydrophilic interactions with decreasing pH.
Rejuvenation in particle filters is necessary to prevent the collapse of the weights when the number of particles is insufficient to properly sample the high-probability regions of the state space. Rejuvenation is often implemented in a heuristic manner by the addition of random noise that widens the support of the ensemble. This work aims at improving canonical rejuvenation methodology by the introduction of additional prior information obtained from climatological samples; the dynamical particles used for importance sampling are augmented with samples obtained from stochastic covariance shrinkage. A localized variant of the proposed method is developed. Numerical experiments with the Lorenz '63 model show that modified filters significantly improve the analyses for low dynamical ensemble sizes. Furthermore, localization experiments with the Lorenz '96 model show that the proposed methodology is extendable to larger systems.
It is shown that for dynamical time scales longer than 10 million years the ensemble of interstellar cold diffuse clouds may be described by equations of mass and momentum conservation where the pressure has the form p/v-squared/ where /v-squared/ is the mass averaged square velocity dispersion. Equilibrium cloud spectra are derived and an equation of state obtained. The equation of state undergoes a phase transition at a critical density close to the local interstellar gas density. At higher densities, more massive clouds are formed and the cloud fluid softens. It is suggested that this mechanism may naturally regulate the interstellar gas density if star formation is very efficient in the cooler high-density systems.