Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Elevated body core temperature in medico-legal investigation of violent death

Pathologically elevated body core temperature, measured at the death scene, is an important finding in medico-legal investigation of violent deaths. An abnormally high rectal temperature at any death scene may point to an underlying pathology, the influence of certain drugs or a hidden cerebral traumatism, and death by suffocation which would remain undetected without further medico-legal investigations. Furthermore, hyperthermia and fever, if unrecognized, may result in an erroneous forensic estimation of time since death in the early postmortem period by the "Henssge method." By a retrospective study of 744 cases, the authors demonstrate that hyperthermia is a finding with an incidence of 10% of all cases of violent death. The main causes are: influence of drugs, malignant tumors, cerebral hypoxia as a result of suffocation, infections, and systemic inflammatory disorders. As a consequence it must be stated, that hyperthermia must be excluded in every medico-legal death scene investigation by a correct measurement of body core temperature and a comparison between the cooling rate of the body and the behavior of early postmortem changes, notably livor and rigor mortis.

Geochemistry of minerals in medium-grained diabase from ODP Holes 137-504B and 140-504B

Samples recovered from Hole 504B during Leg 140 include a number of medium-grained, holocrystalline diabases that appear to represent the cores of thick dikes. The plagioclase and pyroxene in these samples occur in a variety of crystal morphologies. Plagioclase occurs as phenocrysts, microphenocrysts, elongate crystals, skeletal crystals, and branching radial clusters. Pyroxene occurs as phenocrysts, microphenocrysts, ophitic crystals, and poikilitic crystals. Plagioclase compositions became progressively poorer in anorthite and MgO and progressively richer in FeO as crystallization proceeded, while the average grain volume decreased and the aspect ratio of individual grains increased. Pyroxene compositions are largely independent of crystal morphology. The diabase dikes recovered from Hole 504B during Leg 140 appear to have crystallized in situ. Crystal compositions and morphologies are consistent with a rapid cooling rate and solidification times for individual dikes on the order of hours or days. The crystallization rate and nucleation rate of plagioclase lagged behind the cooling rate so that the degree of undercooling progressively increased as crystallization proceeded. Plagioclase crystal morphologies indicate much greater degrees of supersaturation than do pyroxene or olivine crystal morphologies. The 504B diabase magmas appear to have been emplaced with abundant preexisting pyroxene and olivine nuclei, but with few preexisting plagioclase nuclei. The suppression of plagioclase nucleation and crystallization relative to that of pyroxene and olivine could provide a mechanism by which the actual fractionation assemblage is more pyroxene-rich and plagioclase-poor than that predicted from thermodynamic models, or that observed in isothermal crystallization experiments.

(Table 1) Geochemistry of sector-zoned crystal at DSDP Hole 45-395A

Variations in crystal morphologies in pillow basalts and probable sheet flows sampled from the region of the East Pacific Rise drilled during Leg 54 are related both to differences in composition and to an extreme range of cooling rate experienced upon extrusion. The basalts range in composition from olivine-rich tholeiites to tholeiitic ferrobasalts, and include some more alkaline basalts. The kinetics of crystal growth in some samples appears to have been influenced by the amount of initial superheating (or supercooling) of the magma, or possibly by differential retention of volatiles. Olivine in quartznormative ferrobasalts apparently formed metastably at high undercooling. Despite these effects, reliable petrographic criteria are established to distinguish the principal rock types described regardless of the crystallinity and grain size. Microphenocrysts formed prior to pillow formation correspond closely to mineral assemblages inferred from normative plots and variation diagrams to control crystal fractionation at various stages. The details of spherulitic and dendritic growth also provide some clues about composition. Petrographic evidence for magma mixing is scant. Only some Siqueiros fracture zone basalts contain zoned plagioclase phenocrysts with glass inclusions similar to those used to infer mixing among Mid-Atlantic Ridge basalts. All basalts from the summit and flanks of the East Pacific Rise are aphyric. One possible petrographic consequence of mixing between olivine tholeiites and ferrobasalts - formation of clinopyroxene phenocrysts - is not evident in any fracture zone or Rise crest basalt. Highly evolved ferrobasalts with liquidus low-Ca clinopyroxene have not been sampled, nor does textural evidence indicate that any basalts sampled are hybrid compositions between such magmas and less fractionated compositions. Evidently the sampled ferrobasalts are close to the most evolved compositions that occur in any abundance on this portion of the East Pacific Rise.

(Table 2) Chemical composition of basalts from the DSDP Leg 74 Walvis Ridge transect

