Organisers

• Carlo Antonio Pignedoli (Swiss Federal Laboratories for Materials Testing and Research (EMPA))
• Daniele Passerone (Swiss Federal Laboratories for Materials Testing and Research (EMPA))

Supports

   CECAM

   EMPA

Description

Metal/ceramics interfaces are object of extensive study, due to their broad technological applications in electronics, catalysis, medicine and engineering.
The definition and prediction of the physical-chemical properties of such a system is a complex problem which has its fundaments in "common" challenges in materials science: the ability to join ceramic parts to metallic alloys is conditioned by the knowledge of the microscopic processes that take place at the interface as a function of the experimental conditions.
In order to properly engineer a metal/ceramic interface, bridging between different approaches in materials science is needed to realize a multiscale approach which combines different
modeling techniques, namely atomistic, thermodynamic and kinetic modeling, as well as continuum mechanics/FEM, and stochastic-kinematic modeling.
The goal is to simulate the interaction between structure and composition of e.g. brazed cBN grains on a micro and nano scale on one hand, and the mechanical performance and long term stability under complex loading conditions on the other.
The problem of metal/ceramic interface is complex and characterized by the ability of a molten alloy to wet the ceramic substrate. The wetting depends on the bonding characteristics of the liquid and the ceramic as well as on the magnitude of interactive forces at the interface.
Moreover, the reactions at molten alloy/ceramic substrate interface involve the formation and subsequent growth of many reaction products whose layering and morphology are the result of the interplay of the thermodynamics and kinetics of the whole system.
A complex metal/ceramic interface is characterized by the coexistence of different phases (e.g. in the case of Ti/cBN, at the interface the TiB TiN and TiB2 phases coexists) so that thermodynamic approaches that allow to compute the different phase equilibria are needed.

How far can we go with ab initio simulations in a realistic description of a metal/ceramic interface?

What information can we get from the electronic structure at the interface in terms of nature of the bonding and its strength?

For this class of systems is it possible to set up a series of ab initio calculations aimed at defining the basic ingredients for a phenomenological (e.g. classical atomistic force fields and mesoscopic coarse graining) large scale approach?

Which kind of physical outcome can we expect from a large scale simulation? Which kind of failures?

At which point an atomistic description has to be abandoned in favor of thermodynamic phase field ap-proaches or continuum models?

Answering these questions should be the starting point for a multiscale approach to the simulation of a complex metal/ceramic interface. During the workshop, experts in different fields of materials science will discuss all these points, from ab initio simulations and large scale molecular dynamics simulations based on empirical potentials up to thermodynamic approaches and continuum models.

State of the art

The realistic simulation of a metal ceramic interface is a challenging problem, a pure ab initio approach would imply unaffordable large size systems. On the other end, ab initio tecniques are mandatory for a correct description of the complex pattern of interactions at the interface.
Large scale systems (several thousands of atoms) are needed to accommodate, without artifacts, all the geometric features (e.g. defects,dislocations,..) that characterize the interface. Only empirical potentials allow for a treatment of large size systems through classical MD based methods or through Monte Carlo methods.
To go beyond the simulation of a model with a geometry fixed by experimental data, it is necessary to combine classical molecular dynamics with methods that allow for enhancing statistical sampling.
We here give a brief summary of the status of the art of computational techniques related to the subject.
The survey is divided in three parts.


Ab initio simulation of interfaces

