Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical studies”

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 181 records · Page 10

Origin of the reverse optical-contrast change of Ga-Sb phase-change materials—An ab initio molecular-dynamics study

A large number of phase-change materials (PCMs) have been developed experimentally; however, only Ge{sub 2}Sb{sub 2}Te{sub 5}-based PCMs have been significantly explored using ab initio molecular-dynamics (AIMD) simulations. We present an AIMD study of the full melt/quench/anneal PC cycle for Ga-Sb materials, namely, the stoichiometric composition, GaSb, and the near-eutectic alloy, Ga{sub 16}Sb{sub 84}. The calculated electronic densities of states and optical reflectivities are compared between the amorphous and crystalline phases for both compositions, and it is shown that the contrasting opto-electronic properties of each crystalline material can be attributed to different structural transformations of Ga and Sb on crystallization from the amorphous state.

36 MATERIALS SCIENCE↗

Ab initio molecular dynamics study of the properties of cerium in liquid sodium at 1000 K temperature

For liquid-sodium-cooled fast nuclear reactor systems, it is crucial to understand the behavior of lanthanides and other potential fission products in liquid sodium or other liquid metal solutions such as liquid cesium-sodium. In this study, we focus on lanthanide behavior in liquid sodium. Using ab initio molecular dynamics, we found that the solubility of cerium in liquid sodium at 1000 K was less than 0.78 at. %, and the diffusion coefficient of cerium in liquid sodium was calculated to be 5.57 × 10{sup −9} m{sup 2}/s. Furthermore, it was found that cerium in small amounts may significantly alter the heat capacity of the liquid sodium system. Our results are consistent with the experimental results for similar materials under similar conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydration of divalent lanthanides, Sm 2+ and Eu 2+ : A molecular dynamics study with polarizable AMOEBA force field

