Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC MODELS”

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 811 records · Page 45

Effect of high-intensity x-ray radiation on Bragg diffraction in silicon and diamond

We simulated the Bragg reflection of high-intensity short x-ray pulses from single-crystal silicon and diamond by coupling tight-binding-molecular dynamics with a simple atomic kinetics model. We found that even when the pulse-averaged Bragg intensity degrades significantly, the reflectivity drops only slightly at the beginning of the pulse until the lattice is disordered by non-thermal melting. These results suggest that Bragg reflectors could produce shortened x-ray pulses through temporal slicing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Lattice thermal conductivity of nanograined half-Heusler solid solutions

We report a phenomenological model of atomic weight, lattice constant, temperature, and grain size to calculate the high-temperature lattice thermal conductivity of nanograined solid solutions. The theoretical treatment developed here is reasonably consistent with the experimental results of n-type MNiSn and p-type MCoSb alloys, where M is the combination of Hf, Zr, and Ti. For disordered half-Heusler alloys with moderated grain sizes, we predict that the reduction in lattice thermal conductivity due to grain boundary scattering is independent of the scattering parameter, which characterizes the phonon scattering cross section of point defects. In addition, the lattice thermal conductivity falls off with temperature as T{sup –1∕2} around the Debye temperature.

36 MATERIALS SCIENCE↗

Defect induced plasticity and failure mechanism of boron nitride nanotubes under tension

The effects of Stone-Wales (SW) and vacancy defects on the failure behavior of boron nitride nanotubes (BNNTs) under tension are investigated using molecular dynamics simulations. The Tersoff-Brenner potential is used to model the atomic interaction and the temperature is maintained close to 300 K. The effect of a SW defect is studied by determining the failure strength and failure mechanism of nanotubes with different radii. In the case of a vacancy defect, the effect of an N-vacancy and a B-vacancy is studied separately. Nanotubes with different chiralities but similar diameter is considered first to evaluate the chirality dependence. The variation of failure strength with the radius is then studied by considering nanotubes of different diameters but same chirality. It is observed that the armchair BNNTs are extremely sensitive to defects, whereas the zigzag configurations are the least sensitive. In the case of pristine BNNTs, both armchair and zigzag nanotubes undergo brittle failure, whereas in the case of defective BNNTs, only the zigzag ones undergo brittle failure. An interesting defect induced plastic behavior is observed in defective armchair BNNTs. For this nanotube, the presence of a defect triggers mechanical relaxation by bond breaking along the closest zigzag helical path, with the defect as the nucleus. This mechanism results in a plastic failure.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Communication: Newton homotopies for sampling stationary points of potential energy landscapes

One of the most challenging and frequently arising problems in many areas of science is to find solutions of a system of multivariate nonlinear equations. There are several numerical methods that can find many (or all if the system is small enough) solutions but they all exhibit characteristic problems. Moreover, traditional methods can break down if the system contains singular solutions. Here, we propose an efficient implementation of Newton homotopies, which can sample a large number of the stationary points of complicated many-body potentials. We demonstrate how the procedure works by applying it to the nearest-neighbor ϕ{sup 4} model and atomic clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural and electronic properties of GaN nanowires with embedded In{sub x}Ga{sub 1−x}N nanodisks

In the present study, the effects of various types of strain and indium concentration on the total energy and optoelectronic properties of GaN nanowires (NWs) with embedded In{sub x}Ga{sub 1−x}N nanodisks (NDs) are examined. In particular, the bi-axial, hydrostatic, and uniaxial strain states of the embedded In{sub x}Ga{sub 1−x}N NDs are investigated for multiple In concentrations. Density functional theory is employed to calculate the band structure of the NWs. The theoretical analysis finds that the supercell-size-dependent characteristics calculated for our 972-atom NW models are very close to the infinite supercell-size limit. It is established that the embedded In{sub x}Ga{sub 1−x}N NDs do not induce deep states in the band gap of the NWs. A bowing parameter of 1.82 eV is derived from our analysis in the quadratic Vegard's formula for the band gaps at the various In concentrations of the investigated In{sub x}Ga{sub 1−x}N NDs in GaN NW structures. It is concluded that up to ∼10% of In, the hydrostatic strain state is competitive with the bi-axial due to the radial absorption of the strain on the surfaces. Above this value, the dominant strain state is the bi-axial one. Thus, hydrostatic and bi-axial strain components coexist in the embedded NDs, and they are of different physical origin. The bi-axial strain comes from growth on lattice mismatched substrates, while the hydrostatic strain originates from the lateral relaxation of the surfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab initio calculations of the optical properties of crystalline and liquid InSb

