Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross Slip”

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 73 records · Page 4

Mitigating Transition in the Fermilab Booster Using a Triple Phase Jump

The PIP-II project will significantly enhance neutrino production for DUNE, Fermilab’s flagship long-baseline neutrino oscillation experiment, by doubling the beam power delivered by the accelerator complex. To achieve this, the total charge injected from the new PIP-II linac into the Booster rapid cycling synchrotron (RCS) will increase from $4.5\times 10^{12}$ to $6.5\times 10^{12}$ protons per pulse. Simultaneously, the Booster’s ramp rate will rise from 15 to 20 Hz, while the injection energy will go from 400 MeV to 800 MeV. The Booster accelerates the beam to 8 GeV and crosses transition at $\gamma_t=5.45$. In current operations, no dedicated $\gamma_t$ jump system is employed; instead, longitudinal emittance growth is controlled via an active quadrupole-mode feedback system. Downstream, the Recycler Ring uses slip-stacking to accumulate beam and increase bunch intensity. However, this process imposes a constraint on the longitudinal emittance at Booster extraction, limiting it to 0.1 eV-s (95\%) to avoid excessive particle loss. Because collective effects scale with intensity, additional measures to mitigate transition crossing may be required to stay below this limit. One potential approach is the so-called triple phase-jump technique. While the method has known limitations, it offers some advantages: it can be implemented using the existing digital low-level RF (LLRF) system, requires no additional magnets or pulsed power supplies, and remains compatible with quadrupole-mode feedback.

Ostiguy, J.-F. [Fermilab] (ORCID:0000000290883681)↗

Mitigating Transition in the Fermilab Booster Using a Triple Phase Jump

The PIP-II project will significantly enhance neutrino production for DUNE, Fermilab s flagship long-baseline neutrino oscillation experiment, by doubling the beam power delivered by the accelerator complex. To achieve this, the total charge injected from the new PIP-II linac into the Booster rapid cycling synchrotron (RCS) will increase from $4.5\times 10^{12}$ to $6.5\times 10^{12}$ protons per pulse. Simultaneously, the Booster s ramp rate will rise from 15 to 20 Hz, while the injection energy will go from 400 MeV to 800 MeV. The Booster accelerates the beam to 8 GeV and crosses transition at $\gamma_t=5.45$. In current operations, no dedicated $\gamma_t$ jump system is employed; instead, longitudinal emittance growth is controlled via an active quadrupole-mode feedback system. Downstream, the Recycler Ring uses slip-stacking to accumulate beam and increase bunch intensity. However, this process imposes a constraint on the longitudinal emittance at Booster extraction, limiting it to 0.1 eV-s (95\%) to avoid excessive particle loss. Because collective effects scale with intensity, additional measures to mitigate transition crossing may be required to stay below this limit. One potential approach is the so-called triple phase-jump technique. While the method has known limitations, it offers some advantages: it can be implemented using the existing digital low-level RF (LLRF) system, requires no additional magnets or pulsed power supplies, and remains compatible with quadrupole-mode feedback.