Five sites were drilled along a transect of the Walvis Ridge. The basement rocks range in age from 69 to 71 m.y., and the deeper sites are slightly younger, in agreement with the sea-floor-spreading magnetic lineations. Geophysical and petrological evidence indicates that the Walvis Ridge was formed at a mid-ocean ridge at anomalously shallow elevations. The basement complex, associated with the relatively smooth acoustic basement in the area, consists of pillowed basalt and massive flows alternating with nannofossil chalk and limestone that contain a significant volcanogenic component. Basalts are quartz tholeiites at the ridge crest and olivine tholeiites downslope. The sediment sections are dominated by carbonate oozes and chalks with volcanogenic material common in the lower parts of the sediment columns. The volcanogenic sediments probably were derived from sources on the Walvis Ridge. Paleodepth estimates based on the benthic fauna are consistent with a normal crustal-cooling rate of subsidence of the Walvis Ridge. The shoalest site in the transect sank below sea level in the late Paleocene, and benthic fauna suggest a rapid sea-level lowering in the mid-Oligocene. Average accumulation rates during the Cenozoic indicate three peaks in the rate of supply of carbonate to the sea floor, that is, early Pliocene, late middle Miocene, and late Paleocene to early Eocene. Carbonate accumulation rates for the rest of the Cenozoic averaged 1 g/cm**2/kyr. Dissolution had a marked effect on sediment accumulation in the deeper sites, particularly during the late Miocene, Oligocene, and middle to late Eocene. Changes in the rates of accumulation as a function of depth demonstrate that the upper part of the water column had a greater degree of undersaturation with respect to carbonate during times of high productivity. Even when the calcium carbonate compensation depth (CCD) was below 4400 m, a significant amount of carbonate was dissolved at the shallower sites. The flora and fauna of the Walvis Ridge are temperate in nature. Warmer-water faunas are found in the uppermost Maastrichtian and lower Eocene sediments, with cooler-water faunas present in the lower Paleocene, Oligocene, and middle Miocene. The boreal elements of the lower Pliocene are replaced by more temperate forms in the middle Pliocene. The Cretaceous-Tertiary boundary was recovered in four sites drilled, with the sediments containing well-preserved nannofossils but poorly preserved foraminifera.

Band gap control via tuning of inversion degree in CdIn2S4 spinel

Based on theoretical arguments, we propose a possible route for controlling the band-gap in the promising photovoltaic material CdIn2S4. Our ab initio calculations show that the experimental degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approximately to the equilibrium value given by the minimum of the theoretical inversion free energy at a typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthesis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic band-gap of the solid, as shown here by screened hybrid functional calculations.

Glassy State and Cryopreservation of Mint Shoot Tips

Vitrification refers to the physical process by which a liquid supercools to very low tem- peratures and finally solidifies into a metastable glass, without undergoing crystallization at a practical cooling rate. Thus, vitrification is an effective freeze-avoidance mechanism and living tissue cryopreservation is, in most cases, relying on it. As a glass is exceedingly vis- cous and stops all chemical reactions that require molecular diffusion, its formation leads to metabolic inactivity and stability over time. To investigate glassy state in cryopreserved plant material, mint shoot tips were submitted to the different stages of a frequently used cryopreservation protocol (droplet-vitrification) and evaluated for water content reduction and sucrose content, as determined by ion chromatography, frozen water fraction and glass transitions occurrence by differential scanning calorimetry, and investigated by low-tempera- ture scanning electron microscopy, as a way to ascertain if their cellular content was vitri- fied. Results show how tissues at intermediate treatment steps develop ice crystals during liquid nitrogen cooling, while specimens whose treatment was completed become vitrified, with no evidence of ice formation. The agreement between calorimetric and microscopic observations was perfect. Besides finding a higher sucrose concentration in tissues at the more advanced protocol steps, this level was also higher in plants precultured at 25/21?C than in plants cultivated at 25?C.

Development of novel melt-processable biopolymer nanocomposites based on poly(L-lactic acid) and WS2 inorganic nanotubes.

