Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Structure-Based Identification of Novel Histone Deacetylase 4 (HDAC4) Inhibitors

Histone deacetylases (HDACs) are important cancer drug targets. Existing FDA-approved drugs target the catalytic pocket of HDACs, which is conserved across subfamilies (classes) of HDAC. However, engineering specificity is an important goal. Herein, we use molecular modeling approaches to identify and target potential novel pockets specific to Class IIA HDAC-HDAC4 at the interface between HDAC4 and the transcriptional corepressor component protein NCoR. These pockets were screened using an ensemble docking approach combined with consensus scoring to identify compounds with a different binding mechanism than the currently known HDAC modulators. Binding was compared in experimental assays between HDAC4 and HDAC3, which belong to a different family of HDACs. HDAC4 was significantly inhibited by compound 88402 but not HDAC3. Two other compounds (67436 and 134199) had IC50 values in the low micromolar range for both HDACs, which is comparable to the known inhibitor of HDAC4, SAHA (Vorinostat). However, both of these compounds were significantly weaker inhibitors of HDAC3 than SAHA and thus more selective, albeit to a limited extent. Five compounds exhibited activity on human breast carcinoma and/or urothelial carcinoma cell lines. The present result suggests potential mechanistic and chemical approaches for developing selective HDAC4 modulators.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Design and simulation of n -type solar cells based on an iodine-doped CdTe absorber using SCAPS-1D

The performance of conventional p-type CdTe solar cells has plateaued in recent years, motivating exploration of n-type absorbers. Here, we evaluate a homojunction solar cell employing iodine-doped CdTe (CdTe:I) as the absorber in a Ti 3 C 2 T x MXene/p-CdTe:As/n-CdTe:I/indium structure through SCAPS-1D simulations and prototype devices. Optimized simulations predict efficiencies above 25% for thin CdTe:I absorbers (~0.7 μm). In contrast, the first prototype achieved only ~1.36% efficiency with V OC = 0.48 V, J SC = 6.45 mA/cm 2 and FF = 43.8%. When the simulation is adjusted to match the actual device structure, and the effective illumination is reduced to account for front-side light loss, the predicted V OC (0.48 V) and JSC (6.74 mA/cm 2 ) closely reproduce the experimental values. This suggests that the device performance is primarily limited by reduced front-side photon transmission and other material non-idealities. These results highlight the promise of iodine-doped n-type CdTe and identify clear pathways for further efficiency improvement.

14 SOLAR ENERGY↗
DOE OSTI2021

CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. Additionally, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. Finally, the capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.

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

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2019

The Structural, Electronic, and Magnetic Properties of Binary Heusler Alloys ZCl{sub 3} (Z = Li, Na, K, Rb) with DO{sub 3}-Type Structure from First Principle Calculations

The structural, electronic, and magnetic properties of binary Heusler alloys ZCl{sub 3} (Z = Li, Na, K, Rb) with DO{sub 3}-type structure have been studied by means of the first principle calculations. The calculated results reveal that all four alloys are half-metallic (HM) materials at the stable ferromagnetic ground state. The total magnetic moments (M{sub t}) of LiCl{sub 3}, NaCl{sub 3}, KCl{sub 3,} and RbCl{sub 3} alloys are all 2.000 μ{sub B}/f.u., following the Slater-Pauling rule M{sub t} = 24 − Z{sub t} for the former two alloys, while satisfying the Slater-Pauling rule M{sub t} = 30 − Z{sub t} for the latter two alloys. The spin-up band gaps of the binary Heusler alloys ZCl{sub 3} (Z = Li, Na, K, Rb) are 5.598, 5.687, 6.744, and 6.745 eV at equilibrium lattice constants 7.076, 7.237, 7.588, and 7.819 Å, respectively. The positive cohesion energy and negative formation energy indicate that they are thermodynamically stable. Moreover, wide HM regions for binary Heusler alloys ZCl{sub 3} (Z = Li, Na, K, Rb) show their HM characters are robust when the lattices are expanded or compressed under uniform and tetragonal strains. Finally, the total charge density, magnetic charge density, and charge density difference of the (110) plane for LiCl{sub 3} alloy are analyzed in detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Pairing a Global Optimization Algorithm with EXAFS to Characterize Lanthanide Structure in Solution

