Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular statics”

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
NASA NTRS1997

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerenes and their Correlation with Three Geometric Parameters: Group Order, Aromaticity, and Size

Calculations were performed on the valence contribution to the static molecular third-order polarizabilities (gamma) of thirty carbon-cage fullerenes (C60, C70, five isomers of C78, and twenty-three isomers of C84). The molecular structures were obtained from B3LYP/STO-3G calculations. The values of the tensor elements and an associated numerical uncertainty were obtained using the finite-field approach and polynomial expansions of orders four to eighteen of polarization versus static electric field data. The latter information was obtained from semiempirical calculations using the AM1 hamiltonian.

Moore, Craig E.↗
DOE OSTI2020

Molecular Statics Analyses of Thermodynamics and Kinetics of Hydrogen Cottrell Atmosphere Formation Around Edge Dislocations in Aluminum

Aluminum alloys are being explored as lightweight structural materials for use in hydrogen-containing environments. In this paper, to understand hydrogen effects on deformation, we perform molecular statics studies of the hydrogen Cottrell atmosphere around edge dislocations in aluminum. First, we calculate the hydrogen binding energies at all interstitial sites in a periodic aluminum crystal containing an edge dislocation dipole. This allows us to use the Boltzmann equation to quantify the hydrogen Cottrell atmosphere. Based on these binding energies, we then construct a continuum model to study the kinetics of the hydrogen Cottrell atmosphere formation. Finally, we compare our results with existing theories and discuss the effects of hydrogen on deformation of aluminum-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Size effects of 109° domain walls in rhombohedral barium titanate single crystals—A molecular statics analysis

Ferroelectric functional materials are of great interest in science and technology due to their electromechanically coupled material properties. Therefore, ferroelectrics, such as barium titanate, are modeled and simulated at the continuum scale as well as at the atomistic scale. Due to recent advancements in related manufacturing technologies the modeling and simulation of smart materials at the nanometer length scale is getting more important not only to predict but also fundamentally understand the complex material behavior of such materials. In this study, we analyze the size effects of 109° nanodomain walls in ferroelectric barium titanate single crystals in the rhombohedral phase using a recently proposed extended molecular statics algorithm. We study the impact of domain thicknesses on the spontaneous polarization, the coercive field, and the lattice constants. Moreover, we discuss how the electromechanical coupling of an applied electric field and the introduced strain in the converse piezoelectric effect is affected by the thickness of nanodomains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1998

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerene and Their Correlation with Three Geometric Properties: Symmetry, Aromaticity, and Size

The static third-order polarizabilities (gamma) of C60, C70, five isomers of C78 and two isomers of C84 were analyzed in terms of three properties, from a geometric point of view: symmetry, aromaticity and size. The polarizability values were based on the finite field approximation using a semiempirical Hamiltonian (AM1) and applied to molecular structures obtained from density functional theory calculations. Symmetry was characterized by the molecular group order. The selection of 6-member rings as aromatic was determined from an analysis of bond lengths. Maximum interatomic distance and surface area were the parameters considered with respect to size. Based on triple linear regression analysis, it was found that the static linear polarizability (alpha) and gamma in these molecules respond differently to geometrical properties: alpha depends almost exclusively on surface area while gamma is affected by a combination of number of aromatic rings, length and group order, in decreasing importance. In the case of alpha, valence electron contributions provide the same information as all-electron estimates. For gamma, the best correlation coefficients are obtained when all-electron estimates are used and when the dependent parameter is ln(gamma) instead of gamma.

Moore, C. E.↗
DOE OSTI2018

An atomistic investigation of the interaction of dislocations with Guinier-Preston zones in Al-Cu alloys