Ostiguy, Jean-Francois [Fermilab] (ORCID:000000029↗

Androgynous { 10 1 ¯ 2 } twin in zinc

Under ambient conditions, Zn is a hexagonal metal with a large c / a ratio of 1.856. Plastic deformation is predominantly accommodated by basal 〈 a 〉 slip and compression twins on the { 10 1 ¯ 2 } planes. Increasing hydrostatic pressure drastically reduces the c / a ratio of Zn and, when a critical threshold of c / a = 3 at about 10 GPa is crossed, the { 10 1 ¯ 2 } twins are predicted to change from compressive to tensile in nature. What happens at the transition point, when c / a = 3 , remains unknown. Here, we strain-cycle a textured polycrystalline sample of pure Zn at uniform hydrostatic pressures ranging between 2 and 17 GPa, over which the c / a ratio crosses the c / a = 3 compressive-tensile transition for { 10 1 ¯ 2 } twins. During deformation, the state of the sample is monitored through x-ray diffraction to extract texture and internal strain evolution. By comparing the experimental results with the predictions of an elastoviscoplastic polycrystal simulation, we confirm the androgynous nature of { 10 1 ¯ 2 } twin response at low and high pressures. When c / a = 3 , polycrystalline Zn does not display any evidence of twinning and its plastic behavior is controlled by mostly basal and pyramidal 〈 c + a 〉 slip activity, with a very small contribution of prismatic 〈 a 〉 slip. Evidence for the activity of other { 10 1 ¯ n } twinning modes, which have been suggested for Zn under high pressure, are not observed. Published by the American Physical Society 2024

Merkel, Sébastien (ORCID:000000032767581X)↗

Theoretical study of the CO 2 –O 2 van der Waals complex: potential energy surface and applications

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO 2 –O 2 ) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b), and extrapolation to the complete basis set (CBS) limit using the cc-pVTZ-F12/cc-pVQZ-F12 bases and the l –3 formula. An analytic representation of the 4D-PES was fitted using the method of interpolating moving least squares (IMLS). These calculations predict that the most stable configuration of CO 2 –O 2 complex corresponds to a planar slipped-parallel structure with a binding energy of V ~ –243 cm –1 . Another isomer is found on the PES, corresponding to a non-planar cross-shaped structure, with V ~ –218 cm –1 . The transition structure connecting the two minima is found at V ~ –211 cm –1 . We also performed comparisons with some CO 2 –X van der Waals complexes. Moreover, we provide a SAPT analysis of this molecular system. Then, we discuss the complexation induced shifts of CO 2 and O 2 . Afterwards, this new 4D-PES is employed to compute the second virial coefficient including temperature dependence. A comparison between quantities obtained in our calculations and those from experiments found close agreement attesting to the high quality of the PES and to the importance of considering a full description of the anisotropic potential for the derivation of thermophysical properties of CO 2 –O 2 mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

En echelon faults reactivated by wastewater disposal near Musreau Lake, Alberta

We use machine-learning and cross-correlation techniques to enhance earthquake detectability by two magnitude units for the earthquake sequence near Musreau Lake, Alberta, which is induced by wastewater disposal. This deep catalogue reveals a series of en echelon ~N–S oriented strike-slip faults that are favourably oriented for reactivation. These faults require only ~0.6 MPa overpressure for triggering to occur. Earthquake activity occurs in bursts, or episodes; episodes restricted to the largest fault tend to have earthquakes starting near the southern end (distant from injectors) and progressing northwards (towards the injectors). While most events are concentrated along these ~N–S oriented faults, we also delineate smaller faults. Together, these findings suggest pore pressure as the triggering mechanism, where a time-dependent increase in pore pressure likely caused these faults to progressively reawaken. Analysis of the ‘next record-breaking event’, a statistical model that forecasts the sequencing of earthquake magnitudes, suggests that the next largest event would be M L ~4.3. The seismically illuminated length of the largest fault indicates potential magnitudes as large as M w 5.3.

58 GEOSCIENCES↗

Dynamic Deformation in Nuclear Graphite and Underlying Mechanisms

Time-dependent deformation in nuclear graphite is influenced by the creation and migration of radiation-induced defects in the reactor environment. This study investigates the role of pre-existing defects such as point defect clusters and Mrozowski cracks in nuclear graphite IG-110. Separate specimens were irradiated with a 2.8 MeV Au 2+ beam with a fluence of 4.38 × 10 14 cm −2 and an 8 MeV C 2+ beam with a fluence of 1.24 × 10 16 cm −2 . Microscopic specimens were either mechanically loaded inside a transmission electron microscope (TEM) or subjected to ex situ indentation-based creep loading. In situ TEM tests showed significant plasticity in regions highly localized around the Mrozowski cracks, resembling slip or ripplocation bands. Slip bands were also seen near regions without pre-existing defects but at very high stresses. Ex situ self-ion irradiation embrittled the specimens and decreased the creep displacement and rate, while heavy ion irradiation resulted in the opposite behavior. We hypothesize that the large-sized gold ions (compared to the carbon atoms) induced interplanar swelling as well as cross-plane channels for increased defect mobility. These findings illustrate the role of pre-existing defects in the dynamic relaxation of stresses during irradiation and the need for more studies into the radiation environment’s impact on the mechanical response of nuclear graphite.

creep deformation↗

Critical Fluorescence of a Transmon at the Schmid Transition

Here, we investigate inelastic microwave photon scattering by a transmon qubit embedded in a high-impedance circuit. The transmon undergoes a charge-localization (Schmid) transition upon the impedance reaching the critical value. Due to the unique transmon level structure, the fluorescence spectrum carries a signature of the transition point. At higher circuit impedance, quasielastic photon scattering may account for the main part of the inelastic scattering cross-section; we find its dependence on the qubit and circuit parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗

