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.

Gas-phase reactions of aryl radicals with 2-butyne : an experimental and theoretical investigation employing the N-methyl-pyridinium-4-yl radical cation

Aromatic radicals form in a variety of reacting gas-phase systems, where their molecular weight growth reactions with unsaturated hydrocarbons are of considerable importance. We have investigated the ion-molecule reaction of the aromatic distonic N-methyl-pyridinium-4-yl (NMP) radical cation with 2-butyne (CH3C CCH3) using ion trap mass spectrometry. Comparison is made to high-level ab initio energy surfaces for the reaction of NMP and for the neutral phenyl radical system. The NMP radical cation reacts rapidly with 2-butyne at ambient temperature, due to the apparent absence of any barrier. The activated vinyl radical adduct predominantly dissociates via loss of a H atom, with lesser amounts of CH3 loss. High-resolution Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry allows us to identify small quantities of the collisionally deactivated reaction adduct. Statistical reaction rate theory calculations (master equation/RRKM theory) on the NMP + 2-butyne system support our experimental findings, and indicate a mechanism that predominantly involves an allylic resonance-stabilized radical formed via H atom shuttling between the aromatic ring and the C-4 side-chain, followed by cyclization and/or low-energy H atom beta-scission reactions. A similar mechanism is demonstrated for the neutral phenyl radical (Ph center dot)+2-butyne reaction, forming products that include 3-methylindene. The collisionally deactivated reaction adduct is predicted to be quenched in the form of a resonance-stabilized methylphenylallyl radical. Experiments using a 2,5-dichloro substituted methyl-pyridiniumyl radical cation revealed that in this case CH3 loss from the 2-butyne adduct is favoured over H atom loss, verifying the key role of ortho H atoms, and the shuttling mechanism, in the reactions of aromatic radicals with alkynes. As well as being useful phenyl radical analogues, pyridiniumyl radical cations may form in the ionosphere of Titan, where they could undergo rapid molecular weight growth reactions to yield polycyclic aromatic nitrogen hydrocarbons (PANHs).

Soil moisture monitoring at the field scale using neutron probe

Measurement of the moisture variation in soils is required for geotechnical design and research because soil properties and behavior can vary as moisture content changes. The neutron probe, which was developed more than 40 years ago, is commonly used to monitor soil moisture variation in the field. This study reports a full-scale field monitoring of soil moisture using a neutron moisture probe for a period of more than 2 years in the Melbourne (Australia) region. On the basis of soil types available in the Melbourne region, 23 sites were chosen for moisture monitoring down to a depth of 1500 mm. The field calibration method was used to develop correlations relating the volumetric moisture content and neutron counts. Observed results showed that the deepest “wetting front” during the wet season was limited to the top 800 to 1000 mm of soil whilst the top soil layer down to about 550mmresponded almost immediately to the rainfall events. At greater depths (550 to 800mmand below 800 mm), the moisture variations were relatively low and displayed predominantly periodic fluctuations. This periodic nature was captured with Fourier analysis to develop a cyclic moisture model on the basis of an analytical solution of a one-dimensional moisture flow equation for homogeneous soils. It is argued that the model developed can be used to predict the soil moisture variations as applicable to buried structures such as pipes.

Validation tool for traction force microscopy

Traction force microscopy (TFM) is commonly used to estimate cells’ traction forces from the deformation that they cause on their substrate. The accuracy of TFM highly depends on the computational methods used to measure the deformation of the substrate and estimate the forces, and also on the specifics of the experimental set-up. Computer simulations can be used to evaluate the effect of both the computational methods and the experimental set-up without the need to perform numerous experiments. Here, we present one such TFM simulator that addresses several limitations of the existing ones. As a proof of principle, we recreate a TFM experimental set-up, and apply a classic 2D TFM algorithm to recover the forces. In summary, our simulator provides a valuable tool to study the performance, refine experimentally, and guide the extraction of biological conclusions from TFM experiments.