The interaction between edge dislocations and Guinier-Preston zones in an Al-Cu alloy was analyzed by means of atomistic simulations. The different thermodynamic functions that determine the features of these obstacles for the dislocation glide were computed using molecular statics, molecular dynamics and the nudged elastic band method. It was found that Guinier-Preston zones are sheared by dislocations and the rate at which dislocations overcome the precipitate is controlled by the activation energy, Δ U , in agreement with the postulates of the harmonic transition state theory. Moreover, the entropic contribution to the Helmholtz activation free energy was in the range 1.3–1.8 k b , which can be associated with the typical vibrational entropy of solids. Finally, an estimation of the initial shear flow stress as a function of temperature was carried out from the thermodynamic data provided by the atomistic simulations. Furthermore, comparison with experimental results showed that the effect of the random precipitate distribution and of the dislocation character and dislocation/precipitation orientation has to be taken into account in the simulations to better reproduce experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Vacancy enhanced formation and phase transition of Cu-rich precipitates in α - iron under neutron irradiation

In this paper, we employed both molecular statics and molecular dynamics simulation methods to investigate the role of vacancies in the formation and phase transition of Cu-rich precipitates in α-iron. The results indicated that vacancies promoted the diffusion of Cu atoms to form Cu-rich precipitates. After Cu-rich precipitates formed, they further trapped vacancies. The supersaturated vacancy concentration in the Cu-rich precipitate induced a shear strain, which triggered the phase transition from bcc to fcc structure by transforming the initial bcc (110) plane into fcc (111) plane. In addition, the formation of the fcc-twin structure and the stacking fault structure in the Cu-rich precipitates was observed in dynamics simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Atomistic Simulation of Interfaces in Materials of Solid State Ionics

The possibilities of describing correctly interfaces of different types in solids within a computer experiment using molecular statics simulation, molecular dynamics simulation, and quantum chemical calculations are discussed. Heterophase boundaries of various types, including grain boundaries and solid electrolyte‒solid electrolyte and ionic conductor‒electrode material interfaces, are considered. Specific microstructural features and mechanisms of the ion transport in real heterophase structures (cationic conductor‒metal anode and anionic conductor‒cathode) existing in solid state ionics devices (such as solid-state batteries and fuel cells) are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Homogeneous hydride formation path in α-Zr: Molecular dynamics simulations with the charge-optimized many-body potential

A formation path for homogeneous γ hydride formation in hcp α-Zr, from solid solution to the ζ and then the γ hydride, was demonstrated using molecular static calculations and molecular dynamic simulations with the charge-optimized many-body (COMB) potential. Hydrogen has limited solubility in α-Zr. Once the solubility limit is exceeded, the stability of solid solution gives way to that of coherent hydride phases such as the ζ hydride by planar precipitation of hydrogen. At finite temperatures, the ζ hydride goes through a partial hcp-fcc transformation via 1/3 <1¯100> slip on the basal plane, and transforms into a mixture of γ hydride and α-Zr. In the ζ hydride, slip on the basal plane is favored thermodynamically with negligible barrier, and is therefore feasible at finite temperatures without mechanical loading. The transformation process involves slips of three equivalent shear partials, in contrast to that proposed in the literature where only a single shear partial was involved. The adoption of multiple slip partials minimizes the macroscopic shape change of embedded hydride clusters and the shear strain accumulation in the matrix, and thus reduces the overall barrier needed for homogeneous γ hydride formation. In conclusion, this formation path requires finite temperatures for hydrogen diffusion without mechanical loading. Therefore, it should be effective at the cladding operating conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

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

Interactions between displacement cascade and dislocation and their influences on Peierls stress in Fe-20Cr-25Ni alloys

The interactions between dislocations and displacement cascades are investigated using molecular dynamics simulations. The primary knock-on atoms (PKAs) are placed at different distances to edge and screw dislocation lines. Specifically, an fcc Fe-20Cr-25Ni system is studied; the alloy is the based alloy for Alloy 709. The number of defects is calculated using Wigner-Seitz analysis in order to illustrate the effect of dislocations on cascades. The displacement cascades in systems containing dislocations tend to generate more total defects compared with bulk systems, but eventually lead to less surviving defects in the materials matrix due to interactions with dislocations. The changes in atomic structures of the dislocations after interacting with cascades are analyzed to understand how the primary damage affects dislocations. The displacement cascade potentially causes dislocation climb in edge dislocations and cross-slip in screw dislocations, which could serve as additional mechanisms contributing to the radiation-induced hardening (RIH) compared with conventional RIH models. To reveal this, Peierls stresses are calculated before and after cascades using molecular static simulations. This study provides critical information of how defect production and dislocations are correlated, especially when total dose is high, which needs to be taken into account in upper scale models, such as mean-field rate theory.

