Molecular Dynamics Simulations of Aluminium Foams under Tension

Publication Date (Web): October 24, 2018. Copyright © 2018 American Chemical Society. Cite this:J. Phys. Chem. C XXXX, XXX, XXX-XXX ...
0 downloads 0 Views 7MB Size
Subscriber access provided by The University of Texas at El Paso (UTEP)

C: Physical Processes in Nanomaterials and Nanostructures

Molecular Dynamics Simulations of Aluminium Foams under Tension: Influence of Oxidation Nina Gunkelmann, Eduardo M. Bringa, and Yudi Rosandi J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b07146 • Publication Date (Web): 24 Oct 2018 Downloaded from http://pubs.acs.org on October 30, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Molecular Dynamics Simulations of Aluminium Foams under Tension: Influence of Oxidation Nina Gunkelmann,∗,† Eduardo M. Bringa,‡ and Yudi Rosandi¶ †Chair of Computational Material Sciences/Engineering, Institute of Applied Mechanics, Technische Universit¨at Clausthal, Arnold-Sommerfeld-Straße 6, 38678 Clausthal-Zellerfeld, Germany ‡CONICET and Facultad de Ingenier´ıa, Universidad de Mendoza, Mendoza, 5500 Argentina ¶Department of Geophysics, Universitas Padjadjaran, Jatinangor, Sumedang 45363, Indonesia E-mail: [email protected] October 24, 2018 Abstract For materials with high surface-to-volume ratio and high oxygen affinity, oxide layers will significantly change the material properties. However, oxidation effects have not been studied for metal nanofoams which have many applications because of their light weight and high stiffness. We use molecular dynamics simulations to show that oxidized Aluminum nanofoams possess significantly improved ductility without reduction in tensile strength. The Al-O interface leads to an increased defect nucleation rate at the oxide interface preventing localized deformation. At the same time, the enthalpy of mixing between Aluminiun and oxygen decreases for increasing O concentrations, reaching a minimum at the stoichiometric ratio of Al2 O3 , resulting in stabilized bonds and increased strength.

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction Metal foams have many technological applications. They are widely used as shock absorbers due to their light weight and high stiffness. 1–3 They are used for industrial and energy-storage applications, 4 and can be applied as radiation shields in space science. 5 There are experimental and simulation studies dealing with macroscopic metal or plastic foams. 6–8 On the other hand, experiments on foams with nanoscale porosity show large yield strength, 9 and nanoporous metals have been widely used for functional nanostructures such as sensors 10–12 and nanophotonic devices. 13 Molecular Dynamics (MD) simulations can often be used to investigate their response. For instance, many authors have reported results on uniaxial deformation and compression of nanofoams at the atomic level. 14–17 The evolution of foamed aluminium at high rate tension was recently investigated by Mayer et al. using a mechanical model based on molecular dynamics simulations. 18 Zhao et al. 19 studied shock induced melting, finding internal jetting from pores. Xiang et al. investigated the shock response of nanoporous Aluminum and revealed two mechanisms for void collapse: the plasticity mechanism and the internal jetting mechanism. 20 Gunkelmann et al. 21 also studied the compaction and plasticity in an Al nanofoam. Simulations focus on pure, single-element foams but, in atmospheres with oxygen content, oxide layers will change the material properties. For materials with high surface-to-volume ratio and high oxygen affinity this effect is especially pronounced. At aluminum surfaces, an oxide layer may form in seconds, even under vacuum conditions. 22 However, the mechanics describing oxidation behavior are still not fully understood. An efficient oxygen reduction catalysis was found for subnanometer Pt alloy nanowires prepared by a new experimental approach. 23 Recently, foaming of aluminum was investigated under oxidizing and non-oxidizing gas atmospheres using X-ray radioscopy to study the structure and distribution of the oxides. 24 Despite the importance of oxidation to determine the properties of materials, only rela2

ACS Paragon Plus Environment

Page 2 of 25

Page 3 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

tively few simulation studies exist on this topic. There are detailed studies of oxidation of Al nanoparticles, starting with the seminal study by Campbell et al., 25 including the study of Hong et al., 26 and more recent studies. 27,28 Recently, Sun et al. analyzed the oxidation of Fe nanoparticles in solution by small- and wide-angle x-ray scattering and molecular dynamics simulations. 29 They observed the formation of voids within the nanoparticles accompanied by mass diffusion into and out of the nanoparticles. Sen et al. showed with MD simulations that oxidation leads to increased ductility for Al nanowires as a consequence of increased dislocation nucleation due to increased activation volume and increased number of dislocation sites. 22 There are even fewer structures on the role of oxidation in porous structures. The shock Hugoniot of full-density and porous CeO2 was investigated in the liquid regime using ab initio and classical molecular dynamics. 30 The study points to the necessity of acquiring data in the porous regime to increase the reliability of equation of state models. In the present paper we use MD simulations to study tension of an Al nanofoam in an oxygen atmosphere. We characterize the stress-strain response of Al foams: pure, covered by an oxide layer, and in an O2 environment, showing that oxidized foams display increased ductility. We also calculate the mixing enthalpy which is minimal for the stoichiometric ratio of Al2 O3 . Our results are relevant for understanding the behavior of nanoporous structures at ambient conditions.

