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 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 OSTI2026

From 2D to 4D: a containerized workflow and browser to explore dynamic chromatin architecture

Background Characterizing the physical organization of the genome is essential for understanding long-range gene regulation, chromatin compartmentalization, and epigenetic accessibility. Hi-C experiments generate two-dimensional (2D) genome-wide contact maps of chromatin interactions by capturing the spatial proximity between genomic loci, which reveal interaction frequencies but lack the spatial resolution needed to interpret the three-dimensional (3D) genome structure(s). Emerging evidence suggests that epigenetic regulation is closely linked to 3D genome architecture, and that structural changes over time (4D) drive key biological processes in development, disease, and environmental response. Thus, integrating 3D structure with functional data is critical for a more complete understanding of genome regulation. Previous work, most notably the 4DHiC chromosome modeling framework, has shown that physical multi-dimensional modeling approaches rooted in polymer physics and molecular dynamics can resolve these structures at biologically meaningful resolutions by integrating temporal Hi-C data with physical constraints to uncover dynamic chromosome reorganization. Thus, molecular dynamics simulations, constrained by Hi-C contact matrices, can resolve fine-scale structural changes and reveal functionally significant transitions in chromatin conformation. Results Herein, we present the 4D Genome Browser Workflow (4DGBWorkflow) and the 4D Genome Browser (4DGB). The algorithm is based on the 4DHiC method, and the containerized tool is an end-to-end workflow that can transform, filter, and view 4D epigenomics and chromatin datasets, allowing non-specialists to apply three-dimensional modeling principles to diverse datasets and experimental conditions. The software executes on a laptop running macOS, Linux or Windows. From input Hi-C files (.hic), the 4DGBWorkflow produces 3D reconstructions of chromosomes, integrates the reconstruction with track data (e.g., epigenetic marks, transcriptome profiles), and provides comparative visualization of the results in a single workflow. Conclusions The 4DGBWorkflow and 4D Genome Browser are open-source tools for comparative analysis and visualization of 4D chromosome datasets, including chromatin architecture and epigenomic signals. Automatic integration of Hi-C data with molecular dynamics democratizes the construction of time resolved 3D genome structures, simplifying complex simulations and data integration schemes.

3D Genome Browser↗
DOE OSTI2020

Photoinduced Dynamics with Constrained Vibrational Motion: FrozeNM Algorithm

Ab initio molecular dynamics (AIMD) simulation, analyzed in terms of vibrational normal modes, is a widely used technique that facilitates understanding of complex structural motions and coupling between electronic and nuclear degrees of freedom. Usually, only a subset of vibrations is directly involved in the process of interest. The impact of these vibrations can be evaluated by performing AIMD simulations by selectively freezing certain motions. In this paper, we present frozen normal mode (FrozeNM), a new algorithm to apply normal-mode constraints in AIMD simulations, as implemented in the nonadiabatic excited state molecular dynamics code. We further illustrate its capacity by analyzing the impact of normal-mode constraints on the photoinduced energy transfer between polyphenylene ethynylene dendrimer building blocks. Our results show that the electronic relaxation can be significantly slowed down by freezing a well-selected small subset of active normal modes characterized by their contributions in the direction of energy transfer. The application of these constraints reduces the nonadiabatic coupling between electronic excited states during the entire dynamical simulations. Furthermore, we validate reduced dimensionality models by freezing all the vibrations, except a few active modes. Altogether, we consider FrozeNM as a useful tool that can be broadly used to underpin the role of vibrational motion in a studied process and to formulate reduced models that describe essential physical phenomena.

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

Ab initio molecular dynamics of solvation effects on reactivity at electrified interfaces

Using ab initio molecular dynamics as implemented in periodic, self-consistent (generalized gradient approximation Perdew–Burke–Ernzerhof) density functional theory, we investigated the mechanism of methanol electrooxidation on Pt(111). In this work, we investigated the role of water solvation and electrode potential on the energetics of the first proton transfer step, methanol electrooxidation to methoxy (CH 3 O) or hydroxymethyl (CH 2 OH). The results show that solvation weakens the adsorption of methoxy to uncharged Pt(111), whereas the binding energies of methanol and hydroxymethyl are not significantly affected. The free energies of activation for breaking the C–H and O–H bonds in methanol were calculated through a Blue Moon Ensemble using constrained ab initio molecular dynamics. Calculated barriers for these elementary steps on unsolvated, uncharged Pt(111) are similar to results for climbing-image nudged elastic band calculations from the literature. Water solvation reduces the barriers for both C–H and O–H bond activation steps with respect to their vapor-phase values, although the effect is more pronounced for C–H bond activation, due to less disruption of the hydrogen bond network. The calculated activation energy barriers show that breaking the C–H bond of methanol is more facile than the O–H bond on solvated negatively biased or uncharged Pt(111). However, with positive bias, O–H bond activation is enhanced, becoming slightly more facile than C–H bond activation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