Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular modelling”

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 973 records · Page 54

Enhanced Oxygen Activation Achieved by Robust Single Chromium Atom-Derived Catalysts in Aerobic Oxidative Desulfurization

Oxidative desulfurization (ODS) plays critical roles in the production of high-quality sulfur impurity-free fuels, especially procedures using molecular oxygen as the sole oxidant. However, the state-of-the-art systems still rely on noble metal nanocatalysts, with deactivation issues caused by coking and sintering still present. In this work, leveraging the merits of single-atom catalysts (SACs) with earth-abundant metal cores and robust nanoporous supports, a series of catalysts composed of homogeneously distributed single chromium atoms anchored on multiwalled carbon nanotubes were fabricated and deployed to catalyze the aerobic ODS transformation. Adopting aromatic sulfur compound (ASC)-containing alkanes as a model system with molecular oxygen as the oxidant, efficient oxygen-to-active O 2 •– transformation and subsequent ASCs-to-sulfone conversion have been achieved, with the former benefiting from the chromium-active sites and the latter arising from the robust, nanoporous, and π-conjugated architecture of the supports. The attractive catalytic performance and cycling stability of the chromium-based SACs make them promising candidates in practical sulfur species removal from liquid fuels, supplying an alternative guidance on the catalyst design toward cost-effective and energy-efficient ODS procedures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Backpropagated frame transformation theory: A reformulation

In this work, the energy-dependent frame transformation theory of Gao and Greene [H. Gao and C. H. Greene, Phys. Rev. A 42, 6946 (1990)] is extended to yield a quantitatively accurate description of the dissociative recombination process. Evidence is presented to show that direct application of the original theory leads to inaccurate cross sections. A major revision, based on an interaction-free backpropagation of the Born-Oppenheimer solutions, markedly improves the frame transformation theory, reducing its average error by orders of magnitude. The original theory and its extension are tested on the previously explored two-dimensional (2D) model that is tailored to describe the singlet ungerade states of molecular hydrogen. The 2D model can be solved exactly (within the numerical accuracy) without implementing the Born-Oppenheimer approximation. These exact results then serve as a benchmark for the frame transformation theory developed in this paper.

74 ATOMIC AND MOLECULAR PHYSICS↗

Building Toward the Future in Chemical and Materials Simulation with Accessible and Intelligently Designed Web Applications

Over the last few decades, significant progress has been made in the development and use of electronic structure and other molecular simulation methods. As these methods become more mature and are able to simulate larger and more complex chemical simulations, the need for improvement in scientific visualization, molecular builders, simplified input to simulation methods, and the development of new approaches and languages to describe simulations, along with workflows to carry them out, becomes more apparent. In this chapter, we describe our recent efforts in developing a prototype open-source computational tool called Arrows that combines NWChem, SQL and NoSQL databases, email, web APIs, and web applications in a way that make molecular and materials modeling accessible to all scientists and engineers. At the same time, because of its simplified input, it provides a framework for expert users to carry out large numbers of calculations and run complex workflows.

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

Multinucleon transfer in the interaction of 977 MeV and 1143 MeV Hg 204 with Pb 208

A previous study of symmetric collisions of massive nuclei has shown that current models of multinucleon transfer (MNT) reactions do not adequately describe the transfer product yields. To gain further insight into this problem, we have measured the yields of MNT products in the interaction of 977 (E/A = 4.79 MeV) and 1143 MeV (E/A = 5.60 MeV) 204 Hg with 208 Pb. We find that the yield of multinucleon transfer products are similar in these two reactions and are substantially lower than those observed in the reaction of 1257 MeV (E/A = 6.16 MeV) 204 Hg+ 198 Pt. We compare our measurements with the predictions of the GRAZING-F, dinuclear systems (DNS), and improved quantum molecular dynamics (ImQMD) models. We report that for the observed isotopes of the elements Au, Hg, Tl, Pb, and Bi, the measured values of the MNT cross sections are orders of magnitude larger than the predicted values. Furthermore, the various models predict the formation of nuclides near the N = 126 shell, which are not observed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Molecular simulation of liquid–vapor coexistence for NaCl: Full-charge vs scaled-charge interaction models

Scaled-charge models have been recently introduced for molecular simulations of electrolyte solutions and molten salts to attempt to implicitly represent polarizability. Although these models have been found to accurately predict electrolyte solution dynamic properties, they have not been tested for coexistence properties, such as the vapor pressure of the melt. In this work, we evaluate the vapor pressure of a scaled-charge sodium chloride (NaCl) force field and compare the results against experiments and a non-polarizable full-charge force field. The scaled-charge force field predicts a higher vapor pressure than found in experiments, due to its overprediction of the liquid-phase chemical potential. Reanalyzing the trajectories generated from the scaled-charge model with full charges improves the estimation of the liquid-phase chemical potential but not the vapor pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A modular and extensible CHARMM-compatible model for all-atom simulation of polypeptoids

