Influence of Longer Range Transfer of Vapor Interface Modified Caging Constraints on the Spatially Heterogeneous Dynamics of Glass-Forming Liquids

Thickness-Dependent Segmental Dynamics in Supported Thin Films: Insights from a Dynamically Correlated Network Model

A large body of experimental studies shows that the local dynamics in supercooled liquids are significantly altered by spatial nanoconfinement. In a previous study, we proposed a concept of a dynamically correlated network (DCN) model, which assumes that segments in a supercooled liquid undergo cooperative rearrangements within a network-like cluster. We further demonstrated that a model modified for freestanding thin films can predict for the glass transition dynamics in atactic polystyrene (PS) films consistent with experimental results. In this study, we adapted the model to apply it to supported thin films by introducing a layer of virtual vacant segments at the free surface and virtual anchoring segments at the liquid/substrate interface. The latter segments, carrying a finite number of virtual segments, reduce mobility at the interface. We evaluated the cooperative cluster size and distribution with respect to temperature and film thickness, along with the average relaxation time and glass transition temperature Tg for supported thin films of PS. The model predicted that the thickness dependence of Tg for PS becomes stronger with increasing time scale, and this result agreed well with experimental data across different timescales from pseudothermodynamic and dynamic measurements. The results provide insights into the origin of the dynamical decoupling between pseudothermodynamic and dynamic glass transition behaviors.

Surface Glass Transition Temperature Region of Diarylethene Films Determined by Nano-Marangoni Effect

For the last two decades, research has addressed whether the glass transition temperature and the molecular motions on the surface of organic films are significantly different from those inside the bulk glasses. It is reported that the surface of the photochromic diarylethene film prepared by vacuum deposition has fluidity and the vacuum deposition of small amount of rubrene molecules induces surface tension fluctuations, generating dents due to the Marangoni flow in nanoscale. The depth of the dents increases in proportion to these radii for the colorless diarylethene film with a bulk glass transition temperature (Tg) close to room temperature. On the other hand, in the colored diarylethene obtained by UV irradiation to the colorless film, the depth becomes constant at a certain level. The Tg distribution in the depth direction is clarified based on an analysis of the dent depth. By approximating the obtained Tg depth distribution with an exponential function, the outermost surface Tg is about 100 K lower than the bulk Tg in the case of photoisomerized diarylethene.

Effect of the nature of the solid substrate on spatially heterogeneous activated dynamics in glass forming supported films

We extend the force-level elastically collective nonlinear Langevin equation theory to treat the spatial gradients of the alpha relaxation time and glass transition temperature, and the corresponding film-averaged quantities, to the geometrically asymmetric case of finite thickness supported films with variable fluid-substrate coupling. The latter typically nonuniversally slows down motion near the solid-liquid interface as modeled via modification of the surface dynamic free energy caging constraints that are spatially transferred into the film and which compete with the accelerated relaxation gradient induced by the vapor interface. Quantitative applications to the foundational hard sphere fluid and a polymer melt are presented. The strength of the effective fluid-substrate coupling has very large consequences for the dynamical gradients and film-averaged quantities in a film thickness and thermodynamic state dependent manner. The interference of the dynamical gradients of opposite nature emanating from the vapor and solid interfaces is determined, including the conditions for the disappearance of a bulk-like region in the film center. The relative importance of surface-induced modification of local caging vs the generic truncation of the long range collective elastic component of the activation barrier is studied. The conditions for the accuracy and failure of a simple superposition approximation for dynamical gradients in thin films are also determined. The emergence of near substrate dead layers, large gradient effects on film-averaged response functions, and a weak non-monotonic evolution of dynamic gradients in thick and cold films are briefly discussed. The connection of our theoretical results to simulations and experiments is briefly discussed, as is the extension to treat more complex glass-forming systems under nanoconfinement.

Interplay between dynamic heterogeneity and interfacial gradients in a model polymer film

