Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fragmentation approximation method”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Microscopic predictions for the production of neutron-rich nuclei in the reaction Yb 176 + Yb 176

Background: Production of neutron-rich nuclei is of vital importance to both understanding nuclear structure far from stability and to informing astrophysical models of the rapid neutron capture process (r-process). Multinucleon transfer (MNT) in heavy-ion collisions offers a possibility to produce neutron-rich nuclei far from stability. Purpose: The 176 Yb + 176 Yb reaction has been suggested as a potential candidate to explore the neutron-rich region surrounding the principal fragments. The current study has been conducted with the goal of providing guidance for future experiments wishing to study this (or similar) system. Methods: Time-dependent Hartree-Fock (TDHF) and its time-dependent random-phase approximation (TDRPA) extension are used to examine both scattering and MNT characteristics in 176 Yb + 176 Yb. TDRPA calculations are performed to compute fluctuations and correlations of the neutron and proton numbers, allowing for estimates of primary fragment production probabilities. Results: Both scattering results from TDHF and transfer results from the TDRPA are presented for different energies, orvientations, and impact parameters. In addition to fragment composition, scattering angles and total kinetic energies, as well as correlations between these observables are presented. Conclusions: 176 Yb + 176 Yb appears to be an interesting probe for the midmass neutron-rich region of the chart of nuclides. Lastly, the predictions of both TDHF and TDRPA are speculative, and will benefit from future experimental results to test the validity of this approach to studying MNT in heavy, symmetric collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1974

The distribution in lunar soil of hydrogen released by pyrolysis

The hydrogen contents of various Apollo 15 rocks and Apollo 15 and 16 soil particle types and sieve fractions have been determined by the pyrolysis method. The Apollo 15 crystalline rocks analyzed release very little hydrogen (approximately 6 micromoles H/g sample). The mineral, glass, and light-colored breccia fragments in the soils contain approximately 5, 30, and 10 micromoles H/g sample, respectively. The dark breccias and glassy agglutinates examined have much higher hydrogen contents (60 and 110 micromoles H/g sample). For soil particle diameters less than 105 microns, hydrogen content increases linearly with particle surface area per particle mass. The calculated grain surface average concentrations of hydrogen in the Apollo 15 and 16 soils fall short of experimentally determined saturation values by a factor of 20. The volume-correlated component of hydrogen in a soil increases with the maturity of the soil.

Desmarais, D. J.↗
NASA NTRS2015

An Empirical Non-TNT Approach to Launch Vehicle Explosion Modeling

In an effort to increase crew survivability from catastrophic explosions of Launch Vehicles (LV), a study was conducted to determine the best method for predicting LV explosion environments in the near field. After reviewing such methods as TNT equivalence, Vapor Cloud Explosion (VCE) theory, and Computational Fluid Dynamics (CFD), it was determined that the best approach for this study was to assemble all available empirical data from full scale launch vehicle explosion tests and accidents. Approximately 25 accidents or full-scale tests were found that had some amount of measured blast wave, thermal, or fragment explosion environment characteristics. Blast wave overpressure was found to be much lower in the near field than predicted by most TNT equivalence methods. Additionally, fragments tended to be larger, fewer, and slower than expected if the driving force was from a high explosive type event. In light of these discoveries, a simple model for cryogenic rocket explosions is presented. Predictions from this model encompass all known applicable full scale launch vehicle explosion data. Finally, a brief description of on-going analysis and testing to further refine the launch vehicle explosion environment is discussed.

Blackwood, James M.↗
NASA NTRS1972

Exposure ages and neutron capture record in lunar samples from Fra Mauro.

Cosmic-ray exposure ages of Apollo 14 rocks and rock fragments obtained by the Kr81-Kr83 method range from 27 to 700 m.y. Rock 14321, collected near the Cone crater rim, is one of the many approximately 27 m.y. old ejecta which were reported at the Third Lunar Science Conference. All the other rocks have considerably higher exposure ages. Isotopic anomalies from neutron capture in gadolinium, bromine, and barium are used to obtain information on the lunar neutron spectrum at various depths below the lunar surface. The flux ratio of resonance and slow (less than 0.3 eV) neutrons is found to be nearly constant in the topmost approximately 100 g/sq cm.