Lu, Zizhe↗
NASA NTRS1993

Diffusion on Cu surfaces

Understanding surface diffusion is essential in understanding surface phenomena, such as crystal growth, thin film growth, corrosion, physisorption, and chemisorption. Because of its importance, various experimental and theoretical efforts have been directed to understand this phenomena. The Field Ion Microscope (FIM) has been the major experimental tool for studying surface diffusion. FIM have been employed by various research groups to study surface diffusion of adatoms. Because of limitations of the FIM, such studies are only limited to a few surfaces: nickel, platinum, aluminum, iridium, tungsten, and rhodium. From the theoretical standpoint, various atomistic simulations are performed to study surface diffusion. In most of these calculations the Embedded Atom Method (EAM) along with the molecular static (MS) simulation are utilized. The EAM is a semi-empirical approach for modeling the interatomic interactions. The MS simulation is a technique for minimizing the total energy of a system of particles with respect to the positions of its particles. One of the objectives of this work is to develop the EAM functions for Cu and use them in conjunction with the molecular static (MS) simulation to study diffusion of a Cu atom on a perfect as well as stepped Cu(100) surfaces. This will provide a test of the validity of the EAM functions on Cu(100) surface and near the stepped environments. In particular, we construct a terrace-ledge-kink (TLK) model and calculate the migration energies of an atom on a terrace, near a ledge site, near a kink site, and going over a descending step. We have also calculated formation energies of an atom on the bare surface, a vacancy in the surface, a stepped surface, and a stepped-kink surface. Our results are compared with the available experimental and theoretical results.

Karimi, Majid↗
DOE OSTI2017

Atomistic Simulation of Interstitial Dislocation Loop Evolution under Applied Stresses in BCC Iron

Evolution of an interstitial 1/2&lang;111&rang; dislocation loop under tensile, shear, and torsion stresses is studied with molecular statics method. Under a tensile stress, the dependence of ultimate tensile strength on size of loop is calculated. The formation of small shear loops around the initial prismatic loop is confirmed as an intermediate state to form the final dislocation network. Under a shear stress, the rotation of a loop is observed not only by a change of the habit plane but also through a transformation between a shear and a prismatic loop. Under torsion, a perfect BCC crystal may undergo a BCC to FCC or BCC to HCP transformation. The present work indicates that a 1/2&lang;111&rang; loop can delay these transformations, resulting in the formation of micro-crack on the surface.

dislocation loop↗
DOE OSTI2019

Phase-field-based calculations of the disregistry fields of static extended dislocations in FCC metals

In the continuum context, the displacements of atoms induced by a dislocation can be approximated by a continuum disregistry field. In this paper, two phase-field (PF)-based approaches and their variants are employed to calculate the disregistry fields of static, extended dislocations of pure edge and pure screw character in two face-centred cubic metals: Au and Al, which have distinct stable stacking fault energy and elastic anisotropy. A new truncated Fourier series form is developed to approximate the generalised stacking fault energy (GSFE) surface, which shows significant improvement over the previously employed Fourier series form. By measuring the intrinsic stacking fault (ISF) width and partial dislocation core size in different ways, the PF-based disregistry fields are quantitatively compared against those predicted by molecular statics. In particular, two new measures for the ISF widths are proposed and shown to overcome drawbacks of the more commonly used standards in the literature. Our calculations also show that continuum formulation of the elastic energy and the GSFE for a homogeneous surface can successfully characterise the core structure. Finally, our comparisons highlight the significance of including the gradient energy in the free energy formulation when an accurate description of the dislocation core structure is desired.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spectrum of embrittling potencies and relation to properties of symmetric-tilt grain boundaries

