Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

Photodiode with nanocrystalline Si/amorphous Si absorber bilayer

This letter reports a near-ultraviolet/visible/near-infrared n(+)-n-i-delta(i)-p photodiode with an absorber comprising a nanocrystalline silicon n layer and a hydrogenated amorphous silicon i layer. Device modeling reveals that the dominant source of reverse dark current is deep defect states in the n layer, and its magnitude is controlled by the i layer thickness. The photodiode with the 900/400 nm thick n-i layers exhibits a reverse dark current density of 3nA/cm(2) at -1V. Donor concentration and diffusion length of holes in the n layer are estimated from the capacitance-voltage characteristics and from the bias dependence of long-wavelength response, respectively. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3660725]

Monte Carlo simulation of interface alloying

An Ising-like model, with interactions ranging up to next-nearest-neighbor pairs, is used to simulate the process of interface alloying. Interactions are chosen to stabilize an intermediate "antiferromagnetic" ordered structure. The dynamics proceeds exclusively by atom-vacancy exchanges. In order to characterize the process, the time evolution of the width of the intermediate ordered region and the diffusion length is studied. Both lengths are found to follow a power-law evolution with exponents depending on the characteristic features of the model.

Spatial correlations and cross sections of clusters in the A + B -> 0 reaction

We consider the distribution of cross sections of clusters and the density-density correlation functions for the A+B¿0 reaction. We solve the reaction-diffusion equations numerically for random initial distributions of reactants. When both reactant species have the same diffusion coefficients the distribution of cross sections and the correlation functions scale with the diffusion length and obey superuniversal laws (independent of dimension). For different diffusion coefficients the correlation functions still scale, but the scaling functions depend on the dimension and on the diffusion coefficients. Furthermore, we display explicitly the peculiarities of the cluster-size distribution in one dimension.

Spatial correlations and cross sections of clusters in the A + B -> 0 reaction

We consider the distribution of cross sections of clusters and the density-density correlation functions for the A+B¿0 reaction. We solve the reaction-diffusion equations numerically for random initial distributions of reactants. When both reactant species have the same diffusion coefficients the distribution of cross sections and the correlation functions scale with the diffusion length and obey superuniversal laws (independent of dimension). For different diffusion coefficients the correlation functions still scale, but the scaling functions depend on the dimension and on the diffusion coefficients. Furthermore, we display explicitly the peculiarities of the cluster-size distribution in one dimension.

Accessory phases as a tool to investigate magmatic underplating in polymetamorphic metasediments : an example from a fossil lower crust : Ivrea Verbano zne, NW Italy

