Multiscale coupling and multiphysics approaches in earth sciences : applications

Geoscientists are confronted with the challenge of assessing nonlinear phenomena that result from multiphysics coupling across multiple scales from the quantum level to the scale of the earth and from femtoseconds to the 4.5 Ga of history of our planet. We neglect in this review electromagnetic modelling of the processes in the Earth’s core, and focus on four types of couplings that underpin fundamental instabilities in the Earth. These are thermal (T), hydraulic (H), mechanical (M) and chemical (C) processes which are driven and controlled by the transfer of heat to the Earth’s surface. Instabilities appear as faults, folds, compaction bands, shear/fault zones, plate boundaries and convective patterns. Convective patterns emerge from buoyancy overcoming viscous drag at a critical Rayleigh number. All other processes emerge from non-conservative thermodynamic forces with a critical critical dissipative source term, which can be characterised by the modified Gruntfest number Gr. These dissipative processes reach a quasi-steady state when, at maximum dissipation, THMC diffusion (Fourier, Darcy, Biot, Fick) balance the source term. The emerging steady state dissipative patterns are defined by the respective diffusion length scales. These length scales provide a fundamental thermodynamic yardstick for measuring instabilities in the Earth. The implementation of a fully coupled THMC multiscale theoretical framework into an applied workflow is still in its early stages. This is largely owing to the four fundamentally different lengths of the THMC diffusion yardsticks spanning micro-metre to tens of kilometres compounded by the additional necessity to consider microstructure information in the formulation of enriched continua for THMC feedback simulations (i.e., micro-structure enriched continuum formulation). Another challenge is to consider the important factor time which implies that the geomaterial often is very far away from initial yield and flowing on a time scale that cannot be accessed in the laboratory. This leads to the requirement of adopting a thermodynamic framework in conjunction with flow theories of plasticity. This framework allows, unlike consistency plasticity, the description of both solid mechanical and fluid dynamic instabilities. In the applications we show the similarity of THMC feedback patterns across scales such as brittle and ductile folds and faults. A particular interesting case is discussed in detail, where out of the fluid dynamic solution, ductile compaction bands appear which are akin and can be confused with their brittle siblings. The main difference is that they require the factor time and also a much lower driving forces to emerge. These low stress solutions cannot be obtained on short laboratory time scales and they are therefore much more likely to appear in nature than in the laboratory. We finish with a multiscale description of a seminal structure in the Swiss Alps, the Glarus thrust, which puzzled geologists for more than 100 years. Along the Glarus thrust, a km-scale package of rocks (nappe) has been pushed 40 km over its footwall as a solid rock body. The thrust itself is a m-wide ductile shear zone, while in turn the centre of the thrust shows a mm-cm wide central slip zone experiencing periodic extreme deformation akin to a stick-slip event. The m-wide creeping zone is consistent with the THM feedback length scale of solid mechanics, while the ultralocalised central slip zones is most likely a fluid dynamic instability.

Alternation of chemoselective control in Stille-Heck and Heck-Stille reaction sequences

Sequential and one-pot Stille–Heck and Heck–Stille reaction processes have been invoked to give divergent access to polycyclic ring systems. Both reaction conditions and substrate structure are important in determining the nature of the reaction products formed. The Heck–Stille reactions have involved a reversal of the usual Heck regioselectivity and both cine- and ipso-substitutions have been observed in the Stille reaction.

Optical coupling of Au nanoparticles on vertical graphenes to maximize SERS response

Gold particle interaction with few-layer graphenes is of interest for the development of numerous optical nanodevices. The results of numerical studies of the coupling of gold nanoparticles with few-layer vertical graphene sheets are presented. The field strengths are computed and the optimum nanoparticle configurations for the formation of SERS hotpots are obtained. The nanoparticles are modeled as 8 nm diameter spheres atop 1.5 nm (5 layers) graphene sheet. The vertical orientation is of particular interest as it is possible to use both sides of the graphene structure and potentially double the number of particles in the system. Our results show that with the addition of an opposing particle a much stronger signal can be obtained as well as the particle separation can be controlled by the number of atomic carbon layers. These results provide further insights and contribute to the development of next-generation plasmonic devices based on nanostructures with hybrid dimensionality.

Residential balconies with road traffic noise : a combined source image and radiosity model