Glass-forming liquids exhibit long-lived, spatially correlated dynamical heterogeneity, in which some nm-scale regions in the fluid relax more slowly than others. In the nanoscale vicinity of an interface, glass-formers also exhibit the emergence of massive interfacial gradients in glass transition temperature Tg and relaxation time τ. Both of these forms of heterogeneity have a major impact on material properties. Nevertheless, their interplay has remained poorly understood. Here, we employ molecular dynamics simulations of polymer thin films in the isoconfigurational ensemble in order to probe how bulk dynamic heterogeneity alters and is altered by the large gradient in dynamics at the surface of a glass-forming liquid. Results indicate that the τ spectrum at the surface is broader than in the bulk despite being shifted to shorter times, and yet it is less spatially correlated. This is distinct from the bulk, where the τ distribution becomes broader and more spatially organized as the mean τ increases. We also find that surface gradients in slow dynamics extend further into the film than those in fast dynamics-a result with implications for how distinct properties are perturbed near an interface. None of these features track locally with changes in the heterogeneity of caging scale, emphasizing the local disconnect between these quantities near interfaces. These results are at odds with conceptions of the surface as reflecting simply a higher "rheological temperature" than the bulk, instead pointing to a complex interplay between bulk dynamic heterogeneity and spatially organized dynamical gradients at interfaces in glass-forming liquids.

Glass formation and dynamics of model polymer films with one versus two active interfaces

Polymers and other glass-forming liquids can exhibit profound alterations in dynamics in the nanoscale vicinity of interfaces, over a range appreciably exceeding that of typical interfacial thermodynamic gradients. The understanding of these dynamical gradients is particularly complicated in systems with internal or external nanoscale dimensions, where a gradient nucleated at one interface can impinge on a second, potentially distinct, interface. To better understand the interactions that govern system dynamics and glass formation in these cases, here we simulate the baseline case of a glass-forming polymer film, over a wide range of thickness, supported on a dynamically neutral substrate that has little effect on nearby dynamics. We compare these results to our prior simulations of freestanding films. Results indicate that dynamical gradients in our simulated systems, as measured based upon translational relaxation, are simply truncated when they impinge on a secondary surface that is locally dynamically neutral. Altered film behavior can be described almost entirely by gradient effects down to the thinnest films probed, with no evidence for finite-size effects sometimes posited to play a role in these systems. Finally, our simulations predict that linear gradient overlap effects in the presence of symmetric dynamically active interfaces yield a non-monotonic variation of the whole free standing film stretching exponent (relaxation time distribution breadth). The maximum relaxation time distribution breadth in simulation is found at a film thickness of 4-5 times the interfacial gradient range. Observation of this maximum in experiment would provide an important validation that the gradient behavior observed in simulation persists to experimental timescales. If validated, observation of this maximum would potentially also enable determination of the dynamic gradient range from experimental mean-film measurements of film dynamics.

Understanding creep suppression mechanisms in polymer nanocomposites through machine learning

While recent efforts have shown how local structure plays an essential role in the dynamic heterogeneity of homogeneous glass-forming materials, systems containing interfaces such as thin films or composite materials remain poorly understood. It is known that interfaces perturb the molecular packing nearby, however, numerous studies show the dynamics are modified over a much larger range. Here, we examine the dynamics in polymer nanocomposites (PNCs) using a combination of simulations and experiments and quantitatively separate the role of polymer packing from other effects on the dynamics, as a function of distance from the nanoparticle surfaces. After showing good qualitative agreement between the simulations and experiments in glassy structure and creep compliance, we use a machine-learned structure indicator, softness, to decompose polymer dynamics in our simulated PNCs into structure-dependent and structure-independent processes. With this decomposition, the free energy barrier for polymer rearrangement can be described as a combination of packing-dependent and packing-independent barriers. We find both barriers are higher near nanoparticles and decrease with applied stress, quantitatively demonstrating that the slow interfacial dynamics is not solely due to polymer packing differences, but also the change of structure-dynamics relationships. Finally, we present how this decomposition can be used to accurately predict strain-time creep curves for PNCs from their static configuration, providing additional insights into the effects of polymer-nanoparticle interfaces on creep suppression in PNCs.

Influence of the Graft Length on Nanocomposite Structure and Interfacial Dynamics

