Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Dynamics”

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.

At least 235 records · Page 13

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗

Silicon material task. Part 3: Low-cost silicon solar array project

The feasibility of a process for carbon reduction of low impurity silica in a plasma heat source was investigated to produce low-cost solar-grade silicon. Theoretical aspects of the reaction chemistry were studied with the aid of a computer program using iterative free energy minimization. These calculations indicate a threshold temperature exists at 2400 K below which no silicon is formed. The computer simulation technique of molecular dynamics was used to study the quenching of product species.

Roques, R. A.↗

Solid-liquid phase transition in argon

Starting from the Lennard-Jones interatomic potential, a modified cell theory has been used to describe the solid-liquid phase transition in argon. The cell-size variations may be evaluated by a self-consistent condition. With the inclusion of cell-size variations, the transition temperature, the solid and liquid densities, and the liquid-phase radial-distribution functions have been calculated. These ab initio results are in satisfactory agreement with molecular-dynamics calculations as well as experimental data on argon.

Tsang, T.↗

Changes in "thermal lens" measure diffusivity

In an extension of "thermal lens" effect to new applications and better resolution, two laser beams combine to rapidly measure thermal diffusivity and other molecular dynamic properties. New double-beam technique handles very small samples unlike classical techniques for measuring diffusivity. It can be used for measurements on samples undergoing stress, making it applicable to data collection for structural engineering.

Gupta, A.↗

A molecular model for ice nucleation and growth, attachment 2

The central field empirical pair potential model is applied to studying the effects of kinks, ledges, and vacancies on the absorption of water molecules from the vapor. Molecular dynamics simulations indicate that cluster and/or surface modes play a primary role in the absorption process, the flexibility of the hydrogen bond serves to decrease the energy required for structural interconversion, and the rapid distribution of added energy in a hydrogen bonded system lead to aggregate stability which greatly exceeds that predicted by static energy calculations.

Hale, B. N.↗

Time correlation functions in classical dense fluids

The velocity autocorrelation functions and memory functions of dense classical fluids may be directly obtained from the static radial distribution function g(r) in an approximate way. Following the Mori projection operator formalism, the memory functions may be related to the fluctuating force correlation. At low densities, these functions may be evaluated by following the trajectories of particle pairs in the interatomic potential. At higher densities, the force correlation functions can be evaluated approximately from particle pair trajectories via the potential of the mean force. The results for argon fluids at various densities and temperatures agree satisfactorily with the molecular dynamics and the Enskog values. The decrease of the diffusion coefficient with density is partly due to the nature of g(r) which reflects the stronger clustering of atoms at higher densities.

Tsang, T.↗

Effect of laser frequency on a collision-induced radiative process

A review is presented of the principles of collision induced radiative processes, followed by an examination of the effects of laser frequencies on these processes. A one-dimensional problem involving two electron states is considered, and it is found that the Hamiltonian of the radiation field is dominated by electric-dipole interaction which couples states of different parity. Transitions are noted to be dependent on collisions, and the complexities of three-dimensional systems are expressed as considerations of the angular momentum of the photon, the necessity of treating different states simultaneously, and the fact that a radiation field destroys rotational invariance. Changing the radiation frequency alters the crossing point and offers opportunities to study the interplay of potential surfaces with molecular dynamics. Experiments on Na+A systems are outlined for several collision energies and various laser frequencies. Multiple crossings were obtained, although the total cross-section, at all energies, decreased at 18,350/cm.

Devries, P. L.↗

Ion-induced nucleation in a binary mixture

Recent ion-induced nucleation studies in a binary mixture of H2O and H2SO4 vapors indicate that conventional classical nucleation theories cannot account for the experimentally observed ion-H2O-H2SO4 microcluster spectra. An empirical-analytical, semi-molecular theory of nucleation of a binary mixture of vapors on ions has been developed. This theory includes molecular approach to ion-induced changes in dielectric properties in the immediate neighborhood of an ion core as well as the size dependence of the surface tension of the microcluster. In addition to gaseous sulfuric acid and water molecules, the effect of ion-hydrates and sulfuric acid hydrates has also been considered. This theory is further simplified by using an empirical correction factor for the electrostatic term. This correction factor has been obtained by comparing the simplified model with the molecular dynamics calculations for selected ion hydrate formation. These two theoretical models are then compared with the classical theory for a variety of experimental conditions.

Smith, A. C.↗

Computer simulation of surface and film processes

A molecular dynamics technique based upon Lennard-Jones type pair interactions is used to investigate time-dependent as well as equilibrium properties. The case study deals with systems containing Si and O atoms. In this case a more involved potential energy function (PEF) is employed and the system is simulated via a Monte-Carlo procedure. This furnishes the equilibrium properties of the system at its interfaces and surfaces as well as in the bulk.