Ab initio calculations of the electronic and optical properties of InSb were performed for both the crystalline and liquid states. Two sets of atomic structure models for liquid InSb at 900 K were obtained by ab initio molecular dynamics simulations. To reduce the effect of structural peculiarities in the liquid models, an averaging of the two sets of the calculated electronic and optical properties corresponding to the two liquid models was performed. The calculated results indicate that, owing to the phase transition from crystal to liquid, the density of states around the Fermi level increases. As a result, the energy band gap opening near the Fermi level disappears. Consequently, the optical properties change from semiconductor to metallic behavior. Namely, owing to the melting of InSb, the interband transition peaks disappear and a Drude-like dispersion is observed in the optical dielectric functions. The optical absorption at a photon energy of 3.06 eV, which is used in Blu-ray Disc systems, increases owing to the melting of InSb. This increase in optical absorption is proposed to result from the increased optical transitions below 2 eV.

36 MATERIALS SCIENCE↗

Molybdenum carbide nanoparticles as catalysts for oil sands upgrading: Dynamics and free-energy profiles

There is no doubt that a huge gap exists in understanding heterogeneous catalysis between a cluster model of a few atoms and a bulk model of periodic slabs. Nanoparticles, which are crucial in heterogeneous catalysis in industry, lie in the middle of the gap. We present here our work on the computational modelling of molybdenum carbide nanoparticles (MCNPs) as the catalysts for the upgrading of oil sands in the in-situ environment, using benzene hydrogenation as a model reaction. With a cluster model, efforts were first made to understand the mechanism of the reaction with a density functional theory (DFT) study on the adsorption of benzene and its hydrogenation product – cyclohexane, as well as the cyclic hydrogenation reaction intermediates on the Mo{sub 2}C(0001) surface. From the thermodynamic data, along with literature information, it was found that the benzene hydrogenation reaction on molybdenum carbide happens most likely through a Langmuir-Hinshelwood mechanism with the gradual lifting up of the benzene molecule. The electron localization function (ELF) was then used to help understand the nature of the interactions between the MCNPs, identifying strong multi-center interactions between the adsorbates and the MCNPs. To enable the treatment of larger nanoparticles, a fast semi-empirical density functional tight-binding (DFTB) method was parameterized. With this method, the potential energy profiles of benzene hydrogenation reactions on different sizes of MCNPs are calculated. The study was then extended to consider a MCNP embedded in solvent (benzene), using a quantum mechanical (DFTB) / molecular mechanical approach. Calculations on the free energies profiles with the umbrella sampling method show that the entropy of the MCNPs and the solvent are essential in understanding the catalytic activity of the transition metal related nanoparticles for solid/liquid heterogeneous catalysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of configurational disorder on the elastic properties of icosahedral boron-rich alloys based on B{sub 6}O, B{sub 13}C{sub 2}, and B{sub 4}C, and their mixing thermodynamics

The elastic properties of alloys between boron suboxide (B{sub 6}O) and boron carbide (B{sub 13}C{sub 2}), denoted by (B{sub 6}O){sub 1−x}(B{sub 13}C{sub 2}){sub x}, as well as boron carbide with variable carbon content, ranging from B{sub 13}C{sub 2} to B{sub 4}C are calculated from first-principles. Furthermore, the mixing thermodynamics of (B{sub 6}O){sub 1−x}(B{sub 13}C{sub 2}){sub x} is studied. A superatom-special quasirandom structure approach is used for modeling different atomic configurations, in which effects of configurational disorder between the carbide and suboxide structural units, as well as between boron and carbon atoms within the units, are taken into account. Elastic properties calculations demonstrate that configurational disorder in B{sub 13}C{sub 2}, where a part of the C atoms in the CBC chains substitute for B atoms in the B{sub 12} icosahedra, drastically increase the Young’s and shear modulus, as compared to an atomically ordered state, B{sub 12}(CBC). These calculated elastic moduli of the disordered state are in excellent agreement with experiments. Configurational disorder between boron and carbon can also explain the experimentally observed almost constant elastic moduli of boron carbide as the carbon content is changed from B{sub 4}C to B{sub 13}C{sub 2}. The elastic moduli of the (B{sub 6}O){sub 1−x}(B{sub 13}C{sub 2}){sub x} system are also practically unchanged with composition if boron-carbon disorder is taken into account. By investigating the mixing thermodynamics of the alloys, in which the Gibbs free energy is determined within the mean-field approximation for the configurational entropy, we outline the pseudo-binary phase diagram of (B{sub 6}O){sub 1−x}(B{sub 13}C{sub 2}){sub x}. The phase diagram reveals the existence of a miscibility gap at all temperatures up to the melting point. Also, the coexistence of B{sub 6}O-rich as well as ordered or disordered B{sub 13}C{sub 2}-rich domains in the material prepared through equilibrium routes is predicted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin waves in the soft layer of exchange-coupled soft/hard bilayers

