Search NASA⌕ Search

SEARCH · Search NASA

Results for “constrained molecular dynamics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

The role of the Heisenberg principle in constrained molecular dynamics model

In this study, we implement the Heisenberg principle into the Constrained Molecular Dynamics model with a similar approach to the Pauli principle using the one-body occupation probability $\bar{f}_i$. Results of the modified and the original model with comparisons to data are given. The binding energies and the radii of light nuclei obtained with the modified model are more consistent with the experimental data than the original one. The collision term and the density distribution are tested through a comparison to p+ 12 C elastic scattering data. Some simulations for fragmentation and superheavy nuclei production are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Simulation of fusion and quasi-fission in nuclear reactions leading to production of superheavy elements using the Constrained Molecular Dynamics model

Fusion dynamics and the onset of quasi-fission in reactions, leading to production of superheavy nuclei are investigated using the constrained molecular dynamics model. Constraints on the parameters of the nuclear equation of state are derived from experimental fusion probabilities. Here, the obtained constraint on the modulus of incompressibility of nuclear matter K 0 = 240 –260 MeV is consistent with the results of previous study using the Boltzmann-Uehling-Uhlenbeck equation and also with constraints derived using the recent neutron star binary collision event GW170817. Unlike the modulus of incompressibility of symmetric nuclear matter, the stiffness of the density-dependence of symmetry energy influences the fusion probability only weakly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Kinetically constrained ring-polymer molecular dynamics for non-adiabatic chemical reactions

We extend ring-polymer molecular dynamics (RPMD) to allow for the direct simulation of general, electronically non-adiabatic chemical processes. The kinetically constrained (KC) RPMD method uses the imaginary-time path-integral representation in the set of nuclear coordinates and electronic states to provide continuous equations of motion that describe the quantized, electronically non-adiabatic dynamics of the system. KC-RPMD preserves the favorable properties of the usual RPMD formulation in the position representation, including rigorous detailed balance, time-reversal symmetry, and invariance of reaction rate calculations to the choice of dividing surface. However, the new method overcomes significant shortcomings of position-representation RPMD by enabling the description of non-adiabatic transitions between states associated with general, many-electron wavefunctions and by accurately describing deep-tunneling processes across asymmetric barriers. We demonstrate that KC-RPMD yields excellent numerical results for a range of model systems, including a simple avoided-crossing reaction and condensed-phase electron-transfer reactions across multiple regimes for the electronic coupling and thermodynamic driving force.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Many-body expansion for light nuclear systems

Here we extend the many-body expansion (MBE), previously applied to hydrogen bonded and molecular systems, to the light nuclear systems 3 H and 3 He by considering the nucleonic degrees of freedom as fundamental in the expansion. The analysis is based on the Pauli nucleonic dynamics (PND) model, a simple antisymmetrized dynamical code, inspired by the sophisticated constrained molecular dynamics (CoMD) model. The total energy of the 2 H nucleus is calculated with this model at -2.312 MeV, which is within 4% of the experimental value of -2.225 MeV. The application of the MBE yields results for the three-body term in the 3 H nucleus that is comparable with previous estimates, while it is reported for the first time for the nucleus of 3 He . The energies of 3 H and 3 He that include the sum of the one- and two-body terms with the model are - 6.97 ± 0.21 MeV and - 6.19 ± 0.21 MeV and the three-body terms, estimated from the difference of the sum of the one- and two-body terms from the experimentally measured energies, are -1.51 ± 0.21 MeV and -1.53 ± 0.21 MeV, respectively. The MBE for the 3 H and 3 He nuclei is qualitatively similar to the one previously reported for the water trimer: the two- and three-body terms are negative with the latter amounting to ~3% of the former. Additionally, the three-body terms correspond to about 16–23% of the total energy of the bound systems, a percentage that is also comparable to the one in the water trimer (17%). In this manner, the MBE analysis can be applied to light nuclear systems following the same protocol as the one that has been previously extensively used for hydrogen bonded molecular systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Coulomb field correction due to virtual $e$ + $e$ - production in heavy ion collisions