The microscopic characterization of interfaces with the help of ab initio simulations has reached a considerable level of accuracy and complexity.
With the increase of computational power and the development of new algorithms, realistic modeling of extended systems has become possible in physics, chemistry, and biology (see, e.g., (Allen 2001)); in materials science, phenomenological and semiempirical calculations have gradually left the stage to accurate electronic structure models based on density functional theory (DFT).
Metal/ceramics interfaces, for example, due to the broad technological applications in electronics, catalysis, medicine, engineering, have been object of extensive study (see, e.g., (Finnis 1995) (Finnis 1996)).
Quantities directly accessible to experiments, like the work of adhesion and the work of separation could be computed for systems like metal/sapphire (Kruse 1996), metal/ceramics in presence of Titanium (Lee 1995; Yamazaki 1996; Kostlmeier 2000; Dudiy 2001; Christensen 2002; Dudiy 2004), and several authors developed a theoretical framework to include entropic effects and partial pressure of other elements in the model (Zhang 2000; Xiao-Gang 2002; Wang 2003; Wei-Xue 2003; Chatterjee 2006; Feng 2007; Neurock 2007), or of different orientation of the interface constituents (Dudiy 2004).
Ab initio simulations combined with microscopic statistical mechanics models allowed to extract quantitative predictions concerning the growth of metallic films on insulating substrates (Fuks 2000, Fuks 2002); more recently, the enhancement of metallic wetting at CaF2 surfaces was explained with a detailed inspection to the interfacial electronic structure (Barzilai 2008).
Another very active field of research concerns the interfaces between semiconductors and oxides and their structural properties (Bongiorno 2000; Bongiorno 2002; Bongiorno 2003; Giustino 2005; Fischer 2006; Devynck 2007); here the additional problem being the presence of an amorphous phase like silica, resulting in a complex modelization.
All these studies usually start from models of the interface in the slab geometry with the supercell approach: a suitable number of elementary crystal cells is explicitly included. The electronic degrees of freedom are part of the description, and the structure of the interface is optimized with respect to the total energy of the system.
The interface formation energies and work of separation can be extracted from a series of analogous calculations, but they usually refer to solid phases with plasticity suppressed.
Ab initio thermodynamics or analogous approaches (Finnis 1998; Batyrev 2000; Wei-Xue 2003) can connect finite temperature quantities (free energies, works of adhesion) to the results of zero temperature DFT calculations.
Moreover, a DFT calculation allows to obtain a detailed description of the chemical bonding at the interface. For example, very useful tools for this analysis are the local density of electronic states projected onto the interface atoms or physical quantities related to the electron density like the electron localization function.
Dynamical properties like surface reactions and interfacial diffusion can be studied with DFT-based methods Examples include water/alpha alumina interface (Hass 1998) and again Si/SiO2 (Bongiorno 2005; Andreoni 2006; Fischer 2006; Godet 2006).


Large length and time scales: atomistic interface simulations.

Although ab initio simulations have reached a considerable level of complexity and can provide realistic predictions, most of the length and time scales pertinent to the experimentally accessible macroscopic quantities remain out of reach for "canonical" computational techniques like structural optimization and ab initio molecular dynamics.
One possibility is to include a continuum layer of description in the hierarchy of computational models for the interface. On the other hand, if one wants to stick to an atomistic description, different approaches can be used to overcome size and time limitations.
In order to increase the available system size (from nanometers to micrometers) and simulation times (from picoseconds to microseconds) a first logical step is the development of empirical potentials based on high level ab initio simulations.
Just to confine our analysis to the field of interface science a good example is the development of an augmented Tersoff potential for silicon oxynitride compounds, applied to their interfaces with silicon ((Billeter 2006) see also (Ackland) for a general discussion on empirical potentials).
Several thousands of atoms, different coordinations, defects and diffusive behavior can provide a model of the interface of realistic quality. It was recently shown that the contact angle of wetting can be obtained from classical molecular dynamics simulations with reasonable precision (Webb 2005).
When a time scale gap has to be overcome to describe a "rare event" that is fundamental for the physics of the system different methodologies are available to determine activation barriers (Elber 1987); (Henkelman 2000a; Henkelman 2000b; Galvan 2008), to overcome configurational barriers (Hansmann 1997; Sugita 1999; Bussi 2006) or to force a dynamical simulation to quickly sample the most relevant portion of the phase space (Montalenti 2001; Voter 2002); for a review of the different methods see (Christen 2008))


Thermodynamic calculations Phase-Field Models and continuum models