The use of tungsten disulphide inorganic nanotubes (INT-WS2) offers the opportunity to produce novel and advanced biopolymer-based nanocomposite materials with excellent nanoparticle dispersion without the need for modifiers or surfactants via conventional melt blending. The study of the non-isothermal melt-crystallization kinetics provides a clear picture of the transformation of poly(L-lactic acid) (PLLA) molecules from the non-ordered to the ordered state. The overall crystallization rate, final crystallinity and subsequent melting behaviour of PLLA were controlled by both the incorporation of INT-WS2 and the variation of the cooling rate. In particular, it was shown that INT-WS2 exhibits much more prominent nucleation activity on the crystallization of PLLA than other specific nucleating agents or nano-sized fillers. These features may be advantageous for the enhancement of mechanical properties and process-ability of PLLA-based materials. PLLA/INT-WS2 nanocomposites can be employed as low cost biodegradable materials for many eco-friendly and medical applications, and the exceptional crystallization behaviour observed opens new perspectives for scale-up and broader applications.

Bio-based polymer nanocomposites based on nylon 11 and WS2 inorganic nanotubes

Tungsten disulphide nanotubes (INT-WS2) have been successfully dispersed in a bio-based polyamide matrix (nylon 11) by conventional melt processing. The effect of INT-WS2 content on the morphology, thermal stability, crystallization behaviour and dynamic mechanical properties is investigated. The results indicate that these inorganic nanotubes can be efficiently incorporated into the bio-based polymer matrix without the need for modifiers or surfactants. Additionally, it is found that the non-isothermal crystallization behaviour of nylon 11/INT-WS2 depends on both the cooling rate and INT-WS2 concentration. In particular, crystallization kinetics results demonstrate that the nucleating activity of INTs plays a dominant role in accelerating the crystallization of nylon 11. This fact leads to the appearance of the more-disordered phase at higher temperature. More significantly, it was shown that these INT-WS2 nanocomposites can facilitate a good processability and cost efficiency, and will be of interest for many eco-friendly and medical applications.

Investigation of Performance Enhancements for a Liquid-Desiccant Solar Air-Conditioner

Thermally driven liquid-desiccant air-conditioners (LDAC) are a proven but still developing technology. LDACs can use a solar thermal system to reduce the operational cost and environmental impact of the system by reducing the amount of fuel (e.g. natural gas, propane, etc.) used to drive the system. LDACs also have a key benefit of being able to store energy in the form of concentrated desiccant storage. TRNSYS simulations were used to evaluate several different methods of improving the thermal and electrical coefficients of performance (COPt and COPe) and the solar fraction (SF) of a LDAC. The study analyzed a typical June to August cooling season in Toronto, Ontario. Utilizing properly sized, high-efficiency pumps increased the COPe to 3.67, an improvement of 55%. A new design, featuring a heat recovery ventilator on the scavenging-airstream and an energy recovery ventilator on the process-airstream, increased the COPt to 0.58, an improvement of 32%. This also improved the SF slightly to 54%, an increase of 8%. A new TRNSYS TYPE was created to model a stratified desiccant storage tank. Different volumes of desiccant were tested with a range of solar array system sizes. The largest storage tank coupled with the largest solar thermal array showed improvements of 64% in SF, increasing the value to 82%. The COPe was also improved by 17% and the COPt by 9%. When combining the heat recovery systems and the desiccant storage systems, the simulation results showed a 78% increase in COPe and 30% increase in COPt. A 77% improvement in SF and a 17% increase in total cooling rate were also predicted by the simulation. The total thermal energy consumed was 10% lower and the electrical consumption was 34% lower. The amount of non-renewable energy needed from the natural gas boiler was 77% lower. Comparisons were also made between LDACs and vapour-compression (VC) systems. Dependent on set-up, LDACs provided higher latent cooling rates and reduced electrical power consumption. Negatively, a thermal input was required for the LDAC systems but not for the VC systems.

Semisolid structure formation and semisolid casting

Semisolid metal forming has now been accepted as a viable technology for production of components with complex shape and high integrity. The advantages of semisolid metal forming can only be achieved when the feedstock material has a non-dendritic semisolid structure. A controlled nucleation method has been developed to produce such structures for semisolid forming. By controlling grain nucleation and growth, fine-grained and non-dendritic microstructures that are suitable for semisolid casting can be generated. The method was applied to hypoeutectic and hypereutectic Al-Si casting alloys, Al wrought alloys and a Mg alloy. Parameters such as pouring temperature, cooling rate and grain refiner addition were controlled to achieve copious nucleation, nuclei survival and dendritic growth suppression during solidification. The influences of the controlling parameters on the formation of semisolid structure were different for each of these alloy groups. The as-cast structures were then partially remelted and isothermally held. Semisolid structures were developed and followed by semisolid casting into a stepped die.

Iron-rich intermetallic phases and their role in casting defect formation in hypoeutectic Al-Si alloys