Both the dispersion state of nanoparticles (NPs) within polymer nanocomposites (PNCs) and the dynamical state of the polymer altered by the presence of the NP/polymer interfaces have a strong impact on the macroscopic properties of PNCs. In particular, mechanical properties are strongly affected by percolation of hard phases, which may be NP networks, dynamically modified polymer regions, or combinations of both. In this article, the impact on dispersion and dynamics of surface modification of the NPs by short monomethoxysilanes with eight carbons in the alkyl part (C₈) is studied. As a function of grafting density and particle content, polymer dynamics is followed by broadband dielectric spectroscopy and analyzed by an interfacial layer model, whereas the particle dispersion is investigated by small-angle X-ray scattering and analyzed by reverse Monte Carlo simulations. NP dispersions are found to be destabilized only at the highest grafting. The interfacial layer formalism allows the clear identification of the volume fraction of interfacial polymer, with its characteristic time. The strongest dynamical slow-down in the polymer is found for unmodified NPs, while grafting weakens this effect progressively. The combination of all three techniques enables a unique measurement of the true thickness of the interfacial layer, which is ca. 5 nm. Finally, the comparison between longer (C₁₈) and shorter (C₈) grafts provides unprecedented insight into the efficacy and tunability of surface modification. It is shown that C₈-grafting allows for a more progressive tuning, which goes beyond a pure mass effect.

How Tuning Interfaces Impacts the Dynamics and Structure of Polymer Nanocomposites Simultaneously

Fundamental understanding of the macroscopic properties of polymer nanocomposites (PNCs) remains difficult due to the complex interplay of microscopic dynamics and structure, namely interfacial layer relaxations and three-dimensional nanoparticle (NP) arrangements. The effect of surface modification by alkyl methoxysilanes at different grafting densities has been studied in PNCs made of poly(2-vinylpyridine) and spherical 20 nm silica NPs. The segmental dynamics has been probed by broadband dielectric spectroscopy and the filler structure by small-angle X-ray scattering and reverse Monte Carlo simulations. By combining the particle configurations with the interfacial layer properties, it is shown how surface modification tunes the attractive polymer-particle interactions: bare NPs slow down the polymer interfacial layer dynamics over a thickness of ca. 5 nm, while grafting screens these interactions. Our analysis of interparticle spacings and segmental dynamics provides unprecedented insights into the effect of surface modification on the main characteristics of PNCs: particle interactions and polymer interfacial layers.

Nonmonotonic Dynamical Correlations beneath the Surface of Glass-Forming Liquids

Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate how distinct static and dynamic length scales govern the dynamics of vitrifying films. In contrast to a monotonically growing static correlation length, the dynamical correlation length that measures the extent of surface-dynamics acceleration into the bulk displays a striking nonmonotonic temperature evolution that is robust also against changes in detailed interatomic interaction. This nonmonotonic change defines a crossover temperature T_{*} that is distinct from the critical temperature T_{c} of mode-coupling theory. We connect this nonmonotonic change to a morphological change of cooperative rearrangement regions of fast particles, and to the point where the decoupling of fast-particle motion from the bulk relaxation is most sensitive to fluctuations. We propose a rigorous definition of this new crossover temperature T_{*} within a recent extension of mode-coupling theory, the stochastic β-relaxation theory.

Screening and collective effects in randomly pinned fluids: a new theoretical framework

We propose a theoretical framework for the dynamics of bulk isotropic hard-sphere systems in the presence of randomly pinned particles and apply this theory to supercooled water to validate it. Structural relaxation is mainly governed by local and non-local activated process. As the pinned fraction grows, a local caging constraint becomes stronger and the long range collective aspect of relaxation is screened by immobile obstacles. Different responses of the local and cooperative motions results in subtle predictions for how the alpha relaxation time varies with pinning and density. Our theoretical analysis for the relaxation time of water with pinned molecules quantitatively well describe previous simulations. In addition, the thermal dependence of relaxation for unpinned bulk water is also consistent with prior computational and experimental data.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Surface diffusion of a glassy discotic organic semiconductor and the surface mobility gradient of molecular glasses

