Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed inference”

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 703 records · Page 39

The MECA Payload as a Dust Analysis Laboratory on the MSP 2001 Lander

In a companion abstract, the "Mars Environmental Compatibility Assessment" (MECA) payload for Mars Surveyor Program 2001 (MSP 2001) is described in terms of its capabilities for addressing exobiology on Mars. Here we describe how the same payload elements perform in terms of gathering data about surface dust on the planet. An understanding of the origin and properties of dust is important to both human exploration and planetary geology. The MECA instrument is specifically designed for soil/dust investigations: it is a multifunctional laboratory equipped to assess particulate properties with wet chemistry, camera imagery, optical microscopy (potentially with LTV fluorescence capability), atomic force microscopy (AFM; potentially with mineral-discrimination capabilities), electrometry, active & passive external materials-test panels, mineral hardness testing, and electrostatic & magnetic materials testing. Additionally, evaluation of soil chemical and physical properties as a function of depth down to about 50 cm will be facilitated by the Lander/MECA robot arm on which the camera (RAC) and electrometer are mounted. Types of data being sought for the dust include: (1) general textural and grain-size characterization of the soil as a whole --for example, is the soil essentially dust with other components or is it a clast-supported material in which dust resides only in the clast interstices, (2) size frequency distribution for dust particles in the range 0.01 to 10.00 microns, (3) particle-shape distribution of the soil components and of the fine dust fraction in particular, (4) soil fabric such as grain clustering into clods, aggregates, and cemented/indurated grain amalgamations, as well as related porosity, cohesiveness, and other mechanical soil properties, (5) cohesive relationship that dust has to certain types of rocks and minerals as a clue to which soil materials may be prime hosts for dust "piggybacking", (6) particle, aggregate, and bulk soil electrostatic properties, (7) particle hardness, (8) particle magnetic properties, (9) bulk dust geochemistry (solubility, reactivity, ionic and mineral species). All of these quantities are needed in order for the human exploration program to make assessments of hazards on Mars, and to better enable the production on earth, of soil/dust simulants that can act as realistic test materials in terms of those properties that render dust a contaminant.Such properties include the small grain size that enables penetration of space-suit joints, mechanical interfaces and bearings, seals, etc., and presents difficulty for filtration systems. Size also plays a critical role in the potential for lung disease in long-term habitats. The properties of grain shape and hardness are important parameters in determining the abrasiveness of dust as it enters mechanical systems, or bombards helmet visors and habitat windows in dust-laden winds. Adhesive electrostatic and magnetic properties of dust will be prime causes of contamination of space suits and equipment. Contamination causes mechanical malfunction, tracking of dirt into habitats, "piggybacking" of toxins on dust into habitats, changes in albedo and efficiency of solar arrays and heat exchangers, and changes in electrical conductivity of suit surfaces and other materials that may have specific safety requirements regarding electrical conductivity. Other potentially hazardous properties of dust include the possibility of high solubility of some component grains (rendering them reactive), and toxicity of some materials --grains of superoxidants and heavy metals (there is always the slim, but not inconceivable possibility of biogenic components such as spores). Because Mars has no active surface aqueous regime, volcanic emissions, meteoritic debris, weathering products, and photochemical products of Mars have nowhere to go except reside in the surface; there are few mechanical or chemical (buffering) processes to remove the accumulation of eons. From a planetology perspective, there are many enigmatic issues relating to dust and the aeolian regime in general. MECA will be able to address many questions in this area. For example, if MECA determines a particular particle size distribution (size and sorting values), it will be possible to make inferences about the origin of the dust - - is it all aeolian, or a more primitive residue of weathering, volcanic emissions, and meteoritic gardening? Trenching with the Lander/MECA robot arm will enable local stratigraphy to be determined in terms of depositional rates, amounts and cyclicity in dust storms and/or local aeolian transport. Grain shape will betray the origin of the dust fragments as being the product of recent or ancient weathering, or the comminution products of aeolian transport --the dust-silt ratio might be a measure of aeolian comminution energy. Additional information is contained in the original.

Marshall, J.↗

Metallographic Cooling Rate of IVA Irons Revisited