The Ivrea and the Strona-Ceneri zones, NW italy and S Switzerland, offer the possibility to study the continental crust of the Southern Alps. Because of its high metamorphic degree and the abundant Permo- Carboniferous mafic intrusions, the Ivrea Zone is classically interpreted an exposed section trough the Permian lower crust. The present work is focused here on metasedimentary slices (septa) intercalated within Permian gabbro (mafic complex). In particular I studied the evolution of accessory phases such as rutile and zircon and the chemistry of the metasediments. The septa build an irregular and discontinuous band that cut obliquely the mafic complex from its deepest part (N) to its roof (S). The chemistry of the metasediments evolves along the band and the chemical evolution can be compared with that observed in the country-rock surrounding the mafic intrusion to the NE and overprinted by a main regional metamorphic event. This suggests that the degree of chemical depletion of the septa was mainly established during the same regional metamorphic event. Moreover it suggests that incorporation of the septa within the gabbro did not modify their original stratigraphie distribution within the crust. It implies that the mafic complex has been emplaced following a dynamic substantially different from the classic model of « gabbro glacier » (Quick et al., 1992; Quick et al., 1994). It is more likely that it has been emplaced by repeated injections of sills at different depths during a protracted period of time. Zircon trace elements and U-Pb ages suggest that regional metamorphism occurred 330-320Ma, the first sills in the deepest part of the Mafic Complex are injected at ~300Ma, the mafic magmas reached higher levels in the crust at 285Ma and the magmatic activity continued locally until 275Ma. The ages of detrital cores in zircons fix the maximal sedimentation age at ~370Ma, this age corresponds therefore with the maximal age of the incorporation of the Ivrea zone within the lower crust. I propose that the Ivrea zone has been accreted to the lower crust during the Hercynian orogeny sensu lato. The analysis of detrital ages suggests that the source terrains for the Ivrea zone and those for the Strona-Ceneri zone have a completely different Palaeozoic history. The systematic analysis of rutile in partially molten metasediments of the Ivrea zone reveals the occurrence of two generations. The two generations are characterized by a different chemistry and textural distribution. A first generation is formed during pro-grade metamorphism in the restitic counterpart. The second generation is formed in the melts during cooling at the same time that part of the first generation re-equilibrate. Re-equilibration of the first generation seems to be spatially controlled by the presence of fluids. Locally the second generation forms overgrowths on the first generation. Considered the different diffusivity of U and Pb in rutile, U heterogeneities have important implication for U-Pb dating of rutile. ID-TIMS and LA-ICPMS dating coupled with a careful textural investigation (SEM) suggest that rutile grains are characterized by multiple path along which Pb diffusion can occur: volume diffusion is an important process, but intragrain and subgrain boundaries provide additional high diffusivity pathways for Pb escape and reduce drastically the effective diffusion length. -- La zone d'Ivrea et la zone de Strona-Ceneri, en Italie nord-occidentale et Suisse méridionale, offrent la possibilité d'étudier la croûte continentale des Alpes du Sud. En raison du haut degré métamorphique et l'abondance d'intrusions mafiques d'âge Permo-Carbonifère [complexe mafique), la zone d'Ivrea est interprétée classiquement comme de la croûte inférieure permienne. Ce travail ce concentre sur des bandes metasédimentaires (septa) incorporées dans les magmas mafiques lors de l'intrusion. Les septa forment une bande irrégulière qui coupe obliquement le complexe mafique du bas (N) vers le haut (S). La chimie des septa évolue du bas vers le haut et l'évolution chimique se rapproche de l'évolution observé dans la roche encaissante l'intrusion affecté par un événement métamorphique régionale. Cette relation suggère que le degré d'appauvrissement chimique des septa a été établit principalement lors de l'événement métamorphique régional. De plus l'incorporation dans les gabbros n'a pas perturbée la distribution stratigraphique originelle des septa. Ces deux observations impliquent que le métamorphisme dans la roche encaissante précède la mise en place du gabbro et que cette dernière ne se fait pas selon le modèle classique (« gabbro glacier » de Quick et al., 1992, 1994), mais se fait plutôt par injections répétées de sills a différentes profondeurs. Les âges U-Pb et les éléments traces des zircons suggèrent que le métamorphisme régionale a eu lieu 330-320Ma, alors que les premiers sills dans la partie profonde du Mafic Complex s'injectent à ~300Ma, le magmatisme mafique atteigne des niveaux supérieurs à 285Ma et continue localement jusqu'à 270Ma. Les âges des coeurs détritiques des zircons permettent de fixer l'âge maximale de sédimentation à ~370Ma ce qui correspond donc à l'âge maximale de l'incorporation de la zone d'Ivrea dans la croûte inférieur. L'analyse systématique des rutiles, nous a permit de montrer l'existence de plusieurs générations qui ont une répartition texturale et une chimie différente. Une génération se forme lors de l'événement UHT dans les restites, une autre génération se forme dans les liquides lors du refroidissement, au même temps qu'une partie de la première génération se rééquilibre au niveau du Zr. Localement la deuxième génération peut former des surcroissances autour de la première génération. Dans ces cas, des fortes différences en uranium entre les deux générations ont des importantes implications pour la datation U-Pb sur rutile. Classiquement les ratios Pb/U dans le rutile sont interprétés comme indiquant l'âges du refroidissement du minéral sous une température à la quelle la diffusion du Pb dans le minéral n'est plus détectable et la diffusion à plus hautes températures est assumée se faire par «volume diffusion» dans le grain (Mezger et al., 1989). Par des datations ID-TIMS (sur grain entier) et LA-ICPMS (in-situ) et une analyse texturale (MEB) approfondie nous montrons que cette supposition est trop simpliste et que le rutile est repartie en sous-domaines. Chacun de ces domaines a ça propre longueur ou chemin de diffusion spécifique. Nous proposons donc une nouvelle approche plus cohérente pour l'interprétation des âges U-Pb sur rutile.