Surface diffusion has been measured in the glass of an organic semiconductor, MTDATA, using the method of surface grating decay. The decay rate was measured as a function of temperature and grating wavelength, and the results indicate that the decay mechanism is viscous flow at high temperatures and surface diffusion at low temperatures. Surface diffusion in MTDATA is enhanced by 4 orders of magnitude relative to bulk diffusion when compared at the glass transition temperature T_g. The result on MTDATA has been analyzed along with the results on other molecular glasses without extensive hydrogen bonds. In total, these systems cover a wide range of molecular geometries from rod-like to quasi-spherical to discotic and their surface diffusion coefficients vary by 9 orders of magnitude. We find that the variation is well explained by the existence of a steep surface mobility gradient and the anchoring of surface molecules at different depths. Quantitative analysis of these results supports a recently proposed double-exponential form for the mobility gradient: log D(T, z) = log D_v(T) + [log D₀ - log D_v(T)]exp(-z/ξ), where D(T, z) is the depth-dependent diffusion coefficient, D_v(T) is the bulk diffusion coefficient, D₀ ≈ 10^-8 m²/s, and ξ ≈ 1.5 nm. Assuming representative bulk diffusion coefficients for these fragile glass formers, the model reproduces the presently known surface diffusion rates within 0.6 decade. Our result provides a general way to predict the surface diffusion rates in molecular glasses.

A β-NMR study of the depth, temperature, and molecular-weight dependence of secondary dynamics in polystyrene: Entropy–enthalpy compensation and dynamic gradients near the free surface

We investigated the depth, temperature, and molecular-weight (MW) dependence of the γ-relaxation in polystyrene glasses using implanted 8Li+ and β-detected nuclear magnetic resonance. Measurements were performed on thin films with MW ranging from 1.1 to 641 kg/mol. The temperature dependence of the average 8Li spin–lattice relaxation time [Formula: see text] was measured near the free surface and in the bulk. Spin–lattice relaxation is caused by phenyl ring flips, which involve transitions between local minima over free-energy barriers with enthalpic and entropic contributions. We used transition state theory to model the temperature dependence of the γ-relaxation, and hence [Formula: see text]. There is no clear correlation of the average entropy of activation [Formula: see text] and enthalpy of activation [Formula: see text] with MW, but there is a clear correlation between [Formula: see text] and [Formula: see text], i.e., entropy–enthalpy compensation. This results in the average Gibbs energy of activation, [Formula: see text], being approximately independent of MW. Measurements of the temperature dependence of [Formula: see text] as a function of depth below the free surface indicate the inherent entropic barrier, i.e., the entropy of activation corresponding to [Formula: see text] = 0, has an exponential dependence on the distance from the free surface before reaching the bulk value. This results in [Formula: see text] near the free surface being lower than the bulk. Combining these observations results in a model where the average fluctuation rate of the γ-relaxation has a “double-exponential” depth dependence. This model can explain the depth dependence of [Formula: see text] in polystyrene films. The characteristic length of enhanced dynamics is ∼6 nm and approximately independent of MW near room temperature.

Confinement Effects on the Spatially Inhomogeneous Dynamics in Metallic Glass Films

This work develops the elastically collective nonlinear Langevin equation theory to investigate, for the first time, the glassy dynamics in capped metallic glass thin films. Finite-size effects on the spatial gradient of structural relaxation time and glass transition temperature (T_g) are calculated at different temperatures and vitrification criteria. Molecular dynamics is significantly slowed down near rough solid surfaces, and the dynamics at location far from the interfaces is sped up. In thick films, the mobility gradient normalized by the bulk value obeys the double-exponential form since interference effects between two surfaces are weak. Reducing the film thickness induces a strong dynamic coupling between two surfaces and flattens the relaxation gradient. The normalized gradient of the glass transition temperature is independent of vitrification time scale criterion and can be fitted by a superposition function as the films are not ultrathin. The local fragility is found to remain unchanged with location. This finding suggests that one can use Angell plots of bulk relaxation time and the T_g spatial gradient to characterize glassy dynamics in metallic glass films. Our computational results agree well with experimental data and simulation.