The magnetic dynamical properties of the soft layer in exchange-coupled soft/hard bilayers have been investigated numerically using a one-dimensional atomic chain model. The frequencies and spatial profiles of spin wave eigenmodes are calculated during the magnetization reversal process of the soft layer. The spin wave modes exhibit a spatially modulated amplitude, which is especially evident for high-order modes. A dynamic pinning effect of surface magnetic moment is observed. The spin wave eigenfrequency decreases linearly with the increase of the magnetic field in the uniformly magnetized state and increases nonlinearly with field when spiral magnetization configuration is formed in the soft layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Conformational landscape of an amyloid intra-cellular domain and Landau-Ginzburg-Wilson paradigm in protein dynamics

The Landau-Ginzburg-Wilson paradigm is proposed as a framework, to investigate the conformational landscape of intrinsically unstructured proteins. A universal Cα-trace Landau free energy is deduced from general symmetry considerations, with the ensuing all-atom structure modeled using publicly available reconstruction programs Pulchra and Scwrl. As an example, the conformational stability of an amyloid precursor protein intra-cellular domain (AICD) is inspected; the reference conformation is the crystallographic structure with code 3DXC in Protein Data Bank (PDB) that describes a heterodimer of AICD and a nuclear multi-domain adaptor protein Fe65. Those conformations of AICD that correspond to local or near-local minima of the Landau free energy are identified. For this, the response of the original 3DXC conformation to variations in the ambient temperature is investigated, using the Glauber algorithm. The conclusion is that in isolation the AICD conformation in 3DXC must be unstable. A family of degenerate conformations that minimise the Landau free energy is identified, and it is proposed that the native state of an isolated AICD is a superposition of these conformations. The results are fully in line with the presumed intrinsically unstructured character of isolated AICD and should provide a basis for a systematic analysis of AICD structure in future NMR experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrostatic and induction effects in the solubility of water in alkanes

Experiments show that at 298 K and 1 atm pressure, the transfer free energy, μ ex , of water from its vapor to liquid normal alkanes C n H 2n+2 (n=5…12) is negative. Earlier it was found that with the united-atom TraPPE model for alkanes and the SPC/E model for water, one had to artificially enhance the attractive alkane-water cross interaction to capture this behavior. Here we revisit the calculation of μ ex using the polarizable AMOEBA and the non-polarizable Charmm General (CGenFF) forcefields. We test both the AMOEBA03 and AMOEBA14 water models; the former has been validated with the AMOEBA alkane model while the latter is a revision of AMOEBA03 to better describe liquid water. We calculate μ ex using the test particle method. With CGenFF, μ ex is positive and the error relative to experiments is about 1.5 k B T. With AMOEBA, μ ex is negative and deviations relative to experiments are between 0.25 k B T (AMOEBA14) and 0.5 k B T (AMOEBA03). Quantum chemical calculations in a continuum solvent suggest that zero point effects may account for some of the deviation. In conclusion, forcefield limitations notwithstanding, electrostatic and induction effects, commonly ignored in consideration of water-alkane interactions, appear to be decisive in the solubility of water in alkanes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Staged Z-pinch modeling of high and low atomic number liners compressing deuterium targets using parameters of the Z pulsed power facility

The staged Z-pinch is a potentially transformative magneto-inertial fusion energy concept where a high atomic number liner implodes on a deuterium or deuterium-tritium target using multi-MA pulsed current. Over the past several years, this concept was studied experimentally on 1 MA facilities with argon or krypton gas puffs compressing a magnetized deuterium target. Consistent thermonuclear neutron yield of 1010 per shot was measured with krypton liners. In this paper, we investigate the fusion performance of deuterium targets of varying density undergoing compression with low (beryllium) and high atomic number liners (silver, tantalum) using parameters of the Z pulsed power facility. Silver and tantalum liners create strong shocks that preheat the target plasma above 100 eV and pile up liner material at the liner-target interface. The increased mass density at the interface creates strong ram pressure just before the pinch stagnation time. The target plasma is heated to 2–4 keV, in contrast to the < 0.5 kV temperature calculated for the beryllium case. The high atomic number liners produce neutron yield orders of magnitude higher than the yield from the low atomic number liner.

Ruskov, E. (ORCID:0000000336724580)↗

Ionization disequilibrium in K- and L-shell ions