There is long standing problem reconciling the chemical evidence that the IVA iron meteorites formed in a core with the diverse cooling rates reported by several researchers. This large inferred range of cooling rates suggests that the IVA irons were distributed at different depths in a parent body with a complex structure when the Widmanstatten pattern formed. On the other hand, some researchers argued that the diverse cooling rates in group IVA result from inaccurate model parameters such as phase diagram, interdiffusion coefficients, and kamacite nucleation and growth mechanisms. In addition, the measured cooling rates may not apply for the same cooling temperature ranges, and the variation in the crystallographic orientations of the Widmanstatten plates on the analysis surface may result in inaccurate measurements of widths needed for the computer simulation models. We have revaluated the major parameters in computer model developed by Hopfe and Goldstein and measured cooling rates for the IVA irons. Such data are useful in evaluating whether these meteorites were part of a single core of a parent body during the formation of the Widmanstatten pattern.

Yang, J.↗

The effect of a non-Maxwellian electron distribution on oxygen and iron ionization balances in the solar corona

Analytic expressions are derived for ionization and recombination rates in a parameterized non-Maxwellian electron velocity distribution with an enhanced high-energy tail. These expressions are then used in investigating the effect of such an enhancement in the high-energy tail of the coronal electron velocity distribution on the oxygen and iron ionization balances, O(+6) - O(+7) and Fe(+11) - Fe(+12). Relative to a Maxwellian of the same mean electron energy, the degree of ionization allowed by such a distribution is found to be either unchanged or slightly decreased for iron but often substantially increased for oxygen. The greater sensitivity of oxygen ionization balance to the high-energy distribution tail derives from the higher oxygen ionization threshold energy. It is noted that the electron temperature inferred from a measurement of the oxygen ionization ratio, O(+6)/O(+7), could indeed overestimate the actual coronal electron temperature by nearly 10 to the 6th K if the coronal electron distribution is incorrectly assumed to be Maxwellian.

Owocki, S. P.↗

Value-of-Information Analysis of a Fracture Prediction Method

A core-based fracture prediction method is used to illustrate a Value of Information (VOI) decision-analysis protocol to inform completion decisions in tight-gas sandstones. The ratio of late host-rock cement to available pore volume or degradation index uses petrographic observations of cement distributions in core (including sidewall cores) to predict whether nearby but unsampled fractures (widths > 0.5 to 1 mm) are sealed (non-conductive) or open (conductive). Measurements from four sandstone plays suggest that the index correctly predicts open vs. sealed fractures with an accuracy in excess of 80%. The value added is calculated using Bayesian inference in which accuracy of the index serves as the likelihood of the prior distribution of open fractures to assess the posterior probability that data represent a useful predictor of producibility. VOI of the prediction method is more than three times the cost to acquire the data. VOI is most sensitive to play-specific geologic and cost parameters including cost to drill, expected revenue from a successful well, cost of completion, cost of acquiring data for the index, and fracture probability distributions. The approach provides a way to value acquiring fracture data and points to a need for zone-specific production data in tight-gas sandstones.

02 PETROLEUM↗

The Escape Velocity Profile of the Milky Way from Gaia DR3

Abstract The escape velocity profile of the Milky Way offers a crucial and independent measurement of its underlying mass distribution and dark matter (DM) properties. Using a sample of stars from the third data release of Gaia with 6D kinematics and strict quality cuts, we obtain an escape velocity profile of the Milky Way from 4 to 11 kpc in Galactocentric radius. To infer the escape velocity in radial bins, we model the tail of the stellar speed distribution with both traditional power-law models and a new functional form that we introduce. While power-law models tend to rely on extrapolation to high speeds, we find our new functional form gives the most faithful representation of the observed distribution. Using this for the escape velocity profile, we constrain the properties of the Milky Way’s DM halo modeled as a Navarro–Frenck–White profile. Combined with constraints from the circular velocity at the solar position, we obtain a concentration and mass of c 200 c DM = 13.9 − 4.3 + 6.2 and M 200 c DM = 0.55 − 0.14 + 0.15 × 10 12 M ⊙ . This corresponds to a total Milky Way mass of M 200 c = 0.64 − 0.14 + 0.15 × 10 12 M ⊙ , which is on the low end of the historic range of the galaxy’s mass, but in line with other recent estimates.