Several significant studies have been made in recent decades toward understanding road traffic noise and its effects on residential balconies. These previous studies have used a variety of techniques such as theoretical models, scale models and measurements on real balconies. The studies have considered either road traffic noise levels within the balcony space or inside an adjacent habitable room or both. Previous theoretical models have used, for example, simplified specular reflection calculations, boundary element methods (BEM), adaptations of CoRTN or the use of Sabine Theory. This paper presents an alternative theoretical model to predict the effects of road traffic noise spatially within the balcony space. The model includes a specular reflection component by calculating up to 10 orders of source images. To account for diffusion effects, a two compartment radiosity component is utilised. The first radiosity compartment is the urban street, represented as a street with building facades on either side. The second radiosity compartment is the balcony space. The model is designed to calculate the predicted road traffic noise levels within the balcony space and is capable of establishing the effect of changing street and balcony geometries. Screening attenuation algorithms are included to determine the effects of solid balcony parapets and balcony ceiling shields.

Fault detection of slow speed rolling element bearing with noise removal techniques

Rolling Element Bearings (REBs) are vital components in rotating machineries for providing rotating motion. In slow speed rotating machines, bearings are normally subjected to heavy static loads and a catastrophic failure can cause enormous disruption to production and human safety. Due to its low operating speed the impact energy generated by the rotating elements on the defective components is not sufficient to produce a detectable vibration response. This is further aggravated by the inability of general measuring instruments to detect and process the weak signals at the initiation of the defect accurately. Furthermore, the weak signals are often corrupted by background noise. This is a serious problem faced by maintenance engineers today and the inability to detect an incipient failure of the machine can significantly increases the risk of functional failure and costly downtime. This paper presents the application of noise removal techniques for enhancing the detection capability for slow speed REB condition monitoring. Blind deconvolution (BD) and adaptive line enhancer (ALE) are compared to evaluate their performance in enhancing the source signal with consequential removal of background noise. In the experimental study, incipient defects were seeded on a number of roller bearings and the signals were acquired using acoustic emission (AE) sensor. Kurtosis and modified peak ratio (mPR) were used to determine the detectability of signal corrupted by noise.

Identification of broad specificity P450CAM variants by primary screening against indole as substrate

High-throughput screening of cytochrome P450CAM libraries, for their ability to oxidise indole to indigo and indirubin, has resulted in the identification of variants with activity towards the structurally unrelated substrate diphenylmethane.

Substrate choice of membrane-type 1 matrix metalloproteinase is dictated by tissue inhibitor of metalloproteinase-2 levels

Although tissue inhibitor of metalloproteinase-2 (TIMP-2) is known to be not only an inhibitor of matrix metalloproteinases (MMP) but also a cofactor for membrane-type 1 MMP (MT1-MMP)-mediated MMP-2 activation, it is still unclear how TIMP-2 regulates MMP-2 activation and cleavage of substrates by MT1-MMP. In the present study we examined the levels of cell-surface MT1-MMP, MMP-2 activation and cleavage of MT1-MMP substrates in 293T cells transfected with the MT1-MMP and TIMP-2 genes. Co-expression of TIMP-2 at an appropriate level increased the level of cell-surface MT1-MMP, both the TIMP-2-bound and free forms, and generated processed MMP-2 with gelatin-degrading activity. In contrast, MT1-MMP substrates testican-1 and syndecan-1 were cleaved by the cells expressing MT1-MMP, which was inhibited by TIMP-2 even at levels that stimulate MMP-2 activation. These results suggest that TIMP-2 environment determines MT1-MMP substrate choice between direct cleavage of its own substrates and MMP-2 activation.

Using the Gini coefficient with BIOLOG substrate utilisation data to provide an alternative quantitative measure for comparing bacterial soil communities

A measure quantifying unequal use of carbon sources, the Gini coefficient (G), has been developed to allow comparisons of the observed functional diversity of bacterial soil communities. This approach was applied to the analysis of substrate utilisation data obtained from using BIOLOG microtiter plates in a study which compared decomposition processes in two contrasting plant substrates in two different soils. The relevance of applying the Gini coefficient as a measure of observed functional diversity, for soil bacterial communities is evaluated against the Shannon index (H) and average well colour development (AWCD), a measure of the total microbial activity. Correlation analysis and analysis of variance of the experimental data show that the Gini coefficient, the Shannon index and AWCD provided similar information when used in isolation. However, analyses based on the Gini coefficient and the Shannon index, when total activity on the microtiter plates was maintained constant (i.e. AWCD as a covariate), indicate that additional information about the distribution of carbon sources being utilised can be obtained. We demonstrate that the Lorenz curve and its measure of inequality, the Gini coefficient, provides not only comparable information to AWCD and the Shannon index but when used together with AWCD encompasses measures of total microbial activity and absorbance inequality across all the carbon sources. This information is especially relevant for comparing the observed functional diversity of soil microbial communities.

