Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC DISPLACEMENTS”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Uncovering the active sites and demonstrating stable catalyst for the cost-effective conversion of ethanol to 1-butanol

We report the recent emergence of a robust renewable ethanol industry has provided a sustainable platform molecule toward the production of value-added chemicals and fuels; what is lacking now are viable conversion processes from ethanol that can displace the current production pathways from non-renewable pathways. Here in the work, we demonstrate the highly selective conversion of ethanol to higher alcohols over low copper loaded MgAl mixed oxide catalysts, with 50% improvement in higher alcohol yields over the current state of the art. At these copper concentrations, atomically dispersed Cu +1 were found to be stable even at highly reductive conditions and highly active towards higher alcohol products (e.g. butanol, hexanol) while suppressing side reaction pathways and leading to extended lifetimes of over 150 hours time on stream. Technoeconomic analysis conducted based on these experimental results demonstrate that this catalytic system is cost-competitive with the conventional process. This marks significant progress in the development of Guerbet coupling of ethanol as a viable renewable process and offers a pathway toward sustainable chemical and fuel production.

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

DNA Strand Displacement Driven Molecular Additive Manufacturing (DSD-MAM)

The goal of this project was to validate two-dimensional molecular printers, initially selfassembled from DNA and then actuated by externally driven cycles of DNA strand displacement, as prototype integrated nanosystems for molecular additive manufacturing. Novel functionalities of these nanomachines were explored during this project, including the following: nanometer-precision positioning mechanisms based on DNA strand displacement with multivalent interactions for discrete stepping or else diffusive capture; integration of independently moving layers of DNA origami to achieve 2D controllable motion; integration of spatial positioning with deposition functionality. The principal importance of this project was to provide an essential step in the development of a new technology for atomically precise manufacturing. Our first generation molecular 2D printer offers several advantages over conventional DNA-origami patterning, such as faster prototyping, faster dynamic rearrangement of patterns, and the ability to respond with feedback. We anticipate that our first-generation molecular printers may inspire future generations of molecular printers with iterative improvements in robustness and throughput. Potential applications of atomically precise manufacturing include the following: photovoltaics; photosynthetic and fuel cells; thermoelectrics and anisotropic heat spreaders; solid-state lighting; molecular electronic and plasmonic circuits; selectively preamble membranes; self-repairing materials with high strength-to-weight and fracture resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗
NASA NTRS2019

Distorted Perovskites for High Voltage Dielectric Capacitors

There are several mechanisms which have been proposed for the existence of colossal dielectric constant in the class of perovskite calcium copper titanate (CaCu3Ti4O12 or CCTO) materials. Researches indicate that existence of twinning parallel to (100) (001) and (010) planes causes planar defects and causes changes in local electronic structure. This change can cause insulating barriers locally which contribute to the large dielectric values irrespective of processing. The combination of insulating barriers, defects and displacements caused by twinning have been attributed to the generation of large dielectric constant in CCTO. To examine some of these arguments some researchers replaced Ca with other elements and evaluated this concept. In this study we present the synthesis and characterization of Ga2/3Cu3Ti4O12-xNx (GCTON) material. This provides both distortion due to atomic size difference and defects due to insertion of nitrogen. The morphology of the compound was determined to show that processing has tremendous effect on the dielectric values. The resistivity of GCTON was several order higher than CCTO and dielectric constant was higher than 10,000.

Gill, Puneet↗
DOE OSTI2024

Discovery of Dielectric Response and Forces in Sub-Nanoscale Objects (Final Technical Report)