Methods We construct foams with initial filling factor of φ0 = 0.5. The construction algorithm takes advantage of the non-simply connected structure of atoms of fixed temperature in a liquid and is described in. 21 In brief, we first create a template for the foam structure by starting with an fcc crystal with periodic boundary condition using a simple LennardJones (LJ) potential. The target crystal is then heated above the melting temperature for

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

several thousand time steps. The atoms do not all have the same temperature because of equilibrium fluctuations. From the melt we delete atoms with temperature above a certain value, until the desired filling factor is reached. We obtain a spongy structure, which serves as a template for our nanofoam. Thus, we create the crystalline Al foam by using it as a spatial filter to remove all superfluous atoms from an fcc Al crystal. Note that some atoms are isolated inside voids. We relax the structure using a short-ranged LJ potential with a strongly increased binding energy to attach them to filaments. In a final step, all non-connected atoms and clusters are removed by using a cluster detector 31 which isolates the largest interconnected structure. This structure is used as as template to create the crystalline Al foam by using it as a spatial filter to remove all superfluous atoms from an fcc Al crystal. Another often applied procedure to generate Au nanofoams makes use of a phase-field model simulating the spinodal decomposition of a binary alloy. 32,33 The authors found that the obtained microstructure is similar to a microstructure obtained using the Berk analysis based on small angle neutron scattering data taken for porous gold. Our generation procedure leads to similar foam structures.

Figure 1: Different simulation configurations. From left to right: Foam in an O2 environment filling the voids of the foam (‘O2’), oxidized Al in vacuum and Al nanofoam. Red: O; blue: Al. The number of oxygen atoms in the initial cubic box was 64000. The foams have a cubic shape with a side length of 13 nm and the total number of Al atoms amounts to 69367. The foam was placed in a cubic box with a number of O atoms varying from 32000 to

4

ACS Paragon Plus Environment

Page 4 of 25

Page 5 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

80000 in increments of 16000 atoms using the software PACKMOL. 34 The oxygen density was varying between 0.23 and 0.60 g/cm3 , several hundred times that of ambient conditions. We deleted all O atoms overlapping with the Al nanofoam. In detail, pairs of O and Al atoms whose distance of separation is smaller than the sum of their radii are searched for, and O atoms are deleted. After 100 ps relaxation in a NPH ensemble an amorphous oxide layer has formed at the Al surface. According to Hong et al. , oxidation of nanoparticles reaches a limiting thickness after approximately 100 ps, for T < 570 K, 26 and this thickness compares well with experimental thickness meassured after several 100 seconds. Sun et al. 29 compared oxidation kinetics from MD at the ns scale, with experimental kinectics at the minute scale. Both of these simulations use relatively high oxygen gas pressure in order to achieve comparable oxidation at available simulation time scales (initial oxygen density of 0.48g/cm3 compared to approximately 1e-4 g/cm3 in experiments). Kinetics depends on a number of factors, including structure, diffusivities and temperature. 35 Regarding mechanical deformation, we are carrying simulations at high strain rate, where it takes 100 ps to reach 10 % deformation (more details see below), and it takes the same time to reach a steady oxide layer thickness. This might compare well with experiments, where both of these processes take minutes. Plasticity would be controlled by dislocation nucleation and propagation, and by bond-breaking, which take a much shorter time than oxidation reaching a steady thickness. From mean square displacement (MSD) calculations, for instance in, 28 O and Al atoms would move more than 2 nm after 100 ps, traversing the filaments in our simulated foam and allowing partial healing during necking. We removed the atoms outside the simulation box of the Al nanofoam to ensure that both the Al foam and the oxidized foam have the same stress state and strain rate upon deformation. The number of O atoms after this step amounts to [15039, 21678, 29391, 31911], for [32000, 48000, 64000, 80000] starting atoms. To evaluate the difference between oxidized Al in vacuum (in the following denoted as ‘O2 layer’) and in an O2 environment

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 25

