A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

Optical properties of ZnO thin film on gamma-LiAlO2 substrate grown by pulsed laser deposition

Optical properties were investigated of ZnO thin films grown oil (100) gamma-LiAlO2 (LAO) substrates by pulsed laser deposition method. C-axis oriented ZnO film was grown oil (100) LAO substrate at the substrate temperature of 550 degrees C. The transmittances of the films were over 85%. Peaks attributed to excitons were seen in the absorption spectra, indicating that the thin films have high crystallinity. Photoluminescence spectra were observed at room temperature; the peak at 550 urn is ascribed to oxygen vacancies in the ZnO films caused by the diffusion of Li from the substrate into the film during deposition. (c) 2005 Elsevier B.V. All rights reserved.

Structural, morphological and optical properties of ZnO thin films grown on gamma-LiAlO2 substrate by pulsed laser deposition

ZnO thin films were deposited on the substrates of (100) gamma-LiAlO2 at 400, 550 and 700 degrees C using pulsed laser deposition (PLD) with the fixed oxygen pressure of 20 Pa, respectively. When the substrate temperature is 400 degrees C, the grain size of the film is less than 1 mu m observed by Leitz microscope and measured by X-ray diffraction (XRD). As the substrate temperature increases to 550 degrees C, highly-preferred c-orientation and high-quality ZnO film can be attained. While the substrate temperature rises to 700 degrees C, more defects appears on the surface of film and the ZnO films become polycrystalline again possibly because more Li of the substrate diffused into the ZnO film at high substrate temperature. The photoluminescence (PL) spectra of ZnO films at room temperature show the blue emission peaks centered at 430 nm. We suggest that the blue emission corresponds to the electron transition from the level of interstitial Zn to the valence band. Meanwhile, the films grown on gamma-LiAlO2 (LAO) exhibit green emission centered at 540 nm, which seemed to be ascribed to excess zinc and/or oxygen vacancy in the ZnO films caused by diffusion of Li. from the substrates into the films during the deposition.

Deposition of ZnO thin films on (100) γ-LiAlO2 substrate

Optical properties for ZnO thin films grown on (100) γ-LiAlO2 (LAO) substrate by pulsed laser deposition method were investigated. The c-axis oriented ZnO films were grown on (100) γ-LiAlO2 substrates at the substrate temperature of 550 Celsius degrees. The transmittance of the films was over 85%. Peaks attributed to excitons were shown in absorption spectra, which indicated that thin films had high crystallinity. Photoluminescence spectra with the maximum peak at 540 nm were observed at room temperature, which seemed to be ascribed to oxygen vacancy in the ZnO films caused by diffusion of Li from the substrates into the films during the deposition.

Evolution of localized blobs of swirling or buoyant fluid with and without an ambient magnetic field.

We investigate the evolution of localized blobs of swirling or buoyant fluid in an infinite, inviscid, electrically conducting fluid. We consider the three cases of a strong imposed magnetic field, a weak imposed magnetic field, and no magnetic field. For a swirling blob in the absence of a magnetic field, we find, in line with others, that the blob bursts radially outward under the action of the centrifugal force, forming a thin annular vortex sheet. A simple model of this process predicts that the vortex sheet thins exponentially fast and that it moves radially outward with constant velocity. These predictions are verified by high-resolution numerical simulations. When an intense magnetic field is applied, this phenomenon is suppressed, with the energy and angular momentum of the blob now diffusing axially along the magnetic field lines, converting the blob into a columnar structure. For modest or weak magnetic fields, there are elements of both types of behavior, with the radial bursting dominating over axial diffusion for weak fields. However, even when the magnetic field is very weak, the flow structure is quite distinct to that of the nonmagnetic case. In particular, a small but finite magnetic field places a lower bound on the thickness of the annular vortex sheet and produces an annulus of counter-rotating fluid that surrounds the vortex core. The behavior of the buoyant blob is similar. In the absence of a magnetic field, it rapidly develops the mushroomlike shape of a thermal, with a thin vortex sheet at the top and sides of the mushroom. Again, a simple model of this process predicts that the vortex sheet at the top of the thermal thins exponentially fast and rises with constant velocity. These predictions are consistent with earlier numerical simulations. Curiously, however, it is shown that the net vertical momentum associated with the blob increases linearly in time, despite the fact that the vertical velocity at the front of the thermal is constant. As with the swirling blob, an imposed magnetic field inhibits the formation of a vortex sheet. A strong magnetic field completely suppresses the phenomenon, replacing it with an axial diffusion of momentum, while a weak magnetic field allows the sheet to form, but places a lower bound on its thickness. The magnetic field does not, however, change the net vertical momentum of the blob, which always increases linearly with time.

Application of electron beams in thermal processing of semiconductor materials and devices