Computer simulations can provide useful information on the stable and metastable equilibrium phases at the metal/ceramic interface: once all possible reaction products are individually analyzed, the overall equilibrium conditions can be evaluated by CALPHAD method (Sundman 1985).
Recently, CALPHAD simulations have been interconnected with ab initio calculations to provide an even more powerful strategy (see, e.g., Kauf-man 2001). An interesting example of the application of CALPHAD in connection with phase field models (see below) for the evaluation of the equilibrium phases at the interface (together with heterogeneous nu-cleation) between immiscible metals can be found in (Toth 2007).
Another possible approach to characterize the the equilibrium phases at the interface relies on Phase Field Models: the sequence of reaction products is divided into separate steps which are individually analyzed and then coupled in a way to satisfy equilibrium and conservation conditions (2008 Moelans).
These models can incorporate physical effects such as interface kinetics and surface tension anisotropy in a natural manner, since the underlying formalism is based on a consistent thermodynamic description of the system (Visintin 1997).
Finally, to obtain a macroscopic description of the mechanical properties of the interface, continuum models have to be employed, cohesive models have been developed for several purposes (see Tvergaard 2003, Schwalbe 2003), mainly for the description of crack propagation and debonding between hard and ductile phases. In (Wegener 1995, 1999) constitutive equations for time dependent plasticity of metals have been derived and, in (Wegener 1996) the interrelationship between phenomenological modeling and physical principles of metals has been pointed out.

The experimental literature on the metal/ceramic interface is extremely wide and goes beyond the aim of this proposal. However we give a reference to the main experimental techniques that are commonly used to characterize metal/ceramic interfaces. We just give a reference to a review on sessile drop for wetting measurements (Eustathopoulos 1999, Eustathopoulos 2005) and we remind that different microscopic techniques are usually employed to characterize the interface composition: e.g. Scanning/Transmission Electron Microscopy, Electron Probe Micro-Analyzer, Focussed Ion Beam, (S)TEM + Energy Dispersive Xray analysis as well as x-ray diffraction and Raman spectroscopy. Moreover, sessile drop techniques allow to follow the reactive wetting in an efficient fashion (Hodaj 2007). The role of anisotropy in the wetting behavior has been recently reviewed by Chatain (Chatain 2008).
Finally, we cite the remarkable efforts to relate high temperature wetting experiments with thermodynamical models, done in the groups of D. Chatain, N. Eustathopoulos and A. Passerone (for an example concerning oxygen adsorption isotherms, see e.g., Ghetta 1996, Arato 2008, Giuranno 2006, Novakovic 2005a, Novakovic 2005b Passerone 2008, Muolo 2008)

Scientific Objectives

The aim of the workshop is to have a transfer of know how from different fields in materials science; starting from a clear definition of an interface problem from the experimental point of view, it will be outlined the necessary connection between different simulation approaches in order to individuate the main problems that can arise in the realization of such a multiscale approach.
In particular, two main results are expected as outcome of the workshop

-definition of the optimal simulation techniques available to investigate metal/ceramic interfaces
-definition of the feasibility of new multiscale approaches for the investigation of mtal/ceramic interfaces
-definition of a strategy of cooperation between experiments and simulation in the engineering of the inter-faces

References

Logo design by Wrinkly Pea Design (www.wrinklypea.com)

(Ackland) www.ph.ed.ac.uk/cmatter/gja/cargese.pdf

(Allen 2001) Allen, F., G. Almasi, et al. (2001). "Blue Gene: A vision for protein science using a petaflop supercomputer." Ibm Systems Journal 40(2): 310-327.

(Andreoni 2006) Andreoni, W., A. Curioni, et al. (2006). "Studying the Effects of Nitrogen and Hafnium Incorporation into the SiO2/Si(100) Interface with Replica-exchange Molecular Dynamics and Density-functional-theory Calculation."

(Arato 2008) Arato E., Eiici E., Costa P., Surface Science 602, 349, (2008)

(Barzilai 2008) Barzilai, S., N. Argaman, et al. (2008). "First-principles modeling of metal layer adsorption on CaF2(111)." Surface Science 602(8): 1517-1524.

(Batyrev 2000) Batyrev, I. G., A. Alavi, et al. (2000). "Equilibrium and adhesion of Nb/sapphire: The effect of oxygen partial pressure." Physical Review B 62(7): 4698-4706.