Time-gated Sc K-shell and Ge L-shell spectra are presented from a range of characterized thermodynamic states spanning ion densities of 1019–1020cm−3 and plasma temperatures around 2000 eV. For the higher densities studied and temperatures from 1000 to 3000 eV, the Sc and Ge x-ray emission spectra are consistent with steady-state calculations from the modern atomic kinetics model SCRAM. At the lower ion densities achieved through plasma expansion, however, the model calculations require a higher plasma temperature to reproduce the observed Ge spectrum. We attribute this to ionization disequilibrium of the Sc because the ionization time scales exceed the hydrodynamic timescale when the inferred temperatures diverge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density of molten oxides measured in an aero-acoustic levitator

Knowing the thermophysical properties of high-temperature melts can aid the design of melt processes and validate atomic structural models, such as those used in studying glass formation. Property measurements on such melts are challenging, however, due to container-related contamination and heterogeneous nucleation. Containerless processing techniques that employ levitation can be used to avoid these obstacles. In that context, we demonstrate here the application of silhouette imaging to measure the density of molten oxides in an aero-acoustic levitation instrument (AAL). The AAL combines gas jet levitation with actively controlled acoustic positioning to enable full optical access to samples ca. 2–4 mm in diameter, which are laser beam heated and melted. The cross sections of molten drops are imaged using a monochromatic light source and narrowband-filtered camera. Melt volume is calculated from the cross sections and used to find density at several temperatures ranging 1530–1920 K, including up to 350 K of supercooling. We report densities for CaAl 2 O 4 , Ca 12 Al 14 O 33 , CaSiO 3 , their Fe 2 O 3 -bearing analogs, and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). These provide important benchmarks of the capabilities, measurement uncertainties, and future outlook for this technique.

36 MATERIALS SCIENCE↗

Structure of the infectious salmon anemia virus receptor complex illustrates a unique binding strategy for attachment

Orthomyxoviruses are an important family of RNA viruses, which include the various influenza viruses. Despite global efforts to eradicate orthomyxoviral pathogens, these infections remain pervasive. One such orthomyxovirus, infectious salmon anemia virus (ISAV), spreads easily throughout farmed and wild salmonids, constituting a significant economic burden. ISAV entry requires the interplay of the virion-attached hemagglutinin-esterase and fusion glycoproteins. Preventing infections will rely on improved understanding of ISAV entry. Here, we present the crystal structures of ISAV hemagglutinin-esterase unbound and complexed with receptor. Several distinctive features observed in ISAV HE are not seen in any other viral glycoprotein. The structures reveal a unique mode of receptor binding that is dependent on the oligomeric assembly of hemagglutinin-esterase. Importantly, ISAV hemagglutinin-esterase receptor engagement does not initiate conformational rearrangements, suggesting a distinct viral entry mechanism. This work improves our understanding of ISAV pathogenesis and expands our knowledge on the overall diversity of viral glycoprotein-mediated entry mechanisms. Lastly, it provides an atomic-resolution model of the primary neutralizing antigen critical for vaccine development.

59 BASIC BIOLOGICAL SCIENCES↗

Tracking the insulator-to-metal phase transition in VO 2 with few-femtosecond extreme UV transient absorption spectroscopy

We present coulomb correlations can manifest in exotic properties in solids, but how these properties can be accessed and ultimately manipulated in real time is not well understood. The insulator-to-metal phase transition in vanadium dioxide (VO 2 ) is a canonical example of such correlations. Here, few-femtosecond extreme UV transient absorption spectroscopy (FXTAS) at the vanadium M 2,3 edge is used to track the insulator-to-metal phase transition in VO 2 . This technique allows observation of the bulk material in real time, follows the photoexcitation process in both the insulating and metallic phases, probes the subsequent relaxation in the metallic phase, and measures the phase-transition dynamics in the insulating phase. An understanding of the VO 2 absorption spectrum in the extreme UV is developed using atomic cluster model calculations, revealing V 3+ /d 2 character of the vanadium center. We find that the insulator-to-metal phase transition occurs on a timescale of 26 ± 6 fs and leaves the system in a long-lived excited state of the metallic phase, driven by a change in orbital occupation. Potential interpretations based on electronic screening effects and lattice dynamics are discussed. A Mott–Hubbard-type mechanism is favored, as the observed timescales and d 2 nature of the vanadium metal centers are inconsistent with a Peierls driving force. In conclusion, the findings provide a combined experimental and theoretical roadmap for using time-resolved extreme UV spectroscopy to investigate nonequilibrium dynamics in strongly correlated materials.

36 MATERIALS SCIENCE↗