Quantum chemical study of the effect of precursor stereochemistry on dissociative chemisorption and surface redox reactions during the atomic layer deposition of the transition metal copper

Using quantum chemical calculations, we investigate surface reactions of copper precursors and diethylzinc as the reducing agent for effective Atomic Layer Deposition (ALD) of Cu. The adsorption of various commonly used Cu(II) precursors is explored. The precursors vary in the electronegativity and conjugation of the ligands and flexibility of the whole molecule. Our study shows that the overall stereochemistry of the precursor governs the adsorption onto its surface. Formation of different Cu(II)/Cu(I)/Cu(0) intermediate complexes from the respective Cu(II) compounds on the surface is also explored. The surface model is a (111) facet of a Cu55 cluster. Cu(I) compounds are found to cover the surface after the precursor pulse, irrespective of the precursor chosen. We provide new information about the surface chemistry of Cu(II) versus Cu(I) compounds. A pair of CuEt intermediates or the dimer Cu2Et2 reacts in order to deposit a new Cu atom and release gaseous butane. In this reaction, two electrons from the Et anions are donated to copper for reduction to metallic form. This indicates that a ligand exchange between the Cu and Zn is important for the success of this transmetalation reaction. The effect of the ligands in the precursor on the electron density before and after adsorption onto the surface has also been computed through population analysis. In the Cu(I) intermediate, charge is delocalized between the Cu precursor and the bare copper surface, indicating metallic bonding as the precursor densifies to the surface.

Computational study of the growth of copper thin films by atomic layer deposition

Copper is the main interconnect material in microelectronic devices, and a 2 nm-thick continuous Cu film seed layer needs to be deposited to produce microelectronic devices with the smallest features and more functionality. Atomic layer deposition (ALD) is the most suitable method to deposit such thin films. However, the reaction mechanism and the surface chemistry of copper ALD remain unclear, which is deterring the development of better precursors and design of new ALD processes. In this thesis, we study the surface chemistries during ALD of copper by means of density functional theory (DFT). To understand the effect of temperature and pressure on the composition of copper with substrates, we used ab initio atomistic thermodynamics to obtain phase diagram of the Cu(111)/SiO2(0001) interface. We found that the interfacial oxide Cu2O phases prefer high oxygen pressure and low temperature while the silicide phases are stable at low oxygen pressure and high temperature for Cu/SiO2 interface, which is in good agreement with experimental observations. Understanding the precursor adsorption on surfaces is important for understanding the surface chemistry and reaction mechanism of the Cu ALD process. Focusing on two common Cu ALD precursors, Cu(dmap)2 and Cu(acac)2, we studied the precursor adsorption on Cu surfaces by means of van der Waals (vdW) inclusive DFT methods. We found that the adsorption energies and adsorption geometries are dependent on the adsorption sites and on the method used to include vdW in the DFT calculation. Both precursor molecules are partially decomposed and the Cu cations are partially reduced in their chemisorbed structure. It is found that clean cleavage of the ligand−metal bond is one of the requirements for selecting precursors for ALD of metals. 2 Bonding between surface and an atom in the ligand which is not coordinated with the Cu may result in impurities in the thin film. To have insight into the reaction mechanism of a full ALD cycle of Cu ALD, we proposed reaction pathways based on activation energies and reaction energies for a range of surface reactions between Cu(dmap)2 and Et2Zn. The butane formation and desorption steps are found to be extremely exothermic, explaining the ALD reaction scheme of original experimental work. Endothermic ligand diffusion and re-ordering steps may result in residual dmap ligands blocking surface sites at the end of the Et2Zn pulse, and in residual Zn being reduced and incorporated as an impurity. This may lead to very slow growth rate, as was the case in the experimental work. By investigating the reduction of CuO to metallic Cu, we elucidated the role of the reducing agent in indirect ALD of Cu. We found that CuO bulk is protected from reduction during vacuum annealing by the CuO surface and that H2 is required in order to reduce that surface, which shows that the strength of reducing agent is important to obtain fully reduced metal thin films during indirect ALD processes. Overall, in this thesis, we studied the surface chemistries and reaction mechanisms of Cu ALD processes and the nucleation of Cu to form a thin film.

A prochelator activated by beta-secretase inhibits Abeta aggregation and suppresses copper-induced reactive oxygen species formation.

The intersection of the amyloid cascade hypothesis and the implication of metal ions in Alzheimer's disease progression has sparked an interest in using metal-binding compounds as potential therapeutic agents. In the present work, we describe a prochelator SWH that is enzymatically activated by beta-secretase to produce a high affinity copper chelator CP. Because beta-secretase is responsible for the amyloidogenic processing of the amyloid precursor protein, this prochelator strategy imparts disease specificity toward copper chelation not possible with general metal chelators. Furthermore, once activated, CP efficiently sequesters copper from amyloid-beta, prevents and disassembles copper-induced amyloid-beta aggregation, and diminishes copper-promoted reactive oxygen species formation.

Coordination of platinum therapeutic agents to met-rich motifs of human copper transport protein1.