Cosmology at high redshift — a probe of fundamental physics

An observational program focused on the high redshift (2<z<6) Universe has the opportunity to dramatically improve over upcoming LSS and CMB surveys on measurements of both the standard cosmological model and its extensions. Using a Fisher matrix formalism that builds upon recent advances in Lagrangian perturbation theory, we forecast constraints for future spectroscopic and 21-cm surveys on the standard cosmological model, curvature, neutrino mass, relativistic species, primordial features, primordial non-Gaussianity, dynamical dark energy, and gravitational slip. We compare these constraints with those achievable by current or near-future surveys such as DESI and Euclid, all under the same forecasting formalism, and compare our formalism with traditional linear methods. Here, our Python code FishLSS — used to calculate the Fisher information of the full shape power spectrum, CMB lensing, the cross-correlation of CMB lensing with galaxies, and combinations thereof — is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Structural Evolution of the Hogback Monocline and Its Tectonic Significance in the San Juan Basin

The San Juan Basin is recognized as a Laramide foreland basin. It is located within the Colorado Plateau, a broad tectonic province characterized by a thick sedimentary sequence that was segmented into smaller sub basins during the Late Cretaceous to Paleogene Laramide orogeny. The Hogback Monocline lies along the northwestern margin of the San Juan Basin and is considered a Laramide-age structure formed in response to compressional stress. In this study, we interpret surface and subsurface datasets to construct a structural geological model and evaluate its tectonic significance. Through seismic data, we identify key fault and fold geometries at depth. The seismic dataset used in this study was reprocessed in depth and constrained with well log velocity data to enhance seismic imaging quality. Additionally, we performed well log correlations to identify formation tops and assess variations in basin infill and thickness geometry. A series of structural cross-sections, constructed using seismic data and a high density of boreholes, are presented to evaluate geometric variations along the structure and its evolution during basin development. Furthermore, kinematic restoration and forward modeling analyses were conducted to validate our structural interpretation. This work suggests that the Hogback Monocline formed through fault-propagation folding and flexural slip affecting the pre-Laramide sedimentary sequence under compressional stresses associated with the Laramide orogeny. This structure is interpreted as a high-angle reverse fault that influenced the geometry of the late basin infill. Additionally, monocline bending along the structure may have been controlled by fault relay systems and, in some cases, influenced by strike-slip faulting.