Astronomy & Astrophysics↗

Toward designing effective exascale scientific computing workflows: experiences and best practices

Many fields within scientific computing have embraced advances in big-data analysis and machine learning, which often requires the deployment of large, distributed and complicated workflows that may combine training neural networks, performing simulations, running inference, and performing database queries and data analysis in asynchronous, parallel and pipelined execution frameworks. Such a shift has brought into focus the need for scalable, efficient workflow management solutions with reproducibility, error and provenance handling, traceability, and checkpoint-restart capabilities, among other needs. Here, we discuss challenges and best-practices for deploying exascale-generation computational science workflows on resources at the Oak Ridge Leadership Computing Facility (OLCF). We present our experiences with large-scale deployment of distributed workflows on the Summit supercomputer, including for bioinformatics and computational biophysics, materials science, and deep learning model optimization. We also present problems and solutions created by working within a Python-centric software base on traditional HPC systems, and discuss steps that will be required before the convergence of HPC, AI, and data science can be fully realized. Our results point to a wealth of exciting new possibilities for harnessing this convergence to tackle new scientific challenges.

Coletti, Mark↗

The heliospheric ambipolar potential inferred from sunward-propagating halo electrons

ABSTRACT We provide evidence that the sunward-propagating half of the solar wind electron halo distribution evolves without scattering in the inner heliosphere. We assume the particles conserve their total energy and magnetic moment, and perform a ‘Liouville mapping’ on electron pitch angle distributions measured by the Parker Solar Probe SPAN-E instrument. Namely, we show that the distributions are consistent with Liouville’s theorem if an appropriate interplanetary potential is chosen. This potential, an outcome of our fitting method, is compared against the radial profiles of proton bulk flow energy. We find that the inferred potential is responsible for nearly 100 per cent of the proton acceleration in the solar wind at heliocentric distances 0.18-0.79 AU. These observations combine to form a coherent physical picture: the same interplanetary potential accounts for the acceleration of the solar wind protons as well as the evolution of the electron halo. In this picture the halo is formed from a sunward-propagating population that originates somewhere in the outer heliosphere by a yet-unknown mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the emulator based halo model

