Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

New approach for precise computation of Lyman-α forest power spectrum with hydrodynamical simulations

Current experiments are providing measurements of the flux power spectrum from the Lyman-α forests observed in quasar spectra with unprecedented accuracy. Their interpretation in terms of cosmological constraints requires specific simulations of at least equivalent precision. In this paper, we present a suite of cosmological N-body simulations with cold dark matter and baryons, specifically aiming at modeling the low-density regions of the inter-galactic medium as probed by the Lyman-α forests at high redshift. The simulations were run using the GADGET-3 code and were designed to match the requirements imposed by the quality of the current SDSS-III/BOSS or forthcoming SDSS-IV/eBOSS data. They are made using either 2 × 768{sup 3} ≅ 1 billion or 2 × 192{sup 3} ≅ 14 million particles, spanning volumes ranging from (25 Mpc h{sup −1}){sup 3} for high-resolution simulations to (100 Mpc h{sup −1}){sup 3} for large-volume ones. Using a splicing technique, the resolution is further enhanced to reach the equivalent of simulations with 2 × 3072{sup 3} ≅ 58 billion particles in a (100 Mpc h{sup −1}){sup 3} box size, i.e. a mean mass per gas particle of 1.2 × 10{sup 5}M{sub ⊙} h{sup −1}. We show that the resulting power spectrum is accurate at the 2% level over the full range from a few Mpc to several tens of Mpc. We explore the effect on the one-dimensional transmitted-flux power spectrum of four cosmological parameters (n{sub s}, σ{sub 8}, Ω{sub m} and H{sub 0}) and two astrophysical parameters (T{sub 0} and γ) that are related to the heating rate of the intergalactic medium. By varying the input parameters around a central model chosen to be in agreement with the latest Planck results, we built a grid of simulations that allows the study of the impact on the flux power spectrum of these six relevant parameters. We improve upon previous studies by not only measuring the effect of each parameter individually, but also probing the impact of the simultaneous variation of each pair of parameters. We thus provide a full second-order expansion, including cross-terms, around our central model. We check the validity of the second-order expansion with independent simulations obtained either with different cosmological parameters or different seeds. Finally, a comparison to the one-dimensional Lyman-α forest power spectrum obtained with BOSS by [1] shows an excellent agreement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Characterization of the Redesigned Extended MITL for HERMES III

Modeling and simulation of the legacy HERMES III Magnetically Insulated Transmission Line (MITL) has been performed using EMPHASIS, an unstructured time-domain electromagnetic (UTDEM) particle-in-cell (PIC) simulation software. This design when used lost roughly half of its current before the pulse reached the load. The cause of the current loss in the MITL was found to be the vacuum impedance changes along the MITL. The MITL was then redesigned to maintain constant impedance and simulated in EMPHASIS once again. Following predicting simulation results, the new MITL was then built, installed, and tested, showing minimal current loss and good agreement with simulation and theoretical results, all of which are reported here. Additionally, an analysis of experimental voltage calculation techniques using cathode and anode currents is performed and compared to simulation results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Estimation of High-Altitude Electromagnetic Pulse Signal Leakage into Power Generation Facilities: Simulations and Measurements

The coupling of early-time high-altitude electromagnetic pulse signals into a power generation facility was investigated. This study encompasses both electromagnetic simulation and on-site experimental efforts that are relevant to shielding effectiveness estimation with a passive technique. As an initial step in understanding the attenuation models involved, simulations were carried out as a function of frequency, building properties, and signal incidence angle. Testing was then performed at the US Department of Energy’s Oak Ridge National Laboratory campus to identify the equipment needed for signal collection and to fine-tune the experimental procedure. Finally, as an application of the models and methods presented, measurements were conducted at a power generation plant to derive the field transfer functions and coupled high-altitude electromagnetic pulse waveforms for various locations within the facility.

42 ENGINEERING↗
DOE OSTI2016

Solution-state structure and affinities of cyclodextrin: Fentanyl complexes by nuclear magnetic resonance spectroscopy and molecular dynamics simulation