filling the voids of the foam (‘O2’) we removed the oxygen atoms inside the holes which were not attached to the surface by using a cluster detector. The number of O atoms for the oxidized sample is [11523, 17319, 20978, 23962] referring to [32000, 48000, 64000, 80000] O atoms in the initial cubic box. The concentration of O atoms for this configuration thus amounts to c =[0.14, 0.20, 0.23, 0.26]. The different structures are shown in Fig. 1. We relax the samples using high-temperature annealing at 80 % of the melting temperature with a Nos´e/Hoover isenthalpic ensemble (NPH) during 100 ps. 36 The samples were then cooled back to 10 K followed by 2 ps in a microcanonical NVE ensemble. The heating and cooling ramps were performed during 2 ps. The simulations are conducted at a temperature of 10 K in an NVT ensemble in order to minimize thermal noise. The open-source MD code LAMMPS 37 is used in this work to perform the simulations. We apply periodic boundary conditions. For our simulations we employ the reactive force field ReaxFF, including charge transfer between aluminum and oxygen molecules, 38,39 with the implementation by Aktulga et al. 40 The charge is equilibrated instantaneously at each molecular dynamics timestep to minimize the Coulomb energy. An elaborated force field potential equation can be derived in terms of different interactions factors as:

Epot = Ebond + Eover + Eangle + Etors +EvdWaals + ECoulomb + Especific .

(1)

Here, the bonding interaction represents ionic and covalent components, while van-der-Waals interactions are attributed to non-bonding terms. Eangle represents the deviation of the bond angle from equilibrium described by a harmonic term, Etors describes the four-body torsional angle strain. Ebond is a continuous function of interatomic distance and describes the energy associated with forming bonds between atoms. Eover describes an energy penalty term preventing the over coordination of atoms. The non-bonding interactions ECoulomb and EvdWaals

6

ACS Paragon Plus Environment

Page 7 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

are electrostatic and dispersive contributions quantifying long-range interaction between neutral atoms, while Especific represents specific energy contributions of the system, capturing properties particular to the system of interest. The potential we use was fit to describe Al2 O3 , but some Alx Oy clusters were included in the fit. 39 Zeng et al. describe simulations with the same potential used in this work, and the presence of sub-oxides in Al nanoparticles compares favorably with experiments. 28 Douglas-Gallardo et al. recently found good agreement between ReaxFF and DFT for different oxidation degrees in Al nanoclusters. 41 The simulations were carried out by performing uniaxial tension tests along the z axis at a strain rate of 109 s−1 for 320 ps.

Results The shear stress describing the onset of plasticity is defined as 1 pshear = − (pzz − ptrans ), 2

(2)

where pij denote the components of the stress tensor, and the transverse stress is defined as 1 ptrans = (pxx + pyy ). 2

(3)

For tensile stresses the shear stress definition above is equal to the von Mises stress if offdiagonal terms are ignored 42 and it is often used for the analysis of MD simulations. 43 In our simulation, non-diagonal stress components are small. Fig. 2 shows the shear stress in the tension experiment versus strain for pure Al, with oxide layer and in an O2 environment. For pure Al foams the shear stress shows a narrow linear elastic regime up to 3 % strain until dislocation emission leads to a maximum in the stress and subsequent plastic relaxation after 11 % strain. We observe increased ductility for Al covered with an oxide layer as well as for foams in an oxygen environment. In contrast to

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the curve for pure Al, for both the foam with the oxide layer and in the oxygen cloud, the shear stress only slowly decreases after the maximum at around 17 % strain and still has a magnitude of 4 GPa at 60 % strain. In addition, the strength is increased in both cases. Note that the strength is higher in the O2 environment because the number of O atoms is higher in this sample. We note that recent simulations of collision of Al oxide nanoparticles show much smaller maximum shear values, around 5 GPa. 27 The magnitude of the shear stress depends on the number of O atoms. Fig. 3 shows the shear stress versus strain for an Al nanofoam covered by amorphous oxide (‘O2 layer’) with different numbers of O atoms, thus different O concentrations. Although the shear stress increases for increasing number of oxygen atoms, the curves do not change qualitatively for c ≥ 0.14. We observe a maximum at around 12 % strain and subsequent plastic relaxation. If we increase the number of O atoms further the concentration will not increase anymore because the surface contacts are saturated. Only for small oxygen concentrations, c = 0.02, the behavior is similar to the pure Al foam showing decreasing shear stress after 7 % strain. For visualization of the atomistic configurations we use the adaptive common-neighbor analysis together with the dislocation extraction algorithm (DXA) 44 within the free software tool OVITO. 45 We display in Fig. 4 the tension behavior of the pure Al foam, the stress induced changes in our foam covered with an oxide layer and in an oxygen environment. The oxide layer has a thickness of 2 nm. Up to 31 % strain the structure of all foams does not change significantly during tension. Then, we observe fracture occurring in the center of the foam. For the pure Al foam, we observe a defined fracture surface. On the contrary, for the oxidized foam we do not have a complete rupture of the material but chains of Al-O atoms stick the material together, even up to 50 % strain. This effect is even more pronounced for the foam in an oxygen environment. Here, the foam appears to be stable up to 41 % strain. At 51 % strain chains of Al-O connect the upper and lower part of the foam similar to the foam covered by an oxide layer. A reason could be that strain energy is absorbed through

8

ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