We present cosmology results from a blinded joint analysis of cosmic shear, galaxy-galaxy weak lensing, and projected galaxy clustering, measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog—a cosmology analysis. We define luminosity-cut, and therefore nearly volume-limited, samples of SDSS galaxies to serve as the tracers of w p and as the lens samples for Δ⁢Σ in three spectroscopic redshift bins spanning the range 0.15± and Δ⁢Σ measurements, we use a single sample of about seven million source galaxies over 416 deg 2 , selected from HSC-Y3 based on having photometric redshifts (photo-z) greater than 0.75. The deep, high-quality HSC-Y3 data enable significant detections of the Δ⁢Σ signals, with integrated signal-to-noise ratio S/N~24 in the range 3≤R/[ℎ –1 Mpc]≤30 over the three lens samples. ξ ± has S/N ~19 in the range 8' ≤$\theta$ ≤5⁢0' and 30' ≤$\theta$ ≤15⁢0' for ξ + and ξ – , respectively. For cosmological parameter inference, we use the dark emulator package, combined with a halo occupation distribution prescription for the relation between galaxies and halos, to model w p and Δ⁢Σ down to quasinonlinear scales, and we estimate cosmological parameters after marginalizing over nuisance parameters. In our baseline analysis we employ an uninformative flat prior of the residual photo-z error, given by Π⁡(Δz ph )=$U$⁡(–1,1), to model a residual bias in the mean redshift of HSC source galaxies. Comparing the relative lensing amplitudes for Δ⁢Σ in the three redshift bins and for ξ ± with the single HSC source galaxy sample allows us to calibrate the photo-z parameter Δ⁢z ph to the precision of σ⁡(Δ⁢z ph )≃0.09. With these methods, we obtain a robust constraint on the cosmological parameters for the flat Λ⁢CDM model: S 8 =σ 8 ⁢(Ω m /0.3) 0.5 =$0.763^{+0.040}_{–0.036}$, or the best-constrained parameter given by $S^{'}_{8}$ = σ 8 ⁢(Ω m /0.3) 0.22 =0.721±0.028, determined with about 4% fractional precision. Based on multidimensional tension metrics, HSC-Y3 data exhibits about 2.5σ tension with the cosmological constraint inferred by Planck for the Λ⁢CDM model, and hints at a nonzero residual photo-z bias implying that the true mean redshift of the HSC galaxies at z ≳0.75 is higher than that implied by the original photo-z estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

A multiprocessing architecture for real-time monitoring

A multitasking architecture for performing real-time monitoring and analysis using knowledge-based problem solving techniques is described. To handle asynchronous inputs and perform in real time, the system consists of three or more distributed processes which run concurrently and communicate via a message passing scheme. The Data Management Process acquires, compresses, and routes the incoming sensor data to other processes. The Inference Process consists of a high performance inference engine that performs a real-time analysis on the state and health of the physical system. The I/O Process receives sensor data from the Data Management Process and status messages and recommendations from the Inference Process, updates its graphical displays in real time, and acts as the interface to the console operator. The distributed architecture has been interfaced to an actual spacecraft (NASA's Hubble Space Telescope) and is able to process the incoming telemetry in real-time (i.e., several hundred data changes per second). The system is being used in two locations for different purposes: (1) in Sunnyville, California at the Space Telescope Test Control Center it is used in the preflight testing of the vehicle; and (2) in Greenbelt, Maryland at NASA/Goddard it is being used on an experimental basis in flight operations for health and safety monitoring.

Schmidt, James L.↗

Vertical profiles of heating derived from IR-based precipitation estimates during FGGE SOP-1

This paper examines a technique for retrieving from geostationary IR data the vertical profiles of heating and cooling due to moist diabatic processes. First, GOES IR imagery is used to estimate precipitation fields which are independent of fields inferred from residuals in heat budget analysis based on the FGGE level III-b data. Vertical distributions of the associated heating are then obtained using thermodynamic data from the level III-b analysis, one-dimensional cloud models, and the satellite-estimated precipitation. The technique was applied to infer heating in the South Pacific convergence zone during a portion of FGEE SOP-1, and the results were compared with heat-budget calculations made using the ECMWF analyses.

Robertson, Franklin R.↗

Behavior of Langmuir Probes in Non-Equilibrium Plasmas

Langmuir probes are diagnostic tools used to determine electron temperature, number density, and plasma potential. Irving Langmuir first used an electrostatic probe in the 1920s to find these characteristics in ionized gases. Single, double, and triple Langmuir probes are commonly used in plasma diagnostics because of their relative simplicity. In the single probe, a swept voltage is applied between the probe tip and circuit common to acquire a waveform showing the collected current as a function of applied voltage. A double Langmuir probe consists of two tips, both inserted into the plasma, with a voltage applied between them. As this voltage is swept, a current-voltage characteristic is measured. In a triple probe three probe tips are electrically coupled to each other with constant non-swept voltages applied between each of the tips. The voltages are selected to represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Triple Langmuir probe measurements have been widely employed for various types of plasmas, including pulsed and time-varying plasmas such as those seen in pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. The typical Langmuir probe analysis for determining electron temperature and number density of the plasma (for a single, double, or triple Langmuir probe) includes an assumption that the plasma is in thermal equilibrium. While the this assumption may be justified for some applications, it is unlikely that it is fully justifiable for pulsed and time-varying plasmas or for the entire time a plasma device is in use. In the present work, we model the responses of Langmuir probes as they are inserted into a range of simple equilibrium and non-equilibrium plasmas. We return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. Comparing the distribution function that was assumed a priori with the plasma density and temperature one would infer when applying standard probe theory to analyze the collected currents serves to illustrate the effect a non- Maxwellian plasma would have on results interpreted using the equilibrium probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗

Examination of nucleon distribution with Bayesian imaging for isobar collisions

Relativistic collision of isobaric systems is found to be valuable in differentiating the nucleon distributions for nuclei with the same mass number. In recent contrast experiment of $^{96}_{44}$Ru + $^{96}_{44}$Ru versus $^{96}_{40}$Zr + $^{96}_{40}$Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV, the ratios of multiplicity distribution, elliptic flow, triangular flow, and radial flow are precisely measured and found to be significantly different from unity, indicating the difference in the shapes of the isobar pair. Here, in this work, we investigate the feasibility of nuclear structure reconstruction from heavy-ion collision observables. We perform Bayesian Inference with employing the Monte-Carlo Glauber model as an estimator of the mapping from nuclear structure to the final state observables and to provide the mock data for reconstruction. By varying combination of observables included in the mock data, we find it plausible to infer Woods–Saxon parameters from the observables. We also observe that single-system multiplicity distribution for the isobar system, rather than their ratio, is crucial to simultaneously determine the nuclear structure for the isobar system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Attraction Indian Buffet Distribution

We propose the attraction Indian buffet distribution (AIBD), a distribution for binary feature matrices influenced by pairwise similarity information. Binary feature matrices are used in Bayesian models to uncover latent variables (i.e., features) that explain observed data. The Indian buffet process (IBP) is a popular exchangeable prior distribution for latent feature matrices. In the presence of additional information, however, the exchangeability assumption is not reasonable or desirable. The AIBD can incorporate pairwise similarity information, yet it preserves many properties of the IBP, including the distribution of the total number of features. Thus, much of the interpretation and intuition that one has for the IBP directly carries over to the AIBD. A temperature parameter controls the degree to which the similarity information affects feature-sharing between observations. Unlike other nonexchangeable distributions for feature allocations, the probability mass function of the AIBD has a tractable normalizing constant, making posterior inference on hyperparameters straight-forward using standard MCMC methods. A novel posterior sampling algorithm is proposed for the IBP and the AIBD. We demonstrate the feasibility of the AIBD as a prior distribution in feature allocation models and compare the performance of competing methods in simulations and an application.

97 MATHEMATICS AND COMPUTING↗

Advancing Diagnostic Model Evaluation to Better Understand Water Shortage Mechanisms in Institutionally Complex River Basins

Abstract Water resources systems models enable valuable inferences on consequential system stressors by representing both the geophysical processes determining the movement of water and the human elements distributing it to its various competing uses. This study contributes a diagnostic evaluation framework that pairs exploratory modeling with global sensitivity analysis to enhance our ability to make inferences on water scarcity vulnerabilities in institutionally complex river basins. Diagnostic evaluation of models representing institutionally complex river basins with many stakeholders poses significant challenges. First, it needs to exploit a large and diverse suite of simulations to capture important human‐natural system interactions as well as institutionally aware behavioral mechanisms. Second, it needs to have performance metrics that are consequential and draw on decision‐relevant model outputs that adequately capture the multisector concerns that emerge from diverse basin stakeholders. We demonstrate the proposed model diagnostic framework by evaluating how potential interactions between changing hydrologic conditions and human demands influence the frequencies and durations of water shortages of varying magnitudes experienced by hundreds of users in a subbasin of the Colorado River. We show that the dominant factors shaping these effects vary both across users and, for an individual user, across percentiles of shortage magnitude. These differences hold even for users sharing diversion locations, demand levels or water right seniority. Our findings underline the importance of detailed institutional representation for such basins, as institutions strongly shape how dominant factors of stakeholder vulnerabilities propagate through the complex network of users.

Hadjimichael, Antonia↗

Textures and clusters

We discuss the properties of galaxy clusters expected in a texture-seeded, CDM-dominated, Omega = 1 universe. Assuming that the textures are spherical, we use the spherical collapse model to compute the cluster velocity dispersion (or temperature) distribution function. For objects of mass 10 exp 11 to 10 exp 15 solar masses, we find v varies as M super gamma with gamma of about 0.25. An unbiased (b = 1) texture model predicts too many high-velocity dispersion clusters. A biased texture model appears to be compatible with cluster properties inferred from optical and X-ray observations. In the texture model, the cluster velocity distribution functio does not evolve rapidly; thus, the model predicts the existence of rich clusters at moderate redshift (about 1-2).

Bartlett, James G.↗

Statistically-informed deep learning for gravitational wave parameter estimation

We introduce deep learning models to estimate the masses of the binary components of black hole mergers, $(m_1,m_2)$, and three astrophysical properties of the post-merger compact remnant, namely, the final spin, $a_\mathrm f$, and the frequency and damping time of the ringdown oscillations of the fundamental $\ell = m = 2$ bar mode, $(\omega_\mathrm R, \omega_\mathrm I)$. Our neural networks combine a modified WaveNet architecture with contrastive learning and normalizing flow. We validate these models against a Gaussian conjugate prior family whose posterior distribution is described by a closed analytical expression. Upon confirming that our models produce statistically consistent results, we used them to estimate the astrophysical parameters $(m_1,m_2, a_\mathrm f, \omega_\mathrm R, \omega_\mathrm I)$ of five binary black holes: GW150914, GW170104, GW170814, GW190521 and GW190630. We use PyCBC Inference to directly compare traditional Bayesian methodologies for parameter estimation with our deep learning based posterior distributions. Our results show that our neural network models predict posterior distributions that encode physical correlations, and that our data-driven median results and 90% confidence intervals are similar to those produced with gravitational wave Bayesian analyses. This methodology requires a single V100 NVIDIA GPU to produce median values and posterior distributions within two milliseconds for each event. Furthermore, this neural network, and a tutorial for its use, are available at the Data and Learning Hub for Science.

79 ASTRONOMY AND ASTROPHYSICS↗

Unusually Long Path Length for a Nearly Scatter-Free Solar Particle Event Observed By Solar Orbiter at 0.43 au

Context. After their acceleration and release at the Sun, solar energetic particles (SEPs) are injected into the interplanetary medium and are bound to the interplanetary magnetic field (IMF) by the Lorentz force. The expansion of the IMF close to the Sun focuses the particle pitch-angle distribution, and scattering counteracts this focusing. Solar Orbiter observed an unusual solar particle event on 9 April 2022 when it was at 0.43 astronomical units (au) from the Sun. Aims. We show that the inferred IMF along which the SEPs traveled was about three times longer than the nominal length of the Parker spiral and provide an explanation for this apparently long path. Methods. We used velocity dispersion analysis (VDA) information to infer the spiral length along which the electrons and ions traveled and infer their solar release times and arrival direction. Results. The path length inferred from VDA is approximately three times longer than the nominal Parker spiral. Nevertheless, the pitch-angle distribution of the particles of this event is highly anisotropic, and the electrons and ions appear to be streaming along the same IMF structures. The angular width of the streaming population is estimated to be approximately 30 degrees. The highly anisotropic ion beam was observed for more than 12 h. This may be due to the low level of fluctuations in the IMF, which in turn is very probably due to this event being inside an interplanetary coronal mass ejection The slow and small rotation in the IMF suggests a flux-rope structure. Small flux dropouts are associated with very small changes in pitch angle, which may be explained by different flux tubes connecting to different locations in the flare region. Conclusions. The unusually long path length along which the electrons and ions have propagated virtually scatter-free together with the short-term flux dropouts offer excellent opportunities to study the transport of SEPs within interplanetary structures. The 9 April 2022 solar particle event offers an especially rich number of unique observations that can be used to limit SEP transport models.

Robert F. Wimmer-Schweingruber↗

Intelligent resolution: Integrating Cryo-EM with AI-driven multi-resolution simulations to observe the severe acute respiratory syndrome coronavirus-2 replication-transcription machinery in action

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) replication transcription complex (RTC) is a multi-domain protein responsible for replicating and transcribing the viral mRNA inside a human cell. Attacking RTC function with pharmaceutical compounds is a pathway to treating COVID-19. Conventional tools, e.g., cryo-electron microscopy and all-atom molecular dynamics (AAMD), do not provide sufficiently high resolution or timescale to capture important dynamics of this molecular machine. Consequently, we develop an innovative workflow that bridges the gap between these resolutions, using mesoscale fluctuating finite element analysis (FFEA) continuum simulations and a hierarchy of AI-methods that continually learn and infer features for maintaining consistency between AAMD and FFEA simulations. We leverage a multi-site distributed workflow manager to orchestrate AI, FFEA, and AAMD jobs, providing optimal resource utilization across HPC centers. Our study provides unprecedented access to study the SARS-CoV-2 RTC machinery, while providing general capability for AI-enabled multi-resolution simulations at scale.

Trifan, Anda↗