Cyclodextrins (CDs) are investigated for their ability to form inclusion complexes with the analgesic fentanyl and three similar molecules: acetylfentanyl, thiofentanyl, and acetylthiofentanyl. Stoichiometry, binding strength, and complex structure are revealed through nuclear magnetic resonance (NMR) techniques and discussed in terms of molecular dynamics (MD) simulations. It was found that β-cyclodextrin is generally capable of forming the strongest complexes with the fentanyl panel. Two-dimensional NMR data and computational chemical calculations are used to derive solution-state structures of the complexes. Binding of the fentanyls to the CDs occurs at the amide phenyl ring, leaving the majority of the molecule solvated by water, an observation common to all four fentanyls. This finding suggests a universal binding behavior, as the vast majority of previously synthesized fentanyl analogues contain this structural moiety. Furthermore, this baseline study serves as the most complete work on CD:fentanyl complexes to date and provides the insights into strategies for producing future generations of designer cyclodextrins capable of stronger and more selective complexation of fentanyl and its analogues.

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

MO-D-PinS Room/Hall E-02: Multiparametric Fuctional Imaging in Radiation Therapy

MRI, with its excellent soft tissue contrast and its ability to provide physiological as well as anatomical information, is becoming increasingly used in radiation therapy for treatment planning, image-guided radiation therapy, and treatment evaluation. This session will explore solutions to integrating MRI into the simulation process. Obstacles for using MRI for simulation include distortions and artifacts, image acquisition speed, complexity of imaging techniques, and lack of electron density information. Partners in Solutions presents vendor representatives who will present their approaches to meeting these challenges and others. An increased awareness of how MRI simulation works will allow physicists to better understand and use this powerful technique. The speakers are all employees who are presenting information about their company’s products.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Quantum Algorithm for High Energy Physics Simulations

Simulating quantum field theories is a flagship application of quantum computing. However, calculating experimentally relevant high energy scattering amplitudes entirely on a quantum computer is prohibitively difficult. It is well known that such high energy scattering processes can be factored into pieces that can be computed using well established perturbative techniques, and pieces which currently have to be simulated using classical Markov chain algorithms. These classical Markov chain simulation approaches work well to capture many of the salient features, but cannot capture all quantum effects. To exploit quantum resources in the most efficient way, we introduce a new paradigm for quantum algorithms in field theories. This approach uses quantum computers only for those parts of the problem which are not computable using existing techniques. In particular, we develop a polynomial time quantum final state shower that accurately models the effects of intermediate spin states similar to those present in high energy electroweak showers with a global evolution variable. The algorithm is explicitly demonstrated for a simplified quantum field theory on a quantum computer.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Structure and dynamics of ethane confined in silica nanopores in the presence of CO 2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C 2 H 6 ) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C 2 H 6 and CO 2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C 2 H 6 , C 2 H 6 :CO 2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C 2 H 6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C 2 H 6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO 2 molecules are more strongly attracted to the pore surface compared to C 2 H 6 . The C 2 H 6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO 2 molecules play in enhancing C 2 H 6 diffusion (“molecular lubrication”) at high CO 2 :C 2 H 6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

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

FIELD IMPLEMENTATION PLAN FOR A WILLISTON BASIN BRINE EXTRACTION AND STORAGE TEST