Rapid and effective thermal processing methods using electron beams are described in this paper. Heating times ranging from a fraction of a second to several seconds and temperatures up to 1400°C are attainable. Applications such as the annealing of ion implanted material, both without significant dopant diffusion and with highly controlled diffusion of impurities, are described. The technique has been used successfully to activate source/drain regions for fine geometry NMOS transistors. It is shown that electron beams can produce localised heating of semiconductor substrates and a resolution of approximately 1 μm has been achieved. Electron beam heating has been applied to improving the crystalline quality of silicon-on sapphire used in CMOS device fabrication. Silicon layers with defect levels approaching bulk material have been obtained. Finally, the combination of isothermal and selective annealing is shown to have application in recrystallisation of polysilicon films on an insulating layer. The approach provides the opportunity of producing a silicon-on-insulator substrate with improved crystalline quality compared to silicon-on-sapphire at a potentially lower cost. It is suggested that rapid heating methods are expected to provide a real alternative to conventional furnace processing of semiconductor devices in the development of fabrication technology. © 1984 Benn electronics Publications Ltd, Luton.

RESIDUAL DEFECTS FOLLOWING RAPID ANNEALING OF BORON IMPLANTED SILICON WITH AND WITHOUT PREAMORPHISATION BY SILICON IMPLANTATION.

Implants of boron into silicon which has been made amorphous by silicon implantation have a shallower depth profile than the same implants into silicon. This results in higher activation and restricted diffusion of the B implants after annealing, and there are also significant differences in the microstructure after annealing compared with B implants into silicon. Rapid isothermal heating with an electron beam and furnace treatments are used to characterize the defect structure as a function of time and temperature. Defects are seen to influence the diffusion of non-substitutional boron.

An analysis of the cooperative mechano-sensitive feedback between intracellular signaling, focal adhesion development, and stress fiber contractility

Cells communicate with their external environment via focal adhesions and generate activation signals that in turn trigger the activity of the intracellular contractile machinery. These signals can be triggered by mechanical loading that gives rise to a cooperative feedback loop among signaling, focal adhesion formation, and cytoskeletal contractility, which in turn equilibrates with the applied mechanical loads. We devise a signaling model that couples stress fiber contractility and mechano-sensitive focal adhesion models to complete this above mentioned feedback loop. The signaling model is based on a biochemical pathway where IP3 molecules are generated when focal adhesions grow. These IP3 molecules diffuse through the cytosol leading to the opening of ion channels that disgorge Ca2+ from the endoplasmic reticulum leading to the activation of the actin/myosin contractile machinery. A simple numerical example is presented where a one-dimensional cell adhered to a rigid substrate is pulled at one end, and the evolution of the stress fiber activation signal, stress fiber concentrations, and focal adhesion distributions are investigated. We demonstrate that while it is sufficient to approximate the activation signal as spatially uniform due to the rapid diffusion of the IP3 through the cytosol, the level of the activation signal is sensitive to the rate of application of the mechanical loads. This suggests that ad hoc signaling models may not be able to capture the mechanical response of cells to a wide range of mechanical loading events. © 2011 American Society of Mechanical Engineers.

A rhetorical approach to IT diffusion: reconceptualizing the ideology-framing relationship in computerization movements

In this paper we propose rhetoric as a valuable yet underdeveloped alternative paradigm for examining IT diffusion. Building on recent developments of computerization movements theory, our rhetorical approach proposes that two central elements of the theory, framing and ideology, rather than being treated as separate can be usefully integrated. We suggest that IT diffusion can be usefully explored through examining the interrelationship of the deep structures underlying ideology and the type and sequence of rhetorical claims underpinning actors’ framing strategies. Our theoretical developments also allow us to better understand competing discourses influencing the diffusion process. These discourses reflect the ideologies and shape the framing strategies of actors in the broader field context. We illuminate our theoretical approach by drawing on the history of the diffusion of free and open source software.

A new observation on large scale disintegration of eggs of Indian major carps during development

While carp themselves may destroy carp eggs by sucking from outside the happa, destruction in the absence of carp also occurs. During observations of the effect of electrical field on developing carp eggs it was seen that the egg membrane ruptured prematurely under certain potential differences. In natural waters and mud potential differences occur in the presence of high concentrations of certain salts compared with the internal concentration. Temperature and pH fluctuations following strong sun or heavy rain accelerates the diffusion of ions along side water through the egg membrane. The egg cells expand rapidly and the embryos were rapidly dropped.

The internal structure of igniting turbulent sprays as revealed by complex chemistry DNS