Capacity of fatty acid profiles and substrate utilization patterns to describe differences in soil microbial communities associated with increased salinity or alkalinity at three locations in South Australia

Fatty acid methyl ester (FAME) profiles, together with Biolog substrate utilization patterns, were used in conjunction with measurements of other soil chemical and microbiological properties to describe differences in soil microbial communities induced by increased salinity and alkalinity in grass/legume pastures at three sites in SE South Australia. Total ester-linked FAMEs (EL-FAMEs) and phospholipid-linked FAMEs (PL-FAMEs), were also compared for their ability to detect differences between the soil microbial communities. The level of salinity and alkalinity in affected areas of the pastures showed seasonal variation, being greater in summer than in winter. At the time of sampling for the chemical and microbiological measurements (winter) only the affected soil at site 1 was significantly saline. The affected soils at all three sites had lower organic C and total N concentrations than the corresponding non-affected soils. At site 1 microbial biomass, CO 2-C respiration and the rate of cellulose decomposition was also lower in the affected soil compared to the non-affected soil. Biomarker fatty acids present in both the EL- and PL-FAME profiles indicated a lower ratio of fungal to bacterial fatty acids in the saline affected soil at site 1. Analysis of Biolog substrate utilization patterns indicated that the bacterial community in the affected soil at site 1 utilized fewer carbon substrates and had lower functional diversity than the corresponding community in the non-affected soil. In contrast, increased alkalinity, of major importance at sites 2 and 3, had no effect on microbial biomass, the rate of cellulose decomposition or functional diversity but was associated with significant differences in the relative amounts of several fatty acids in the PL-FAME profiles indicative of a shift towards a bacterial dominated community. Despite differences in the number and relative amounts of fatty acids detected, principal component analysis of the EL- and PL-FAME profiles were equally capable of separating the affected and non-affected soils at all three sites. Redundancy analysis of the FAME data showed that organic C, microbial biomass, electrical conductivity and bicarbonate-extractable P were significantly correlated with variation in the EL-FAME profiles, whereas pH, electrical conductivity, NH 4-N, CO 2-C respiration and the microbial quotient were significantly correlated with variation in the PL-FAME profiles. Redundancy analysis of the Biolog data indicated that cation exchange capacity and bicarbonate-extractable K were significantly correlated with the variation in Biolog substrate utilization patterns.

Hydrodynamic theory of liquid slippage on a solid substrate near a moving contact line

In this Letter a hydrodynamic theory of liquid slippage on a solid substrate near a moving contact line is proposed. A family of spatially varying slip lengths in the Navier slip law recovers the results of past formulations for slip in continuum theories and molecular dynamics simulations and is consistent with well-established experimental observations of complete wetting. This formulation gives a general approach for continuum hydrodynamic theories. New fluid flow behaviors are also predicted yet to be seen in experiment. © 2013 American Physical Society.

Effect of Coulomb blockade, gold resistance, and thermal expansion on the electrical resistance of ultrathin gold films

Controlling the electrical resistance of granular thin films is of great importance for many applications, yet a full understanding of electron transport in such films remains a major challenge. We have studied experimentally and by model calculations the temperature dependence of the electrical resistance of ultrathin gold films at temperatures between 2 K and 300 K. Using sputter deposition, the film morphology was varied from a discontinuous film of weakly coupled meandering islands to a continuous film of strongly coupled coalesced islands. In the weak-coupling regime, we compare the regular island array model, the cotunneling model, and the conduction percolation model with our experimental data. We show that the tunnel barriers and the Coulomb blockade energies are important at low temperatures and that the thermal expansion of the substrate and the island resistance affect the resistance at high temperatures. At low temperatures our experimental data show evidence for a transition from electron cotunneling to sequential tunneling but the data can also be interpreted in terms of conduction percolation. The resistivity and temperature coefficient of resistance of the meandering gold islands are found to resemble those of gold nanowires. We derive a simple expression for the temperature at which the resistance changes from non-metal-like behavior into metal-like behavior. In the case of strong island coupling, the total resistance is solely determined by the Ohmic island resistance.