(Billeter 2006) Billeter, S. R., A. Curioni, et al. (2006). "Ab initio derived augmented Tersoff potential for silicon oxynitride compounds and their interfaces with silicon." Physical Review B 73(15).

(Bongiorno 2000) Bongiorno, A. and A. Pasquarello (2000). "Validity of the bond-energy picture for the energetics at Si-SiO2 interfaces." Physical Review B 62(24): 16326-16329.

(Bongiorno 2002) Bongiorno, A. and A. Pasquarello (2002). "Atomic structure at the Si(001)-SiO2 interface: from the interpretation of Si 2p core-level shifts to a model structure." Materials Science and Engineering B-Solid State Materials for Advanced Technology 96(2): 102-106.

(Bongiorno 2003) Bongiorno, A., A. Pasquarello, et al. (2003). "Transition structure at the Si(100)-SiO2 interface." Physical Review Letters 90(18).

(Bongiorno 2005) Bongiorno, A. and A. Pasquarello (2005). "Atomic-scale modeling of kinetic processes occurring during silicon oxidation." Journal of Physics-Condensed Matter 17(21): S2051-S2063.

(Bussi 2006) G. Bussi, F.L. Gervasio, et al. "Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics" Journal od the American Chemical Society 128, 13435 (2006)

(Chatain 2008) D. Chatain (2008), "Anisotropy of wetting", Annual Rev. of Materials Research, v.38, 45-70.

(Chatterjee 2006) Chatterjee, A. (2006). "Comparative behavior of Pd and Ag deposition phenomenon on clean alpha-Al2O3(001) surface - A first principle study." Molecular Simulation 32(2): 155-162.

(Christen 2008) Christen, M. And W. F. Van Gunsteren (2008). "On searching in, sampling of, and dynamically moving through conformational space of biomolecular systems: A review." Journal of Computational Chemistry 29(2): 157-166.

(Christensen 2002) Christensen, M., S. Dudiy, et al. (2002). "First-principles simulations of metal-ceramic interface adhesion: Co/WC versus Co/TiC." Physical Review B 65(4).

(Devynck 2007) Devynck, F., F. Giustino, et al. (2007). "Structural and electronic properties of an abrupt 4H-SiC(0001)/SiO2 interface model: Classical molecular dynamics simulations and density functional calculations." Physical Review B 76(7).

(Dudiy 2001) Dudiy, S. V. and B. I. Lundqvist (2001). "First-principles density-functional study of metal-carbonitride interface adhesion: Co/TiC(001) and Co/TiN(001)." Physical Review B 6404(4).

(Dudiy 2004) Dudiy, S. V. and B. I. Lundqvist (2004). "Wetting of TiC and TiN by metals." Physical Review B 69(12).

(Elber 1987) Elber, R. and M. Karplus (1987). "A Method for Determining Reaction Paths in Large Molecules - Application to Myoglobin." Chemical Physics Letters 139(5): 375-380.

(Eustathopoulos 1999) Eustathopoulos N., Nicholas M., Drevet B. Wettability at high temperature. Oxford, Pergamon, (1999).

(Eustathopoulos 2005) Eustathopoulos, N., N. Sobczak, A. Passerone, K. Nogi, "Measurement of contact angle and work of adhesion at high temperature", J. Mat. Sci 40, 2271 (2005), 2271-2280.

(Feng 2007) Feng, J. W., W. Q. Zhang, et al. (2007). "Atomistic thermodynamic simulation of Ag/Al2O3 interfaces under O-2 pressure." Journal of Inorganic Materials 22(1): 119-122.

(Finnis 1995) Finnis, M. W., C. Kruse, et al. (1995). "Ab initio calculations of metal/ceramic interfaces: What have we learned, what can we learn?" Nanostructured Materials 6(1-4): 145-155.

(Finnis 1996) Finnis, M. W. (1996). "The theory of metal-ceramic interfaces." Journal of Physics-Condensed Matter 8(32): 5811-5836.

(Finnis 1998) Finnis, M. W. (1998). "Accessing the excess: An atomistic approach to excesses at planar defects and dislocations in ordered compounds." Physica Status Solidi a-Applied Research 166(1): 397-416.