Lugmair, G. W.↗
NASA NTRS1995

Fission foil detector calibrations with high energy protons

Fission foil detectors (FFD's) are passive devices composed of heavy metal foils in contact with muscovite mica films. The heavy metal nuclei have significant cross sections for fission when irradiated with neutrons and protons. Each isotope is characterized by threshold energies for the fission reactions and particular energy-dependent cross sections. In the FFD's, fission fragments produced by the reactions are emitted from the foils and create latent particle tracks in the adjacent mica films. When the films are processed surface tracks are formed which can be optically counted. The track densities are indications of the fluences and spectra of neutrons and/or protons. In the past, detection efficiencies have been calculated using the low energy neutron calibrated dosimeters and published fission cross sections for neutrons and protons. The problem is that the addition of a large kinetic energy to the (n,nucleus) or (p,nucleus) reaction could increase the energies and ranges of emitted fission fragments and increase the detector sensitivity as compared with lower energy neutron calibrations. High energy calibrations are the only method of resolving the uncertainties in detector efficiencies. At high energies, either proton or neutron calibrations are sufficient since the cross section data show that the proton and neutron fission cross sections are approximately equal. High energy proton beams have been utilized (1.8 and 4.9 GeV, 80 and 140 MeV) for measuring the tracks of fission fragments emitted backward and forward.

Benton, E. V.↗
NASA NTRS1984

An abrasion-ablation model description of galactic heavy-ion fragmentation

The fragmentation of high-energy galactic heavy ions by nuclear interactions with arbitrary target nuclei is described within the context of a simple abrasion-ablation fragmentation model. The abrasion part of the theory utilizes a quantum-mechanical formalism based upon an optical model potential approximation to the exact nucleus-nucleus multiple-scattering series. Nuclear charge distributions of the excited prefragments are calculated using either a hypergeometric distribution or a method based upon the zero-point oscillations of the giant dipole resonance. The excitation energy of the prefragment is estimated from the geometric clean-cut abrasion-ablation model. The decay probabilities for the various particle emission channels, in the ablation stage of the fragmentation, are obtained from the EVAP-4 Monte Carlo computer program. Elemental production cross sections for 1.88-GeV/nucleon iron colliding with carbon, silver, and lead targets are calculated and compared with experimental data and with the predictions from the semiempirical relations of Silberberg and Tsao.

Townsend, L. W.↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Transport methods and interactions for space radiations

This report presents a brief history leading to the involvement of the Langley Research Center of the National Aeronautics and Space Administration (NASA) in space-radiation physics and protection. Indeed, a relatively complete summary of technical capability as of the summer of 1990 is given. The Boltzmann equations for coupled ionic and neutronic fields are presented and inversion techniques for the Boltzmann operator are discussed. Errors generated by the straight ahead approximation are derived and are shown to be negligible for most problems of space-radiation protection. A decoupling of projectile propagation from the target fields greatly simplifies the Boltzmann equations and allows an analytic solution of the target fragment transport. Analytic and numerical methods of solving the projectile transport equations are discussed. The nuclear physics underlying the coefficients in the Boltzmann equation is discussed. A coupled-channel optical model is found as a consequence of the loose binding of nuclear matter and closure of the nuclear states in high-energy reactions. Transport solutions with the developed data base are used with laboratory experiments to validate both the transport code and the data base. Numerical benchmarks and comparison with Monte Carlo calculations are also used for code validation.

Wilson, John W.↗
DOE OSTI2021

Projector-Free Capped-Fragment Scheme within Density Functional Embedding Theory for Covalent and Ionic Compounds

Quantum-mechanics-(QM)-based simulations now routinely aid in understanding and even discovering new chemistries involving molecules and materials exhibiting desired functionalities. Ab initio correlated wavefunction (CW) theories systematically improve QM methods, with many exhibiting high accuracy. However, execution of CW methods requires expensive computations that typically scale poorly with system size. Divide-and-conquer approaches partition large systems into smaller fragments; a lower level of theory treats fragment interactions while a preferred higher level of theory describes the important fragment. These methods offer ways to incorporate CWs into chemical simulations of large systems, e.g., biomolecules, surfaces, large inorganic clusters, bulk crystals, etc. Here we propose a partitioning protocol that utilizes capping atoms to saturate severed covalent bonds at fragment interfaces and density functional embedding theory (DFET) to describe fragment interactions. The capping groups in each fragment provide an ad-hoc potential that approximates the effects of the environment. An embedding potential optimized via DFET then serves as an augmentation of the capping group to simulate the effects of the environment. We concurrently use an auxiliary fragment (a separate system comprised of only the combined capping groups) to account for, and thereby correct, the electron density contributions of all the capping groups added to all of the fragments. This method depends only on the capped-subsystem and auxiliary-fragment electron densities, forgoing, as with the original DFET developed for metallic systems, orbital-based projector approaches that determine a nonlocal action of the embedding potential onto the fragment electrons. By using an auxiliary fragment, the method maintains a purely electron-density-dependent embedding potential, substantially lessening the cost and leading to simpler implementation. In this work, we demonstrate the utility of our capped-DFET and ensuing capped embedded CW method in two contrasting systems, namely, an organic molecule and an ionic metal oxide cluster.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1978

The Yamato-74659 ureilite and some new findings on the Yamato achondritic pyroxenes

A cube-like fragment with an approximately 4 mm edge was taken from the Yamato-74659 meteorite for the mineralogical and crystallographic examination discussed in the present paper. An analysis of powdered fragments by a standard wet chemical method showed that it matched to weathered ureilite. The chemical composition was found to differ from other ureilites in that the FeO content (8.83 wt%) is the lowest among known ureilites and it is richer in SiO2 contents (42.91 wt%).

Takeda, H.↗
DOE OSTI2025

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗
NASA NTRS1986

Optical model calculations of heavy-ion target fragmentation

The fragmentation of target nuclei by relativistic protons and heavy ions is described within the context of a simple abrasion-ablation-final-state interaction model. Abrasion is described by a quantum mechanical formalism utilizing an optical model potential approximation. Nuclear charge distributions of the excited prefragments are calculated by both a hypergeometric distribution and a method based upon the zero-point oscillations of the giant dipole resonance. Excitation energies are estimated from the excess surface energy resulting from the abrasion process and the additional energy deposited by frictional spectator interactions of the abraded nucleons. The ablation probabilities are obtained from the EVA-3 computer program. Isotope production cross sections for the spallation of copper targets by relativistic protons and for the fragmenting of carbon targets by relativistic carbon, neon, and iron projectiles are calculated and compared with available experimental data.

Townsend, L. W.↗
NASA NTRS2017

A New Type of Foreign Clast in A Polymict Ureilite: A CAI or AL-Rich Chondrule

Introduction: Polymict ureilites are breccias interpreted to represent regolith formed on a ureilitic asteroid [1-3]. They consist of approximately 90-95% clasts of various ureilite types (olivine-pyroxene rocks with Fo 75-95), a few % indigenous feldspathic clasts, and a few % foreign clasts [4-20]. The foreign clasts are diverse, including fragments of H, L, LL and R chondrites, angrites, other achondrites, and dark clasts similar to CC [6,7,9-19]. We report a new type of foreign clast in polymict ureilite DaG 999. Methods: Clast 8 in Dar al Gani (DaG) 999/1 (Museum fur Naturkunde) was discovered during a survey of feldspathic clasts in polymict ureilites [19,20]. It was studied by BEI, EMPA, and X-ray mapping on the JEOL 8530F electron microprobe at ARES, JSC. Petrography and Mineral Compositions: Clast 8 is sub-rounded to irregular in shape, approximately 85 micrometers in diameter, and consists of approximately 68% pyroxene and 32% mesostasis (by area). Part of the pyroxene (top half of clast in Fig. 1a and 2) shows a coarse dendritic morphology; the rest appears massive. Mesostasis may be glassy and contains fine needles/grains of pyroxene. The pyroxene has very high CaO (23.5 wt.%) and Al2O3 (19.7 wt.%), with the formula: (Ca(0.91)Mg(0.63)Fe(0.01)Al(sup VI) (0.38)Cr(0.01)Ti(0.05)1.99 Si2O6. The bulk mesostasis also has very high Al2O3 (approximately 26 wt.%). A bulk composition for the clast was obtained by combining modal abundances with phase compositions (Table 1, Fig. 3). Discussion: The pyroxene in clast 8 has a Ca-Al-(Ti)- rich (fassaitic) composition that is clearly distinct from compositions of pyroxenes in main group ureilites [22] or indigenous feldspathic clasts in polymict ureilites [4-8]. It also has significantly higher Al than fassaite in angrites (up to approximately 12 wt.% [23]), which occur as xenoliths in polymict ureilites. Ca-Al-Ti rich pyroxenes are most commonly found in CAIs, Al-rich chondrules and other types of refractory inclusions in chondrites [21,24-31]. However, the clast 8 pyroxene matches only the most Al-Ca-rich of these, e.g., pyroxenes in type B CAIs in CV3 chondrites [25,30,31], a pyroxene-hibonite spherule and a pyroxene-anorthitespinel fragment from unique CC Acfer 094 [29], and one Al-rich chondrule from Chainpur (LL3.4) [21]. The mineralogy of clast 8 is not consistent with the mineral assemblages of any of these objects (since it lacks hibonite, spinel and/or anorthite), which suggests that it is unrepresentatively sectioned or is a fragment of a more mineralogically diverse object. Its bulk composition (Table 1; Fig. 3) is similar to bulk compositions of some Al-rich chondrules, as well as those of Type C CAIs (which plot in the sp+An+L field in Fig. 3), although it is enriched in silica relative to type C CAIs [e.g., 31]. This suggests a more likely affinity to Al-rich chondrules, although most Al-rich chondrules have less Al-Ca-rich pyroxene [21,26,27]. These bulk compositional comparisons may not be definitive, however, if the clast is unrepresentatively sampled. One of eleven Al-rich chondrules from UOCs described by [21] has textural and compositional characteristics that make it a possible progenitor type for clast 8. This chondrule (Chainpur 1251-14-2) is anorthiteporphyritic, with an interstitial dendritic intergrowth of pyroxene (similar in composition to that in clast 8) and plagioclase [21]. Clast 8 is conceivably a fragment from the interstitial area of such an object. The occurrence of glassy mesostasis (in clast 8) rather than plagioclase may not be a significant difference; it could result from a difference only in cooling rate. Al-rich chondrules with glassy mesostasis are rare, and known occurrences are Ca-poor [26], unlike clast 8. Polymict ureilites are known to contain xenoliths of various chondrites (including OC, R and CC) as well as individual ferromagnesian and silica-pyroxene chondrules probably derived from OC or RC [6,9,15,16,18]. This is the first report of an individual chondritic refractory inclusion as a xenolith in a polymict ureilite. An RC-like sample from anomalous polymict ureilite Almahata Sitta contains CAIs, but they are spinel-rich and not similar to clast 8 [13,14]. Further studies of this clast (which, unfortunately, may not be possible), or the discovery of additional (more representative?) materials of this type would be needed to determine the exact nature of this xenolith and the type of chondrite from which it is derived.

Goodrich, C. A.↗
DOE OSTI2024

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.

Kurzak, Jakub↗
NASA NTRS1978

User's guide to computer programs JET 5A and CIVM-JET 5B to calculate the large elastic-plastic dynamically-induced deformations of multilayer partial and/or complete structural rings

These structural ring deflections lie essentially in one plane and, hence, are called two-dimensional (2-d). The structural rings may be complete or partial; the former may be regarded as representing a fragment containment ring while the latter may be viewed as a 2-d fragment-deflector structure. These two types of rings may be either free or supported in various ways (pinned-fixed, locally clamped, elastic-foundation supported, mounting-bracket supported, etc.). The initial geometry of each ring may be circular or arbitrarily curved; uniform-thickness or variable-thickness rings may be analyzed. Strain-hardening and strain-rate effects of initially-isotropic material are taken into account. An approximate analysis utilizing kinetic energy and momentum conservation relations is used to predict the after-impact velocities of each fragment and of the impact-affected region of the ring; this procedure is termed the collision-imparted velocity method (CIVM) and is used in the CIVM-JET 5 B program. This imparted-velocity information is used in conjunction with a finite-element structural response computation code to predict the transient, large-deflection, elastic-plastic responses of the ring. Similarly, the equations of motion of each fragment are solved in small steps in time. Provisions are made in the CIVM-JET 5B code to analyze structural ring response to impact attack by from 1 to 3 fragments, each with its own size, mass, translational velocity components, and rotational velocity. The effects of friction between each fragment and the impacted ring are included.

Wu, R. W. H.↗
DOE OSTI2023

Breakups are complicated: an efficient representation of collisional breakup in the superdroplet method

Abstract. A key constraint of particle-based methods for modeling cloud microphysics is the conservation of total particle number, which is required for computational tractability. The process of collisional breakup poses a particular challenge to this framework, as breakup events often produce many droplet fragments of varying sizes, which would require creating new particles in the system. This work introduces a representation of collisional breakup in the so-called “superdroplet” method which conserves the total number of superdroplets in the system. This representation extends an existing stochastic collisional-coalescence scheme and samples from a fragment size distribution in an additional Monte Carlo step. This method is demonstrated in a set of idealized box model and single-column warm-rain simulations. We further discuss the effects of the breakup dynamic and fragment size distribution on the particle size distribution, hydrometeor population, and microphysical process rates. Box model experiments serve to characterize the impacts of properties such as coalescence efficiency and fragmentation function on the relative roles of collisional breakup and coalescence. The results demonstrate that this representation of collisional breakup can produce a stationary particle size distribution, in which breakup and coalescence rates are approximately equal, and that it recovers expected behavior such as a reduction in precipitate-sized particles in the column model. The breakup algorithm presented here contributes to an open-source pythonic implementation of the superdroplet method, PySDM, which will facilitate future research using particle-based microphysics.

de Jong, Emily (ORCID:0000000253104554)↗
DOE OSTI2025

Azimuthal anisotropies of charged particles with high transverse momentum in Pb + Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV with the ATLAS detector

A measurement is presented of elliptic (𝑣 2 ) and triangular (𝑣 3 ) azimuthal anisotropy coefficients for charged particles produced in Pb + Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV using a dataset corresponding to an integrated luminosity of 0.44nb −1 collected with the ATLAS detector at the LHC in 2018. The values of 𝑣 2 and 𝑣 3 are measured for charged particles over a wide range of transverse momentum (𝑝 T ), 1–400 GeV, and Pb + Pb collision centrality, 0–60%, using the scalar-product and multiparticle cumulant methods. These methods are sensitive to event-by-event fluctuations and nonflow effects in the measurements of azimuthal anisotropies. Positive values of 𝑣 2 are observed up to a 𝑝 T of approximately 100 GeV from both methods across all centrality intervals. Positive values of 𝑣 3 are observed up to approximately 25 GeV using both methods, though the application of the three-subevent technique to the multiparticle cumulant method leads to significant changes at the highest 𝑝 T . At high 𝑝 T (𝑝 T ⪆ 10 GeV), charged particles are dominantly from jet fragmentation. These jets, and hence the measurements presented here, are sensitive to the path-length dependence of parton energy loss in the quark-gluon plasma produced in Pb + Pb collisions.

Collective flow↗
DOE OSTI2020

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