Rapid decay of vacancy islands at step edges on Ag(111): step orientation dependence

Previous work has established that vacancy islands or pits fill much more quickly when they are in contact with a step edge, such that the common boundary is a double step. The present work focuses on the effect of the orientation of that step, with two possibilities existing for a face centered cubic (111) surface: A- and B-type steps. We find that the following features can depend on the orientation: (1) the shapes of islands while they shrink; (2) whether the island remains attached to the step edge; and (3) the rate of filling. The first two effects can be explained by the different rates of adatom diffusion along the A- and B-steps that define the pit, enhanced by the different filling rates. The third observation--the difference in the filling rate itself--is explained within the context of the concerted exchange mechanism at the double step. This process is facile at all regular sites along B-steps, but only at kink sites along A-steps, which explains the different rates. We also observe that oxygen can greatly accelerate the decay process, although it has no apparent effect on an isolated vacancy island (i.e. an island that is not in contact with a step).

Polymer length distributions for catalytic polymerization within mesoporous materials: non-Markovian behavior associated with partial extrusion

We analyze a model for polymerization at catalytic sites distributed within parallel linear pores of a mesoporous material. Polymerization occurs primarily by reaction of monomers diffusing into the pores with the ends of polymers near the pore openings. Monomers and polymers undergo single-file diffusion within the pores. Model behavior, including the polymer length distribution, is determined by kinetic Monte Carlo simulation of a suitable atomistic-level lattice model. While the polymers remain within the pore, their length distribution during growth can be described qualitatively by a Markovian rate equation treatment. However, once they become partially extruded, the distribution is shown to exhibit non-Markovian scaling behavior. This feature is attributed to the long-tail in the "return-time distribution" for the protruding end of the partially extruded polymer to return to the pore, such return being necessary for further reaction and growth. The detailed form of the scaled length distribution is elucidated by application of continuous-time random walk theory.

The changing face of emotion: age-related patterns of amygdala activation to salient faces

The present study investigated age-related differences in the amygdala and other nodes of face-processing networks in response to facial expression and familiarity. fMRI data were analyzed from 31 children (3.5-8.5 years) and 14 young adults (18-33 years) who viewed pictures of familiar (mothers) and unfamiliar emotional faces. Results showed that amygdala activation for faces over a scrambled image baseline increased with age. Children, but not adults, showed greater amygdala activation to happy than angry faces; in addition, amygdala activation for angry faces increased with age. In keeping with growing evidence of a positivity bias in young children, our data suggest that children find happy faces to be more salient or meaningful than angry faces. Both children and adults showed preferential activation to mothers' over strangers' faces in a region of rostral anterior cingulate cortex associated with self-evaluation, suggesting that some nodes in frontal evaluative networks are active early in development. This study presents novel data on neural correlates of face processing in childhood and indicates that preferential amygdala activation for emotional expressions changes with age.

Group specific optimisation of fMRI processing steps for child and adult data

Motion is a major issue in functional magnetic resonance imaging (fMRI) dataseries and causes artifacts or increased overall noise obscuring signals of interest. It is particularly important to be able to control for and correct these artifacts when dealing with child data. We analysed the data from 35 children (4-8 years old) and 13 adults (18-30 years old) during an emotional face paradigm. The children were split into low and high motion groups on the basis of having less or more than an estimated maximal movement of one voxel (3.75 mm) and one degree of rotation in any motion direction between any pair of scans in the run. Several different preprocessing steps were evaluated for their ability to correct for the excess motion using agnostic canonical variates analysis (aCVA) in the NPAIRS (Nonparametric, Prediction, Activation, Influence, Reproducibility, re-Sampling) framework. The adult dataset was reasonably stable whereas the motion-prone child datasets benefited greatly from motion parameter regression (MPR). Motion parameter regression had a strong beneficial impact on all datasets, a result that was largely unaffected by other preprocessing choices; however, motion correction on its own did not have as much impact. The low motion child group subjected to MPR had reproducibility values at par with those of the adult group, but needed almost twice as many subjects to achieve this result, indicating weaker responses in young children. The aCVA showed greater sensitivity to the task response pattern than the mixed effects general linear model (mGLM) in the expected face processing regions, although the mGLM showed more responses in some other areas. This work illustrates that preprocessing choices must be made in a group-specific fashion to optimise fMRI results.