Grain boundaries in metallic materials can exist in a wide range of stable and metastable structures. In addition, the properties of a grain boundary may be altered through solute segregation. In this work, we present a formulation that combines the spectrum of embrittling potencies associated with solute segregation with site-occupancy statistics. As a prototype problem, we illustrate the relation between segregation and embrittlement in the case of S segregation to grain boundaries in Ni. To obtain a population of site segregation energies, we perform molecular statics calculations on 378 different symmetric-tilt grain boundaries and their free surface equivalents, using an embedded-atom method interatomic potential developed specifically for studying embrittlement. Our results show that it is important to consider both the energies associated with embrittlement and the probability of occupancy to describe the general embrittling nature of a grain boundary. When analyzed in isolation, certain grain boundaries show large embrittling potencies; however, that effect is diminished when the probability of S segregation to that grain boundary is considered within a polycrystal. We propose a new quantity, the embrittling estimator, which not only categorizes grain boundaries as embrittling or strengthening, but also considers site occupancy probabilities, so that the embrittlement behavior of grain boundaries within a network of grain boundaries can be compared. Finally, we examine the relationship between embrittlement behavior and innate grain boundary properties, such as the free volume, and find statistical evidence that the complex nature of embrittlement cannot be explained by linear correlations with excess volumes or energies. Ultimately, this combined approach provides a theoretical tool to assist grain boundary engineering of metastable alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Density functional theory based generalized effective fragment potential method

We present a generalized Kohn-Sham (KS) density functional theory (DFT) based effective fragment potential (EFP2-DFT) method for the treatment of solvent effects. Similar to the original Hartree-Fock (HF) based potential with fitted parameters for water (EFP1) and the generalized HF based potential (EFP2-HF), EFP2-DFT includes electrostatic, exchange-repulsion, polarization, and dispersion potentials, which are generated for a chosen DFT functional for a given isolated molecule. The method does not have fitted parameters, except for implicit parameters within a chosen functional and the dispersion correction to the potential. The electrostatic potential is modeled with a multipolar expansion at each atomic center and bond midpoint using Stone's distributed multipolar analysis. The exchange-repulsion potential between two fragments is composed of the overlap and kinetic energy integrals and the nondiagonal KS matrices in the localized molecular orbital basis. The polarization potential is derived from the static molecular polarizability. The dispersion potential includes the intermolecular D3 dispersion correction of Grimme et al. [J. Chem. Phys. 132, 154104 (2010)]. The potential generated from the CAMB3LYP functional has mean unsigned errors (MUEs) with respect to results from coupled cluster singles, doubles, and perturbative triples with a complete basis set limit (CCSD(T)/CBS) extrapolation, of 1.7, 2.2, 2.0, and 0.5 kcal/mol, for the S22, water-benzene clusters, water clusters, and n-alkane dimers benchmark sets, respectively. The corresponding EFP2-HF errors for the respective benchmarks are 2.41, 3.1, 1.8, and 2.5 kcal/mol. Thus, the new EFP2-DFT-D3 method with the CAMB3LYP functional provides comparable or improved results at lower computational cost and, therefore, extends the range of applicability of EFP2 to larger system sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Computer simulation of surface and film processes

All the investigations which were performed employed in one way or another a computer simulation technique based on atomistic level considerations. In general, three types of simulation methods were used for modeling systems with discrete particles that interact via well defined potential functions: molecular dynamics (a general method for solving the classical equations of motion of a model system); Monte Carlo (the use of Markov chain ensemble averaging technique to model equilibrium properties of a system); and molecular statics (provides properties of a system at T = 0 K). The effects of three-body forces on the vibrational frequencies of triatomic cluster were investigated. The multilayer relaxation phenomena for low index planes of an fcc crystal was analyzed also as a function of the three-body interactions. Various surface properties for Si and SiC system were calculated. Results obtained from static simulation calculations for slip formation were presented. The more elaborate molecular dynamics calculations on the propagation of cracks in two-dimensional systems were outlined.

Tiller, W. A.↗