The Energy & Environmental Research Center (EERC) successfully completed all technical work of Phase I, including development of a field implementation plan (FIP) for a brine extraction and storage test (BEST) in the North Dakota portion of the Williston Basin. This implementation plan was commissioned by the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL) as a proxy for managing formation pressure plumes and measuring/monitoring the movement of differential pressure and CO2 plumes in the subsurface for future saline CO2 storage projects. BEST comprises the demonstration and validation of active reservoir management (ARM) strategies and extracted brine treatment technologies. Two prospective commercial brine injection sites were evaluated for BEST to satisfy DOE’s goals. Ultimately, an active saltwater disposal (SWD) site, Johnsons Corner, was selected because it possesses an ideal combination of key factors making it uniquely suited to host BEST. This site is located in western North Dakota and operated by Nuverra Environmental Solutions (Nuverra), a national leader in brine handling, treatment, and injection. An integrated management approach was used to incorporate local and regional geologic characterization activities with geologic and simulation models, inform a monitoring, verification, and accounting (MVA) plan, and to conduct a risk assessment. This approach was used to design a FIP for an ARM schema and an extracted brine treatment technology test bed facility. The FIP leverages an existing pressure plume generated by two commercial SWD wells. These wells, in conjunction with a new brine extraction well, will be used to conduct the ARM schema. Results of these tests will be quantified based on their impact on the performance of the existing SWD wells and the surrounding reservoir system. Extracted brine will be injected into an underlying deep saline formation through a new injection well. The locations of proposed extraction and injection wells were selected during the Phase I efforts. These wells will be permitted as North Dakota Administrative Code Underground Injection Control Class II wells and will yield additional characterization data which will further refine the FIP in Phase II. An array of surface and downhole monitoring techniques will validate ARM performance against predictive simulation results. Infrastructure will be constructed to manage extracted fluids at the surface and provide brine to a treatment test bed facility. Treatment of extracted brine can provide a means of reducing extracted brine disposal volumes, an alternate source of water, and/or salable products for beneficial use. A test bed facility will be constructed to provide a means of demonstrating these technologies on a wide range of brine concentrations. Screening criteria based on a techno-economic and life cycle assessment were developed to select high-salinity brine treatment technologies for extended duration treatment (30–60 days) in Phase II. A detailed cost assessment determined total implementation costs for BEST of $19,901,065 million (DOE share $15,680,505). These costs are inclusive of all necessary equipment, infrastructure construction, operations and project closeout costs required to implement BEST. An ideal combination of key factors makes the Johnsons Corner site uniquely suited to be the BEST demonstration.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Augmented Kalman filter with a reduced mechanical model to estimate tower loads on a land-based wind turbine: a step towards digital-twin simulations

This article presents an application of the Kalman filtering technique to estimate loads on a wind turbine. The approach combines a mechanical model and a set of measurements to estimate signals that are not available in the measurements, such as wind speed, thrust, tower position, and tower loads. The model is severalfold faster than real time and is intended to be run online, for instance, to evaluate real-time fatigue life consumption of a field turbine using a digital twin, perform condition monitoring, or assess loads for dedicated control strategies. The mechanical model is built using a Rayleigh–Ritz approach and a set of joint coordinates. We present a general method and illustrate it using a 2-degrees-of-freedom (DOF) model of a wind turbine and using rotor speed, generator torque, pitch, and tower-top acceleration as measurement signals. The different components of the model are tested individually. The overall method is evaluated by computing the errors in estimated tower-bottom-equivalent moment from a set of simulations. From this preliminary study, it appears that the tower-bottom-equivalent moment is obtained with about 10 % accuracy. The limitation of the model and the required steps forward are discussed.

17 WIND ENERGY↗
DOE OSTI2025

Accurate and efficient predictions of keyhole dynamics in laser materials processing using machine learning-aided simulations

The keyhole phenomenon has been widely observed in laser materials processing, including laser welding, remelting, cladding, drilling, and additive manufacturing. Keyhole-induced defects, primarily pores, dramatically affect the performance of final products, impeding the broad use of these laser-based technologies. The formation of these pores is typically associated with the dynamic behavior of the keyhole. So far, the accurate characterization and prediction of keyhole features, particularly keyhole depth, as a function of time, has been a challenging task. In situ characterization of keyhole dynamic behavior using the synchrotron X-ray technique is informative but complicated and expensive. Current simulations are generally hindered by their poor accuracy and generalization abilities in predicting keyhole depths due to the lack of accurate laser absorptance data. In this study, we develop a machine learning-aided simulation method that accurately predicts keyhole dynamics, especially in keyhole depth fluctuations, over a wide range of processing parameters. In two case studies involving titanium and aluminum alloys, we achieve keyhole depth prediction with a mean absolute percentage error of 10 %, surpassing those simulated using the ray-tracing method with an error margin of 30 %, while also reducing computational time. This exceptional fidelity and efficiency empower our model to serve as a cost-effective alternative to synchrotron experiments. Our machine learning-aided simulation method is affordable and readily deployable for a large variety of materials, opening new doors to eliminate or reduce defects for a wide range of laser materials processing techniques.

Computational fluid dynamics↗
DOE OSTI2013