breaking of bonds between highly coordinated atoms (see Fig. 7, below). This finding was also observed by Gao et al. for silica and alumina coated Silicon nanowires. 46 From Fig. 5 we see that many Shockley partial dislocations nucleate at the Al-O interface for oxidized foams; 73 % of the detected dislocations are Shockley partials. The reason is the small size of the filaments which cannot accommodate a full stacking fault ribbon. The disordered structure at the Al-O interface eases dislocation nucleation so that they mainly originate at the interface. Further dislocation pictures showing nucleation at the interfaces and stacking faults at different strains for the sample covered by the oxide layer can be seen in the Supplemental Material. The dislocation evolution is in agreement with the results by Sen et al. 22 where it was shown that the oxide shell decreases the aluminum dislocation nucleation stress by increasing the activation volume and the number of nucleation sites in Al nanowires. By increasing the defect nucleation rate at the interface, localized strains at the glide planes of dislocations are reduced and toughness is increased. In contrast, for pure Al dislocations do not only nucleate at the surface of the foam but are homogeneously distributed inside the filaments. Here, stair-rod dislocations cross the filaments. Although the total amount of dislocations is comparable for both samples the dislocation analysis is more reliable for pure Al. The reason is that the detector cannot sufficiently analyze the structure at the Al-O interface resulting in a smaller defect surface for the oxidized foam. Slices of the oxidized sample at different strains colored by their lattice types can be found in the Supplemental material showing that the majority of the foam is detected as disordered (up to 60 %). Our results for Al nanofoam are similar to the results for a single Al nanowire by Sen et al. 22 This is not unexpected, since a nanofoam has often been considered as a collection of connected nanowires. 47 This approach, however, does not take into account plasticity at nanowire junctions, which could control the overall mechanical response. 48 In Fig. 6 we color the atoms of the pure and oxidized foam by the shear strain invariant EvM which can be computed from the components of the atomic strain tensor Ei . For details

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of the calculation of the atomic strain tensor please refer to. 49 r EvM =

 1 2 + E2 + E2 (Exx − Eyy )2 + (Exx − Ezz )2 + (Ezz − Eyy )2 + 6 · Exy xz yz . 2

We observe that strain is localized in the center of the specimen for the Al foam. In contrast, the shear strain appears less concentrated for the foam covered by the oxide layer. Here, the material has deformed by dislocations nucleated at the Al-O interface. Shear localization is related to bond-breaking, and Al-O bonds have an important role modifying ductility. To quantify this, Fig. 7 shows the number of bonds in the sample. The pure Al sample displays a large localized fracture starting at 20 % strain, leaving only a few connecting thin wires, shown in Fig. 4. This is the reason why there is a significant number of surviving Al-Al bonds (about 99 %), and the number of bonds stabilizes because the upper and bottom part are no longer connected. For the samples with O, O mobility leads to bond-healing, and there is creation of new Al-O bonds which delay fracture, as already noted by Sen et al. for a single nanowire. 22 Note that the creation of new Al-O bonds originates mostly from rearrangement of the AlOx layer while the number of molecular O2 remains almost constant. Bond breaking is not localized in a single region, but spread. For the sample with O2 atmosphere, there is nearly 6 % of Al-Al bond breaking, compensated by 6 % Al-O bond formation. For this sample, the O-O bonds also increase, unlike what happens for the sample without the O2 atmosphere, where Al-O bond formation is also more modest. The underlying mechanism leading to increased ductility without reduction in strength can therefore be understood from two observations. First, oxygen absorbs strain energy through breaking of highly coordinated bonds as explained in Fig. 7. Due to the strong interaction between O and Al atoms, bond breaking is not localized but distributed across a larger volume. Second, localized deformation is prevented by increased defect nucleation rate from the Al-O interface (Fig. 5). To find out the correct concentration of O in the oxidized Al target, we measure the

10

ACS Paragon Plus Environment

Page 10 of 25

Page 11 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

mixing enthalpy of the system with varying concentration. We take a mixture of Al and O, ˚3 . To obtain the desired concentration initialized by Alumina structure of size 33x29x26 A we randomly remove Al or O atoms from the target. The system is then relaxed using an NpT ensemble to 1 K. Upon annealing the sample is amorphized due to the breaking of bonds of the corundum structure. The enthalpy of mixing ∆H is measured using a modified formula by Delannay et al. 50 The plot of ∆H is shown in Fig. 8. The figure shows that O atoms only do not mix with Al at very small concentration. At concentrations > 0.15 oxygen atoms are always bounded to Al atoms. The minimum shows the correct stoichiometric ratio of Alumina emphasizing the reliability of the potential used in this study. The increasing ductility, shown in Fig. 3, corresponds to more negative values of the heat of mixing. At these concentrations, it is shown that the bonds between O and Al are preferable and play an important role in modifying ductility as seen in Fig. 7.