A parametric study of spark ignition in a uniform monodisperse turbulent spray is performed with complex chemistry three-dimensional Direct Numerical Simulations in order to improve the understanding of the structure of the ignition kernel. The heat produced by the kernel increases with the amount of fuel evaporated inside the spark volume. Moreover, the heat sink by evaporation is initially higher than the heat release and can have a negative effect on ignition. With the sprays investigated, heat release occurs over a large range of mixture fractions, being high within the nominal flammability limits and finite but low below the lean flammability limit. The burning of very lean regions is attributed to the diffusion of heat and species from regions of high heat release, and from the spark, to lean regions. Two modes of spray ignition are reported. With a relatively dilute spray, nominally flammable material exists only near the droplets. Reaction zones are created locally near the droplets and have a non-premixed character. They spread from droplet to droplet through a very lean interdroplet spacing. With a dense spray, the hot spark region is rich due to substantial evaporation but the cold region remains lean. In between, a large surface of flammable material is generated by evaporation. Ignition occurs there and a large reaction zone propagates from the rich burned region to the cold lean region. This flame is wrinkled due to the stratified mixture fraction field and evaporative cooling. In the dilute spray, the reaction front curvature pdf contains high values associated with single droplet combustion, while in the dense spray, the curvature is lower and closer to the curvature associated with gaseous fuel ignition kernels. © 2011 The Combustion Institute.

Effects of preferential transport in turbulent bluff-body-stabilized lean premixed CH 4/air flames

Preferential species diffusion is known to have important effects on local flame structure in turbulent premixed flames, and differential diffusion of heat and mass can have significant effects on both local flame structure and global flame parameters, such as turbulent flame speed. However, models for turbulent premixed combustion normally assume that atomic mass fractions are conserved from reactants to fully burnt products. Experiments reported here indicate that this basic assumption may be incorrect for an important class of turbulent flames. Measurements of major species and temperature in the near field of turbulent, bluff-body stabilized, lean premixed methane-air flames (Le=0.98) reveal significant departures from expected conditional mean compositional structure in the combustion products as well as within the flame. Net increases exceeding 10% in the equivalence ratio and the carbon-to-hydrogen atom ratio are observed across the turbulent flame brush. Corresponding measurements across an unstrained laminar flame at similar equivalence ratio are in close agreement with calculations performed using Chemkin with the GRI 3.0 mechanism and multi-component transport, confirming accuracy of experimental techniques. Results suggest that the large effects observed in the turbulent bluff-body burner are cause by preferential transport of H 2 and H 2O through the preheat zone ahead of CO 2 and CO, followed by convective transport downstream and away from the local flame brush. This preferential transport effect increases with increasing velocity of reactants past the bluff body and is apparently amplified by the presence of a strong recirculation zone where excess CO 2 is accumulated. © 2011 The Combustion Institute.

Steady-state modelling of the universal exhaust gas oxygen (UEGO) sensor

The universal exhaust gas oxygen (UEGO) sensor is a well-established device which was developed for the measurement of relative air fuel ratio in internal combustion engines. There is, however, little information available which allows for the prediction of the UEGO's behaviour when exposed to arbitrary gas mixtures, pressures and temperatures. Here we present a steady-state model for the sensor, based on a solution of the Stefan-Maxwell equation, and which includes a momentum balance. The response of the sensor is dominated by a diffusion barrier, which controls the rate of diffusion of gas species between the exhaust and a cavity. Determination of the diffusion barrier characteristics, especially the mean pore size, porosity and tortuosity, is essential for the purposes of modelling, and a measurement technique based on identification of the sensor pressure giving zero temperature sensitivity is shown to be a convenient method of achieving this. The model, suitably calibrated, is shown to make good predictions of sensor behaviour for large variations of pressure, temperature and gas composition. © 2012 IOP Publishing Ltd.

The behaviour of laminar stratified methane/air flames in counterflow

Flames propagating through a mixture with a gradient of equivalence ratio have been previously demonstrated to travel faster or slower than their equivalent premixed flames. The present study aims to numerically investigate the response of strained laminar methane-air flames to such gradients. The flames are simulated in a counterflow configuration where a premixed reactant stream at equivalence ratio φR opposes a hot equilibrium stream at equivalence ratio φP. Premixed and stratified flames are compared with respect to the equivalence ratio φ* and the corresponding gradient ∇φ* at the point of peak heat release rate, for three strain rates, a=50, 300 and 500s-1 and a range of φ*. The effect of different stratification levels is also investigated by varying the ratio of φP to φR, Θ. Results indicate that, as long as flames stabilize within the diffusion layer and Θ>1, increased heat release rate Q is seen throughout the progress variable space in comparison to the premixed state. In contrast, an attenuation of heat release rate is seen for Θ<1. The enhancement (or attenuation) of heat release varies monotonically with Θ. The effect of stratification on flame behavior becomes more pronounced as the strain rate increases. The present study reveals the mechanisms for the propagation of quasi-steady stratified flames under lean and rich conditions: stratified flames are primarily dominated by the diffusion of heat under lean conditions, and diffusion of H2 under rich conditions. Thanks to species and thermal support, stratified flames continue to burn beyond the premixed lean and rich flammability limits. Further investigation on the unsteady response of flames to the fluctuating equivalence ratio implies that the steady results represent the unsteady response well, as long as φ* and ∇φ* are similar in both steady and unsteady cases. © 2013 The Combustion Institute.