Ensemble-average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. Here, to accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO 3 ) 3 dissolved in acetonitrile with a terpyridine ligand (terpyNO 2 ), and Nd(NO 3 ) 3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO 3 ) 3 (terpyNO 2 ) and Nd(NO 3 )(acetonitrile) 3 , with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.

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

Local Order in Liquid Gallium–Indium Alloys

Liquid metals such as eutectic Ga–In alloys have low melting points and low toxicity and are used in catalysis and micro-robotics. This study investigates the local atomic structure of liquid gallium-indium alloys by a combination of density measurements, diffraction data, and Monte-Carlo simulation via the empirical potential structure refinement approach. A high-Q shoulder observed in liquid Ga is related to structural rearrangements in the second coordination shell. Structure analysis found coordination environments close to a random distribution for eutectic Ga–In alloy, while electronic effects appear to dominate the mixing enthalpy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine-learned interatomic potentials by active learning: amorphous and liquid hafnium dioxide

We propose an active learning scheme for automatically sampling a minimum number of uncorrelated configurations for fitting the Gaussian Approximation Potential (GAP). Our active learning scheme consists of an unsupervised machine learning (ML) scheme coupled with a Bayesian optimization technique that evaluates the GAP model. We apply this scheme to a Hafnium dioxide (HfO 2 ) dataset generated from a "melt-quench" ab initio molecular dynamics (AIMD) protocol. Our results show that the active learning scheme, with no prior knowledge of the dataset, is able to extract a configuration that reaches the required energy fit tolerance. Further, molecular dynamics (MD) simulations performed using this active learned GAP model on 6144 atom systems of amorphous and liquid state elucidate the structural properties of HfO 2 with near ab initio precision and quench rates (i.e., 1.0 K/ps) not accessible via AIMD. The melt and amorphous X-ray structural factors generated from our simulation are in good agreement with experiment. In addition, the calculated diffusion constants are in good agreement with previous ab initio studies.

36 MATERIALS SCIENCE↗
DOE OSTI2017

The Implications of 3D Thermal Structure on 1D Atmospheric Retrieval

Using the atmospheric structure from a 3D global radiation-hydrodynamic simulation of HD 189733b and the open-source Bayesian Atmospheric Radiative Transfer (BART) code, we investigate the difference between the secondary-eclipse temperature structure produced with a 3D simulation and the best-fit 1D retrieved model. Synthetic data are generated by integrating the 3D models over the Spitzer , the Hubble Space Telescope ( HST ), and the James Web Space Telescope ( JWST ) bandpasses, covering the wavelength range between 1 and 11 μ m where most spectroscopically active species have pronounced features. Using the data from different observing instruments, we present detailed comparisons between the temperature–pressure profiles recovered by BART and those from the 3D simulations. We calculate several averages of the 3D thermal structure and explore which particular thermal profile matches the retrieved temperature structure. We implement two temperature parameterizations that are commonly used in retrieval to investigate different thermal profile shapes. To assess which part of the thermal structure is best constrained by the data, we generate contribution functions for our theoretical model and each of our retrieved models. Our conclusions are strongly affected by the spectral resolution of the instruments included, their wavelength coverage, and the number of data points combined. We also see some limitations in each of the temperature parametrizations, as they are not able to fully match the complex curvatures that are usually produced in hydrodynamic simulations. The results show that our 1D retrieval is recovering a temperature and pressure profile that most closely matches the arithmetic average of the 3D thermal structure. When we use a higher resolution, more data points, and a parametrized temperature profile that allows more flexibility in the middle part of the atmosphere, we find a better match between the retrieved temperature and pressure profile and the arithmetic average. The Spitzer and HST simulated observations sample deep parts of the planetary atmosphere and provide fewer constraints on the temperature and pressure profile, while the JWST observations sample the middle part of the atmosphere, providing a good match with the middle and most complex part of the arithmetic average of the 3D temperature structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Active Sites in the Dealuminated Beta Zeolite-Supported Cobalt Catalyst for Non-Oxidative Ethane Dehydrogenation