A kinetic Monte Carlo study on the role of defects and detachment in the formation and growth of In chains on Si(100)

Deposition on a Si(100) surface and subsequent self-assembly of In atoms into one-dimensional (1D) atomic chains at room temperature is investigated via kinetic Monte Carlo simulation of a suitable atomistic model. Model development is guided by recent experimental observations in which 1D In chains nucleate effectively exclusively at C-type defects, although In atoms can detach from chains. We find that a monotonically decreasing form of the scaled island size distribution (ISD) is consistent with a high defect density which facilitates persistent chain nucleation even at relatively high coverages. The predominance of heterogeneous nucleation may be attributed to several factors including low surface diffusion barriers, a high defect density, and relatively weak In-In binding.

Neural correlates of personally familiar faces: parents, partner and own faces

Investigations of the neural correlates of face recognition have typically used old/new paradigms where subjects learn to recognize new faces or identify famous faces. Familiar faces, however, include one's own face, partner's and parents' faces. Using event-related fMRI, we examined the neural correlates of these personally familiar faces. Ten participants were presented with photographs of own, partner, parents, famous and unfamiliar faces and responded to a distinct target. Whole brain, two regions of interest (fusiform gyrus and cingulate gyrus), and multiple linear regression analyses were conducted. Compared with baseline, all familiar faces activated the fusiform gyrus; own faces also activated occipital regions and the precuneus; partner faces activated similar areas, but in addition, the parahippocampal gyrus, middle superior temporal gyri and middle frontal gyrus. Compared with unfamiliar faces, only personally familiar faces activated the cingulate gyrus and the extent of activation varied with face category. Partner faces also activated the insula, amygdala and thalamus. Regions of interest analyses and laterality indices showed anatomical distinctions of processing the personally familiar faces within the fusiform and cingulate gyri. Famous faces were right lateralized whereas personally familiar faces, particularly partner and own faces, elicited bilateral activations. Regression analyses show experiential predictors modulated with neural activity related to own and partner faces. Thus, personally familiar faces activated the core visual areas and extended frontal regions, related to semantic and person knowledge and the extent and areas of activation varied with face type.

Statistical mechanical modeling of catalytic polymerization within surface-functionalized mesoporous materials

A discrete lattice model is developed to describe diffusion-mediated polymerization occurring within mesopores, where reaction is enhanced at catalytic sites distributed within the interior of the pores. Diffusive transport of monomers and polymers is one-dimensional, diffusion coefficients for the latter decreasing with polymer length. Kinetic Monte Carlo simulation is utilized to analyze model behavior focusing on a "clogging" regime, where the amount of polymer within the pores grows. We characterize the evolution of the overall and mean length of polymers, the mean number of polymers, as well as the polymer spatial and length distributions.

Structure and growth of height-selected Ag islands on fivefold i-AlPdMn quasicrystalline surfaces: STM analysis and step dynamics modeling

The development and local structure of height-selected 3-layer Ag islands on fivefold surfaces of icosahedral Al-Pd-Mn quasicrystals is characterized by STM for Ag deposition at 365 K. Heterogeneous nucleation of pseudomorphic single layer high islands is followed by rapid formation of 2nd and 3rd layers and subsequent lateral spreading, where each of these 3 layers consists of a family of nonfcc structures. The behavior is elucidated by step dynamics modeling incorporating strain buildup for larger islands, enhanced binding in higher layers, and height selection due to quantum size effects.

Formation of complex wedding-cake morphologies during homoepitaxial film growth of Ag on Ag(111): atomistic, step-dynamics, and continuum modeling