Conclusions Using MD simulations we studied the response under uniaxial tension of pure Al nanofoams and Al nanofoams containing different oxygen concentration. Our results show that oxygen helps to increase the ductility without decreasing the tensile strength of Al nanofoams. The total amount of dislocations under strain is comparable for pure Al and oxidized foams, but for Al-O Shockley partials nucleate at the interface bettween the metal and the oxide. The enthalpy of mixing is positive for low oxygen concentrations and its minimum corresponds to the stoichiometric relation of corundum leading to increased stability of the bonds. These results are relevant for designing nanofoams at ambient conditions. The increased toughness could lead to an improved long-term stability of nanodevices.

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

12

O2 O2 layer Al

10 shear stress (GPa)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 25

8 6 4 2 0 -2 0

0.1

0.2

0.3

0.4

0.5

0.6

strain

Figure 2: Shear stress versus strain for pure Al, an Al nanofoam covered by amorphous oxide (‘O2 layer’) and an Al foam in an oxygen environment (‘O2’). The number of oxygen atoms in the initial cubic box was 64000. The blue arrows mark the linear elastic regime up to 3 % strain and the inversion point for pure Al at 11 % strain where the shear stress starts decreasing. The green arrow shows the maximum shear stress for O2 at 0.17 strain.

12

ACS Paragon Plus Environment

Page 13 of 25

12

c=0.02 c=0.14 c=0.20 c=0.23 c=0.26

10 shear stress (GPa)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

8 6 4 2 0 -2 0

0.1

0.2

0.3

0.4

0.5

0.6

strain

Figure 3: Shear stress versus strain for an Al nanofoam covered by amorphous oxide (‘O2 layer’) with different concentrations c of O atoms. The arrow denotes the maximum of the shear stress. The mechanical behavior of the foam under tension is somewhat similar to the mechanical behavior of a single nanowire, but the sudden failure of a single wire, 22 is replaced by extended ductility, because there is a collection of wires with irregular cross-sections and junctions. It would be interesting to study an oxidation scenario under compression, since there is a large tension-compression asymmetry in foams, 17 which might be modified by oxygen on the foam surface, affecting surface stress. 51 Given that temperature affects significantly oxygen mobility, 28 ductility of the oxidized foam might be enhanced even further with respect to the pure foam as temperature is increased. Al oxide, as well as other metallic oxide materials display reduced ductility when compared to their pure metallic counterparts. Oxide layers fracture and lead to failure, 52 as shown for instance in early work on thin Al films. 53 There must be a critical oxide thickness where there is a transition from this cracking behavior observed at the micron scale, to the assisted-ductility behavior observed at the nano-scale. Our oxide layer is amorphous, while

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4: Snapshots of the foams at different strains. Top: pure Al nanofoam, middle: Al nanofoam covered by amorphous oxide (‘O2 layer’), bottom: foam in an oxygen environment (‘O2’). Red: O; blue: Al. The number of oxygen atoms in the initial cubic box was 64000.

14

ACS Paragon Plus Environment

Page 14 of 25

Page 15 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5: Dislocation pictures for the foam covered by an oxide layer (left) and the Al foam (right) at 30 % strain. White: defect surface; dislocations: green: Shockley partials 16 h112i, gray: stair-rod dislocations 16 h110i, red: other dislocations.

Figure 6: Snapshots for the foam covered by an oxide layer (left) and the Al foam (right) at 30 % strain. Atoms are colored by shear strain.

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1.01

1.08

Al Al-Al O-O Al-O

Al Al-Al O-O Al-O

1.06 Relative number of bonds

1 Relative number of bonds

0.99 0.98 0.97 0.96 0.95

1.04 1.02 1 0.98 0.96

0.94

0.94 0

0.2

0.4

0.6

0.8

1

0

strain

0.2

0.4

0.6

0.8

1

strain

Figure 7: Relative number of bonds versus strain for the pure Al foam and the sample covered by an oxide layer (top) and in an oxygen environment (bottom). Bonds between O-O (rcut =0.4 nm), Al-O (rcut =0.25 nm) and Al-Al (v=0.35 nm). Numbers of bonds have been normalized to the number of bonds at 0 % strain. Under tension, bond breaking is expected and would lead to a decreasing number of bonds, as observed in all samples for Al-Al bonds. However, O mobility leads to formation of an increasing number of Al-O bonds, which allow extended ductility of samples with O, compared to the pure Al sample.

5 0 -5 ∆H (kcal/mol)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 25

-10 -15 -20 -25 -30 -35 0

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 concentration

1

Figure 8: Mixing enthalpy versus oxygen concentration, calculated using ReaxFF. The minimum at concentration of 0.6 corresponds to the stoichiometric ratio of Alumina .

16

ACS Paragon Plus Environment

Page 17 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