The chemistry of divalent lanthanides, Ln 2+ , is a growing sub-field of heavy element chemistry owing to new synthetic approaches. However, some theoretical aspects of these unusual cations are currently underdeveloped, especially as they relate to their dynamic properties in solution. In this study, we address the hydration of two of the classical Ln 2+ cations, Sm 2+ and Eu 2+ , using atomic multipole optimized energetic for biomolecular applications (AMOEBA) force fields. These cations have not been parameterized to date with AMOEBA, and few studies are available because of their instability with respect to oxidation in aqueous media. Coordination numbers (CN's) of 8.2 and 8.1 respectively for Sm 2+ and Eu 2+ , and 8.8 for both Sm 3+ and Eu 3+ have been obtained and are in good agreement with the few available AIMD and X-ray absorption fine structures studies. The decreased CN of Ln 2+ compared with Ln 3+ arises from progressive water exchange events that indicates the gradual stabilization of 8-coordinate structures with respect to 9-coordinate geometries. Moreover, the effects of the chloride counter anions on the coordination of Ln 2+ cations have been studied at different chloride concentrations in this work. Lastly, water exchange times of Ln 2+ cations have been calculated to provide a comprehensive understanding of the behavior of Eu 2+ and Sm 2+ in aqueous chloride media.

AMOEBA force field parameterization↗

Phase transition in bilayer MoS 2 under tensile loading: a molecular dynamics study

Abstract Molybdenum disulfide (MoS 2 ), especially single-layer MoS 2 , has been experimentally and computationally discovered to exist in several different polymorphs exhibiting various electronic and mechanical properties. The morphology of MoS 2 can be tuned through strain engineering. Molecular dynamics simulations are conducted to systematically study the phase transition of single-layer MoS 2 and bilayer MoS 2 under the uniaxial tensile condition at room temperature. The roles of edge and S-line vacancy are investigated. Phase transitions are always triggered near the edge and vacancy sites. The initiation of the metastable T″ phase can release the tensile stress in the lattice, followed by I4/mmm phase initiation, regardless of the edge conditions. The growth of the I4/mmm phase can cause the local buckling of the MoS 2 plane. With a tilted S-line vacancy, I4/mmm phase is first initiated to reduce the local shear stress accumulated near the vacancy line. Overall, the phase transition mechanism of single layer and bilayer MoS 2 under the uniaxial tensile loading is provided, which guides the future strain engineering of MoS 2 in nanoelectronics applications.

Materials Science↗

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A non-adiabatic dynamics study of octatetraene: The radiationless conversion from S{sub 2} to S{sub 1}

Simulation of the excited state dynamics of all-trans-1,3,5,7-octatetraene has been performed to investigate the ultrafast radiationless S{sub 2} → S{sub 1} internal conversion process. Multireference configuration interaction with single excitation method has been employed to optimize the equilibrium structure of the excited states, as well as the S{sub 2}/S{sub 1} conical intersection, and to investigate the non-adiabatic molecular dynamics of the S{sub 2}/S{sub 1} state transition. At the conical intersection, the molecule is found to be distorted from the original planar trans structure to a nearly perpendicular conformation around C{sub 3}−C{sub 4} bond, with the torsion angle being about 107°. Such structural change can result in mutual approaching of states S{sub 2} and S{sub 1} in energy, and drastically increase the nonadiabatic coupling between the two states by destroying the inter-state symmetry prohibition in the electronic wavefunctions. Surface-hopping molecular dynamics simulations are performed to describe the non-adiabatic process. Upon the Franck-Condon excitation to the S{sub 2} state, the molecule quickly twists its C{sub 3}−C{sub 4} bond and approaches the conical intersection region, where it can undergo efficient internal conversion to S{sub 1}. The decay time constant (τ) of S{sub 2} state is estimated to be around 251 fs by fitting the occupation number of average fraction of trajectories using an exponential damping function. This value is reasonably consistent with previous experimental measurements of around 300–400 fs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Studies of Photo- and Electron-Induced Reactions on Nanostructured Surfaces

Our current research program examines the photon- and electron-initiated reaction mechanisms, half-collision dynamics, and other nonequilibrium-excited dynamics effects, occurring with excitation of adsorbates on well-characterized metal-oxide of metalchalcogenide surfaces, both in single-crystal and in nanocrystalline forms. The program has made extensive use of STM methods, collaborations with the Nanoprobing SPLEEM system at Brookhaven National Lab, and XPS methods. As our current program terminates, it brings to an end a sustained and major effort to understand reactions and electronic structure on materials of interest for photochemistry on precision surfaces. The last two close-out years of our program allowed us to complete our STM surface studies, interact with scientists at BNL on programs useful for two dimensional (2D) materials, and initiate a new effort in ultrafast two-photon photoemission spectroscopy. Our work has then shown how nanoscale stressed surface regions on sold crystal surfaces (typically TiO 2 ) lead to enhanced reactivity of adsorbate molecules – a result of much importance for nanocatalysis. In addition we have developed several new surface and near-surface probing methods – these, also, are of general use to the community. Finally due to the pervasive interest in DOE labs in these materials, we have exploited our earlier synchrotron-based understanding of the electronic structure in (2D) materials to begin a study of surface chemical reactions in these same materials and their corresponding molecular dynamics. These efforts are described briefly below and in the references or DOE publications presented in the publications section of this report. We are pleased to have been able to contribute so many publications in the last two years of this many year DOE program. We also wish to thank DOE for the sustained funding that has made this program such a success.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A molecular dynamics study of intramolecular proton transfer reaction of malonaldehyde in solutions based upon mixed quantum-classical approximation. I. Proton transfer reaction in water

In order to investigate proton transfer reaction in solution, mixed quantum-classical molecular dynamics calculations have been carried out based on our previously proposed quantum equation of motion for the reacting system [A. Yamada and S. Okazaki, J. Chem. Phys. 128, 044507 (2008)]. Surface hopping method was applied to describe forces acting on the solvent classical degrees of freedom. In a series of our studies, quantum and solvent effects on the reaction dynamics in solutions have been analysed in detail. Here, we report our mixed quantum-classical molecular dynamics calculations for intramolecular proton transfer of malonaldehyde in water. Thermally activated proton transfer process, i.e., vibrational excitation in the reactant state followed by transition to the product state and vibrational relaxation in the product state, as well as tunneling reaction can be described by solving the equation of motion. Zero point energy is, of course, included, too. The quantum simulation in water has been compared with the fully classical one and the wave packet calculation in vacuum. The calculated quantum reaction rate in water was 0.70 ps{sup −1}, which is about 2.5 times faster than that in vacuum, 0.27 ps{sup −1}. This indicates that the solvent water accelerates the reaction. Further, the quantum calculation resulted in the reaction rate about 2 times faster than the fully classical calculation, which indicates that quantum effect enhances the reaction rate, too. Contribution from three reaction mechanisms, i.e., tunneling, thermal activation, and barrier vanishing reactions, is 33:46:21 in the mixed quantum-classical calculations. This clearly shows that the tunneling effect is important in the reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A seven-degree-of-freedom, time-dependent quantum dynamics study on the energy efficiency in surmounting the central energy barrier of the OH + CH{sub 3} → O + CH{sub 4} reaction

A time-dependent, quantum reaction dynamics calculation with seven degrees of freedom was carried out to study the energy efficiency in surmounting the approximate center energy barrier of OH + CH{sub 3}. The calculation shows the OH vibration excitations greatly enhance the reactivity, whereas the vibrational excitations of CH{sub 3} and the rotational excitations hinder the reactivity. On the basis of equal amount of total energy, although this reaction has a slight early barrier, it is the OH vibrational energy that is the dominate force in promoting the reactivity, not the translational energy. The studies on both the forward O + CH{sub 4} and reverse OH + CH{sub 3} reactions demonstrate, for these central barrier reactions, a small change of the barrier location can significantly change the energy efficacy roles on the reactivity. The calculated rate constants agree with the experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A discrete dislocation dynamics study of precipitate bypass mechanisms in nickel-based superalloys

Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.

36 MATERIALS SCIENCE↗

Molecular dynamics study of Si(100)-oxidation: SiO and Si emissions from Si/SiO{sub 2} interfaces and their incorporation into SiO{sub 2}

Dynamics of Si(100)-oxidation processes at the Si/SiO{sub 2} interface and in the SiO{sub 2} region are investigated focusing on SiO and Si emissions from the interface and the following incorporation into the SiO{sub 2} and/or substrate. Classical molecular dynamics (MD) simulations with variable charge interatomic potentials are performed to clarify these atomic processes. By incorporating oxygen atoms, two-folded Si atoms are formed after structural relaxation at the interface and are emitted as SiO molecules into SiO{sub 2}. The energy barrier of the SiO emission is estimated to be 1.20 eV on the basis of the enthalpy change in an MD simulation. The emitted SiO molecule is incorporated into the SiO{sub 2} network through a Si-O rebonding process with generating an oxygen vacancy. The energy barrier of the SiO incorporation is estimated to be 0.79–0.81 eV. The elementary process of oxygen vacancy diffusion leading to the complete SiO incorporation is also simulated, and the energy barriers are found to be relatively small, 0.71–0.79 eV. The energy changes of Si emissions into the substrate and SiO{sub 2} are estimated to be 2.97–7.81 eV, which are larger than the energy barrier of the SiO emission. This result suggests that, at the ideally flat Si/SiO{sub 2} interface, the SiO emission into the SiO{sub 2} region occurs prior to the Si emission, which is consistent with previous theoretical and experimental studies. The above mentioned typical atomic processes are successfully extracted from some (or one) of MD simulations among many trials in which a statistical procedure is partly employed. Our results give a unified understanding of Si oxidation processes from an atomistic point of view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unraveling Photodimerization of Cyclohexasilane from Molecular Dynamics Studies

Photoinduced reactions of a pair of cyclohexasilane (CHS) monomers are explored by time-dependent excited-state molecular dynamics (TDESMD) calculations. In TDESMD trajectories, one observes vivid reaction events including dimerization and fragmentation. A general reaction pathway is identified as (i) ring-opening formation of a dimer, (ii) rearrangement induced by bond breaking, and (iii) decomposition through the elimination of small fragments. The identified pathway supports the chemistry proposed for the fabrication of silicon-based materials using CHS as a precursor. Additionally, we find dimers have smaller HOMO–LUMO gaps and exhibit a red shift and line-width broadening in the computed photoluminescence spectra compared with a pair of CHS monomers.

36 MATERIALS SCIENCE↗

Molecular dynamics studies of grain boundary mobility and anisotropy in BCC γ-uranium

Grain morphologies such as grain size and aspect ratio in uranium-based metallic fuels are important microstructural features that can impact various fuel performance properties such as fission-gas-induced swelling, thermal transport, high burnup structure formation, and radiation resistance. Accurate prediction of the fuel grain morphologies requires knowledge of critical grain growth parameters such as grain boundary (GB) mobility and anisotropy. In this work, molecular dynamics simulations were performed to study the GB mobility and its anisotropy in pure body-centered-cubic (BCC) γ uranium. Nine GBs with different combinations of misorientation angles (20°, 30°, 45°) and rotation axes (<100>, <110>, <111>), as well as an additional <111> 38.2° GB were studied using three interatomic potentials. Here, it is found that the GB mobility anisotropy has complex trends, depending on both rotation axis and misorientation. However, in general the <110> rotation axis has the fastest GB mobility at the same misorientation. The results of this work can be used as not only baselines for future studies of GB mobility in uranium-based alloys such as uranium-molybdenum (U-Mo) fuels, but also input for mesoscale modeling of grain growth in uranium-based alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A molecular dynamics study of domain switching in BiFeO 3 nanofilm under DC electric field

The performance of ferroelectric materials is deeply rooted in the domain switching process triggered by an applied electric field. To understand the effect of controlling factors on the domain switching characteristics, nanoscale studies are essential. In this work, an atomistic-scale computational framework is developed to study the domain switching process of bismuth ferrite (BiFeO 3 ) nanofilms. Through a series of molecular dynamics simulations, the effects of duration and magnitude of an electric field, domain size and nanofilm thickness on the domain switching behavior are investigated. The evolution of domains under different magnitudes and durations of the applied electric field is studied, and the results of unswitched and switched cases are summarized in a switching map. Analysis of various nanofilms reveals that by increasing the film thickness, the switching time is increased. Additionally, the nanofilms with different domain sizes and configurations are examined, and the minimum domain size for complete reversal switching is obtained. Our results are helpful to predict and tune the switching process in ferroelectric materials.

36 MATERIALS SCIENCE↗

Thermal transpiration: A molecular dynamics study

Thermal transpiration is a phenomenon where fluid molecules move from the cold end towards the hot end of a channel under the influence of longitudinal temperature gradient alone. Although the phenomenon of thermal transpiration is observed at rarefied gas conditions in macro systems, the phenomenon can occur at atmospheric pressure if the characteristic dimensions of the channel is less than 100 nm. The flow through these nanosized channels is characterized by the free molecular flow regimes and continuum theory is inadequate to describe the flow. Thus a non-continuum method like molecular dynamics (MD) is necessary to study such phenomenon. In the present work, MD simulations were carried out to investigate the occurance of thermal transpiration in copper and platinum nanochannels at atmospheric pressure conditions. The mean pressure of argon gas confined inside the nano channels was maintained around 1 bar. The channel height is maintained at 2nm. The argon atoms interact with each other and with the wall atoms through the Lennard-Jones potential. The wall atoms are modelled using an EAM potential. Further, separate simulations were carried out where a Harmonic potential is used for the atom-atom interaction in the platinum channel. A thermally insulating wall was introduced between the low and high temperature regions and those wall atoms interact with fluid atoms through a repulsive potential. A reduced cut off radius were used to achieve this. Thermal creep is induced by applying a temperature gradient along the channel wall. It was found that flow developed in the direction of the increasing temperature gradient of the wall. An increase in the volumetric flux was observed as the length of the cold and the hot regions of the wall were increased. The effect of temperature gradient and the wall-fluid interaction strength on the flow parameters have been studied to understand the phenomenon better.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Influence of molecular shape on self-diffusion under severe confinement: A molecular dynamics study

In this study, we have investigated the effect of molecular shape and charge asymmetry on the translation dynamics of confined hydrocarbon molecules having different shapes but similar kinetic diameters, inside ZSM-5 pores using molecular dynamics simulations. The mean square displacement of propane, acetonitrile, acetaldehyde, and acetone in ZSM-5 exhibit two different regimes – ballistic and diffusive/sub-diffusive. All the molecules except propane exhibit sub-diffusive motion at time scales greater than 1 ps. The intermediate scattering functions reveal that there is a considerable rotational-translational coupling in the motion of all the molecules, due to the strong geometrical restriction imposed by ZSM-5. Overall the difference in shape and asymmetry in charge imposes severe restriction inside the ZSM-5 channels for all the molecules to different extents. Additionally, the behavior of molecules confined in ZSM-5 in the present study, quantified wherever possible, is compared to their behavior in bulk or in other porous media reported in literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electro-impulse deicing - Structural dynamic studies, icing tunnel tests and applications

This paper contains a discussion of the comparison of analytical and experimental results of the dynamic response of a flat rectangular plate subjected to electro-impulse type deicing forces. Early attempts in this correlation have been hampered by the complex leading edge geometries of the airfoils tested to date. The lack of a suitable analytical model for a typical leading edge structure has prompted these preliminary investigations of simple geometries, beginning with a flat plate and to be followed by a semi-cylindrical leading edge representation. The comparisons reported herein are thus limited to rectangular plate models. Project plans and icing tunnel results are given.

Bernhart, W. D.↗

A tight-binding molecular dynamics study of the noble metals Cu, Ag and Au

Here we have used the Naval Research Lab (NRL) tight-binding (TB) method to study the electronic and mechanical properties of the noble metals. In order to perform molecular dynamics simulations, we used new TB parameters that work well at smaller interatomic distances. The TB parameters were fitted to the fcc, bcc and sc periodic structures and were demonstrated to be transferable and robust for calculating additional dynamical properties which they had not been fitted to. We calculated the phonon frequencies and density of states at finite temperature and we also performed simulations to determine the temperature dependence of the coefficient of thermal expansion and the mean squared displacement. The energy for vacancy formation as well as energy for fcc-based, bcc-based clusters and icosahedral clusters of different sizes were also calculated. The results compared very well with experimental observations and independent (non-fitted) first-principles density functional calculations.

36 MATERIALS SCIENCE↗