An atomistic lattice-gas model is developed which successfully describes all key features of the complex mounded morphologies which develop during deposition of Ag films on Ag(111) surfaces. We focus on this homoepitaxial thin film growth process below 200 K. The unstable multilayer growth mode derives from the presence of a large Ehrlich-Schwoebel step-edge barrier, for which we characterize both the step-orientation dependence and the magnitude. Step-dynamics modeling is applied to further characterize and elucidate the evolution of the vertical profiles of these wedding-cake-like mounds. Suitable coarse-graining of these step-dynamics equations leads to instructive continuum formulations for mound evolution.

Schloegl's second model for autocatalysis with particle diffusion: Lattice-gas realization exhibiting generic two-phase coexistence

We analyze a discontinuous nonequilibrium phase transition between an active (or reactive) state and a poisoned (or extinguished) state occurring in a stochastic lattice-gas realization of Schloegl's second model for autocatalysis. This realization, also known as the quadratic contact process, involves spontaneous annihilation, autocatalytic creation, and diffusion of particles on a square lattice, where creation at empty sites requires a suitable nearby pair of particles. The poisoned state exists for all annihilation rates p>0 and is an absorbing particle-free "vacuum" state. The populated active steady state exists only for p below a critical value, p(e). If p(f) denotes the critical value below which a finite population can survive, then we show that p(f)<p(e). This strict inequality contrasts a postulate of Durrett, and is a direct consequence of the occurrence of coexisting stable active and poisoned states for a finite range p(f)<or=p<or=p(e) (which shrinks with increasing diffusivity). This so-called generic two-phase coexistence markedly contrasts behavior in thermodynamic systems. However, one still finds metastability and nucleation phenomena similar to those in discontinuous equilibrium transitions.

Kinetics of facile bilayer island formation at low temperature: Ag/NiAl(110)

Facile nucleation and growth of bilayer Ag(110) islands on NiAl(110) is observed by STM for Ag deposition at temperatures as low as 127 K. Density functional theory analysis for supported Ag films determines adatom adsorption energies (which favor bilayer islands), interaction energies, and diffusion barriers. Analysis of an atomistic lattice-gas model incorporating these energies elucidates the role of strongly anisotropic interactions in enabling the upward mass transport needed for bilayer island formation.

Simultaneous imaging of human cone mosaic with adaptive optics enhanced scanning laser ophthalmoscopy and high-speed transversal scanning optical coherence tomography

We describe a novel instrument capable of acquiring, simultaneously, adaptive optics enhanced scanning laser ophthalmoscopy and optical coherence tomography (OCT) images of the human cone mosaic in vivo. The OCT system is based on transversal scanning of the sample with a line scan rate of 14 kHz, approximately 20 times faster than a previously reported instrument. We demonstrate the capability of this instrument with the measurement of the human cone spacing in perifoveal retina.

Generic two-phase coexistence, relaxation kinetics, and interface propagation in the quadratic contact process: simulation studies

The quadratic contact process is formulated as an adsorption-desorption model on a two-dimensional square lattice. It involves random adsorption at empty sites and correlated desorption requiring diagonally adjacent pairs of empty neighbors. We assess the model behavior utilizing kinetic Monte Carlo simulations. One finds generic two-phase coexistence between a low-coverage active steady state and a completely covered or "poisoned" absorbing steady state; i.e., both states are stable over a finite range of adsorption rates or "pressures." This behavior is in marked contrast to that for equilibrium phase separation. For spatially homogeneous systems, we provide a comprehensive characterization of the kinetics of relaxation to the steady states. We analyze rapid poisoning for higher pressures above an effective spinodal point terminating a metastable active state, nucleation-mediated poisoning in the metastable region, the dynamics of poisoned droplets within the two-phase coexistence region, and behavior reminiscent of bootstrap percolation dynamics for lower pressures. For spatially inhomogeneous systems, we analyze the propagation of planar interfaces between active and absorbing states, fully characterizing an orientation dependence which underlies the generic two-phase coexistence.