Nature of dynamic gradients, glass formation, and collective effects in ultrathin freestanding films

Molecular, polymeric, colloidal, and other classes of liquids can exhibit very large, spatially heterogeneous alterations of their dynamics and glass transition temperature when confined to nanoscale domains. Considerable progress has been made in understanding the related problem of near-interface relaxation and diffusion in thick films. However, the origin of "nanoconfinement effects" on the glassy dynamics of thin films, where gradients from different interfaces interact and genuine collective finite size effects may emerge, remains a longstanding open question. Here, we combine molecular dynamics simulations, probing 5 decades of relaxation, and the Elastically Cooperative Nonlinear Langevin Equation (ECNLE) theory, addressing 14 decades in timescale, to establish a microscopic and mechanistic understanding of the key features of altered dynamics in freestanding films spanning the full range from ultrathin to thick films. Simulations and theory are in qualitative and near-quantitative agreement without use of any adjustable parameters. For films of intermediate thickness, the dynamical behavior is well predicted to leading order using a simple linear superposition of thick-film exponential barrier gradients, including a remarkable suppression and flattening of various dynamical gradients in thin films. However, in sufficiently thin films the superposition approximation breaks down due to the emergence of genuine finite size confinement effects. ECNLE theory extended to treat thin films captures the phenomenology found in simulation, without invocation of any critical-like phenomena, on the basis of interface-nucleated gradients of local caging constraints, combined with interfacial and finite size-induced alterations of the collective elastic component of the structural relaxation process.

Polymers under nanoconfinement: where are we now in understanding local property changes?

Polymers are increasingly being used in applications with nanostructured morphologies where almost all polymer molecules are within a few tens to hundreds of nanometers from some interface. From nearly three decades of study on polymers in simplified nanoconfined systems such as thin films, we have come to understand property changes in these systems as arising from interfacial effects where local dynamical perturbations are propagated deeper into the material. This review provides a summary of local glass transition temperature T_g changes near interfaces, comparing across different types of interfaces: free surface, substrate, liquid, and polymer-polymer. Local versus film-average properties in thin films are discussed, making connections to other related property changes, while highlighting several historically important studies. By experimental necessity, most studies are on high enough molecule weight chains to be well entangled, although aspects that connect to lower molecule weight materials are described. Emphasis is made to identify observations and open questions that have yet to be fully understood such as the evidence of long-ranged interfacial effects, finite domain size, interfacial breadth, and chain connectivity.

Energy barriers and cooperative motion at the surface of freestanding glassy polystyrene films

We investigate the near-surface relaxation of freestanding atactic polystyrene films with molecular dynamics simulations. As in previous coarse-grained simulations, relaxation times for backbone segments and phenyl rings are linked to their bulk relaxation times via a power-law coupling relation. Variation of the coupling exponent with distance from the surface is consistent with depth-dependent activation barriers. We also quantify a reduction in dynamical heterogeneity at the interface, which can be interpreted in the framework of cooperative models for glassy dynamics.

Theory of microstructure-dependent glassy shear elasticity and dynamic localization in melt polymer nanocomposites

We present an integrated theoretical study of the structure, thermodynamic properties, dynamic localization, and glassy shear modulus of melt polymer nanocomposites (PNCs) that spans the three microstructural regimes of entropic depletion induced nanoparticle (NP) clustering, discrete adsorbed layer driven NP dispersion, and polymer-mediated bridging network. The evolution of equilibrium and dynamic properties with NP loading, total packing fraction, and strength of interfacial attraction is systematically studied based on a minimalist model. Structural predictions of polymer reference interaction site model integral equation theory are employed to establish the rich behavior of the interfacial cohesive force density, surface excess, and a measure of free volume as a function of PNC variables. The glassy dynamic shear modulus is predicted to be softened, reinforced, or hardly changed relative to the pure polymer melt depending on system parameters, as a result of the competing and qualitatively different influences of interfacial cohesion (physical bonding), free volume, and entropic depletion on dynamic localization and shear elasticity. The localization of polymer segments is the dominant factor in determining bulk PNC softening and reinforcement effects for moderate to strong interfacial attractions, respectively. While in the athermal entropy-dominated regime, the primary origin of mechanical reinforcement is the stress stored in the aggregated NP subsystem. The PNC shear modulus is often qualitatively correlated with the segment localization length but with notable exceptions. The present work provides the foundation for developing a theory of segmental relaxation, T_g changes, and collective NP dynamics in PNCs based on a self-consistent treatment of the cooperative activated motions of segments and NPs.

