Surface Kinetics and Stability of Faceted Crystal Growth
The applicability of Gibbs-Thompson Law and other fundamental rules to faceted crystal growth will be discussed. Tutorial lecture will include fundamentals of interface kinetics.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The applicability of Gibbs-Thompson Law and other fundamental rules to faceted crystal growth will be discussed. Tutorial lecture will include fundamentals of interface kinetics.
The optical absorbance of the vapor phase over various In-Se compositions between 33.3-60.99 at.% Se and 673-1418 K was measured and used to obtain the partial pressures of Se2(g) and In2Se(g). The results are in agreement with silica Bourdon gauge measurements for compositions between 50-61 at.%, but significantly higher than those from Knudsen cell and simultaneous Knudsen-torsion cell measurements. It is found that 60.99 at.% Se lies outside the sesquiselenide homogeneity range and 59.98 at.% Se lies inside and is the congruently melting composition. The Gibbs energy of formation of the liquid from its pure liquid elements between 1000-1300 K is essentially independent of temperature and falls between -36 to -38 kJ per g atomic weight for 50 and 56% Se at 1200 and 1300 K.
We have pursued the use of biological crosslinking molecules of several types to make colloidal materials at relatively low volume fraction of colloidal particles. The objective is to make binary alloys of colloidal particles, made of two different colloidal particles coated with complementary biological lock-and-key binding molecules, which assemble due to the biological specificity. The long-term goal is to use low affinity lock-and-key biological interactions, so that the can anneal to form crystalline states. We have used a variety of different surface chemistries in order to make colloidal materials. Our first system involved using selectin-carbohydrate (sialyl-Lewis) interactions; this chemistry is derived from immune system. This chemical interaction is of relatively low affinity, with timescales for dissociation of several seconds. Furthermore, the adhesion mediated by these molecules can be reversed by the chelation of calcium atoms; thus assembled structures can be disassembled reversibly. Our second system employed avidin-biotin chemistry. This well-studied system is of high affinity, and is generally irreversible on a laboratory time-scale. Thus, we would expect selectin-carbohydrate interactions at high molecular density and avidin-biotin interactions to give kinetically-trapped structures; however, at low densities, we would expect significant differences in the structure and dynamics of the two materials, owing to their very different release rates. We have also begun to use a third chemistry - DNA hybridization. By attaching single stranded DNA oligonucleotide chains to beads, we can drive the assembly of colloidal materials by hybridization of complementary DNA chains. It is well known that DNA adenosine-thymine (A-T) and guanine-cytosine (G-C) bases hybridize pairwise with a Gibbs free energy change of 1.7 kcal/mol per base; thus, the energy of the assembly can be modulated by altering the number of complementary bases in the DNA chains. Using these different crosslinking molecules, we have assembled colloidal materials from different-sized colloidal particles, A and B. In the first sets of experiment, we used high densities of adhesion molecules, and 0.96 micron (A) and 5.5 micron (B) diameter particles. The high density of adhesion molecules means that the structures are kinetically trapped in nonequilibrium configurations. The structure of the suspension can be varied by changing the number ratio of the two types of colloidal particles, NA and NB, where A is the smaller particle. With carbohydrate-selectin or avidin-biotin interactions, large NA/NB leads to the formation of colloidal micelles, with the large center B particle surrounded by many smaller A particles. As the ratio NA/NB decreases, the structures become more extended, approaching the formation of macro-Rouse polymers - extended linear chains where A beads are connected with intervening small B linkers.
In this paper we present a new form of the collocation method that allows one to find very accurate solutions to time marching problems without the unwelcome appearance of Gibb's phenomenon oscillations. The basic method is applicable to any partial differential equation whose solution is a continuous, albeit possibly rapidly varying function. Discontinuous functions are dealt with by replacing the function in a small neighborhood of the discontinuity with a spline that smoothly connects the function segments on either side of the discontinuity. This will be demonstrated when the solution to the inviscid Burgers equation is discussed.
Development of hydrogen in sealed silica glass ampoules during annealing at elevated temperatures was investigated. The dependence of hydrogen pressure in the ampoules as a function of time, for different temperatures and ampoule parameters was measured. The process was modeled assuming chemical solution of hydrogen according to the reaction: silica + H2 = H- Si= + H-O-Si=. The equilibrium constant of the reaction was determined by fitting the theoretical curves to the experimental data. The Gibbs function for this reaction was estimated at deltaG = -25.8 + 54T.
A numerical mathematical model capable of accurately describing the evolution of the shape of the solid/liquid interface in the proximity of a foreign particle is presented in this paper. The model accounts for the influence of the temperature gradient and the Gibbs-Thomson and disjoining pressure effects. It shows that for the systems characterized by k(sub P) < k(sub L) the disjoining pressure causes the interface curvature to change its sign in the close-contact particle/interface region. It also shows that the increase of the temperature gradient diminishes the effect of the disjoining pressure. Calculated critical solidification velocities for the pushing/engulfment transition are compared with experimental measurements performed in microgravity conditions.
The crystals are built of molecules of protein, nucleic acid and their complexes, like viruses, approx. 5x10(exp 3)+ 3x10(exp 6) Da in weight and 2 + 20 nm in effective diameter. This size strongly exceeds action range of molecular forces and makes a big difference with inorganic crystals. Intermolecular contacts form patches on the biomacromolecular surface. Each patch may occupy only a small percent of the whole surface and vary from polymorph to polymorph of the same protein. Thus, under different conditions (pH, solution chemistry, temperature, any area on the macromolecular surface may form a contact. The crystal Young moduli, E approx. equals 0.1 + 0.5 GPa are more than 10 times lower than that of inorganics and the biomolecules themselves. Water within biocrystals (30-70%) is unable to flow unless typical deformation time is longer than approx. 10(exp -5)s. This explains the discrepancy between light scattering and static measurements of E. Nucleation and Growth requires typically concentrations exceeding the equilibrium ones up to 100 times - because of the new size scale results in 10 - 10(exp 3) times lower kinetic coefficients than that needed for inorganic solution growth. All phenomena observed in the latter occur with protein crystallization and are even better studied by AFM. Crystals are typically facetted. Among unexpected findings of general significance are - net molecular exchange flux at kinks is much lower than that expected from supersaturation, steps with low (< approx. 10(exp -2)) kink density at steps follow Gibbs-Thomson law only at very low supersaturations, step segment growth rate may be independent of step energy. Crystal perfection is a must of biocrystallization to achieve the major goal to find 3-D atomic structure of biomacromolecules by x-ray diffraction. Poor diffraction resolution (> 3Angstrom) makes crystallization a bottleneck for structural biology. All defects typical of small molecule crystals are found in biocrystals, but the defects responsible for poor resolution are not identified. Conformational changes are one of them. Biocrystallization in microgravity reportedly results in 20% cases of better crystals. The mechanism of how lack of convection can do this is still not clear. Lower supersaturation, self-purification &om preferentially trapped homologous impurities and step bunching are viable hypotheses.
The distribution of ceramic particles in a metal matrix composite material depends primarily on the interaction of the particles with the solid/liquid interface during the solidification process. A numerical model that describes the evolution of the shape of the solid/liquid interface in the proximity of a foreign particle will presented in this paper. The model accounts for the influence of the temperature gradient and the Gibbs-Thomson and disjoining pressure effects. It shows that for the systems characterized by k(sub p) < k(sub L) the disjoining pressure causes the interface curvature to change its sign in the close-contact particle/interface region. It also shows that the increase of the temperature gradient diminishes the effect of the disjoining pressure. The analysis of the numerical results obtained for a large range of processing conditions and materials parameters has led to the development of an analytical solution for the critical velocity of pushing/engulfinent transition. The theoretical results will be discussed and compared with the experimental measurements performed under microgravity conditions.
The two most common approaches used to formulate thermodynamic properties of pure substances are fundamental (or characteristic) equations of state (Helmholtz and Gibbs functions) and a piecemeal approach that is described in Adebiyi and Russell (1992). This paper neither presents a different method to formulate thermodynamic properties of pure substances nor validates the aforementioned approaches. Rather its purpose is to present a method to generate property tables from existing property packages and a method to facilitate the accurate interpretation of fluid thermodynamic property data from those tables. There are two parts to this paper. The first part of the paper shows how efficient and usable property tables were generated, with the minimum number of data points, using an aerospace industry standard property package. The second part describes an innovative interpolation technique that has been developed to properly obtain thermodynamic properties near the saturated liquid and saturated vapor lines.
The availability and proper utilization of fluid properties is of fundamental importance in the process of mathematical modeling of propulsion systems. Real fluid properties provide the bridge between the realm of pure analytiis and empirical reality. The two most common approaches used to formulate thermodynamic properties of pure substances are fundamental (or characteristic) equations of state (Helmholtz and Gibbs functions) and a piecemeal approach that is described, for example, in Adebiyi and Russell (1992). This paper neither presents a different method to formulate thermodynamic properties of pure substances nor validates the aforementioned approaches. Rather its purpose is to present a method to be used to facilitate the accurate interpretation of fluid thermodynamic property data generated by existing property packages. There are two parts to this paper. The first part of the paper shows how efficient and usable property tables were generated, with the minimum number of data points, using an aerospace industry standard property package (based on fundamental equations of state approach). The second part describes an innovative interpolation technique that has been developed to properly obtain thermodynamic properties near the saturated liquid and saturated vapor lines.
This study utilizes Saccharomyces cerevisiae to study genetic responses to suspension culture. The suspension culture system used in this study is the high-aspect-ratio vessel, one type of the rotating wall vessel, that provides a high rate of gas exchange necessary for rapidly dividing cells. Cells were grown in the high-aspect-ratio vessel, and DNA microarray and metabolic analyses were used to determine the resulting changes in yeast gene expression. A significant number of genes were found to be up- or downregulated by at least twofold as a result of rotational growth. By using Gibbs promoter alignment, clusters of genes were examined for promoter elements mediating these genetic changes. Candidate binding motifs similar to the Rap1p binding site and the stress-responsive element were identified in the promoter regions of differentially regulated genes. This study shows that, as in higher order organisms, S. cerevisiae changes gene expression in response to rotational culture and also provides clues for investigations into the signaling pathways involved in gravitational response.
The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.
Selected topics that may be of interest for both crystal-structure and crystal-growth communities are overviewed. The growth of protein crystals, along with that of some other compounds, is one of the topics, and recent insights into related phenomena are considered as examples of applications of general principles. The relationship between crystal growth shape and structure is reviewed and an attempt to introduce semiquantitative characterization of binding for proteins is made. The concept of kinks for complex structures is briefly discussed. Even at sufficiently low supersaturations, the fluctuation of steps may not be sufficient to implement the Gibbs-Thomson law if the kink density is low enough. Subsurface ordering of liquids and growth of rough interfaces from melts is discussed. Crystals growing in microgravity from solution should be more perfect if they preferentially trap stress-inducing impurities, thus creating an impurity-depleted zone around themselves. Evidently, such a zone is developed only around the crystals growing in the absence of convection. Under terrestrial conditions, the self-purified depleted zone is destroyed by convection, the crystal traps more impurity and grows stressed. The stress relief causes mosaicity. In systems containing stress-inducing but poorly trapped impurities, the crystals grown in the absence of convection should be worse than those of their terrestrial counterparts.
Analysis of the chemical disequilibrium provided by the mixing of hydrothermal fluids and seawater in present-day systems indicates that organic synthesis from CO2 or carbonic acid is thermodynamically favoured in the conditions in which hyperthermophilic microorganisms are known to live. These organisms lower the Gibbs free energy of the chemical mixture by synthesizing many of the components of their cells. Primary productivity is enormous in hydrothermal systems because it depends only on catalysis of thermodynamically favourable, exergonic reactions. It follows that hydrothermal systems may be the most favourable environments for life on Earth. This fact makes hydrothermal systems logical candidates for the location of the emergence of life, a speculation that is supported by genetic evidence that modern hyperthermophilic organisms are closer to a common ancestor than any other forms of life. The presence of hydrothermal systems on the early Earth would correspond to the presence of liquid water. Evidence that hydrothermal systems existed early in the history of Mars raises the possibility that life may have emerged on Mars as well. Redox reactions between water and rock establish the potential for organic synthesis in and around hydrothermal systems. Therefore, the single most important parameter for modelling the geochemical emergence of life on the early Earth or Mars is the composition of the rock which hosts the hydrothermal system.
The basic kinetic and thermodynamic concepts of crystal growth will be revisited in view of recent AFM and interferometric findings. These concepts are as follows: 1) The Kossel crystal model that allows only one kink type on the crystal surface. The modern theory is developed overwhelmingly for the Kessel model; 2) Presumption that intensive step fluctuations maintain kink density sufficiently high to allow applicability of Gibbs-Thomson law; 3) Common experience that unlimited step bunching (morphological instability) during layer growth from solutions and supercooled melts always takes place if the step flow direction coincides with that of the fluid.
1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.
Nanotechnology has created a demand for new fabrication methods with an emphasis on simple, low-cost techniques. Directional solidification of eutectics (DSE) is an unconventional approach in comparison to low-temperature biomimetic approaches. A technical challenge for DSE is producing microstructural architectures on the nanometer scale. In both processes, the driving force is the minimization of Gibb's free energy. Selfassembly by biomimetic approaches depends on weak interaction forces between organic molecules to define the architectural structure. The architectural structure for solidification depends on strong chemical bonding between atoms. Constituents partition into atomic-level arrangements at the liquid-solid interface to form polyphase structures, and this atomic-level arrangement at the liquid-solid interface is controlled by atomic diffusion and total undercooling due to composition (diffusion), kinetics, and curvature of the boundary phases. Judicious selection of the materials system and control of the total undercooling are the keys to producing structures on the nanometer scale. The silicon-titanium silicide (Si-TiSi2) eutectic forms a rod structure under isothermal cooling conditions. At the NASA Glenn Research Center, directional solidification was employed along with a thermal gradient to promote uniform rods oriented with the thermal gradient. The preceding photomicrograph shows the typical transverse microstructure of a solidified Si-TiSi2 eutectic composition. The dark and light gray regions are Si and TiSi2, respectively. Preferred rod orientation along the thermal gradient was poor. The ordered TiSi2 rods have a narrow distribution in diameter of 2 to 3 m, as shown. The rod diameter showed a weak dependence on process conditions. Anisotropic etch behavior between different phases provides the opportunity to fabricate structures with high aspect ratios. The photomicrographs show the resulting microstructure after a wet chemical etch and a dry plasma etch. The wet chemical etches the silicon away, exposing the TiSi2 rods, whereas plasma etching preferentially etches the Si-TiSi2 interface to form a crater. The porous architectures are applicable to fabricating microdevices or creating templates for part fabrication. The porous rod structure can serve as a platform for fabricating microplasma devices for propulsion or microheat exchangers and for fabricating microfilters for miniatured chemical reactors. Although more work is required, self-assembly from DSE can have a role in microdevice fabrication.
High temperature high vacuum electrostatic levitation was combined with DSC experiments to determine the specific heat C(sub p) of the undercooled Zr(sub 41.2)Ti(sub 13.8)Cu(sub 12.5)Ni(sub 10.0)Be(sub 22.5) liquid as a function of temperature. The containerless approach made it possible to undercool the melt to the glass transition temperature without inducing nucleation. Because the cooling process was purely radiative, noncontact temperature measurement techniques could be used to determine the specific heat to total hemispherical emissivity ratio, C(sub p)/epsilon(sub t), for the undercooled liquid region. Using C(sub p) values which were independently obtained by DSC, epsilon(sub t) could be determined. With knowledge of C(sub p) of the undercooled liquid it was possible to determine other thermodynamic properties such as Gibbs free energy and entropy as a function of undercooling.