Platinum therapeutic agents are widely used in the treatment of several forms of cancer. Various mechanisms for the transport of the drugs have been proposed including passive diffusion across the cellular membrane and active transport via proteins. The copper transport protein Ctr1 is responsible for high affinity copper uptake but has also been implicated in the transport of cisplatin into cells. Human hCtr1 contains two methionine-rich Mets motifs on its extracellular N-terminus that are potential platinum-binding sites: the first one encompasses residues 7-14 with amino acid sequence Met-Gly-Met-Ser-Tyr-Met-Asp-Ser and the second one spans residues 39-46 with sequence Met-Met-Met-Met-Pro-Met-Thr-Phe. In these studies, we use liquid chromatography and mass spectrometry to compare the binding interactions between cisplatin, carboplatin and oxaliplatin with synthetic peptides corresponding to hCtr1 Mets motifs. The interactions of cisplatin and carboplatin with Met-rich motifs that contain three or more methionines result in removal of the carrier ligands of both platinum complexes. In contrast, oxaliplatin retains its cyclohexyldiamine ligand upon platinum coordination to the peptide.

Development of Stimulus-Responsive Ligands for the Modulation of Copper and Iron Coordination

The ability to manipulate the coordination chemistry of metal ions has significant ramifications for the study and treatment of metal-related health concerns, including iron overload, UV skin damage, and microbial infection among many other conditions. To address this concern, chelating agents that change their metal binding characteristics in response to external stimuli have been synthesized and characterized by several spectroscopic and chromatographic analytical methods. The primary stimuli of interest for this work are light and hydrogen peroxide.

Herein we report the previously unrecognized photochemistry of aroylhydrazone metal chelator ((E)-N′-[1-(2-hydroxyphenyl)ethyliden]isonicotinoylhydrazide) (HAPI) and its relation to HAPI metal binding properties. Based on promising initial results, a series of HAPI analogues was prepared to probe the structure-function relationships of aroylhydrazone photochemistry. These efforts elucidate the tunable nature of several aroylhydrazone photoswitching properties.

Ongoing efforts in this laboratory seek to develop compounds called prochelators that exhibit a switch from low to high metal binding affinity upon activation by a stimulus of interest. In this context, we present new strategies to install multiple desired functions into a single structure. The prochelator 2-((E)-1-(2-isonicotinoylhydrazono)ethyl)phenyl (E)-3-(2,4-dihydroxyphenyl)acrylate (PC-HAPI) is masked with a photolabile trans-cinnamic acid protecting group that releases umbelliferone, a UV-absorbing, antioxidant coumarin along with a chelating agent upon UV irradiation. In addition to the antioxidant effects of the coumarin, the released chelator (HAPI) inhibits metal-catalyzed production of damaging reactive oxygen species. Finally a peroxide-sensitive prochelator quinolin-8-yl (Z)-3-(4-hydroxy-2-((4-(4,4,5,5-tetramethyl-1,3,2-dioxaborolan-2-yl)benzyl)oxy)phenyl)acrylate (BCQ) has been prepared using a novel synthetic route for functionalized cis-cinnamate esters. BCQ uses a novel masking strategy to trigger a 90-fold increase in fluorescence emission, along with the release of a desired chelator, in the presence of hydrogen peroxide.

Copper signaling axis as a target for prostate cancer therapeutics.

Previously published reports indicate that serum copper levels are elevated in patients with prostate cancer and that increased copper uptake can be used as a means to image prostate tumors. It is unclear, however, to what extent copper is required for prostate cancer cell function as we observed only modest effects of chelation strategies on the growth of these cells in vitro. With the goal of exploiting prostate cancer cell proclivity for copper uptake, we developed a "conditional lethal" screen to identify compounds whose cytotoxic actions were manifested in a copper-dependent manner. Emerging from this screen was a series of dithiocarbamates, which, when complexed with copper, induced reactive oxygen species-dependent apoptosis of malignant, but not normal, prostate cells. One of the dithiocarbamates identified, disulfiram (DSF), is an FDA-approved drug that has previously yielded disappointing results in clinical trials in patients with recurrent prostate cancer. Similarly, in our studies, DSF alone had a minimal effect on the growth of prostate cancer tumors when propagated as xenografts. However, when DSF was coadministered with copper, a very dramatic inhibition of tumor growth in models of hormone-sensitive and of castrate-resistant disease was observed. Furthermore, we determined that prostate cancer cells express high levels of CTR1, the primary copper transporter, and additional chaperones that are required to maintain intracellular copper homeostasis. The expression levels of most of these proteins are increased further upon treatment of androgen receptor (AR)-positive prostate cancer cell lines with androgens. Not surprisingly, robust CTR1-dependent uptake of copper into prostate cancer cells was observed, an activity that was accentuated by activation of AR. Given these data linking AR to intracellular copper uptake, we believe that dithiocarbamate/copper complexes are likely to be effective for the treatment of patients with prostate cancer whose disease is resistant to classical androgen ablation therapies.

Three-dimensional free surface modelling in an unstructured mesh environment for metal processing applications