Iron is the most common and detrimental impurity in aluminum casting alloys and has long been associated with an increase in casting defects. While the negative effects of iron are clear, the mechanism involved is not fully understood. It is generally believed to be associated with the formation of Fe-rich intermetallic phases. Many factors, including alloy composition, melt superheating, Sr modification, cooling, rate, and oxide bifilms, could play a role. In the present investigation, the interactions between iron and each individual element commonly present in aluminum casting alloys, were investigated using a combination of thermal analysis and interrupted quenching tests. The Fe-rich intermetallic phases were characterized using optical microscope, scanning electron microscope, and electron probe microanalysis (EPMA), and the results were compared with the predictions by Thermocalc. It was found that increasing the iron content changes the precipitation sequence of the beta phase, leading to the precipitation of coarse binary beta platelets at a higher temperature. In contrast, manganese, silicon, and strontium appear to suppress the coarse binary beta platelets, and Mn further promotes the formation of a more compact and less harmful a phase. They are therefore expected to reduce the negative effects of the phase. While reported in the literature, no effect of P on the amount of beta platelets was observed. Finally, attempts are made to correlate the Fe-rich intermetallic phases to the formation of casting defects. The role of the beta phase as a nucleation site for eutectic Si and the role of the oxide bifilms and AIP as a heterogeneous substrate of Fe intermetallics are also discussed.

Modelling of transient heat transfer in annealing furnaces

This is a study of heat transfer in a lift-off furnace which is employed in the batch annealing of a stack of coils of steel strip. The objective of the project is to investigate the various factors which govern the furnace design and the heat transfer resistances, so as to reduce the time of the annealing cycle, and hence minimize the operating costs. The work involved mathematical modelling of patterns of gas flow and modes of heat transfer. These models are: Heat conduction and its conjectures in the steel coils;Convective heat transfer in the plates separating the coils in the stack and in other parts of the furnace; and Radiative and convective heat transfer in the furnace by using the long furnace model. An important part of the project is the development of numerical methods and computations to solve the transient models. A limited number of temperature measurements was available from experiments on a test coil in an industrial furnace. The mathematical model agreed well with these data. The model has been used to show the following characteristics of annealing furnaces, and to suggest further developments which would lead to significant savings: - The location of the limiting temperature in a coil is nearer to the hollow core than to the outer periphery. - Thermal expansion of the steel tends to open the coils, reduces their thermal conductivity in the radial direction, and hence prolongs the annealing cycle. Increasing the tension in the coils and/or heating from the core would overcome this heat transfer resistance. - The shape and dimensions of the convective channels in the plates have significant effect on heat convection in the stack. An optimal design of a channel is shown to be of a width-to-height ratio equal to 9. - Increasing the cooling rate, by using a fluidized bed instead of the normal shell and tube exchanger, would shorten the cooling time by about 15%, but increase the temperature differential in the stack. - For a specific charge weight, a stack of different-sized coils will have a shorter annealing cycle than one of equally-sized coils, provided that production constraints allow the stacking order to be optimal. - Recycle of hot flue gases to the firing zone of the furnace would produce a. decrease in the thermal efficiency up to 30% but decreases the heating time by about 26%.

Fatigue crack growth behaviour in molten-metal processed SiC particle-reinforced aluminium alloys

Fatigue crack initiation and subsequent short crack growth behaviour of 2014-5wt%SiC aluminium alloy composites has been examined in 4-point bend loading using smooth bar specimens. The growth rates of long fatigue cracks have also been measured at different stress ratios using pre-cracked specimens. The distributions of SiC particles and of coarse constituent particles in the matrix (which arise as a result of the molten-metal processing and relatively slow cooling rate) have been investigated. Preferential crack initiation sites were found to be SiC-matrix interfaces, SiC particles associated with constituent particles and the coarse constituent particles themselves. For microstructurally short cracks the dispersed SiC particles also act as temporary crack arresters. In the long crack growth tests, higher fatigue crack growth rates were obtained than for monolithic alloys. This effect is attributed to the contribution of void formation, due to the decohesion of SiC particles, to the fatigue crack growth process in the composite. Above crack depths of about 200 μm 'short' crack growth rates were in good agreement with the long crack data, showing a Pris exponent, m = 4 in both cases. For the long crack and short crack growth tests little effect of specimen orientation and grain size was observed on fatigue crack growth rates, but, specimen orientation affected the toughness. No effect of stress ratio in the range R = 0.2-0.5 was seen for long crack data in the Paris region.