Surface diffusion in glasses of rod-like molecules posaconazole and itraconazole: effect of interfacial molecular alignment and bulk penetration

The method of surface grating decay has been used to measure surface diffusion in the glasses of two rod-like molecules posaconazole (POS) and itraconazole (ITZ). Although structurally similar antifungal medicines, ITZ forms liquid-crystalline phases while POS does not. Surface diffusion in these systems is significantly slower than in the glasses of quasi-spherical molecules of similar volume when compared at the glass transition temperature T_g. Between the two systems, ITZ has slower surface diffusion. These results are explained on the basis of the near-vertical orientation of the rod-like molecules at the surface and their deep penetration into the bulk where mobility is low. For molecular glasses without extensive hydrogen bonds, we find that the surface diffusion coefficient at T_g decreases smoothly with the penetration depth of surface molecules and the trend has the double-exponential form for the surface mobility gradient observed in simulations. This supports the view that these molecular glasses have a similar mobility vs. depth profile and their different surface diffusion rates arise simply from the different depths at which molecules are anchored. Our results also provide support for a previously observed correlation between the rate of surface diffusion and the fragility of the bulk liquid.

Theory of Spatial Gradients of Relaxation, Vitrification Temperature and Fragility of Glass-Forming Polymer Liquids Near Solid Substrates

We employ a new force-level statistical mechanical theory to predict spatial gradients of the structural relaxation time and T_g of polymer liquids near microscopically rough and smooth hard surfaces and contrast the results with vapor interface systems. Repulsive rough (smooth) surfaces induce large slowing down (modest speeding up) of the relaxation time compared to the bulk. Nevertheless, a remarkable degree of universality of distinctive dynamical behaviors is predicted for different polymer chemistries and all interfaces, including a double exponential form of the alpha time gradient, power law decoupling of the relaxation time from its bulk value with exponential spatial variation of the exponent, exponential spatial gradient of T_g, weak dependence of normalized T_g gradients on vitrification criterion, and near linear growth with cooling of the slowed down layer thickness near a rough hard interface. The results appear consistent with simulations and experiments, and multiple testable predictions are made.

Progress towards a phenomenological picture and theoretical understanding of glassy dynamics and vitrification near interfaces and under nanoconfinement

The nature of alterations to dynamics and vitrification in the nanoscale vicinity of interfaces-commonly referred to as "nanoconfinement" effects on the glass transition-has been an open question for a quarter century. We first analyze experimental and simulation results over the last decade to construct an overall phenomenological picture. Key features include the following: after a metrology- and chemistry-dependent onset, near-interface relaxation times obey a fractional power law decoupling relation with bulk relaxation; relaxation times vary in a double-exponential manner with distance from the interface, with an intrinsic dynamical length scale appearing to saturate at low temperatures; the activation barrier and vitrification temperature T_g approach bulk behavior in a spatially exponential manner; and all these behaviors depend quantitatively on the nature of the interface. We demonstrate that the thickness dependence of film-averaged T_g for individual systems provides a poor basis for discrimination between different theories, and thus we assess their merits based on the above dynamical gradient properties. Entropy-based theories appear to exhibit significant inconsistencies with the phenomenology. Diverse free-volume-motivated theories vary in their agreement with observations, with approaches invoking cooperative motion exhibiting the most promise. The elastically cooperative nonlinear Langevin equation theory appears to capture the largest portion of the phenomenology, although important aspects remain to be addressed. A full theoretical understanding requires improved confrontation with simulations and experiments that probe spatially heterogeneous dynamics within the accessible 1-ps to 1-year time window, minimal use of adjustable parameters, and recognition of the rich quantitative dependence on chemistry and interface.