(Fischer 2006) Fischer, D., A. Curioni, et al. (2006). "The structure of the SiO2/Si(100) interface from a restraint-free search using computer simulations." Applied Physics Letters 88(1).

(Fuks 2000) (Fuks 2000) Fuks, D., S. Dorfman, et al. (2000). "Theoretical analysis of the growth mode for thin metallic films on oxide substrates." Physical Review Letters 85(20): 4333-4336.

(Fuks 2002) Fuks, D., S. Dorfman, et al. (2002). "Theory of the growth mode for a thin metallic film on an insulating substrate." Surface Science 499(1): 24-40.

(Galvan 2008) Galvan, I. F. and M. J. Field (2008). "Improving the efficiency of the NEB reaction path finding algorithm." Journal of Computational Chemistry 29(1): 139-143.

(Ghetta 1996) Ghetta, V., J. Fouletier, D. Chatain (1996). "Oxygen adsorption isotherms at the surfaces of liquid Cu and Au-Cu alloys and their interfaces with Al2O3 detected by wetting experiments." Acta Materialia 44(5): 1927-1936.

(Giuranno 2006) Giuranno D., Ricci E., Arato E., Costa P. Acta Mat., 54, 2625, (2006).

(Giustino 2005) Giustino, F., A. uongiorno, et al. (2005). "Atomistic models of the Si(100)-SiO2 interface: structural, electronic and dielectric properties." Journal of Physics-Condensed Matter 17(21): S2065-S2074.

(Godet 2006) Godet, J. and A. Pasquarello (2006). "Proton diffusion mechanism in amorphous SiO2." Physical Review Letters 97(15).

(Hansmann 1997) U.H.E. Hansmann, "Parallel tempering algorithm for conformational studies of biological molecules "Chemical Physics Letters 281, 140 (1997)

(Hass 1998) Hass, K. C., W. F. Schneider, et al. (1998). "The chemistry of water on alumina surfaces: Reaction dynamics from first principles." Science 282(5387): 265-268.

(Henkelman 2000a) Henkelman, G. and H. Jonsson (2000a). "Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points." Journal of Chemical Physics 113(22): 9978-9985

(Henkelman 2000b) Henkelman, G., B. P. Uberuaga, et al. (2000b). "A climbing image nudged elastic band method for finding saddle points and minimum energy paths." Journal of Chemical Physics 113(22): 9901-9904.

(Hodaj 2007) Hodaj, F., O. Dezellus, et al. (2007). "Diffusion-limited reactive wetting: effect of interfacial reaction behind the advancing triple line." Journal of Materials Science 42(19): 8071-8082.

(Kaufman 2001) Kaufman, L. (2001). "Computational thermodynamics and materials design." Calphad-Computer Coupling of Phase Diagrams and Thermochemistry 25(2): 141-161.

(Kostlmeier 2000) Kostlmeier, S. and C. Elsasser (2000). "Density functional study of the 'titanium effect' at metal-ceramic interfaces." Journal of Physics-Condensed Matter 127): 1209-1222.

(Kruse 1996) Kruse, C., M. W. Finnis, et al.(1996). "First-principles study of the atomistic and electronic structure of the niobium-alpha-alumina0001) interface." Philosophical Magazine Letters 736): 377-383.

(Lee 1995) Lee, Y. M. and J. Yu(1995). "Interface Strengths in a Zro2/Ag-Cu-Al-Ti System." Ceramics International 216): 395-398.

(Moelans 2008) N. Moelans, B. Blanpain, P. Wollants, "An introduction to phase-field modeling of microstructure evolution", Computer Coupling of Phase Diagrams and Thermochemistry, pp. 268-294

(Montalenti 2001) Montalenti, F., M. R. Sorensen, et al. (2001). "Closing the gap between experiment and theory: Crystal growth by temperature accelerated dynamics." Physical Review Letters 871212).

(Muolo 2008) Muolo M. L., Valenza F., Passerone A., Passerone D. "Oxygen influence on ceramic wettability by liquid metals: Ag/a-Al2O3 - Experiments and modelling."
Materials Science and Engineering A-Structural Materials Properties Microstructure and Processing, vol. 495 (1-2) pp. 153 - 158. 2008.