Dispersed metal species in siliceous zeolites have been actively studied for non-oxidative dehydrogenation of ethane (NDE). Fundamental insights into the dynamics of metal species in zeolites under reaction conditions have rarely been explored. Herein, we report an atomic level understanding of the dynamics and activity of cobalt (Co) sites in dealuminated Beta zeolite (DeAl-BEA) for NDE during induction and reaction conditions with extensive characterization techniques such as diffuse reflectance UV–vis, solid state nuclear magnetic resonance and X-ray photoelectron, X-ray diffraction along with in situ Fourier transform infrared and X-ray absorption spectroscopy. For a catalyst with 0.5 mass % Co loading, tetrahedral Co 2+ mononuclear sites, di-coordinated to the zeolite framework and with two silanol groups in vicinity (i.e., (≡SiO) 2 Co(HO–Si≡) 2 ), form upon exposure to hydrogen during induction and persist through the NDE reaction. Increasing the Co loading to 3.0 mass % yielded Co sites with similar electronic and coordination structures but slightly elongated Co–O bonds. Upon cooling to room temperature, the Co sites persisted in the same coordination environment, though the disappearance of a feature in the Co K-edge near-edge region revealed changes in the active site’s electronic structure coinciding with modest shifts in bond lengths. The electronic structure and activity of (≡SiO) 2 Co(HO–Si≡) 2 sites were studied comparatively to a few other hypothetical Co 2+ coordination structures, using electronic structure calculations and microkinetic simulations. The simulations showed that NDE is controlled by β-hydride elimination following C–H bond activation and that Co-sites possessing flexibility because of neighboring silanol defects are more active. Interestingly, dinuclear Co–O–Co sites (i.e., (≡SiO)Co(HO–Si≡) 2 –O–(HO–Si≡) 2 Co(≡SiO)) were more active than the mononuclear (≡SiO) 2 Co(HO–Si≡) 2 sites because of favorable hydrogen bonding with the vicinal silanol groups. In conclusion, the present study bridges the gap between the knowledge acquired by ex-situ characterizations and the active sites under the reaction conditions in alkane dehydrogenation chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Molecular Structures of Liquid and Glassy Nifedipine and Felodipine and Their Incorporation into PVP

Background: Amorphous drug formulations are commonly used to improve the solubility and bioavailability of poorly soluble molecular pharmaceuticals, yet less is known about their molecular conformations and local bonding interactions than their crystalline phases. Methods: High-energy X-ray diffraction structure factor measurements have been made on liquid and glassy nifedipine (NIF), felodipine (FEL), NIF 1:3 polyvinylpyrrolidone (PVP), and FEL 1:3 PVP wt.% mixtures. The corresponding X-ray pair distribution functions have been interpreted using empirical potential structure refinement using different models and density functional theory conformer calculations. Results: In both NIF and FEL, the NH···O inter-molecular hydrogen bonds between the pyridyl nitrogen and ester carbonyls are found to be considerably weaker than those observed in the crystalline polymorphs. For nifedipine, it is proposed that either inter-molecular NH…ON nitro bonds are present and/or a fraction (<20%) of conformational changes, with the aryl ring flipped, occur in the liquid state. For felodipine, the models indicate significant disorder associated with the methyl and ethyl side chains in the liquid state, with the main peak intensity at 3.0 Å arising from intra-molecular Cl-Cl atom pairs. When nifedipine molecules are incorporated into PVP, our models show they possess stronger NH···O bonds to the PVP polymer than felodipine molecules, which have stronger affinity for bonding to the polymer than to other felodipine molecules. Conclusions: The amorphous forms of both NIF and FEL show much weaker hydrogen bonding than found in their crystalline phases. Liquid NIF also exhibits configurations which are not observed in the crystal phases.

Amorphous Solid Dispersion↗
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Simulation of the Formation of a Cascade of Displacements and Transient Ionization Processes in Silicon Semiconductor Structures under Neutron Exposure

The formation of a disordered defect region in bulk silicon is simulated using the molecular-dynamics method for various energies of a primary recoil atom. Variations in the volume and number of radiation-induced defects in a cluster during its formation are calculated. The generation rates of nonequilibrium carriers and amplitude-temporal dependences of pulses of ionization currents in test Schottky diodes with hyperhigh frequencies are found theoretically.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Learning to predict the cosmological structure formation