In the casting of metals, tundish flow, welding, converters, and other metal processing applications, the behaviour of the fluid surface is important. In aluminium alloys, for example, oxides formed on the surface may be drawn into the body of the melt where they act as faults in the solidified product affecting cast quality. For this reason, accurate description of wave behaviour, air entrapment, and other effects need to be modelled, in the presence of heat transfer and possibly phase change. The authors have developed a single-phase algorithm for modelling this problem. The Scalar Equation Algorithm (SEA) (see Refs. 1 and 2), enables the transport of the property discontinuity representing the free surface through a fixed grid. An extension of this method to unstructured mesh codes is presented here, together with validation. The new method employs a TVD flux limiter in conjunction with a ray-tracing algorithm, to ensure a sharp bound interface. Applications of the method are in the filling and emptying of mould cavities, with heat transfer and phase change.

Multiphysics modeling and its application to the casting process

A modeling strategy is presented to solve the governing equations of fluid flow, temperature (with solidification), and stress in an integrated manner. These equations are discretized using finite volume methods on unstructured grids, which provide the capability to represent complex domains. Both the cell-centered and vertex-based forms of the finite volume discretization procedure are explained, and the overall integrated solution procedure using these techniques with suitable solvers is detailed. Two industrial processes, based on the casting of metals, are used to demonstrate the capabilities of the resultant modeling framework. This manufacturing process requires a high degree of coupling between the governing physical equations to accurately predict potential defects. Comparisons between model predictions and experimental observations are given.

Modeling of ingot development during the start-up phase of direct chill casting

Direct chill (DC) casting is a core primary process in the production of aluminum ingots. However, its operational optimization is still under investigation with regard to a number of features, one of which is the issue of curvature at the base of the ingot. Analysis of these features requires a computational model of the process that accounts for the fluid flow, heat transfer, solidification phase change, and thermomechanical analysis. This article describes an integrated approach to the modeling of all the preceding phenomena and their interactions.

Modelling the performance of lead-Free solder interconnects for copper bumped flip-chip devices

Traditionally, before flip chips can be assembled the dies have to be attached with solder bumps. This process involves the deposition of metal layers on the Al pads on the dies and this is called the under bump metallurgy (UBM). In an alternative process, however, Copper (Cu) columns can be used to replace solder bumps and the UBM process may be omitted altogether. After the bumping process, the bumped dies can be assembled on to the printed circuit board (PCB) by using either solder or conductive adhesives. In this work, the reliability issues of flip chips with Cu column bumped dies have been studied. The flip chip lifetime associated with the solder fatigue failure has been modeled for a range of geometric parameters. The relative importance of these parameters is given and solder volume has been identified as the most important design parameter for long-term reliability. Another important problem that has been studied in this work is the dissolution of protection metals on the pad and Cu column in the reflow process. For small solder joints the amount of Cu which dissolves into the molten solder after the protection layers have worn out may significantly affect solder joint properties.

Macro-segregation of multi-component alloys in shape casting simulation

The objective of this work is to present a new scheme for temperature-solute coupling in a solidification model, where the temperature and concentration fields simultaneously satisfy the macro-scale transport equations and, in the mushy region, meet the constraints imposed by the thermodynamics and the local scale processes. A step-by-step explanation of the macrosegregation algorithm, implemented in the finite volume unstructured mesh multi-physics modelling code PHYSICA, is initially presented and then the proposed scheme is validated against experimental results obtained by Krane for binary and a ternary alloys.

A comprehensive copper stockpile leach model: background and model formulation

The first phase in the sign, development and implementation of a comprehensive computational model of a copper stockpile leach process is presented. The model accounts for transport phenomena through the stockpile, reaction kinetics for the important mineral species, oxgen and bacterial effects on the leach reactions, plus heat, energy and acid balances for the overall leach process. The paper describes the formulation of the leach process model and its implementation in PHYSICA+, a computational fluid dynamic (CFD) software environment. The model draws on a number of phenomena to represent the competing physical and chemical features active in the process model. The phenomena are essentially represented by a three-phased (solid liquid gas) multi-component transport system; novel algorithms and procedures are required to solve the model equations, including a methodology for dealing with multiple chemical species with different reaction rates in ore represented by multiple particle size fractions. Some initial validation results and application simulations are shown to illustrate the potential of the model.

A comprehensive copper stockpile leach model: background and parameter sensitivity

The design and development of a comprehensive computational model of a copper stockpile leach process is summarized. The computational fluid dynamic software framework PHYSICA+ and various phenomena were used to model transport phenomena, mineral reaction kinetics, bacterial effects, and heat, energy and acid balances for the overall leach process. In this paper, the performance of the model is investigated, in particular its sensitvity to particle size and ore permeability. A combination of literature and laboratory sources was used to parameterize the model. The simulation results from the leach model are compared with closely controlled column pilot scale tests. The main performance characteristics (e.g. copper recovery rate) predicted by the model compare reasonably well with the experimental data and clearly reflect the qualitiative behavior of the process in many respects. The model is used to provide a measure of the sensitivity of ore permeability on leach behavior, and simulation results are examined for several different particle size distributions.