Preparation and characterization of composites from starch with sugarcane bagasse nanofibres

This paper reports on the results of using unbleached sugar cane bagasse nanofibres (average diameter 26.5 nm; aspect ratio 247 assuming a dry fibre density of 1,500 kg/m3) to improve the physico-chemical properties of starch-based films. The addition of bagasse nanofibres (2.5 to 20 wt%) to modified potato starch (i.e. soluble starch) reduced the moisture uptake by up to 17 % at 58 % relative humidity. The film’s tensile strength and Young’s modulus increased by up to 100 % (3.1 to 6.2 MPa) and 300 % (66.3 to 198.3 MPa) respectively with 10 and 20 wt% fibre addition. However, the strain at yield dropped by 50 % for the film containing 10 wt% fibre. Models for composite materials were used to account for the strong interactions between the nanofibres and the starch matrix. The storage and loss moduli as well as the glass transition temperature (Tg) obtained from dynamic mechanical thermal analysis, were increased with the starch-nanofibre films indicating decreased starch chain mobility due to the interacting effect of the nanofibres. Evidence of the existence of strong interactions between the starch matrix and the nanofibres was revealed from detailed Fourier transform infra-red and scanning electron microscopic evaluation.

Polarized infrared absorption spectrum of matrix-isolated methylperoxyl radicals, CH3OO (X)over-tilde 2A ''

We have used a tandem pair of supersonic nozzles to produce clean samples of CH3OO radicals in cryogenic matrices. One hyperthermal nozzle decomposes azomethane (CH3NNCH3) to generate intense pulses of CH3 radicals, While the second nozzle alternately fires a burst Of O-2/Ar at the 20 K matrix. The CH3/O-2/20 K argon radical sandwich acts to produce target methylperoxyl radicals: CH3 + O-2 --> CH3OO. The absorption spectra of the radicals are monitored with a Fourier transform infrared spectrometer. We report 10 of the 12 fundamental infrared bands of the methylperoxyl radical CH3OO, (X) over tilde (2)A", in an argon matrix at 20 K. The experimental frequencies (cm(-1)) and polarizations follow: the a' modes are 3032, 2957, 1448, 1410, 1180, 1109, 90, 492, while the a" modes are 3024 and 1434. We cannot detect the asymmetric CH3 rocking mode, nu(11), nor the torsion, nu(12). The infrared spectra of (CH3OO)-O-18-O-18, (CH3OO)-C-13, and CD3OO have been measured as well in order to determine the isotopic shifts. The experimental frequencies, {nu}, for the methylperoxyl radicals are compared to harmonic frequencies, {omega}, resulting from a UB3LYP/6-311G(d,p) electronic structure calculation. Linear dichroism spectra were measured with photooriented radical samples in order to establish the experimental polarizations of most vibrational bands. The methylperoxyl radical matrix frequencies listed above are within +/-2% of the gas-phase vibrational frequencies. A final set of vibrational frequencies for the H radical, are recommended. See also http://ellison.colorado.edu/methylperoxyl.

Insight into morphology and structure of different particle sized kaolinites with same origin

The particle size, morphology, crystallinity order and structural defects of four kaolinite samples are characterized by the techniques including particle size analysis, scanning electron microscopy (SEM), X-ray diffraction (XRD), Raman spectroscopy, Fourier transform infrared spectroscopy (FTIR) and magic angle spinning nuclear magnetic resonance spectroscopy (MAS NMR). The particle size of four kaolinite samples gradually increases. Four samples all belong to the ordered kaolinite and show a decrease in structural order with the increase of kaolinite particle size. The changes of structural defect are proved by the increase of the band splitting in Raman spectroscopy, the decrease of the intensity of absorption bands in infrared spectroscopy, and the decrease of equivalent silicon atom and the increase of nonequivalent aluminum atom in MAS NMR spectroscopy. The differences in morphology and structural defect are attributed to the broken bonds of Al–O–Si, Al–O–Al and Si–O–Si and the Al substitution for Si in tetrahedral sheets.