The correction to the Coulomb energy due to virtual production of $e$ + $e$ - pairs, which is on the order of one percent of the Coulomb energy at nuclear scales is discussed. The effects of including a pair-production term in the semi-empirical mass formula and the correction to the Coulomb barrier for a handful of nuclear collisions using the Bass and Coulomb potentials are studied. Furthermore, with an eye toward future work using Constrained Molecular Dynamics (CoMD) model, we also calculate the correction to the Coulomb energy and force between protons after folding with a Gaussian spatial distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Correlation between time and angular alignment in molecular dynamics simulations of heavy ion collisions

Neutron-proton equilibration is a process which has been used to study the density dependence of the symmetry energy term in the nuclear equation-of-state. Here, this study utilizes constrained molecular dynamics (CoMD) simulations of 70 Zn + 70 Zn with collision energies of 35 and 45 MeV/nucleon. An algorithm is used which searches through CoMD events and identifies the PLF* after it separates from the target and determines its lifetime, Δt. It also determines the fragments that the PLF* breaks apart into and determines their angular alignment. This technique gives an opportunity to explore how the average alignment of dynamically produced fragments, $\langle$α $\rangle$ dyn , evolves with PLF* lifetime. An approximately linear relationship was determined with d$\langle$α $\rangle$ dyn /dΔt = 0.98 ± 0.08 rad / zs and 1.06 ± 0.09 rad / zs for the 35 and 45 MeV / nucleon, respectively, indicating a correlation with magnitude consistent with classically determined values which were used for prior experimental studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Neutron-rich rare-isotope production from projectile fission of heavy nuclei near 20 MeV/nucleon beam energy

We investigate the possibilities of producing neutron-rich nuclides in projectile fission of heavy beams in the energy range of 20 MeV/nucleon expected from low-energy facilities. We report our efforts to theoretically describe the reaction mechanism of projectile fission following a multinucleon transfer collision at this energy range. Our calculations are mainly based on a two-step approach: The dynamical stage of the collision is described with either the phenomenological deep-inelastic transfer model (DIT) or with the microscopic constrained molecular dynamics model (CoMD). The de-excitation or fission of the hot heavy projectile fragments is performed with the statistical multifragmentation model (SMM). We compared our model calculations with our previous experimental projectile-fission data of 238 U (20 MeV/nucleon) + 208 Pb and 197 Au (20 MeV/nucleon) + 197 Au and found an overall reasonable agreement. Furthermore, our study suggests that projectile fission following peripheral heavy-ion collisions at this energy range offers an effective route to access very neutron-rich rare isotopes toward and beyond the astrophysical r-process path.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Investigating the time dependence of neutron-proton equilibration using molecular dynamics simulations

Neutron-proton equilibration has previously been studied experimentally in dynamically deformed nuclei in heavy ion collisions [Phys. Rev. Lett. 118, 062501 (2017); Phys. Rev. C 95, 044604 (2017)] using the NIMROD detector array at the Cyclotron Institute at Texas A&M University. Results indicated the composition of the two heaviest fragments originating from the excited projectile-like fragment evolved exponential with respect to their orientation angle to be more similar. Constrained Molecular Dynamics and Antisymmetrized Molecular Dynamics simulations were performed for different formulations of the density dependence of the asymmetry energy term of the nuclear equation of state. The simulations are compared to experimental results, which indicate a better agreement with a softer interaction in the nuclear equation of state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

First-Principles Calculation of Water p K a Using the Newly Developed SCAN Functional

Acid/base chemistry is an intriguing topic that still constitutes a challenge for computational chemistry. While estimating the acid dissociation constant (or pK a ) could shed light on many chemistry processes, especially in the fields of biochemistry and geochemistry, evaluating the relative stability between protonated and nonprotonated species is often very difficult. Indeed, a prerequisite for calculating the pK a of any molecule is an accurate description of the energetics of water dissociation. In this study, we applied constrained molecular dynamics simulations, a noncanonical sampling technique, to investigate the water deprotonation process by selecting the OH distance as the reaction coordinate. The calculation is based on density functional theory and the newly developed SCAN functional, which has shown excellent performance in describing water structure. This first benchmark of SCAN on a chemical reaction shows that this functional accurately models the energetics of proton transfer reactions in an aqueous environment. After taking Coulomb long-range corrections and nuclear quantum effects into account, the estimated water pK a is only 1.0 pK a unit different from the target experimental value. Our results show that the combination of SCAN and constrained MD successfully reproduces the chemistry of water and constitutes a good framework for calculating the free energy of chemical reactions of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Thermalization of nuclear matter in heavy-ion collisions at Fermi energies