(Neurock 2007) Neurock, M., R. Meyer, et al. (2007). "An ab initio analysis of adsorption and diffusion of silver atoms on alumina surfaces." Surface Science|Surface Science 6011): 134-45.

(Novakovic 2005a) R. Novakovic, T. Tanaka, M.L. Muolo, J. Lee, A. Passerone, "Bulk and surface properties of liquid Ag- X (X = Ti, Hf) compound forming alloys", Surface Science, 591 (2005) 56-69.

(Novakovic 2005b) R. Novakovic, E. Ricci, D. Giuranno, A. Passerone, "Surface and transport properties of Ag-Cu liquid alloys", Surface Science, 576 [1-3] (2005) 175-187.

(Passerone 2008) Passerone A., Muolo M.L., Valenza F. , Monteverde F. , Sobczak N., "Wetting and interfacial phenomena in Ni-HfB2 systems." Acta Materialia, doi: 10.1016/j.actamat.2008.09.16

(Sugita 1999) Y. Sugita, Y. Okamoto, "Replica-exchange molecular dynamics method for protein folding " Chemical Physics Letters 314, 141 (1999)

(Sundman 1985) Sundman, B., B. Jansson, et al. 1985). "The Thermo-Calc Databank System." Calphad-Computer Coupling of Phase Diagrams and Thermochemistry 92): 153-190.

(Toth 2007) Toth, G. I. and L. Granasy (2007). "Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: I. Transitions in the one-phase liquid region." Journal of Chemical Physics|Journal of Chemical Physics: 74709/1-11.

(Visintin 1997) Visintin A. (1997). "Two-scale model of phase transition", Physica D 106, pp.66-80.

(Voter 2002) Voter, A. F., F. Montalenti, et al. (2002). "Extending the time scale in atomistic simulation of materials." Annual Review of Materials Research 32: 321-346.

(Wang 2003) Wang, X. G., J. R. Smith, et al. (2003). "Adhesion of copper and alumina from first principles." Journal of the American Ceramic Society 864): 696-700.

(Webb 2005) Webb, E. B., J. J. Hoyt, et al. (2005). "High temperature-wetting: Insights from atomistic simulations." Current Opinion in Solid State & Materials Science 94-5): 174-180.

(Wegener 1995) Wegener, K. (1995). “A constitutive model for time-dependent plasticity of metals.” Proceedings of Plasticity ’95, Dynamic Plasticity and Structural Behaviors, (Hrsg.: S. Tanimura, A.S. Khan) Gordon and Breach Publishers. pp. 725-728.

(Wegener 1996) Wegener, K. (1996). “Metallphysik und Stoffmodellierung.” Technische Mechanik 16: 23-38.

(Wegener 1999) Wegener, K. (1999). “Constitutive modeling for simulation in forming technology” In: Proceedings of the 6th ICTP, Advanced Technology of Plasticity, Vol 1. (ed. M. Geiger), Springer, Berlin.

(Wei-Xue 2003) Wei-Xue, L., C. Stampfl, et al. (2003). "Insights into the function of silver as an oxidation catalyst by ab initio atomistic thermodynamics." Physical Review B Condensed Matter and Materials Physics)|Physical Review B Condensed Matter and Materials Physics) 6816): 165412-1-15.

(Xiao-Gang 2002) Xiao-Gang, W., J. R. Smith, et al. (2002). "Effect of hydrogen on Al/sub 2/O/sub 3//Cu interfacial structure and adhesion." Physical Review B Condensed Matter and Materials Physics)|Physical Review B Condensed Matter and Materials Physics) 667): 073411/1-4.

(Yamazaki 1996) Yamazaki, T. and A. Suzumura (1996). "Monitoring of reactive metal's diffusion and reaction at Ag-Cu-Ti/ZrO2 brazed interface." Materials Transactions Jim 375): 1103-1108.

(Zhang 2000) Zhang, W. and J. R. Smith (2000). "Nonstoichiometric interfaces and Al2O3 adhesion with Al and Ag." Physical Review Letters 8515): 3225-3228.