Intercalation of dodecylamine into kaolinite and its layering structure investigated by molecular dynamics simulation

Dodecylamine was successfully intercalated into the layer space of kaolinite by utilizing the methanol treated kaolinite–dimethyl sulfoxide (DMSO) intercalation complex as an intermediate. The basal spacing of kaolinite, measured by X-ray diffraction (XRD), increased from 0.72 nm to 4.29 nm after the intercalation of dodecylamine. Also, the significant variation observed in the Fourier Transform Infrared Spectroscopy (FTIR) spectra of kaolinite when intercalated with dodecylamine verified the feasibility of intercalation of dodecylamine into kaolinite. Isothermal-isobaric (NPT) molecular dynamics simulation with the use of Dreiding force field was performed to probe into the layering behavior and structure of nanoconfined dodecylamine in the kaolinite gallery. The concentration profiles of the nitrogen atom, methyl group and methylene group of intercalated dodecylamine molecules in the direction perpendicular to the kaolinite basal surface indicated that the alkyl chains within the interlayer space of kaolinite exhibited an obvious layering structure. However, the unified bilayer, pseudo-trilayer, or paraffin-type arrangements of alkyl chains deduced based on their chain length combined with the measured basal spacing of organoclays were not found in this study. The alkyl chains aggregated to a mixture of ordered paraffin-type-like structure and disordered gauche conformation in the middle interlayer space of kaolinite, and some alkyl chains arranged in two bilayer structures, in which one was close to the silica tetrahedron surface, and the other was close to the alumina octahedron surface with their alkyl chains parallel to the kaolinite basal surface.

Reinforced insulation properties of epoxy resin/SiO2 nanocomposites by atmospheric pressure plasma modification

Nanocomposite dielectrics hold a promising future for the next generation of insulation materials because of their excellent physical, chemical, and dielectric properties. In the presented study, we investigate the use of plasma processing technology to further enhance the dielectric performance of epoxy resin/SiO2 nanocomposite materials. The SiO2 nanoparticles are treated with atmospheric-pressure non-equilibrium plasma prior to being added into the epoxy resin host. Fourier transform infrared spectroscopy (FTIR) results reveal the effects of the plasma process on the surface functional groups of the treated nanoparticles. Scanning electron microscopy (SEM) results show that the plasma treatment appreciably improves the dispersion uniformity of nanoparticles in the host polymer. With respect to insulation performance, the epoxy/plasma-treated SiO2 specimen shows a 29% longer endurance time than the epoxy/untreated SiO2 nanocomposite under electrical aging. The Weibull plots of the dielectric breakdown field intensity suggest that the breakdown strength of the nanocomposite with the plasma pre-treatment on the nanoparticles is improved by 23.3%.

Plasma functionalization of SiO2 nanoparticles for the synthesis of polymer nano-dielectrics

In this study, we improve the insulation performance of polymeric nano-dielectrics by using plasma pre-treatment on the filled nanoparticles. Non-equilibrium atmospheric-pressure plasma is employed to modify a commercial type of silane-coated SiO2 nanoparticles. The treated nanoparticles and the synthesized epoxy-based nanocomposites are characterized using scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FTIR), and X-ray photoelectron spectroscopy (XPS). The plasma-treated SiO2 nanoparticles can disperse uniformly and form strong covalent bonds with the molecules of the polymer matrix. Moreover, the electrical insulation properties of the synthesized nanocomposites are investigated. Results show that the nanocomposites with plasma-treated SiO2 nanoparticles obtain improved dielectric breakdown strength and extended endurance under intense electrical ageing process.

Si quantum dots embedded in an amorphous SiC matrix : nanophase control by non-equilibrium plasma hydrogenation