Increased size selectivity of Si quantum dots on SiC at low substrate temperatures : an ion-assisted self-organization approach

A simple, effective, and innovative approach based on ion-assisted self-organization is proposed to synthesize size-selected Si quantum dots (QDs) on SiC substrates at low substrate temperatures. Using hybrid numerical simulations, the formation of Si QDs through a self-organization approach is investigated by taking into account two distinct cases of Si QD formation using the ionization energy approximation theory, which considers ionized in-fluxes containing Si3+ and Si1+ ions in the presence of a microscopic nonuniform electric field induced by a variable surface bias. The results show that the highest percentage of the surface coverage by 1 and 2 nm size-selected QDs was achieved using a bias of -20 V and ions in the lowest charge state, namely, Si1+ ions in a low substrate temperature range (227-327 °C). As low substrate temperatures (≤500 °C) are desirable from a technological point of view, because (i) low-temperature deposition techniques are compatible with current thin-film Si-based solar cell fabrication and (ii) high processing temperatures can frequently cause damage to other components in electronic devices and destroy the tandem structure of Si QD-based third-generation solar cells, our results are highly relevant to the development of the third-generation all-Si tandem photovoltaic solar cells.

Distribution of ion-current density on substrate between carbon nanotubes grown from low-temperature plasma

Using Monte Carlo simulation technique, we have calculated the distribution of ion current extracted from low-temperature plasmas and deposited onto the substrate covered with a nanotube array. We have shown that a free-standing carbon nanotube is enclosed in a circular bead of the ion current, whereas in square and hexagonal nanotube patterns, the ion current is mainly concentrated along the lines connecting the nearest nanotubes. In a very dense array (with the distance between nanotubes/nanotube-height ratio less than 0.05), the ions do not penetrate to the substrate surface and deposit on side surfaces of the nanotubes.

Carbon saturation of arrays of Ni catalyst nanoparticles of different size and pattern uniformity on a silicon substrate

The kinetics of saturation of Ni catalyst nanoparticle patterns of the three different degrees of order, used as a model for the growth of carbon nanotips on Si, is investigated numerically using a complex model that involves surface diffusion and ion motion equations. It is revealed that Ni catalyst patterns of different degrees of order, with Ni nanoparticle sizes up to 12.5 nm, exhibit different kinetics of saturation with carbon on the Si surface. It is shown that in the cases examined (surface coverage in the range of 1-50%, highly disordered Ni patterns) the relative pattern saturation factor calculated as the ratio of average incubation times for the processes conducted in the neutral and ionized gas environments reaches 14 and 3.4 for Ni nanoparticles of 2.5 and 12.5 nm, respectively. In the highly ordered Ni patterns, the relative pattern saturation factor reaches 3 for nanoparticles of 2.5 nm and 2.1 for nanoparticles of 12.5 nm. Thus, more simultaneous saturation of Ni catalyst nanoparticles of sizes in the range up to 12.5 nm, deposited on the Si substrate, can be achieved in the low-temperature plasma environment than with the neutral gas-based process.

Nanoparticle manipulation in the near-substrate areas of low-temperature, high-density rf plasmas

Manipulation of a single nanoparticle in the near-substrate areas of high-density plasmas of low-temperature glow discharges is studied. It is shown that the nanoparticles can be efficiently manipulated by the thermophoretic force controlled by external heating of the substrate stage. Particle deposition onto or repulsion from nanostructured carbon surfaces critically depends on the values of the neutral gas temperature gradient in the near-substrate areas, which is directly measured in situ in different heating regimes by originally developed temperature gradient probe. The measured values of the near-surface temperature gradient are used in the numerical model of nanoparticle dynamics in a variable-length presheath. Specific conditions enabling the nanoparticle to overcome the repulsive potential and deposit on the substrate during the discharge operation are investigated. The results are relevant to fabrication of various nanostructured films employing structural incorporation of the plasma-grown nanoparticles, in particular, to nanoparticle deposition in the plasma-enhanced chemical-vapor deposition of carbon nanostructures in hydrocarbon-based plasmas.