Closely spaced nanomagnets by dual e-beam exposure for low-energy nanomagnet logic

The effect of nanomagnet spacing on required clock field has been studied by micromagnetic simulation for supermalloy (Ni{sub 79}Fe{sub 16}Mo{sub 5}) dots with dimensions 90 Multiplication-Sign 60 Multiplication-Sign 20 nm{sup 3} and 120 Multiplication-Sign 60 Multiplication-Sign 20 nm{sup 3}. Reduction of the inter-magnet spacing for both dimensions has resulted in reduction of the required clock field in the simulation. A dual e-beam exposure technique has been developed to allow fabrication of ultra dense features using conventional poly(methylmethacrylate) e-beam resist. Nanomagnet logic (NML) datalines of supermalloy dots with {approx}10 nm and {approx}15 nm spacing have been fabricated using dual e-beam exposure with a 3{sigma} overlay accuracy of {approx}4 nm. Fabricated NML datalines have been characterized using magnetic force microscopy for various clock fields. Datalines of both spacing have shown proper NML functionality with a clock field as low as 60 mT.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase Space Reconstruction from Accelerator Beam Measurements Using Neural Networks and Differentiable Simulations

Characterizing the phase space distribution of particle beams in accelerators is a central part of accelerator understanding and performance optimization. However, conventional reconstruction-based techniques either use simplifying assumptions or require specialized diagnostics to infer high-dimensional (> $2D$) beam properties. In this Letter, we introduce a general-purpose algorithm that combines neural networks with differentiable particle tracking to efficiently reconstruct high-dimensional phase space distributions without using specialized beam diagnostics or beam manipulations. Furthermore, we demonstrate that our algorithm accurately reconstructs detailed 4D phase space distributions with corresponding confidence intervals in both simulation and experiment using a single focusing quadrupole and diagnostic screen. This technique allows for the measurement of multiple correlated phase spaces simultaneously, which will enable simplified 6D phase space distribution reconstructions in the future.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Reducing Southern Ocean Shortwave Radiation Errors in the ERA5 Reanalysis with Machine Learning and 25 Years of Surface Observations

Earth system models struggle to simulate clouds and their radiative effects over the Southern Ocean, partly due to a lack of measurements and targeted cloud microphysics knowledge. We have evaluated biases of downwelling shortwave radiation in the ERA5 climate reanalysis using 25 years (1995–2019) of summertime surface measurements, collected on the Research and Supply Vessel (RSV) Aurora Australis, the Research Vessel (R/V) Investigator, and at Macquarie Island. During October–March daylight hours, the ERA5 simulation of SW down exhibited large errors (mean bias = 54 W m -2 , mean absolute error = 82 W m -2 , root-mean-square error = 132 W m -2 , and R 2 = 0.71). To determine whether we could improve these statistics, we bypassed ERA5’s radiative transfer model for SW down with machine learning–based models using a number of ERA5’s gridscale meteorological variables as predictors. These models were trained and tested with the surface measurements of SW down using a 10-fold shuffle split. An extreme gradient boosting (XGBoost) and a random forest–based model setup had the best performance relative to ERA5, both with a near complete reduction of the mean bias error, a decrease in the mean absolute error and root-mean-square error by 25% ± 3%, and an increase in the R 2 value of 5% ± 1% over the 10 splits. Large improvements occurred at higher latitudes and cyclone cold sectors, where ERA5 performed most poorly. We further interpret our methods using Shapley additive explanations. Our results indicate that data-driven techniques could have an important role in simulating surface radiation fluxes and in improving reanalysis products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Relativistic Modeling Capabilities in Perseus Extended MHD Simulation Code for HED Plasmas