Matter evolved under the influence of gravity from minuscule density fluctuations. Nonperturbative structure formed hierarchically over all scales and developed non-Gaussian features in the Universe, known as the cosmic web. To fully understand the structure formation of the Universe is one of the holy grails of modern astrophysics. Astrophysicists survey large volumes of the Universe and use a large ensemble of computer simulations to compare with the observed data to extract the full information of our own Universe. However, to evolve billions of particles over billions of years, even with the simplest physics, is a daunting task. We build a deep neural network, the Deep Density Displacement Model ( D 3 M ), which learns from a set of prerun numerical simulations, to predict the nonlinear large-scale structure of the Universe with the Zel’dovich Approximation (ZA), an analytical approximation based on perturbation theory, as the input. Our extensive analysis demonstrates that D 3 M outperforms the second-order perturbation theory (2LPT), the commonly used fast-approximate simulation method, in predicting cosmic structure in the nonlinear regime. We also show that D 3 M is able to accurately extrapolate far beyond its training data and predict structure formation for significantly different cosmological parameters. Our study proves that deep learning is a practical and accurate alternative to approximate 3D simulations of the gravitational structure formation of the Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Nanoscale structure of the oil-water interface

X-ray reflectivity (XR) and atomistic molecular dynamics (MD) simulations, carried out to determine the structure of the oil-water interface, provide new insight into the simplest liquid-liquid interface. For several oils (hexane, dodecane, and hexadecane) the XR shows very good agreement with a monotonic interface-normal electron density profile (EDP) broadened only by capillary waves. Similar agreement is also found for an EDP including a sub-Å thick electron depletion layer separating the oil and the water. As a result, the XR and MD derived depletions are much smaller than reported for the interface between solid-supported hydrophobic monolayers and water.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Explicit Filtering and Reconstruction to Reduce Grid Dependence in Convective Boundary Layer Simulations Using WRF-LES

As model grid resolutions move from the mesoscale to the microscale, turbulent structures represented in atmospheric boundary layer simulations change dramatically. At intermediate resolutions, the so-called gray zone, turbulent motions are not resolved accurately, posing a challenge to numerical simulations. The representation of turbulence is also highly sensitive to the choice of closure model. Here, we examine explicit filtering and reconstruction in the gray zone as a technique to better represent atmospheric turbulence. The convective boundary layer is simulated using the Weather Research and Forecasting (WRF) Model with horizontal resolutions ranging from 25 m to 1 km. Four large-eddy simulation (LES) turbulence models are considered: the Smagorinsky model, the TKE-1.5 model, and two versions of the dynamic reconstruction model (DRM). The models are evaluated by their ability to produce consistent mean potential temperature profiles, heat and momentum fluxes, velocity fields, and turbulent kinetic energy spectra as the grids become coarser. The DRM, a mixed model that uses an explicit filtering and reconstruction technique to account for resolvable subfilter-scale (RSFS) stresses, performs very well at resolutions of 500 m and 1 km without any special tuning, whereas the Smagorinsky and TKE-1.5 models produce heavily grid-dependent results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Local structure orders and glass forming ability of Ni-Nb liquids

In this paper, ab initio molecular dynamics simulations were performed to study the atomistic structures of Ni 48 Nb 52 , Ni 59.5 Nb 40.5 and Ni 75 Nb 25 liquids near their melting temperatures. It is found that the total coordination number only changes slightly across the compositions. However, there is significant difference in the chemical ordering. The Cluster Alignment method was used to identify the dominant local structure motifs. It is found that distorted and perfect icosahedral clusters are the most popular motifs around Ni atoms, while Z14, Z15, and Z16 Frank-Kasper polyhedra are dominant around Nb atoms. Among the three compositions, the eutectic composition (Ni 59.5 Nb 40.5 ) has the largest population of these dominant motifs and the highest degree of five-fold local symmetry. Chemical order in the distorted icosahedral motif is studied and the occupation probability of Nb in the distorted region is found to be much different from other sites. The network formed by distorted and perfect icosahedral clusters and that by the interconnection of the perfect icosahedra, Z14, Z15, and Z16 motifs are the strongest at the eutectic composition. Finally, the high percentage of perfect and distorted icosahedral clusters and their strong network structures in this system (especially around the eutectic composition) correlate well with the excellent glass forming ability in this system.

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

Enhancement of nuclear reaction rates in asymmetric binary ionic mixtures

Here, using orbital-free molecular dynamics simulations we study the structure and dynamics of increasingly asymmetric mixtures such as hydrogen–carbon, hydrogen–aluminium, hydrogen–copper, and hydrogen–silver. We show that, whereas the heavy component structure is close to an effective one-component plasma (OCP), the light component appears more structured than the corresponding OCP. This effect is related to the crossover towards a Lorentz-type diffusion triggered by strongly coupled, highly charged heavy ions, and witnessed by the change of temperature scaling laws of diffusion. Finally, this over-correlation translates into an enhancement of nuclear reaction rates much higher than its classical OCP counterpart.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