This is the final technical report for DOE Award DE-SC0005132 entitled “Discovery of Dielectric Response and Forces in Sub-Nanoscale Objects.” This award originated on August 15, 2010 with Principle Investigator Prof. Phil Batson of Rutgers, The State University of New Jersey. Prof. Batson transitioned to Professor Emeritus status on January 31, 2023. At that time PI status was transferred to Prof. Robert Bartynski of the Department of Physics and Astronomy at Rutgers, The State University of New Jersey, and Director of the Rutgers Laboratory for Surface Modification. The central theme of the research performed under this award is develop and refine a NION aberration-corrected Scanning Transmission Electron Microscope to attain atomic spatial resolution and sub 10meV energy resolution in electron energy loss spectroscopy (EELS) performed on the transmitted electron beam. These capabilities allow examination of the excitation properties (primarily plasmonic and vibrational [ie., phonons] of nanoscale objects when excited with a highly-localized (~ Angstrom-scale) high energy (~ 60 keV) electron beam. Direct excitations (i.e., the response to electrons impinging directly on the sample) as well as excitations in response to the dynamic electric (and magnetic!) fields of an aloof electron beam (i.e., a beam that is close to but displaced from the target) have been studied. Our experimental work has benefited greatly from close collaborations with theoretical colleagues, but at Rutgers and from around the world, providing a much more complete understanding of the observed phenomena and suggesting avenues for further study and practical applications. This report summarizes the technical and scientific achievements accomplished during the entire award period. More extensive details are available in the Progress Reports that have already been filed with the Department of Energy. The report is divided into six sections. The first four sections focus on exploring, developing and understanding several unique capabilities and phenomena discovered and accessible owing to the high spatial- and energy-resolution we have been able to achieve. The latter two sections focus on applications of the STEM’s novel capabilities to study atomic- and nanometer-scale properties of solids, as well as broader applications to advanced and bio-materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Thickness and temperature dependence of the atomic-scale structure of SrRuO 3 thin films

The temperature-dependent layer-resolved structure of 3 to 44 unit cell thick SrRuO 3 (SRO) films grown on Nb-doped SrTiO 3 substrates is investigated using a combination of high-resolution synchrotron x-ray diffraction and high-resolution electron microscopy to understand the role that structural distortions play in suppressing ferromagnetism in ultra-thin SRO films. The oxygen octahedral tilts and rotations and Sr displacements characteristic of the bulk orthorhombic phase are found to be strongly dependent on temperature, the film thickness, and the distance away from the film–substrate interface. For thicknesses, t, above the critical thickness for ferromagnetism (t > 3 uc), the orthorhombic distortions decrease with increasing temperature above T C . Below T C , the structure of the films remains constant due to the magneto-structural coupling observed in bulk SRO. The orthorhombic distortions are found to be suppressed in the 2–3 interfacial layers due to structural coupling with the SrTiO 3 substrate and correlate with the critical thickness for ferromagnetism in uncapped SRO films.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Super-Resolved Single-Molecule Tracking Studies of Rhodamine B Accumulation on Fresh and Aged Polyethylene Terephthalate

It is well-known that toxic organic micropollutants (OMs) accumulate on the surfaces of microplastics. However, much remains to be learned about the exact molecular level mechanisms of OM accumulation and how these evolve as the plastics age. In this work, super-resolved single-molecule tracking (SMT) is used for the first time to investigate the accumulation of Rhodamine B (RhB) dye on fresh and artificially aged polyethylene terephthalate (PET) surfaces. PET thin films serve as models for microplastics, while RhB serves as a proxy for the OMs they accumulate. Artificial aging of the films is accomplished by exposing them in a UV-ozone chamber. Water contact angle, spectroscopic ellipsometry, and carbonyl index measurements reveal a gradual decrease in film hydrophobicity, thickness, and carbonyl content with age. Atomic force microscopy (AFM) data reveal an increase in surface roughness and confirm that the films remain largely intact and continuous across the aging times explored. In SMT experiments, wide-field fluorescence videos acquired from the water/PET interface under 7.5 pM RhB reveal both mobile and immobile dye molecules. Measurements of the frame-to-frame displacements of the dye show that diffusion occurs by a desorption-mediated mechanism and that the diffusion rate varies with PET film age. The surface density of mobile dye molecules decreases with increasing PET age, while the population of immobile molecules becomes relatively larger, suggesting an age-dependent transformation of the mechanism(s) by which the dye is accumulated. SMT data reveal that both mobile and immobile molecules repeatedly adsorb over the same surface sites, consistent with the emergence of nanoscale PET surface heterogeneity also revealed by AFM. Estimates of the adsorption coefficients are obtained using a nearest-neighbor analysis, giving values from 9.9 × 10 5 to 2.1 × 10 6 M –1 for immobile molecules and from 1.8 × 10 5 to 2.5 × 10 5 M –1 for mobile molecules on fresh and 5 min aged PET, respectively. Here, anomalous age-dependent variations in the velocity of molecular motion on the PET surface and in the population of immobile molecules are shown to correlate with changes in the strength of RhB adsorption.

Adsorption↗
NASA NTRS1978

An atlas of Copernicus ultraviolet spectra of Wolf-Rayet stars

An atlas of Copernicus UV scans is presented, and line identifications are tabulated, for the Wolf-Rayet stars Gamma-2 Vel (WC 8 + O7), HD 50896 (= EZ CMa; WN 5), and HD 92740 (WN 7). The atlas covers the wavelength ranges from 946.8 to 3182 A for Gamma-2 Vel, from 1012 to 1294 A for HD 50896, and from 1051 to 1243 A for HD 92740. The wavelengths include corrections for components of satellite velocity, earth velocity, and stellar heliocentric velocity; each spectral feature is classified as interstellar, photospheric, emission, UV-displaced P Cygni line absorption, or P Cygni line emission. UV-edge velocities of the P Cygni profiles are estimated, P Cygni profile types are discussed, and the results are compared with Copernicus scans of OB stars exhibiting UV P Cygni profiles. It is noted that: (1) the line-strength ratio of molecular hydrogen to atomic species appears to be substantially greater in the scans of the WN stars than in the Gamma-2 Vel scans; (2) some of the P Cygni profiles in Gamma-2 Vel differ significantly from the corresponding profiles in OB stars; and (3) there may be a slight inverse correlation between ejection velocities and excitation potentials in Gamma-2 Vel.

Johnson, H. M.↗
DOE OSTI2026

Laser interference structuring of Cu for adhesive joining of HVAC equipment

Due to potential energy and cost savings benefits, adhesive joining has been recently considered for heating, ventilation, air conditioning, and refrigeration (HVAC&R) systems. HVAC adhesive bonding requires cost-effective and adequate surface preparations of adherents. Here, this article investigates the use of abrasion, traditional single-beam laser, and a laser-interference technique as surface preparations of copper surfaces. Surface morphology is characterized using scanning electron microscopy (SEM) and atomic force microscopy (AFM). Effective submicrometer peak-to-valley structuring with a periodicity of ∼2.7μ⁢m was demonstrated for laser-interference processing. Single-lap shear tests were conducted for 89 joints made with bondline thicknesses of 0.15 and 0.3 mm. Data on process variables and measured variables included open-time, bond length, maximum load, displacement at failure, shear lap strength, and failure mode. A statistical analysis was conducted on each lot to determine the lower limit with 95% confidence intervals for displacements at failure and shear lap strengths. A comparison is presented between the properties of laser-structured joints with respect to those prepared by abrasion, which is considered the baseline surface preparation technique. Based on this comparison, one single-beam laser technique and two laser-interference techniques were shown to exhibit vastly superior performance over the joints made with abraded specimens.

Sabau, Adrian S. [Oak Ridge National Laboratory (O↗
DOE OSTI2020

Giant Domain Wall Conductivity in Self‐Assembled BiFeO 3 Nanocrystals

Abstract Ever‐increasing demand on electronic devices with ultrahigh‐density non‐volatile data storage has attracted great interest in novel ferroelectric memories based on conductive ferroelectric domain walls. Embedded in an insulating material, ferroelectric domain walls have the capability of being (re)created, displaced, erased, and altered in their spatial configurations and electronic characteristics. However, the domain wall conductivities are in most cases not yet sufficiently high to ensure the current density required to drive read‐out circuits operating at high speeds. In this work, a giant domain wall current (>10 µA) of a single charged domain wall is obtained through conductive atomic force microscopy with a bias field of 4 V. This is achieved in self‐assembled BiFeO 3 nanocrystals grown by sol‐gel method on Nb‐doped SrTiO 3 substrates. Local configurations of domains and domain wall types are studied using vector piezoresponse force microscopy and high‐resolution transmission electronic microscopy. The enhancement of the wall current is shown to be due to the formation of conducting pathways of charged defects accumulated along domain walls and traversing the nanocrystals. The diverse domain walls can be manipulated by electric field in a perpendicular architecture. The perpendicular array structure of BiFeO 3 nanocrystals should have great potentials for developing perpendicular nanoelectronic prototypes.

Liu, Lisha↗
DOE OSTI2025

Precise Linker Length and Dynamic Bond Exchange Control Penetrant Diffusion in Dense Vitrimers

Polymer networks with dynamic covalent bonds have been investigated for their self-healing ability, recyclability, and potential as more sustainable materials. Recent results have indicated that in some cases, bond exchange can enhance the transport of penetrants in dense networks, pointing to their potential for separations of membranes. Here, imine dynamic bonds in ethylene oxide (EO) networks with precise linker lengths were synthesized to investigate the transport of N,N′-bis(2,5-di-tert-butylphenyl)-3,4,9,10-perylenedicarboximide (BTBP), a large, anisotropic dye molecule. Networks with mesh sizes smaller than, comparable to, and greater than the size of the penetrant axes were investigated to probe the effects of bond exchange and network confinement on transport. Mesh sizes, which ranged from 0.5 to 1.62 nm, were determined from shear rheology, glass transitions by calorimetry, and probe diffusion coefficients by fluorescence recovery after photobleaching. Permanent networks with identical EO chain lengths were prepared as control samples, and up to a 3 orders of magnitude increase in diffusion coefficient is observed in the dynamic systems for short linkers containing 13 backbone atoms. The longest linkers with 71 backbone atoms show no difference between the permanent and dynamic networks. Linkers shorter than 11 backbone atoms, corresponding to a mesh size smaller than the penetrant small axis, diffusion is no longer observable on the experimental time scale, indicating a sharp cutoff attributed to the precise linkers and narrow mesh size distribution. The dynamic imine exchange time scales were compared to the diffusive hopping times of penetrants and indicate that exchange can occur during a diffusive displacement. Furthermore, these findings provide insights into the factors affecting penetrant transport in dense polymers and inspire the development of next-generation selective polymer membranes.

Diffusion↗
DOE OSTI2021

Multiphonon diffuse scattering in solids from first principles: Application to layered crystals and two-dimensional materials

Time-resolved diffuse scattering experiments have gained increasing attention due to their potential to reveal nonequilibrium dynamics of crystal lattice vibrations with full momentum resolution. Although progress has been made in interpreting experimental data on the basis of one-phonon scattering, understanding the role of individual phonons can be sometimes hindered by multiphonon excitations. In Ref. [M. Zacharias, H. Seiler, F. Caruso, D. Zahn, F. Giustino, P. C. Kelires, and R. Ernstorfer, Phys. Rev. Lett. 127, 207401 (2021)], we have introduced a rigorous approach for the calculation of the all-phonon inelastic scattering intensity of solids from first-principles. In the present work, we describe our implementation in detail and show that multiphonon interactions are captured efficiently by exploiting translational and time-reversal symmetries of the crystal. We demonstrate its predictive power by calculating the scattering patterns of monolayer molybdenum disulfide (MoS 2 ), bulk MoS 2 , and black phosphorus (bP), and we obtain excellent agreement with our measurements of thermal electron diffuse scattering. Remarkably, our results show that multiphonon excitations dominate in bP across multiple Brillouin zones, while in MoS 2 they play a less pronounced role. We expand our analysis for each system and examine the effect of individual atomic and interatomic vibrational motion on the diffuse scattering signals. We further demonstrate the high-throughput capability of our approach by reporting all-phonon scattering maps of two-dimensional MoSe 2 , WSe 2 , WS 2 , graphene, and CdI 2 , rationalizing in each case the effect of multiphonon processes. As a side point, we show that the special displacement method reproduces the thermally distorted configuration that generates precisely the all-phonon diffuse pattern. In conclusion, the present methodology opens the way for systematic calculations of the scattering intensity in crystals and the accurate interpretation of static and time-resolved inelastic scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

Energy Minimization of Molecular Features Observed on the (110) Face of Lysozyme Crystals

Molecular dynamics and energy minimization have been carried out using the program XPLOR to check the plausibility of a model lysozyme crystal surface. The molecular features of the (110) face of lysozyme were observed using atomic force microscopy (AFM). A model of the crystal surface was constructed using the PDB file 193L, and was used to simulate an AFM image. Molecule translations, van der Waals radii, and assumed AFM tip shape were adjusted to maximize the correlation coefficient between the experimental and simulated images. The highest degree of 0 correlation (0.92) was obtained with the molecules displaced over 6 A from their positions within the bulk of the crystal. The quality of this starting model, the extent of energy minimization, and the correlation coefficient between the final model and the experimental data will be discussed.

Perozzo, Mary A.↗
DOE OSTI2024

Effect of impurities on phase transformation and precipitation in a low-carbon steel

Impurity elements have been added to a commercial low carbon steel grade to simulate the levels that could arise as a result of increased scrap recycling during steel production. In this study, continuous cooling transformation (CCT) diagrams were constructed for the steels with varying levels of impurities, and it is shown that impurities suppress the phase transformation across a wide range of cooling rates. It was found that a step was formed in the start temperature curve, separating the reconstructive and displacive transformations. The influence of impurities on both the reconstructive transformation and displacive transformation are discussed. Additionally, Cu precipitates were observed using scanning transmission electron microscopy (STEM) in the highest impurity-containing steel after slow cooling (0.05°C/s), fast cooling (5°C/s) and interrupted cooling. It was found that the precipitation kinetics is in the following order: cementite within the secondary phase and cementite at secondary phase-ferrite interface> ferrite grain boundaries> ferrite grain matrix. Atom-probe tomography (APT) revealed Cu precipitates formed on the surface of cementite lamellae, but not within it. This work offers insights for the phase transformation control and precipitation regulation during the thermomechanical processing of low carbon steels containing impurity elements due to scrap recycling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗
DOE OSTI2021

Site 2 of the Yersinia pestis substrate-binding protein YfeA is a dynamic surface metal-binding site

The substrate-binding protein YfeA (also known as YPO2439 or y1897) is a polyspecific metal-binding protein that is crucial for nutrient acquisition and virulence in Yersinia pestis, the causative microbe of plague. YfeA folds into a monomeric c-clamp like other substrate-binding proteins and has two metal-binding sites (sites 1 and 2). Site 2 is a bidentate surface site capable of binding Zn and Mn atoms and is a unique feature of YfeA. Occasionally, the site 2 residues of two YfeA molecules will cooperate with the histidine tag of a third YfeA molecule in coordinating the same metal and lead to metal-dependent crystallographic packing. Here, three crystal structures of YfeA are presented at 1.85, 2.05 and 2.25 Å resolution. A comparison of the structures reveals that the metal can be displaced at five different locations ranging from ~4 to ~16 Å away from the canonical site 2. These observations reveal different configurations of site 2 that enable cooperative metal binding and demonstrate how site 2 is dynamic and freely available for inter-protein metal coordination.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗
NASA NTRS2008

Multiscale Modeling of Intergranular Fracture in Aluminum: Constitutive Relation For Interface Debonding

Intergranular fracture is a dominant mode of failure in ultrafine grained materials. In the present study, the atomistic mechanisms of grain-boundary debonding during intergranular fracture in aluminum are modeled using a coupled molecular dynamics finite element simulation. Using a statistical mechanics approach, a cohesive-zone law in the form of a traction-displacement constitutive relationship, characterizing the load transfer across the plane of a growing edge crack, is extracted from atomistic simulations and then recast in a form suitable for inclusion within a continuum finite element model. The cohesive-zone law derived by the presented technique is free of finite size effects and is statistically representative for describing the interfacial debonding of a grain boundary (GB) interface examined at atomic length scales. By incorporating the cohesive-zone law in cohesive-zone finite elements, the debonding of a GB interface can be simulated in a coupled continuum-atomistic model, in which a crack starts in the continuum environment, smoothly penetrates the continuum-atomistic interface, and continues its propagation in the atomistic environment. This study is a step towards relating atomistically derived decohesion laws to macroscopic predictions of fracture and constructing multiscale models for nanocrystalline and ultrafine grained materials.

Yamakov, V.↗