Monte Carlo simulation of interface alloying

An Ising-like model, with interactions ranging up to next-nearest-neighbor pairs, is used to simulate the process of interface alloying. Interactions are chosen to stabilize an intermediate "antiferromagnetic" ordered structure. The dynamics proceeds exclusively by atom-vacancy exchanges. In order to characterize the process, the time evolution of the width of the intermediate ordered region and the diffusion length is studied. Both lengths are found to follow a power-law evolution with exponents depending on the characteristic features of the model.

Influence of nanowire density on the shape and optical properties of ternary InGaAs nanowires

We have synthesized ternary InGaAs nanowires on (111)B GaAs surfaces by metal-organic chemical vapor deposition. Au colloidal nanoparticles were employed to catalyze nanowire growth. We observed the strong influence of nanowire density on nanowire height, tapering, and base shape specific to the nanowires with high In composition. This dependency was attributed to the large difference of diffusion length on (111)B surfaces between In and Ga reaction species, with In being the more mobile species. Energy dispersive X-ray spectroscopy analysis together with high-resolution electron microscopy study of individual InGaAs nanowires shows large In/Ga compositional variation along the nanowire supporting the present diffusion model. Photoluminescence spectra exhibit a red shift with decreasing nanowire density due to the higher degree of In incorporation in more sparsely distributed InGaAs nanowires.

Microwave properties of an inhomogeneous optically illuminated plasma in a microstrip gap

The optical illumination of a microstrip gap on a thick semiconductor substrate creates an inhomogeneous electron-hole plasma in the gap region. This allows the study of the propagation mechanism through the plasma region. This paper uses a multilayer plasma model to explain the origin of high losses in such structures. Measured results are shown up to 50 GHz and show good agreement with the simulated multilayer model. The model also allows the estimation of certain key parameters of the plasma, such as carrier density and diffusion length, which are difficult to measure by direct means. The detailed model validation performed here will enable the design of more complex microwave structures based on this architecture. While this paper focuses on monocrystalline silicon as the substrate, the model is easily adaptable to other semiconductor materials such as GaAs.

Understanding the Reaction Mechanism of Nanocomposite Thermites

Nanocomposite energetics are a relatively new class of materials that combine nanoscale fuels and oxidizers to allow for the rapid release of large amounts of energy. In thermite systems (metal fuel with metal oxide oxidizer), the use of nanomaterials has been illustrated to increase reactivity by multiple orders of magnitude as a result of the higher specific surface area and smaller diffusion length scales. However, the highly dynamic and nanoscale processes intrinsic to these materials, as well as heating rate dependencies, have limited our understanding of the underlying processes that control reaction and propagation. For my dissertation, I have employed a variety of experimental approaches that have allowed me to probe these processes at heating rates representative of free combustion with the goal of understanding the fundamental mechanisms. Dynamic transmission electron microscopy (DTEM) was used to study the in situ morphological change that occurs in nanocomposite thermite materials subjected to rapid (10^11 K/s) heating. Aluminum nanoparticle (Al-NP) aggregates were found to lose their nanostructure through coalescence in as little as 10 ns, which is much faster than any other timescale of combustion. Further study of nanoscale reaction with CuO determined that a condensed phase interfacial reaction could occur within 0.5-5 µs in a manner consistent with bulk reaction, which supports that this mechanism plays a dominant role in the overall reaction process. Ta nanocomposites were also studied to determine if a high melting point (3280 K) affects the loss of nanostructure and rate of reaction. The condensed phase reaction pathway was further explored using reactive multilayers sputter deposited onto thin Pt wires to allow for temperature jump (T-Jump) heating at rates of ~5x10^5 K/s. High speed video and a time of flight mass spectrometry (TOFMS) were used to observe ignition temperature and speciation as a function of bilayer thickness. The ignition process was modeled and a low activation energy for effective diffusivity was determined. T-Jump TOFMS along with constant volume combustion cell studies were also used to determine the effect of gas release in nanoparticle systems by comparing the reaction properties of CuO and Cu2O.