Tiller, W. A.↗

Models of Rapid Solidification

Laser annealing studies provide much information on various consequences of rapid solidification, including the trapping of impurities in the crystal, the generation of vacancies and twins, and on the fundamental limits to the speed of the crystal-melt interface. Some results obtained by molecular dynamics methods of the solidification of a Lennard-Jones liquid are reviewed. An indication of the relationship of interface speed to undercooling for certain materials can be derived from this model. Ising model simulations of impurity trapping in silicon are compared with some of the laser annealing results. The consequences of interface segregation and atomic strain are discussed.

Gilmer, G. H.↗

Quantum noise and the threshold of hearing

It is argued that the sensitivity of the ear reaches a limit imposed by the uncertainty principle. This is possible only if the receptor cell holds the detector elements in a special nonequilibrium state which has the same noise characteristics as a ground (T = 0 K) state. To accomplish this 'active cooling' the molecular dynamics of the system must maintain quantum mechanical coherence over the time scale of the measurement.

Bialek, W.↗

Algorithmic commonalities in the parallel environment

The ultimate aim of this project was to analyze procedures from substantially different application areas to discover what is either common or peculiar in the process of conversion to the Massively Parallel Processor (MPP). Three areas were identified: molecular dynamic simulation, production systems (rule systems), and various graphics and vision algorithms. To date, only selected graphics procedures have been investigated. They are the most readily available, and produce the most visible results. These include simple polygon patch rendering, raycasting against a constructive solid geometric model, and stochastic or fractal based textured surface algorithms. Only the simplest of conversion strategies, mapping a major loop to the array, has been investigated so far. It is not entirely satisfactory.

Mcanulty, Michael A.↗

Surface potential of the water liquid-vapor interface

An analysis of an extended molecular dynamics calculation of the surface potential (SP) of the water liquid-vapor interface is presented. The SP predicted by the TIP4P model is -(130 + or - 50) mV. This value is of reasonable magnitude but of opposite sign to the expectations based on laboratory experiments. The electrostatic potential shows a nonmonotonic variation with depth into the liquid.

Wilson, Michael A.↗

A proposed non-intrusive method for finding coefficients of slip and molecular reflectivity in microgravity

A proposed experimental program to look at a series of vapor transport properties measured along solid and liquid surfaces is described. The research objectives proposed are: (1) with accuracy otherwise unobtainable on ground, to determine the coefficient of slip measured between gases and the surfaces of liquids and solids; (2) for the first time, to classify and tabulate dominant surface effects found for a variety of solids, particularly those crystalized by vapor transport; and (3) to extend understanding of settling rates predicted for cosmic dust and condensed vapor falling through planetary atmospheres. The method used to obtain these objectives, has aided, to an order of magnitude, understanding of various liquid-gas interfaces such as oil and water. But to date, no similar characterization has proved successful for solids or liquids of uncertain densities. Likewise, no data exist in either ground-based research or as part of a microgravity program that, when collected with the high accuracy expected in low gravity, could definitely settle outstanding questions in kinetic theory, molecular dynamics, and cosmic physics.

Noever, D. A.↗

Perception of numerical methods in rarefied gasdynamics

The relationships between various numerical methods applied to problems in rarefied gasdynamics are discussed, with emphasis on conflicting viewpoints and computational requirements associated with physical simulation versus the numerical solution of the Boltzmann equation. The basic differences between the molecular dynamics and direct simulation methods are shown to affect their applicability to dense and rarefied flows. Methods for the probabilistic selection of representative collision in the direct simulation Monte Carlo method are reviewed. A method combining the most desirable features of the earlier methods is presented.

Bird, G. A.↗

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of the evolution of life. Using computational and modeling techniques and guided by analogies from present day systems, we seek to understand the emergence of the genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. Modeling of the ancestral aminoacyl-tRNA-synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. The minimal structural requirements for the catalysis of tRNA aminoacylation are being explored. A formation of an aminoacyl adenylate was studied in the framework of ab initio molecular orbital theory. The role of individual residues in the vicinity of the TyrRS active site was examined, and the effect of all possible amino acids substitutions near the active site was examined. A formation of aminoacyl tRNA was studied by the molecular modeling system SYBYL with the high resolution crystallographic structures of the present day tRNA, aRS's complexes. The ultimate goal is to propose a simple RNA segment that is small enough to be build in the primordial chemical environment but maintains the specificity and catalytic activity of the contemporary RNA enzyme. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation were studied. The structural characteristics of the target helices, particularly a possible role for the G-T pair, is also studied by a molecular dynamics (MD) simulation technique.

Rein, Robert↗