macroscopic layers are polycrystalline and fracture is aided by grain boundaries and crystal plasticity in the oxide. Further research is needed to understand the transition and the critical size at which cracking start occurring. Most studies of mechanical degradation of Al samples focus on the influence of hydrogen on embrittlement and failure. H accumulates at the interface between Al and the oxide layer, leading to blistering, 54,55 especially under irradiation conditions. 56 Future simulation studies might consider the early stages of H-bubble formation at the Al-AlOx interface, which under tension would provide additional pathways for H diffusion and bubble growth. Finally, we note that metallic nanofoams have been shown to be radiation resistant, due to their free surfaces acting as defect sink. 47 Indentation of irradiated Au nanofoam leads to hardness increase under irradiation due to defect accumulation 57 and to topology changes. 58 Oxidation would modify foam evolution under irradiation, and competing effects might not lead to hardening. However, it is likely that oxidized nanofoams would provide additional defect sinks at the oxide-metal interface, likely increasing radiation tolerance.

Supporting Information • Additional dislocation pictures for the foam covered by an oxide layer at different strains. • Slices of the oxidized sample at different strains.

Acknowledgement This work has been supported by the Simulation Science Center Clausthal-G¨ottingen. EMB thanks support from grants PICT2014-0696 and SeCTyP-UNCuyo. The authors acknowledge the North-German Supercomputing Alliance (HLRN) for providing HPC resources that have contributed to the research results reported in this paper.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References (1) Santosa, S.; Wierzbicki, T. Crash behavior of box columns filled with aluminum honeycomb or foam. Computers and Structures 1998, 68, 343–367. (2) Villanueva, G. R.; Cantwell, W. J. The high velocity impact response of composite and FML-reinforced sandwich structures. Composites Science and Technology 2004, 64, 35 – 54. (3) Bourne, N. Materials in mechanical extremes; Cambridge Univ. Press, 2013. (4) Kim, S.; Lee, C.-W. A review on manufacturing and application of open-cell metal Foam. Procedia Materials Science 2014, 4, 305 – 309, 8th International Conference on Porous Metals and Metallic Foams. (5) Thibeault, S. A.; Kang, J. H.; Sauti, G.; Park, C.; Fay, C. C.; King, G. C. Nanomaterials for radiation shielding. MRS Bulletin 2015, 40, 836–841. (6) Boade, R. R. Compression of porous copper by shock waves. J. Appl. Phys. 1968, 39, 5693–5702. (7) Tan, P. J.; Harrigan, J. J.; Reid, S. R. Inertia effects in uniaxial dynamic compression of a closed cell aluminium alloy foam. Mater. Sci. Tech. 2002, 18, 480. (8) Dezulian, R. et al. Hugoniot data of plastic foams obtained from laser-driven shocks. Phys. Rev. E 2006, 73, 047401. (9) Biener, J.; Hodge, A. M.; Hamza, A. V.; Hsiung, L. M.; Satcher, J. H. Nanoporous Au: A high yield strength material. J. Appl. Phys. 2005, 97, 024301. (10) Erlebacher, J.; Aziz, M. J.; Karma, A.; Dimitrov, N.; Sieradzki, K. Evolution of nanoporosity in dealloying. Nature 2001, 410, 450.

18

ACS Paragon Plus Environment

Page 18 of 25