Quadratic contact process: phase separation with interface-orientation-dependent equistability

The quadratic contact process is implemented on a square lattice as a model with random adsorption and correlated desorption requiring empty pairs of diagonal neighbors. The model exhibits a discontinuous phase transition between an active state and an absorbing state, but equistability between these states depends on the orientation of the separating interface. Correspondingly, for a generalized class of models, we find phase coexistence over a finite region of their two-dimensional parameter space. This is in stark contrast to behavior in equilibrium systems.

Chemical diffusion of CO in mixed CO+O adlayers and reaction-front propagation in CO oxidation on Pd(100)

Within the framework of a realistic atomistic lattice-gas model, we present the theoretical formulation and simulation procedures for precise analysis of the chemical diffusion flux of highly mobile CO within a nonuniform interacting mixed CO + O adlayer on a Pd(100) surface. The approach applies in both regimes of relatively immobile unequilibrated and fairly mobile near-equilibrated O adlayer distributions. Spatiotemporal behavior in surface reactions is controlled by chemical diffusion in mixed adlayers. Thus, we naturally integrate the above analysis with a previously developed multiscale modeling strategy to describe mesoscale reaction front propagation in CO oxidation on Pd(100). This treatment avoids using a simplified prescription of chemical diffusion and reaction kinetics as in traditional mean-field reaction-diffusion equation approaches.

Atomistic lattice-gas modeling of CO oxidation on Pd(100): Temperature-programed spectroscopy and steady-state behavior

We have developed an atomistic lattice-gas model for the catalytic oxidation of CO on single-crystal Pd(100) surfaces under ultrahigh vacuum conditions. This model necessarily incorporates an detailed description of adlayer ordering and adsorption-desorption kinetics both for CO on Pd(100), and for oxygen on Pd(100). Relevant energetic parameters are determined by comparing model predictions with experiment, together with some guidance from density functional theory calculations. The latter also facilitates description of the interaction and reaction of adsorbed CO and oxygen. Kinetic Monte Carlo simulations of this reaction model are performed to predict temperature-programmed reaction spectra, as well as steady-state bifurcation behavior.

Theoretical analysis of mound slope selection during unstable multilayer growth

A "step dynamics" model is developed for mound formation during multilayer homoepitaxy. Downward funneling of atoms deposited at step edges is incorporated and controls mound slope selection. Behavior of the selected slope differs from that predicted by phenomenological continuum treatments where the lateral mass current vanishes identically. Instead, this current is shown to vary periodically and vanish only on average. An exact coarse-grained continuum formulation with appropriate boundary conditions is derived and recovers step dynamics results.

Critical behavior in an atomistic model for a bistable surface reaction: CO oxidation with rapid CO diffusion

We analyze critical behavior associated with the loss of bistability for an atomistic model for CO oxidation on surfaces in the limit of infinite diffusion of CO. The model includes infinite nearest-neighbor repulsions between adsorbed immobile O. We use a "hybrid" treatment incorporating a lattice-gas description of the O adlayer, but tracking just the number of adsorbed CO (which are randomly distributed on non-O sites). The critical exponents obtained from a finite-size-scaling analysis on LxL site surfaces with periodic boundary conditions show that the "hybrid" reaction model belongs to the mean-field universality class, despite strong spatial correlations in the O adlayer. We also quantify finite-size effects in the global bifurcation diagram, revealing a significant shift of the bistable region with decreasing system size. Our study elucidates fluctuation effects observed in experiments of CO oxidation on the nanoscale facets of metal-field emitter tips.

Passive and active oxidation of Si(100) by atomic oxygen: a theoretical study of possible reaction mechanisms