This project is submitted in response to DE-FOA-0001153, and will be carried out by Dr. Charles Seyler and Dr. Nathaniel Hamlin as post doc under the direction of Professor Seyler. The objectives are to enable PERSEUS to self-consistently model relativistic high-energy-density (HED) laboratory and astrophysical phenomena. Electrons and ions in a plasma are generally modeled either as fluids using a two-fluid code, like PERSEUS, or as particles using a Particle-In- Cell (PIC) code. Nonrelativistic PERSEUS solves the two-fluid equations, formulated in terms of a generalized Ohm’s law (GOL), so as to model about nine orders of magnitude in density variation using a local semi-implicit method. Relativistic PERSEUS preserves this structure and therefore retains these advantageous properties, enabling it to model a broader range of relativistic HED phenomena than previous codes. We have also overcome a major technical challenge associated with solving a relativistic system of equations relating conserved quantities (momentum, energy, etc.) to primitive variables (density, velocity, and pressure). We have thus far developed a relativistic version of PERSEUS and demonstrated that it recovers expected nonrelativistic results. We are currently using relativistic PERSEUS to model two HED phenomena, namely laser-plasma interactions and X-pinches, both of which show spectroscopic evidence of relativistic electrons. In an X-pinch, a section of wire carrying sufficiently large current forms a plasma column that pinches down to a point-like source, generating X-rays for use in radiographic imaging. Laser-plasma interactions have been simulated using PIC codes, and we use their results as a benchmark for comparison. However, neither PIC codes nor existing two-fluid schemes can simulate X-pinches, due to the broad range of densities which have thus far made such a simulation computationally prohibitive. In both laser-plasma interaction and Xpinch simulations, we have observed relativistic phenomena that PERSEUS could not previously model. This includes relativistic channeling of a laser into an over-dense deuterium gas (i.e. induced transparency), consistent with PIC simulation results, along with energetic electrons and ions following an X-pinch and likely accompanied by strong X-ray emission, consistent with experimental observations of energetic X-rays. This proposal includes our latest simulation results. We will improve computational efficiency with techniques that allow the use of a less refined numerical grid. This will enable extensions to three-dimensional geometries. We will improve modeling of laser-plasma interactions and X-pinches, including radiative transport processes, in order to improve comparisons with PIC simulations and experimental diagnostics. We will also model power loss through electrode surface plasmas in magnetically insulated transmission lines (MITLs). We will model relativistic astrophysical phenomena in both astrophysical and HED regimes, which could shed light on the validity of laboratory astrophysics for certain problems. These phenomena are currently modeled using PIC codes and existing relativistic two-fluid codes, neither of which can model the range of length scales of which relativistic PERSEUS is capable. There is strong potential for this project to have a broad impact on both the astrophysical and laboratory plasma communities. In the laboratory community, our work is applicable to those concerned with laser-plasma interactions such as laser-driven fusion, and those in the dense Z pinch community concerned with X-pinches and electrode surface plasmas. An improved understanding of the mechanisms of hard (i.e. high-energy) X-ray generation will likely make important contributions toward the use of X-pinches for radiography through the harnessing or mitigation of hard X-rays, depending on the application. The modeling of electrode surface plasmas, including vacuum electron flow, will improve our understanding of power loss to various loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

The Arrow of Time in the Collapse of Collisionless Self-gravitating Systems: Non-validity of the Vlasov–Poisson Equation during Violent Relaxation

The collapse of a collisionless self-gravitating system, with the fast achievement of a quasi-stationary state, is driven by violent relaxation, with a typical particle interacting with the time-changing collective potential. It is traditionally assumed that this evolution is governed by the Vlasov–Poisson equation, in which case entropy must be conserved. We run N-body simulations of isolated self-gravitating systems, using three simulation codes, NBODY-6 (direct summation without softening), NBODY-2 (direct summation with softening), and GADGET-2 (tree code with softening), for different numbers of particles and initial conditions. At each snapshot, we estimate the Shannon entropy of the distribution function with three different techniques: Kernel, Nearest Neighbor, and EnBiD. For all simulation codes and estimators, the entropy evolution converges to the same limit as N increases. During violent relaxation, the entropy has a fast increase followed by damping oscillations, indicating that violent relaxation must be described by a kinetic equation other than the Vlasov–Poisson equation, even for N as large as that of astronomical structures. This indicates that violent relaxation cannot be described by a time-reversible equation, shedding some light on the so-called “fundamental paradox of stellar dynamics.” The long-term evolution is well-described by the orbit-averaged Fokker–Planck model, with Coulomb logarithm values in the expected range 10--12. By means of NBODY-2, we also study the dependence of the two-body relaxation timescale on the softening length. The approach presented in the current work can potentially provide a general method for testing any kinetic equation intended to describe the macroscopic evolution of N-body systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Atomic-Scale Structure and Stress Release Mechanism in Core–Shell Nanoparticles