Page 19 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(11) Biener, J.; Hodge, A. M.; Hayes, J. R.; Volkert, C. A.; Zepeda-Ruiz, L. A.; Hamza, A. V.; Abraham, F. F. Size effects on the mechanical behavior of nanoporous Au. Nano Lett. 2006, 6, 2379. (12) Hopkins, P. E.; Norris, P. M.; Phinney, L. M.; Policastro, S. A.; Kelly, R. G. Thermal conductivity in nanoporous gold films during electron-phonon nonequilibrium. J. Nanomater. 2008, 2008, 418050. (13) Gibson, L. J. Mechanical behavior of metallic foams. Annu. Rev. Mater. Sci. 2000, 30, 191–227. (14) Sun, X.-Y.; Xu, G.-K.; Li, X.; Feng, X.-Q.; Gao, H. Mechanical properties and scaling laws of nanoporous gold. Journal of Applied Physics 2013, 113, 023505. (15) Ruestes, C. J.; Bringa, E. M.; Stukowski, A.; Rodriguez-Nieva, J. F.; Tang, Y.; Meyers, M. A. Plastic deformation of a porous bcc metal containing nanometer sized voids. Comput. Mater. Sci. 2014, 88, 92 – 102. (16) Ngˆo, B.-N. D.; Stukowski, A.; Mameka, N.; Markmann, J.; Albe, K.; Weissm¨ uller, J. Anomalous compliance and early yielding of nanoporous gold. Acta Materialia 2015, 93, 144 – 155. (17) Farkas, D.; Caro, A.; Bringa, E. M.; Crowson, D. Mechanical response of nanoporous gold. Acta Materialia 2013, 61, 3249 – 3256. (18) Mayer, P. N.; Mayer, A. E. Evolution of foamed aluminum melt at high rate tension: A mechanical model based on atomistic simulations. J. Appl. Phys. 2018, 124, 035901. (19) Zhao, F. P.; An, Q.; Li, B.; Wu, H. A.; Goddard, W. A.; Luo, S. N. Shock response of a model structured nanofoam of Cu. J. Appl. Phys. 2013, 113, 063516. (20) Xiang, M.; Cui, J.-Z.; Yang, Y.; Liao, Y.; Wang, K.; Chen, Y.; Chen, J. Shock responses

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of nanoporous aluminum by molecular dynamics simulations. Int. J. Plast. 2017, 97, 24–45. (21) Gunkelmann, N.; Rosandi, Y.; Ruestes, C. J.; Bringa, E. M.; Urbassek, H. M. Compaction and plasticity in nanofoams induced by shock waves: A molecular dynamics study. Comput. Mater. Sci. 2016, 119, 27–32. (22) Sen, F. G.; Alpas, A. T.; van Duin, A. C. T.; Qi, Y. Oxidation-assisted ductility of aluminium nanowires. Nature Comm. 2014, 5, 3959. (23) Jiang, K.; Zhao, D.; Guo, S.; Zhang, X.; Zhu, X.; Guo, J.; Lu, G.; Huang, X. Efficient oxygen reduction catalysis by subnanometer Pt alloy nanowires. Science Advances 2017, 3 . (24) Dudka, A.; Garc´ıa-Moreno, F.; Wanderka, N.; Banhart, J. Structure and distribution of oxides in aluminium foam. Acta Mater. 2008, 56, 3990–4001. (25) Campbell, T.; Kalia, R. K.; Nakano, A.; Vashishta, P.; Ogata, S.; Rodgers, S. Dynamics of oxidation of aluminum nanoclusters using variable charge molecular-dynamics simulations on parallel computers. Phys. Rev. Lett. 1999, 82, 4866–4869. (26) Hong, S.; van Duin, A. C. Molecular dynamics simulations of the oxidation of aluminum nanoparticles using the ReaxFF Reactive Force Field. J. Phys. Chem. C 2015, 119, 17876–17886. (27) Ma, B.; Zhao, F.; Cheng, X.; Miao, F.; Zhang, J. The mechanical and thermal responses of colliding oxide-coated aluminum nanoparticles. J. Appl. Phys. 2017, 121, 145108. (28) Zeng, H.; Cheng, X.; Zhang, C.; Lu, Z. Responses of core-Shell Al/Al2O3 nanoparticles to heating: ReaxFF molecular dynamics simulations. J. Phys. Chem. C 2018, 122, 9191–9197.

20

ACS Paragon Plus Environment

Page 20 of 25

Page 21 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(29) Sun, Y.; Zuo, X.; Sankaranarayanan, S. K. R. S.; Peng, S.; Narayanan, B.; Kamath, G. Quantitative 3D evolution of colloidal nanoparticle oxidation in solution. Science 2017, 356, 303–307. (30) Weck, P. F.; Cochrane, K. R.; Root, S.; Lane, J. M. D.; Shulenburger, L.; Carpenter, J. H.; Sjostrom, T.; Mattsson, T. R.; Vogler, T. J. Shock compression of strongly correlated oxides: A liquid-regime equation of state for cerium(IV) oxide. Phys. Rev. B 2018, 97, 125106. (31) Stoddard, S. D. Identifying clusters in computer experiments on systems of particles. J. Comput. Phys. 1978, 27, 291–293. (32) Crowson, D. A.; Farkas, D.; Corcoran, S. G. Geometric relaxation of nanoporous metals: The role of surface relaxation. Scr. Mater. 2007, 56, 919. (33) Crowson, D. A.; Farkas, D.; Corcoran, S. G. Mechanical stability of nanoporous metals with small ligament sizes. Scr. Mater. 2009, 61, 497. (34) Mart´ınez, L.; Andrade, R.; Birgin, E. G.; Mart´ınez, J. M. PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations. J. Comp. Chem. 2009, 30, 2157–2164. (35) Saceleanu, F.; Atashin, S.; Wen, J. Z. Investigation of the effects of phase transformations in micro and nano aluminum powders on kinetics of oxidation using thermogravimetric analysis. Phys. Chem. Chem. Phys. 2017, 19, 18996–19009. (36) Gunkelmann, N.; Bringa, E. M.; Kang, K.; Ackland, G. J.; Ruestes, C. J.; Urbassek, H. M. Polycrystalline iron under compression: Plasticity and phase transitions. Phys. Rev. B 2012, 86, 144111. (37) Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19, http://lammps.sandia.gov/ . 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(38) van Duin, A. C. T.; Dasgupta, S.; Lorant, F.; Goddard III, W. A. ReaxFF: A reactive force field for hydrocarbons. J. Phys. Chem. A 2001, 105, 9396–9409. (39) Zhang, Q.; C ¸ aˇg ın, T.; van Duin, A.; Goddard, W. A.; Qi, Y.; Hector, L. G. Adhesion and nonwetting-wetting transition in the Al/α − Al2 O3 interface. Phys. Rev. B 2004, 69, 045423. (40) Aktulga, H. M.; Fogarty, J. C.; Pandit, S. A.; Grama, A. Y. Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques. Parallel Comput. 2012, 38, 245 – 259. (41) Douglas-Gallardo, O. A.; Soldano, G. J.; Mariscal, M. M.; S´anchez, C. G. Effects of oxidation on the plasmonic properties of aluminum nanoclusters. Nanoscale 2017, 9, 17471–17480. (42) von Mises, R. Mechanik der festen K¨orper im plastisch-deformablen Zustand. Nachrichten von der K¨oniglichen Gesellschaft der Wissenschaften zu G¨ottingen. Mathematisch-physikalische Klasse 1913, 1913, 582. (43) Tramontina, D.; Erhart, P.; Germann, T.; Hawreliak, J.; Higginbotham, A.; Park, N.; Ravelo, R.; Stukowski, A.; Suggit, M.; Tang, Y.; Wark, J.; Bringa, E. Molecular dynamics simulations of shock-induced plasticity in tantalum. High Energy Density Physics 2014, 10, 9 – 15. (44) Stukowski, A. Structure identification methods for atomistic simulations of crystalline materials. Model. Simul. Mater. Sci. Eng. 2012, 20, 045021. (45) Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO - the Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012, http://www.ovito.org/ .