Reyes, Martin [New Mexico Bureau o fGeology and Mi↗

Chemo-Mechanics of α-V 2 O 5 During Lithiation and Implications for Rechargeable Battery Cathodes

Chemo-mechanical degradation of layered oxide electrodes is strongly influenced by crystallographic anisotropy, local stress evolution, and ion insertion, yet the intrinsic mechanical response of layered materials remains incompletely understood. Indeed, most prior studies have focused on polycrystalline materials but single crystals enable direct observation of coupling between anisotropic ion diffusion and mechanical response. This study aims to determine how crystallographic anisotropy and lithiation affect deformation, fracture, and mechanical properties in single-crystal V 2 O 5 , and compares this behavior with polycrystalline counterparts. Polycrystalline V 2 O 5 thin films and single-crystal α-V 2 O 5 were studied using nanoindentation, scanning electron microscopy, focused ion beam cross-sectioning, and Raman spectroscopy. Single crystals were tested in pristine and chemically lithiated states, including experiments in which crystals were first plastically deformed via nanoindentation and subsequently lithiated. Polycrystalline films exhibited significantly higher hardness and elastic modulus than single crystals. Single crystals indented normal to the exposed (001) basal plane exhibited pronounced anisotropic deformation, including directional slip, crystallographically-guided cracking, anisotropic crack propagation, interlayer separation, and shear localization. Lithiation caused substantial softening, reduced hardness and modulus, and suppressed displacement bursts during nanoindentation, while previously indented regions showed crack formation and growth upon lithiation. Mechanical behavior of α-V 2 O 5 is strongly governed by crystallographic anisotropy and further altered by lithiation, with pre-existing deformation serving as a strong driver of fracture during ion insertion. These findings illuminate the coupling among ion insertion, deformation, and fracture in layered oxides and provide a basis for understanding and mitigating mechanical failure in electrochemical energy-storage materials.

Anisotropy↗

Full-Cycle Simulations of the Fermilab Booster

The Proton Improvement Plan phase II (PIP-II) project currently under construction at FNAL will replace the existing 400 MeV normal conducting linac with a new 800 MeV superconducting linac. The beam power in the downstream rapid-cycling Booster synchrotron will be doubled by raising the machine cycle frequency from 15 to 20 Hz and by increasing the injected beam intensity by a factor 1.5. This has to be accomplished without raising uncontrolled losses beyond the administrative limit of 500 W. In addition, slip-stacking efficiency in the Recycler, the next machine in the accelerator chain, sets an upper limit on the longitudinal emittance of the beam delivered by the Booster. As part of an effort to better understand potential losses and emittance blow-up in the Booster, we have been conducting full cycle 6D simulations using the code PyORBIT. The simulations include space charge, wall impedance effects and transition crossing. In this paper, we discuss our experience with the code and present representative results for possible operational scenarios.

43 PARTICLE ACCELERATORS↗

Full Cycle Simulations of The Fermilab Booster

The Proton Improvement Plan phase II (PIP-II) project currently under construction at FNAL will replace the existing 400 MeV normal conducting linac with a new 800 MeV superconducting linac. The beam power in the downstream rapid-cycling Booster synchrotron will be doubled by raising the machine cycle frequency from 15 to 20 Hz and by increasing the injected beam intensity by a factor 1.5. This has to be accomplished without raising uncontrolled losses beyond the administrative limit of 500 W. In addition, slip-stacking efficiency in the Recycler, the next machine in the accelerator chain, sets an upper limit on the longitudinal emittance of the beam delivered by the Booster. As part of an effort to better understand potential losses and emittance blow-up in the Booster, we have been conducting full cycle 6D simulations using the code PyORBIT. The simulations include space charge, wall impedance effects and transition crossing. In this paper, we discuss our experience with the code and present representative results for possible operational scenarios.

43 PARTICLE ACCELERATORS↗

Operating in a deep underground facility improves the locking of gradiometric fluxonium qubits at the sweet spots

We demonstrate flux-bias locking and operation of a gradiometric fluxonium artificial atom using two symmetric granular aluminum (grAl) loops to implement the superinductor. The gradiometric fluxonium shows two orders of magnitude suppression of sensitivity to homogeneous magnetic fields, which can be an asset for hybrid quantum systems requiring strong magnetic field biasing. By cooling down the device in an external magnetic field while crossing the metal-to-superconductor transition, the gradiometric fluxonium can be locked either at 0 or Φ 0 /2 effective flux bias, corresponding to an even or odd number of trapped fluxons, respectively. At mK temperatures, the fluxon parity prepared during initialization survives to magnetic field bias exceeding 100 Φ 0 ⁠. However, even for states biased in the vicinity of 1 Φ 0 ⁠, we observe unexpectedly short fluxon lifetimes of a few hours, which cannot be explained by thermal or quantum phase slips. When operating in a deep-underground cryostat of the Gran Sasso laboratory, the fluxon lifetimes increase to days, indicating that ionizing events activate phase slips in the grAl superinductor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗

Strain-rate dependent deformation mechanisms in single-layered Cu, Mo, and multilayer Cu/Mo thin films

Here, strain-rate sensitivity and rate-dependent hardness, over a range of 10 -2 to 10 2 s -1 , of sputter-deposited single-layered Cu, Mo, and 5 nm Cu/ 5 nm Mo, and 100 nm Cu/ 100 nm Mo multilayer films with a total film thickness of 5 μm were measured using nanoindentation. The plastic zone underneath the nanoindents was characterized via cross-sectional transmission electron microscopy (XTEM). The multilayer films exhibited enhanced hardness but slightly reduced strain-rate sensitivity with decreasing layer thickness from 100 nm to 5 nm. Only the 5 nm Cu/ 5 nm Mo multilayer film exhibited shear bands underneath the nanoindents, and the size of the shear bands increased with increasing strain rate. In contrast, the 100 nm Cu/ 100 nm Mo multilayer film exhibited material pile-up around the indents and significant nanotwinning within Cu grains. The effect of strain rate and layer thickness on the hardness and strain rate sensitivity of the multilayer thin films is interpreted using a modified confined layer slip (CLS) model. The reduced rate sensitivity at 5 nm as compared to 100 nm correlates with abundant growth nanotwins in the Cu grains in 100 nm and formation of shear bands in 5 nm multilayers. In single layer films, a substructure with a high density of dislocations was observed consistent with the plastic strain gradient in the indent plastic zone. No evidence of deformation twins was noted in any of the samples.

36 MATERIALS SCIENCE↗