Peptoids (N-substituted glycines) are a class of sequence-defined synthetic peptidomimetic polymers with applications including drug delivery, catalysis, and biomimicry. Classical molecular simulations have been used to predict and understand the conformational dynamics of single chains and their self-assembly into morphologies including sheets, tubes, spheres, and fibrils. The CGenFF-NTOID model based on the CHARMM General Force Field has demonstrated success in accurate all-atom molecular modeling of peptoid structure and thermodynamics. Extension of this force field to new peptoid side chains has historically required reparameterization of side chain bonded interactions against ab initio data. This fitting protocol improves the accuracy of the force field but is also burdensome and precludes modular extensibility of the model to arbitrary peptoid sequences. In this work, we develop and demonstrate a Modular Side Chain CGenFF-NTOID (MoSiC-CGenFF-NTOID) as an extension of CGenFF-NTOID employing a modular decomposition of the peptoid backbone and side chain parameterizations, wherein arbitrary side chains within the large family of substituted methyl groups (i.e., –CH 3 , –CH 2 R, –CHRR', and –CRR'R") are directly ported from CGenFF. We validate this approach against ab initio calculations and experimental data to develop a MoSiC-CGenFF-NTOID model for all 20 natural amino acid side chains along with 13 commonly used synthetic side chains and present an extensible paradigm to efficiently determine whether a novel side chain can be directly incorporated into the model or whether refitting of the CGenFF parameters is warranted. We make the model freely available to the community along with a tool to perform automated initial structure generation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular clouds and galactic spiral structure

Galactic CO line emission at 115 GHz was surveyed in order to study the distribution of molecular clouds in the inner galaxy. Comparison of this survey with similar H1 data reveals a detailed correlation with the most intense 21 cm features. To each of the classical 21 cm H1 spiral arms of the inner galaxy there corresponds a CO molecular arm which is generally more clearly defined and of higher contrast. A simple model is devised for the galactic distribution of molecular clouds. The modeling results suggest that molecular clouds are essentially transient objects, existing for 15 to 40 million years after their formation in a spiral arm, and are largely confined to spiral features about 300 pc wide.

Dame, T. M.↗

Computational approaches to model X-ray photon correlation spectroscopy from molecular dynamics

X-ray photon correlation spectroscopy (XPCS) allows for the resolution of dynamic processes within a material across a wide range of length and time scales. X-ray speckle visibility spectroscopy (XSVS) is a related method that uses a single diffraction pattern to probe ultrafast dynamics. Interpretation of the XPCS and XSVS data in terms of underlying physical processes is necessary to establish the connection between the macroscopic responses and the microstructural dynamics. To aid the interpretation of the XPCS and XSVS data, we present a computational framework to model these experiments by computing the X-ray scattering intensity directly from the atomic positions obtained from molecular dynamics simulations. We compare the efficiency and accuracy of two alternative computational methods: the direct method computing the intensity at each diffraction vector separately, and a method based on fast Fourier transform that computes the intensities at all diffraction vectors at once. Furthermore, the computed X-ray speckle patterns capture the density fluctuations over a range of length and time scales and are shown to reproduce the known properties and relations of experimental XPCS and XSVS for liquids.

36 MATERIALS SCIENCE↗

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Freely jointed chain models with extensible links

We report analytical relations for the mechanical response of single polymer chains are valuable for modeling purposes, on both the molecular and the continuum scale. These relations can be obtained using statistical thermodynamics and an idealized single-chain model, such as the freely jointed chain model. To include bond stretching, the rigid links in the freely jointed chain model can be made extensible, but this almost always renders the model analytically intractable. Here, an asymptotically correct statistical thermodynamic theory is used to develop analytic approximations for the single-chain mechanical response of this model. The accuracy of these approximations is demonstrated using several link potential energy functions. This approach can be applied to other single-chain models, and to molecular stretching in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An adaptive granular representative volume element model with an evolutionary periodic boundary for hierarchical multiscale analysis

The hierarchical multiscale analysis normally utilizes a microscopic representative volume element (RVE) model to capture path/history-dependent macroscopic responses instead of using phenomenological constitutive models. However, for problems involving large deformation, the current RVE model used in geomechanics may lose representative properties due to the progressive distortion of the RVE box, unless a proper reinitialization is applied. This work develops an adaptive RVE model in conjunction with an evolutionary periodic boundary (EPB) algorithm for hierarchical multiscale analysis of granular materials undergoing large deformation based on a recent RVE model proposed for coupling molecular dynamics and the material point method. The proposed adaptive RVE model avoids the reinitialization of the RVE box that even undergoes extremely large shear deformation; meanwhile, it accounts for the deformation history of the RVE model and treats the interaction between boundary particles and other image particles in a more efficient way. Numerical examples with extremely large deformation are used to illustrate the adaptive granular RVE model enhanced by the proposed EPB algorithm. Furthermore, some key features of this new methodology are further discussed for clarification.

42 ENGINEERING↗

Challenges in the detection of gases in exoplanet atmospheres

Claims of detections of gases in exoplanet atmospheres often rely on comparisons between models including and excluding specific chemical species. However, the space of molecular combinations available for model construction is vast and highly degenerate. Only a limited subset of these combinations is typically explored for any given detection. As a result, apparent detections of trace gases risk being artefacts of incomplete modelling rather than robust identification of atmospheric constituents, especially in the low-signal-to-noise regime. Here, using the sub-Neptune K2-18 b as a case study, we show that recent biosignature claims vanish when the model space is expanded, with numerous alternatives providing equally good or better fits. We demonstrate that the significance of a claimed detection relies on the choice of models being compared, and that model preference does not in itself imply the presence of a specific gas. In conclusion, we recommend treating model comparisons instead as relative adequacy tests, which should be supported by theoretical predictions and complementary metrics of statistical significance to attribute a signal to a particular gas.

Astronomy and AstroPhysics↗

Current status and future outlook for patient‐derived cancer models from a rare cancer research perspective

Abstract Rare cancers are a group of approximately 200 malignancies with extremely low incidences and with a wide variety of genotypes and phenotypes. Collectively, they are more common than any single malignancy. However, given the small numbers of individuals diagnosed with rare cancers, it is difficult to design clinical trials with sufficient patient numbers. Therefore, few effective anticancer drugs have been developed, and evidence‐based medicine is not always feasible for rare cancers. Consequently, their clinical outcomes are generally poorer. Cancer research requires adequate models that faithfully recapitulate molecular features and reproduce treatment responses of the original tumors. Such models allow us to focus on more efficacious drugs in the clinical studies. For rare cancers, patient‐derived cancer models are particularly important because the enrollment of sufficient patients is rarely attainable within a reasonable period of time. However, extremely few models are available for rare cancers. For example, cell lines and xenografts are available for only a limited number of histological subtypes of sarcomas; therefore, most sarcoma research is performed without such models, and a lack of adequate cancer models causes a lag in therapeutic development. The establishment of novel rare cancer models will dramatically facilitate rare cancer research and treatment development in the near future. This review focuses on the status of patient‐derived rare cancer models and discusses their pivotal problems and possibilities, using sarcomas as a representative rare cancer type. Multi‐institutional collaboration will help address the scarcity of patient‐derived rare cancer models.

Kondo, Tadashi↗

Studying light–matter interactions and energy transfer at the nanoscale with a trapped–ion quantum computer

A grand challenge in materials sciences is to control energy transduction, transfer and dissipation at the nanoscale. Meeting this challenge requires progress in nanoscale fabrication and synthesis, but equally importantly, also progress in our understanding of how nanoscale structure influences these phenomena. Although detailed quantum mechanical models exist that capture the structure and physics of materials and molecular species, solving these models to calculate the relevant emergent phenomena, which often occur over vastly different time and energy scales, has so far proven difficult. In this project, we studied how to model energy transfer and light-matter interactions using trapped-ion quantum computers. Our studies focused on emulation of critical processes in natural systems, such as the well-known photosynthetic pigment-protein complexes responsible for capturing and transducing light to electrons. In addition, we implemented proof-of-principle experiments demonstrating that such calculations are possible on trapped-ion quantum computing hardware. These studies and experiments pave the way to investigate energy transport and transduction in a detail not accessible with classical computing technologies. We expect that our studies will allow to extract important guiding principles on how to improve the efficiency of nano-scale light-harvesting and transduction devices.

36 MATERIALS SCIENCE↗

Engineering models of the Venus atmosphere.

Engineering models of the Venus atmosphere are presented for use in spacecraft design and mission planning. The models are developed from (1) data measured by the U.S. Mariner 5 and USSR Venera spacecraft, (2) earth-based measurements, and (3) theory which best satisfies the observations. The models include a most probable model and models for extremes of molecular mass and solar activity. In addition to the basic engineering models, other atmospheric characteristics and parameters affecting the models are discussed including planetary radius, atmospheric composition and clouds.

Noll, R. B.↗

Evolution of models for evolution

The paper discusses models of evolution and their values, and some critiques of these models and the value of these critiques. A model is investigated which involves the formation of prebiotic protein from amino acids. Its formation by four theoretical critiques that suggest alternative environmental conditions is discussed. Experiments are reviewed so as to illustrate the experimental testing of the possible weaknesses of a model for a single molecular evolutionary phase and to suggest some necessary changes in the model.

Rohlfing, D. L.↗

Mass and energy balance of the cold Io torus

A new model of the cold Io torus is described. Ions and energy are injected into the system by independent processes so that the mass balance is isolated from the energy balance. The primary source of energy is local ionization and acceleration of hot pickup ions resulting from charge exchange between thermal ions and an extended cloud of Iogenic sulfur and oxygen atoms. The primary energy loss mechanism of the plasma is collisionally excited line emission at optical wavelengths. The primary ion source is radial diffusion inward from the hot torus on a time scale of 140-710 days. The primary ion loss mechanism is a novel two-step enhanced recombination mechanism involving charge exchange between thermal ions and an extended cloud of neutral SO2 molecules, followed by rapid dissociative recombination of the resultant molecular ion. The model provides a self-consistent solution which reconciles a number of diverse observations with known physical processes.

Moreno, M. A.↗