Dysprosium transport in Nd-Fe-B pellets

The addition of heavy rare earth (RE) elements to Nd2Fe14B based magnets to form (Nd,Dy)2Fe14B is known to increase the coercivity and high temperature performance required for hybrid vehicle electric motors and other extreme temperature applications. Attempts to conserve heavy rare earth elements for high temperature (RE)2Fe14B based magnets have led to the development of a grain boundary diffusion process for bulk magnets. This process relies on transport of a heavy rare earth, such as Dy, into a bulk Nd2Fe14B magnet along pores, a low volume fraction of eutectic liquid along grain boundary grain triple junctions and grain boundaries. This enriches the grain surfaces in Dy through the thickness of the bulk magnet, leading to larger increases coercivity with a smaller Dy concentration than can be achieved with homogeneous alloys. Attempts to carry out the same process during sintering require significant control of Dy transport efficiency. The macroscopic transport of Dy in Nd2.7Fe14B1.4 based powder packs is studied using a 'layered' pellet, where Nd2.7Fe14B1.4powder is an interlayer and Dy source as a center layer. The sintering of this layered pellet provided evidence for very large effective diffusion lengths aided by Dy rich liquid flow through connected porosity. Approaches to controlling Dy transportation include decreasing the liquid phase transport capability of the powder pack by increasing the melting point of the Dy source and the decreasing amount of RE rich liquid in the powder packs. The solid-liquid reaction is studied in which melt spun Nd2.7Fe14B1.4 ribbons are PVD coated with Dy-Fe eutectic composition and then thermally treated. The resulting microstructure from the reaction between Dy-Fe eutectic coating and Nd2.7Fe14B1.4 ribbon is interpreted as support for a proposed dissolution/reprecipitation process between solid and liquid phases. The estimate the diffusion coefficient and the effective diffusion length of Dy sources in Nd2.7Fe14B1.4 layered pellets and melt spun ribbons were obtained from the calculation of Fick's second law combined with EDS results from the experiment. The results indicate that the effective diffusion coefficient of Dy in the layered pellets is higher than the diffusion in ribbons due to its higher porosity than ribbons.

Differential diffusion: Often a finite-mising length effect

Many instances of differential diffusion, i e, different species having different turbulent diffusion coefficients in the same flow, can be explained as a finite mixing length effect. That is, in a simple mixing length scenario, the turbulent diffusion coefficient has the form 1 ( m )2 m m c l K w l OL   =  +    where, wm is the mixing velocity, lm the mixing length and Lc the overall distribution scale for a particular species. The first term represents the familiar gradient diffusion while the second term becomes important when lm/Lc is finite. This second term shows that different species will have different diffusion coefficients if they have different overall distribution scales. Such different Lcs may come about due to different boundary conditions and different intrinsic properties (molecular diffusivity, settling velocity etc) for different species. For momentum transfer in turbulent oscillatory boundary layers the second term is imaginary and explains observed phase leads of shear stresses ahead of velocity gradients.

Transient response of diffusion cell containing a porous solid

Transient response of an adsorbing or non-adsorbing tracer injected as step or square pulse input in a diffusion cell with two flowing streams across the pellet is theoretically investigated in this paper. Exact solutions and the asymptotic solutions in the time domain and in three different limits are obtained by using an integral transform technique and a singular perturbation technique, respectively. Parametric dependence of the concentrations in the top and bottom chambers can be revealed by investigating the asymptotic solutions, which are far simpler than their exact counterpart. In the time domain investigation, it is found that the bottom-chamber concentration is very sensitive to the value of the macropore effective diffusivity. Therefore this concentration could be used to extract diffusivity by fitting in the time domain. The bottom-chamber concentration is also sensitive to flow rate, pellet length chamber volume and the type of input (step and square input).