According to the article, core–shell nanoparticles find applications in catalysts, sensors, and theranostics. The full internal 3D atomic structure, however, cannot be resolved by current imaging and diffraction techniques. We analyzed the atomic positions and stress-release mechanism in a cubic Au–Pd core–shell nanoparticle in approximately 1000 times higher resolution than current experimental techniques using large-scale molecular dynamics simulation to overcome these limitations. The core–shell nanocube of 73 nm size was modeled similarly to solution synthesis by random epitaxial deposition of a 4 nm thick shell of Pd atoms onto a Au core of 65 nm side length using reliable interatomic potentials. The internal structure reveals specific deformations and stress relaxation mechanisms that are caused by the +4.8% lattice mismatch of gold relative to palladium and differential confinement of extended particle facets, edges, and corners by one, two, or three Au–Pd interfaces, respectively. The three-dimensional lattice strain causes long-range, arc-like bending of atomic rows along the faces and edges of the particle, especially near the Au–Pd interface, a bulging deformation of the Pd shell, and stacking faults in the Pd shell at the corners of the particle. The strain pattern and mechanism of stress release were further characterized by profiles of the atomic layer spacing in the principal crystallographic directions. Accordingly, strain in the Pd shell is several times larger in the extended facets than near the edges and corners of the nanoparticle, which likely affects adsorption, optical, and electrochemical properties. The findings are consistent with available experimental data, including 3D reconstructions of the same cubic nanoparticle by coherent diffractive imaging (CDI) and may be verified by more powerful experimental techniques in the future. The stress release mechanisms are representative for cubic core–shell nanoparticles with fcc structure and can be explored for different shapes by the same methods.

Chemistry↗
DOE OSTI2014

SU-D-18C-01: A Novel 4D-MRI Technology Based On K-Space Retrospective Sorting

Purpose: Current 4D-MRI techniques lack sufficient temporal/spatial resolution and consistent tumor contrast. To overcome these limitations, this study presents the development and initial evaluation of an entirely new framework of 4D-MRI based on k-space retrospective sorting. Methods: An important challenge of the proposed technique is to determine the number of repeated scans(NR) required to obtain sufficient k-space data for 4D-MRI. To do that, simulations using 29 cancer patients' respiratory profiles were performed to derive the relationship between data acquisition completeness(Cp) and NR, also relationship between NR(Cp=95%) and the following factors: total slice(NS), respiratory phase bin length(Lb), frame rate(fr), resolution(R) and image acquisition starting-phase(P0). To evaluate our technique, a computer simulation study on a 4D digital human phantom (XCAT) were conducted with regular breathing (fr=0.5Hz; R=256×256). A 2D echo planer imaging(EPI) MRI sequence were assumed to acquire raw k-space data, with respiratory signal and acquisition time for each k-space data line recorded simultaneously. K-space data was re-sorted based on respiratory phases. To evaluate 4D-MRI image quality, tumor trajectories were measured and compared with the input signal. Mean relative amplitude difference(D) and cross-correlation coefficient(CC) are calculated. Finally, phase-sharing sliding window technique was applied to investigate the feasibility of generating ultra-fast 4D-MRI. Result: Cp increased with NR(Cp=100*[1-exp(-0.19*NR)], when NS=30, Lb=100%/6). NR(Cp=95%) was inversely-proportional to Lb (r=0.97), but independent of other factors. 4D-MRI on XCAT demonstrated highly accurate motion information (D=0.67%, CC=0.996) with much less artifacts than those on image-based sorting 4D-MRI. Ultra-fast 4D-MRI with an apparent temporal resolution of 10 frames/second was reconstructed using the phase-sharing sliding window technique. Conclusions: A novel 4D-MRI technology based on k-space sorting has been successfully developed and evaluated on the digital phantom. Framework established can be applied to a variety of MR sequences, showing great promises to develop the optimal 4D-MRI technique for many radiation therapy applications. NIH (1R21CA165384-01A1)

60 APPLIED LIFE SCIENCES↗