Reaction mechanisms for oxidation of the Si(100) surface by atomic oxygen were studied with high-level quantum mechanical methods in combination with a hybrid QM/MM (Quantum mechanics/Molecular Mechanics) method. Consistent with previous experimental and theoretical results, three structures, "back-bond", "on-dimer", and "dimer-bridge", are found to be the most stable initial surface products for O adsorption (and in the formation of SiO(2) films, i.e., passive oxidation). All of these structures have significant diradical character. In particular, the "dimer-bridge" is a singlet diradical. Although the ground state of the separated reactants, O+Si(100), is a triplet, once the O atom makes a chemical bond with the surface, the singlet potential energy surface is the ground state. With mild activation energy, these three surface products can be interconverted, illustrating the possibility of the thermal redistribution among the initial surface products. Two channels for SiO desorption (leading to etching, i.e., active oxidation) have been found, both of which start from the back-bond structure. These are referred to as the silicon-first (SF) and oxygen-first (OF) mechanisms. Both mechanisms require an 89.8 kcal/mol desorption barrier, in good agreement with the experimental estimates of 80-90 kcal/mol. "Secondary etching" channels occurring after initial etching may account for other lower experimental desorption barriers. The calculated 52.2 kcal/mol desorption barrier for one such secondary etching channel suggests that the great variation in reported experimental barriers for active oxidation may be due to these different active oxidation channels.

Catalytic CO oxidation on nanoscale Pt facets: effect of interfacet CO diffusion on bifurcation and fluctuation behavior

We present lattice-gas modeling of the steady-state behavior in CO oxidation on the facets of nanoscale metal clusters, with coupling via interfacet CO diffusion. The model incorporates the key aspects of the reaction process, such as rapid CO mobility within each facet and strong nearest-neighbor repulsion between adsorbed O. The former justifies our use of a "hybrid" simulation approach treating the CO coverage as a mean-field parameter. For an isolated facet, there is one bistable region where the system can exist in either a reactive state (with high oxygen coverage) or a (nearly CO-poisoned) inactive state. Diffusion between two facets is shown to induce complex multistability in the steady states of the system. The bifurcation diagram exhibits two regions with bistabilities due to the difference between adsorption properties of the facets. We explore the role of enhanced fluctuations in the proximity of a cusp bifurcation point associated with one facet in producing transitions between stable states on that facet, as well as their influence on fluctuations on the other facet. The results are expected to shed more light on the reaction kinetics for supported catalysts.

Metabolism of tirapazamine by multiple reductases in the nucleus

Tirapazamine (TPZ, 3-amino-1,2,4-benzotriazine 1,4-di-N-oxide, SR4233, Tirazone), a bioreductive drug currently in clinical trials, is selectively toxic to hypoxic cells commonly found in solid tumors. The toxicity results from the intracellular metabolism of TPZ to a highly toxic radical. When oxygen levels are low, the TPZ radical reacts with cellular molecules, producing DNA damage and cell death. The much lower toxicity towards aerobic cells results from the back-oxidation of the TPZ radical by oxygen. A major unresolved aspect of the mechanism of TPZ is the identity of the reductase(s) in the cell responsible for activating the drug to its toxic form. We have studied both the metabolism of the drug using HPLC and the formation of the TPZ radical with a fluorescence assay using dihydrorhodamine 123. We also measured DNA double- and single-strand breaks produced by TPZ, using the comet assay. We demonstrated that multiple reductases in the nucleus metabolize TPZ under hypoxia. Using the cofactor dependence of the reductases for metabolizing TPZ and of the DNA damage caused by TPZ, we show that DNA single-strand breaks after TPZ metabolism are probably caused by the most abundant source of reductase in the nucleus. DNA double-strand breaks, on the other hand, are formed by TPZ metabolism by an unknown nuclear reductase that requires only NADPH for its activity. This study is the first to characterize multiple nuclear reductases capable of activating TPZ.

Evolution of two-dimensional wormlike nanoclusters on metal surfaces

A pinch-off phenomenon is discovered in the evolution of 2D wormlike nanoclusters formed in homoepitaxial adlayers. This feature is shown to distinguish mass transport via periphery diffusion from other mechanisms. Continuum modeling of such evolution accurately describes experimental observations, particularly if one incorporates the anisotropy in step-edge line tension.