Nanophase nc-Si/a-SiC films that contain Si quantum dots (QDs) embedded in an amorphous SiC matrix were deposited on single-crystal silicon substrates using inductively coupled plasma-assisted chemical vapor deposition from the reactive silane and methane precursor gases diluted with hydrogen at a substrate temperature of 200 °C. The effect of the hydrogen dilution ratio X (X is defined as the flow rate ratio of hydrogen-to-silane plus methane gases), ranging from 0 to 10.0, on the morphological, structural, and compositional properties of the deposited films, is extensively and systematically studied by scanning electron microscopy, high-resolution transmission electron microscopy, X-ray diffraction, Raman spectroscopy, Fourier-transform infrared absorption spectroscopy, and X-ray photoelectron spectroscopy. Effective nanophase segregation at a low hydrogen dilution ratio of 4.0 leads to the formation of highly uniform Si QDs embedded in the amorphous SiC matrix. It is also shown that with the increase of X, the crystallinity degree and the crystallite size increase while the carbon content and the growth rate decrease. The obtained experimental results are explained in terms of the effect of hydrogen dilution on the nucleation and growth processes of the Si QDs in the high-density plasmas. These results are highly relevant to the development of next-generation photovoltaic solar cells, light-emitting diodes, thin-film transistors, and other applications.

Low-temperature plasma-assisted growth of optically transparent, highly oriented nanocrystalline AlN

Optically transparent, highly oriented nanocrystalline AlN(002) films have been synthesized using a hybrid plasma enhanced chemical vapor deposition and plasma-assisted radio frequency (rf) magnetron sputtering process in reactive Ar+ N2 and Ar+ N2 + H2 gas mixtures at a low Si(111)/glass substrate temperature of 350 °C. The process conditions, such as the sputtering pressure, rf power, substrate temperature, and N2 concentration were optimized to achieve the desired structural, compositional, and optical characteristics. X-ray diffractometry reveals the formation of highly c -oriented AlN films at a sputtering pressure of 0.8 Pa. Field emission scanning electron microscopy suggests the uniform distribution of AlN grains over large surface areas and also the existence of highly oriented in the (002) direction columnar structures of a typical length ∼100-500 nm with an aspect ratio of ∼7-15. X-ray photoelectron and energy dispersive x-ray spectroscopy suggest that films deposited at a rf power of 400 W feature a chemically pure and near stoichiometric AlN. The bonding states of the AlN films have been confirmed by Raman and Fourier transform infrared spectroscopy showing strong E2 (high) and E1 transverse optical phonon modes. Hydrogenated AlN films feature an excellent optical transmittance of ∼80% in the visible region of the spectrum, promising for advanced optical applications.

Rapid, low-temperature synthesis of nc-Si in high-density, non-equilibrium plasmas : enabling nanocrystallinity at very low hydrogen dilution

Nanocrystalline silicon thin films were deposited on single-crystal silicon and glass substrates simultaneously by inductively coupled plasma-assisted chemical vapor deposition from the reactive silane reactant gas diluted with hydrogen at a substrate temperature of 200 °C. The effect of hydrogen dilution ratio X (X is defined as the flow rate ratio of hydrogen to silane gas), ranging from 1 to 20, on the structural and optical properties of the deposited films, is extensively investigated by Raman spectroscopy, X-ray diffraction, Fourier transform infrared absorption spectroscopy, UV/VIS spectroscopy, and scanning electron microscopy. Our experimental results reveal that, with the increase of the hydrogen dilution ratio X, the deposition rate Rd and hydrogen content CH are reduced while the crystalline fraction Fc, mean grain size δ and optical bandgap ETauc are increased. In comparison with other plasma enhanced chemical vapor deposition methods of nanocrystalline silicon films where a very high hydrogen dilution ratio X is routinely required (e.g. X > 16), we have achieved nanocrystalline silicon films at a very low hydrogen dilution ratio of 1, featuring a high deposition rate of 1.57 nm/s, a high crystalline fraction of 67.1%, a very low hydrogen content of 4.4 at.%, an optical bandgap of 1.89 eV, and an almost vertically aligned columnar structure with a mean grain size of approximately 19 nm. We have also shown that a sufficient amount of atomic hydrogen on the growth surface essential for the formation of nanocrystalline silicon is obtained through highly-effective dissociation of silane and hydrogen molecules in the high-density inductively coupled plasmas. © 2009 The Royal Society of Chemistry.