We analyze the time evolution of the kinetic properties of nuclear matter produced in heavy-ion collisions at Fermi energies. The collision system is simulated using Constrained Molecular Dynamics (CoMD) transport calculations whose output is the isospin, position, and momentum of the nucleons. Focusing on central 35 A·MeV 40 Ca+ 40 Ca collisions we utilize this information to extract localized momentum distributions in volume elements of 8 fm 3 and time steps of 5 fm/$c$. We then parameterize the single-particle momentum distributions with thermally motivated fit functions in the local rest frame of each cell. While the transverse-momentum distributions are well reproduced by thermal ones, the longitudinal ones carry a marked imprint of the initial nuclear motion which we capture by introducing a centroid motion into our fit functions. In particular, we find that Fermi distributions yield significantly better fits than Boltzmann ones, a consequence of the Pauli blocking implemented in CoMD. From the fits we extract the time dependence of the thermodynamic and collective properties of the excited nuclear medium. We find that the transverse temperature gradually rises to about 6 MeV, which is accompanied by a dissipation of the initial centroid motion of the incoming nuclei which vanishes at about 100 fm/c after initial impact. In conclusion, we are therefore able to track the transition of beam energy into random kinetic energy for nucleons, suggesting a three-dimensional equilibration of energy in the late stages of the collision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Insights into lithium ion deposition on lithium metal surfaces

Lithium metal is among the most promising anodes for the next generation of batteries due to its high theoretical energy density and high capacity. Challenges such as extreme reactivity and lithium dendrite formation have kept lithium metal anodes away from practical applications. However, the underlying mechanisms of Li ion deposition from the electrolyte solution onto the anode surface are still poorly understood due to their inherent complexity. In this work, density functional theory calculations and thermodynamic integration via constrained molecular dynamics simulations are conducted to study the electron and ion transfer between lithium metal slab and the electrolyte in absence of an external field. Here, we explore the effect of the solvent chemistry and structure, distance of the solvated complex from the surface, anion–cation separation, and concentration of Li-salts on the deposition of lithium ions from the electrolyte phase onto the surface. Ethylene carbonate (EC), 1,2-dimethoxyethane (DME), 1,3-dioxolane (DOL), and mixtures of them are used as solvents. These species compete with the salt anion and the Li cation for electron transfer from the surface. It is found that the structure and properties of the solvation shell around the lithium cation has a great influence on the ability of the cation to diffuse as well as on its surrounding electron environment. DME molecules allow easier motion of the lithium ion compared with EC and DOL molecules. The slow growth approach allows the study of energy barriers for the ion diffusion and desolvation during the deposition pathway. This method helps elucidating the underlying mechanisms on lithium-ion deposition and provides a better understanding of the early stages of Li nucleation.

25 ENERGY STORAGE↗
DOE OSTI2024

Multinucleon transfer in 48 Ti + 48 Ti collisions at 11.5 MeV/nucleon

In this work, experimental data of projectile fragment distributions from a preliminary experimental study of the reaction of 11.5 MeV/nucleon 48 Ti on 48 Ti analyzed with the MARS recoil separator at the Cyclotron Institute of Texas A&M University are presented. Production cross sections, momentum distributions and excitation-energy distributions are extracted and compared with calculations with the Deep Inelastic Transfer (DIT) model and the Constrained Molecular Dynamics model (CoMD). Despite the limited extent of the present dataset, an overall agreement of the models with the data is found pointing at the prevailing nucleon-exchange character of the collisions at this energy.