22

ACS Paragon Plus Environment

Page 22 of 25

Page 23 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(46) Gao, A.; Mukherjee, S.; Srivastava, I.; Daly, M.; Singh, C. V. Atomistic origins of ductility enhancement in metal oxide coated silicon nanowires for Li-Ion Battery Anodes. Adv. Mater. Interfaces 2017, 4, 1700920. (47) Bringa, E. M.; Monk, J. D.; Caro, A.; Misra, A.; Zepeda-Ruiz, L.; Duchaineau, M.; Abraham, F.; Nastasi, M.; Picraux, S. T.; Wang, Y. Q.; et al., Are nanoporous materials radiation resistant? Nano Lett. 2012, 12, 3351. (48) Zepeda-Ruiz, L. A.; Martinez, E.; Caro, M.; Fu, E. G.; Caro, A. Deformation mechanisms of irradiated metallic nanofoams. Appl. Phys. Letters 2013, 103, 031909. (49) Stukowski, A. In Multiscale Materials Modeling for Nanomechanics; Weinberger, C. R., Tucker, G. J., Eds.; Springer International Publishing: Cham, 2016; pp 317–336. (50) Delannay, F.; Froyen, L.; Deruyttere, A. The wetting of solids by molten metals and its relation to the preparation of metal-matrix composites composites. J. Mater. Sci. 1987, 22, 1–16. (51) Nagl, M. M.; Evans, W. T. The mechanical failure of oxide scales under tensile or compressive load. J. Mater. Sci. 1993, 28, 6247–6260. (52) Osgerby, S.; Dyson, B. F. Effects of oxygen on creep performance: Mechanisms and predictive modelling. Mater. Sci. Tech. 1990, 6, 2–8. (53) Edeleanu, C.; Law, T. J. Fracture of oxide films on strained aluminium. Philos. Mag. 1962, 7, 573–581. (54) Li, M.; Xie, D.-G.; Ma, E.; Li, J.; Zhang, X.; Shan, Z. Effect of hydrogen on the integrity of aluminium-oxide interface at elevated temperatures. Nat. Commun. 2017, 8, 14564. (55) Xie, D.-G.; Wang, Z.-J.; Sun, J.; Li, J.; Ma, E.; Shan, Z. In situ study of the initiation of hydrogen bubbles at the aluminium metal/oxide interface. Nature Mater. 2015, 14, 899. 23

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(56) Sznajder, M.; Geppert, U.; Dudek, M. R. Hydrogen blistering under extreme radiation conditions. npj Materials Degradation 2018, 2, 2397–2106. (57) Caro, M.; Mook, W. M.; Fu, E. G.; Wang, Y. Q.; Sheehan, C.; Martinez, E.; Baldwin, J. K.; Caro, A. Radiation induced effects on mechanical properties of nanoporous gold foams. Appl. Phys. Lett. 2014, 104, –. (58) Ruestes, C. J.; Anders, C.; Bringa, E. M.; Urbassek, H. M. Nanoindentation tests of heavy-ion-irradiated Au foams - molecular dynamics simulation. J. Appl. Phys. 2018, 123, 225903.

24

ACS Paragon Plus Environment

Page 24 of 25

Page 25 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC graphic

25

ACS Paragon Plus Environment