Structural evolution of nanocrystalline silicon thin films synthesized in high-density, low-temperature reactive plasmas

Silicon thin films with a variable content of nanocrystalline phase were deposited on single-crystal silicon and glass substrates by inductively coupled plasma-assisted chemical vapor deposition using a silane precursor without any hydrogen dilution in the low substrate temperature range from 100 to 300 °C. The structural and optical properties of the deposited films are systematically investigated by Raman spectroscopy, x-ray diffraction, Fourier transform infrared absorption spectroscopy, UV/vis spectroscopy, scanning electron microscopy and high-resolution transmission electron microscopy. It is shown that the structure of the silicon thin films evolves from the purely amorphous phase to the nanocrystalline phase when the substrate temperature is increased from 100 to 150 °C. It is found that the variations of the crystalline fraction fc, bonded hydrogen content CH, optical bandgap ETauc, film microstructure and growth rate Rd are closely related to the substrate temperature. In particular, at a substrate temperature of 300 °C, the nanocrystalline Si thin films of our interest feature a high growth rate of 1.63nms-1, a low hydrogen content of 4.0at.%, a high crystalline fraction of 69.1%, a low optical bandgap of 1.55eV and an almost vertically aligned columnar structure with a mean grain size of approximately 10nm. It is also shown that the low-temperature synthesis of nanocrystalline Si thin films without any hydrogen dilution is attributed to the outstanding dissociation ability of the high-density inductively coupled plasmas and effective plasma-surface interactions during the growth process. Our results offer a highly effective yet simple and environmentally friendly technique to synthesize high-quality nanocrystalline Si films, vitally needed for the development of new-generation solar cells and other emerging nanotechnologies.

High-rate, low-temperature synthesis of composition controlled hydrogenated amorphous silicon carbide films in low-frequency inductively coupled plasmas

It is commonly believed that in order to synthesize high-quality hydrogenated amorphous silicon carbide (a-Si1-xCx : H) films at competitive deposition rates it is necessary to operate plasma discharges at high power regimes and with heavy hydrogen dilution. Here we report on the fabrication of hydrogenated amorphous silicon carbide films with different carbon contents x (ranging from 0.09 to 0.71) at high deposition rates using inductively coupled plasma (ICP) chemical vapour deposition with no hydrogen dilution and at relatively low power densities (∼0.025 W cm -3) as compared with existing reports. The film growth rate R d peaks at x = 0.09 and x = 0.71, and equals 18 nm min-1 and 17 nm min-1, respectively, which is higher than other existing reports on the fabrication of a-Si1-xCx : H films. The extra carbon atoms for carbon-rich a-Si1-xCx : H samples are incorporated via diamond-like sp3 C-C bonding as deduced by Fourier transform infrared absorption and Raman spectroscopy analyses. The specimens feature a large optical band gap, with the maximum of 3.74 eV obtained at x = 0.71. All the a-Si1-xCx : H samples exhibit low-temperature (77 K) photoluminescence (PL), whereas only the carbon-rich a-Si1-xCx : H samples (x ≥ 0.55) exhibit room-temperature (300 K) PL. Such behaviour is explained by the static disorder model. High film quality in our work can be attributed to the high efficiency of the custom-designed ICP reactor to create reactive radical species required for the film growth. This technique can be used for a broader range of material systems where precise compositional control is required. © 2008 IOP Publishing Ltd.