Souliotis, Georgios↗
DOE OSTI2022

Investigating the possible existence of hyper-heavy nuclei in a neutron-star environment

The synthesis of hyper-heavy elements is investigated under conditions simulating neutron star environment. The constrained molecular dynamics approach is used to simulate low energy collisions of extremely n-rich nuclei. A new type of the fusion barrier due to a “neutron wind” is observed when the effect of neutron star environment (screening of Coulomb interaction) is introduced implicitly. When introducing also a background of surrounding nuclei, the nuclear fusion becomes possible down to temperatures of 10 8 K and synthesis of extremely heavy and n-rich nuclei appears feasible. A possible existence of hyper-heavy nuclei in a neutron star environment could provide a mechanism of extra coherent neutrino scattering or an additional mechanism, resulting in x-ray burst or a gravitational wave signal and, thus, becoming another crucial process adding new information to the suggested models on neutron star evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Kinetically-constrained ring-polymer molecular dynamics for non-adiabatic chemistries involving solvent and donor–acceptor dynamical effects

We investigate the performance of the recently developed kinetically-constrained ring polymer molecular dynamics (KC-RPMD) method for the description of model condensed-phase electron transfer (ET) reactions in which solvent and donor–acceptor dynamics play an important role. Comparison of KC-RPMD with results from Golden-Rule rate theories and numerically exact quantum dynamics calculations demonstrates that KC-RPMD accurately captures the combination of electronic- and nuclear-dynamical effects throughout the Marcus (intermediate solvent friction) and Zusman (large solvent friction) regimes of ET. It is also demonstrated that KC-RPMD accurately describes systems in which the magnitude of the diabatic coupling depends on the position of a dynamical donor–acceptor mode. In addition to these successes, however, we present an unsurprising failure of KC-RPMD to capture the enhancement of the ET rate in the low solvent friction regime associated with nuclear coherence effects. In this analysis, we re-visit several aspects of the original KC-RPMD formulation, including the form of the kinetic constraint and the choice of the mass of the auxiliary electronic variable. In particular, we introduce a Langevin bath for the auxiliary electronic variable to correct for its unphysically low coupling with the nuclear degrees of freedom. This work demonstrates that the KC-RPMD method is well suited for the direct simulation of non-adiabatic donor–acceptor chemistries associated with many complex systems, including those for which solvent dynamics plays an important role in the reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Overcoming potential energy distortions in constrained internal coordinate molecular dynamics simulations

The Internal Coordinate Molecular Dynamics (ICMD) method is an attractive molecular dynamics (MD) method for studying the dynamics of bonded systems such as proteins and polymers. It offers a simple venue for coarsening the dynamics model of a system at multiple hierarchical levels. For example, large scale protein dynamics can be studied using torsional dynamics, where large domains or helical structures can be treated as rigid bodies and the loops connecting them as flexible torsions. ICMD with such a dynamic model of the protein, combined with enhanced conformational sampling method such as temperature replica exchange, allows the sampling of large scale domain motion involving high energy barrier transitions. Once these large scale conformational transitions are sampled, all-torsion, or even all-atom, MD simulations can be carried out for the low energy conformations sampled via coarse grained ICMD to calculate the energetics of distinct conformations. Such hierarchical MD simulations can be carried out with standard all-atom forcefields without the need for compromising on the accuracy of the forces. Using constraints to treat bond lengths and bond angles as rigid can, however, distort the potential energy landscape of the system and reduce the number of dihedral transitions as well as conformational sampling. We present here a two-part solution to overcome such distortions of the potential energy landscape with ICMD models. To alleviate the intrinsic distortion that stems from the reduced phase space in torsional MD, we use the Fixman compensating potential. To additionally alleviate the extrinsic distortion that arises from the coupling between the dihedral angles and bond angles within a force field, we propose a hybrid ICMD method that allows the selective relaxing of bond angles. This hybrid ICMD method bridges the gap between all-atom MD and torsional MD. We demonstrate with examples that these methods together offer a solution to eliminate the potential energy distortions encountered in constrained ICMD simulations of peptide molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