Monoatomic layer removal mechanism in chemical mechanical polishing process: A molecular dynamics study

NASA Astrophysics Data System (ADS)

Si, Lina; Guo, Dan; Luo, Jianbin; Lu, Xinchun

2010-03-01

Molecular dynamics simulations of nanoscratching processes were used to study the atomic-scale removal mechanism of single crystalline silicon in chemical mechanical polishing (CMP) process and particular attention was paid to the effect of scratching depth. The simulation results under a scratching depth of 1 nm showed that a thick layer of silicon material was removed by chip formation and an amorphous layer was formed on the silicon surface after nanoscratching. By contrast, the simulation results with a depth of 0.1 nm indicated that just one monoatomic layer of workpiece was removed and a well ordered crystalline surface was obtained, which is quite consistent with previous CMP experimental results. Therefore, monoatomic layer removal mechanism was presented, by which it is considered that during CMP process the material was removed by one monoatomic layer after another, and the mechanism could provide a reasonable understanding on how the high precision surface was obtained. Also, the effects of the silica particle size and scratching velocity on the removal mechanism were investigated; the wear regimes and interatomic forces between silica particle and workpiece were studied to account for the different removal mechanisms with indentation depths of 0.1 and 1 nm.

Numerical Optimization of Density Functional Tight Binding Models: Application to Molecules Containing Carbon, Hydrogen, Nitrogen, and Oxygen

DOE PAGES

Krishnapriyan, A.; Yang, P.; Niklasson, A. M. N.; ...

2017-10-17

New parametrizations for semiempirical density functional tight binding (DFTB) theory have been developed by the numerical optimization of adjustable parameters to minimize errors in the atomization energy and interatomic forces with respect to ab initio calculated data. Initial guesses for the radial dependences of the Slater- Koster bond integrals and overlap integrals were obtained from minimum basis density functional theory calculations. The radial dependences of the pair potentials and the bond and overlap integrals were represented by simple analytic functions. The adjustable parameters in these functions were optimized by simulated annealing and steepest descent algorithms to minimize the value ofmore » an objective function that quantifies the error between the DFTB model and ab initio calculated data. The accuracy and transferability of the resulting DFTB models for the C, H, N, and O system were assessed by comparing the predicted atomization energies and equilibrium molecular geometries of small molecules that were not included in the training data from DFTB to ab initio data. The DFTB models provide accurate predictions of the properties of hydrocarbons and more complex molecules containing C, H, N, and O.« less

Effect of orientation of prismatic dislocation loops on interaction with free surfaces in BCC iron

NASA Astrophysics Data System (ADS)

Fikar, Jan; Gröger, Roman; Schäublin, Robin

2017-12-01

The prismatic loops appear in metals as a result of high-energy irradiation. Understanding their formation and interaction is important for quantification of irradiation-induced deterioration of mechanical properties. Characterization of dislocation loops in thin foils is commonly made using transmission electron microscopy (TEM), but the results are inevitably influenced by the proximity of free surfaces. The prismatic loops are attracted to free surfaces by image forces. Depending on the type, shape, size, orientation and depth of the loop in the foil, they can escape to the free surface creating denuded loop-free zones and thus invalidating TEM observations. In our previous studies we described a simple general method to determine the critical depth and the critical stress to move prismatic dislocation loops. The critical depths can be further used to correct measurements of the loop density by TEM. Here, we use this procedure to compare 〈100〉 loops and 1/2 〈111〉 loops in body-centered cubic (BCC) iron. The influences of the interatomic potential and the loop orientation are studied in detail. The difference between interstitial and vacancy type loop is also investigated.

Reexamination of the interaction of atoms with a LiF(001) surface

NASA Astrophysics Data System (ADS)

Miraglia, J. E.; Gravielle, M. S.

2017-02-01

Pairwise additive potentials for multielectronic atoms interacting with a LiF(001) surface are revisited by including an improved description of the electron density associated with the different lattice sites, as well as nonlocal electron density contributions. Within this model, the electron distribution around each ionic site of the crystal is described by means of a so-called "onion" approach that accounts for the influence of the Madelung potential. From such densities, binary interatomic potentials are then derived by using well-known nonlocal functionals. Rumpling and long-range contributions due to projectile polarization and van der Waals forces are also included. We apply this pairwise additive approximation to evaluate the interaction potential between closed-shell (He, Ne, Ar, Kr, and Xe) and open-shell (N, S, and Cl) atoms and the LiF surface, analyzing the relative importance of the different contributions. The performance of the proposed potentials is assessed by contrasting angular positions of rainbow and supernumerary rainbow maxima produced by fast grazing incidence with available experimental data. One important result of our model is that both van der Waals contributions and thermal lattice vibrations play a negligible role for normal energies in the eV range.

Effects of functional group mass variance on vibrational properties and thermal transport in graphene

NASA Astrophysics Data System (ADS)

Lindsay, L.; Kuang, Y.

2017-03-01

Intrinsic thermal resistivity critically depends on features of phonon dispersions dictated by harmonic interatomic forces and masses. Here we present the effects of functional group mass variance on vibrational properties and thermal conductivity (κ ) of functionalized graphene from first-principles calculations. We use graphane, a buckled graphene backbone with covalently bonded hydrogen atoms on both sides, as the base material and vary the mass of the hydrogen atoms to simulate the effect of mass variance from other functional groups. We find nonmonotonic behavior of κ with increasing mass of the functional group and an unusual crossover from acoustic-dominated to optic-dominated thermal transport behavior. We connect this crossover to changes in the phonon dispersion with varying mass which suppress acoustic phonon velocities, but also give unusually high velocity optic modes. Further, we show that out-of-plane acoustic vibrations contribute significantly more to thermal transport than in-plane acoustic modes despite breaking of a reflection-symmetry-based scattering selection rule responsible for their large contributions in graphene. This work demonstrates the potential for manipulation and engineering of thermal transport properties in two-dimensional materials toward targeted applications.

QMD and classical MD simulation of alpha boron and boron-carbide behavior under pressure

NASA Astrophysics Data System (ADS)

Yanilkin, Alexey; Korotaev, Pavel; Kuksin, Alexey; Pokatashkin, Pavel

2015-06-01

Boron and some boron-rich compounds are super-hard and light-weighted material with a wide range of different applications. Nevertheless, the question of its behavior under pressure is still open. In the present work we study the equation of state (EOS), stability and deformation of α-B and B4C under high pressure within quantum and classical molecular dynamics (QMD and MD). Based on QMD results the finite temperature EOSs are revealed. CBC chain bending, amorphization and recrystallization of B4C are investigated under hydrostatic, uniform and uniaxial compression. The influence of nonhydrostatic loading is discussed. Angular dependent interatomic potentials are derived by force-matching method. The properties of α-B and B4C, obtained by classical potential, are verified. Structure, bulk modulus, pressure-volume relation, Gruneisen and thermal expansion coefficients are in good agreement with both ab initio and experimental data. These potentials are used to study shock wave propagation in a single crystal of α-B and B4C. Two mechanisms of shear deformation are observed: stacking fault formation and local amorphization. The crystallographic orientations of defects are in a good agreement with experiments.

Effects of functional group mass variance on vibrational properties and thermal transport in graphene

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lindsay, L.; Kuang, Y.

Intrinsic thermal resistivity critically depends on features of phonon dispersions dictated by harmonic interatomic forces and masses. We present the effects of functional group mass variance on vibrational properties and thermal conductivity (κ ) of functionalized graphene from first principles calculations. We also use graphane, a buckled graphene backbone with covalently bonded Hydrogen atoms on both sides, as the base material and vary the mass of the Hydrogen atoms to simulate the effect of mass variance from other functional groups. We find non-monotonic behavior of κ with increasing mass of the functional group and an unusual cross-over from acoustic-dominated tomore » optic-dominated thermal transport behavior. We connect this cross-over to changes in the phonon dispersion with varying mass which suppress acoustic phonon velocities, but also give unusually high velocity optic modes. Further, we show that out-of-plane acoustic vibrations contribute significantly more to thermal transport than in-plane acoustic modes despite breaking of a reflection symmetry based scattering selection rule responsible for their large contributions in graphene. Our work demonstrates the potential for manipulation and engineering of thermal transport properties in two dimensional materials toward targeted applications.« less

DOE Office of Scientific and Technical Information (OSTI.GOV)

Quashie, Edwin E.; Saha, Bidhan C.; Correa, Alfredo A.

Here, we present an ab initio study of the electronic stopping power of protons in copper over a wide range of proton velocities v = 0.02–10a.u. where we take into account nonlinear effects. Time-dependent density functional theory coupled with molecular dynamics is used to study electronic excitations produced by energetic protons. A plane-wave pseudopotential scheme is employed to solve the time-dependent Kohn-Sham equations for a moving ion in a periodic crystal. The electronic excitations and the band structure determine the stopping power of the material and alter the interatomic forces for both channeling and off-channeling trajectories. Our off-channeling results aremore » in quantitative agreement with experiments, and at low velocity they unveil a crossover region of superlinear velocity dependence (with a power of ~1.5) in the velocity range v = 0.07–0.3a.u., which we associate to the copper crystalline electronic band structure. The results are rationalized by simple band models connecting two separate regimes. We find that the limit of electronic stopping v → 0 is not as simple as phenomenological models suggest and it is plagued by band-structure effects.« less

Interatomic Coulombic Decay Effects in Theoretical DNA Recombination Systems Involving Protein Interaction Sites

NASA Astrophysics Data System (ADS)

Vargas, E. L.; Rivas, D. A.; Duot, A. C.; Hovey, R. T.; Andrianarijaona, V. M.

2015-03-01

DNA replication is the basis for all biological reproduction. A strand of DNA will ``unzip'' and bind with a complimentary strand, creating two identical strands. In this study, we are considering how this process is affected by Interatomic Coulombic Decay (ICD), specifically how ICD affects the individual coding proteins' ability to hold together. ICD mainly deals with how the electron returns to its original state after excitation and how this affects its immediate atomic environment, sometimes affecting the connectivity between interaction sites on proteins involved in the DNA coding process. Biological heredity is fundamentally controlled by DNA and its replication therefore it affects every living thing. The small nature of the proteins (within the range of nanometers) makes it a good candidate for research of this scale. Understanding how ICD affects DNA molecules can give us invaluable insight into the human genetic code and the processes behind cell mutations that can lead to cancer. Authors wish to give special thanks to Pacific Union College Student Senate in Angwin, California, for their financial support.

Physical Realization of von Neumann Lattices in Rotating Bose Gases with Dipole Interatomic Interactions.

PubMed

Cheng, Szu-Cheng; Jheng, Shih-Da

2016-08-22

This paper reports a novel type of vortex lattice, referred to as a bubble crystal, which was discovered in rapidly rotating Bose gases with long-range interactions. Bubble crystals differ from vortex lattices which possess a single quantum flux per unit cell, while atoms in bubble crystals are clustered periodically and surrounded by vortices. No existing model is able to describe the vortex structure of bubble crystals; however, we identified a mathematical lattice, which is a subset of coherent states and exists periodically in the physical space. This lattice is called a von Neumann lattice, and when it possesses a single vortex per unit cell, it presents the same geometrical structure as an Abrikosov lattice. In this report, we extend the von Neumann lattice to one with an integral number of flux quanta per unit cell and demonstrate that von Neumann lattices well reproduce the translational properties of bubble crystals. Numerical simulations confirm that, as a generalized vortex, a von Neumann lattice can be physically realized using vortex lattices in rapidly rotating Bose gases with dipole interatomic interactions.

Effect of alloying on screw dislocation structure in Mo: atomistic modelling approach with ab-initio parametrization

NASA Astrophysics Data System (ADS)

Gornostyrev, Yu. N.

2005-03-01

The plastic deformation in bcc metals is realized by the motion of screw dislocations with a complex star-like non-planar core. In this case, the direct investigation of the solute effect by first principles electronic structure calculations is a challenging problem for which we follow a combined approach that includes atomistic dislocation modelling with ab-initio parametrization of interatomic interactions. The screw dislocation core structure in Mo alloys is described within the model of atomic row displacements along a dislocation line with the interatomic row potential estimated from total energy full-potential linear muffin-tin orbital (FLMTO) calculations with the generalized gradient approximation (GGA) for the exchange-correlation potential. We demonstrate (1) that the solute effect on the dislocation structure is different for ``hard'' and ``easy'' cores and (2) that the softener addition in a ``hard'' core gives rise to a structural transformation into a configuration with a lower energy through an intermediate state. The softener solute is shown to disturb locally the three-fold symmetry of the dislocation core and the dislocation structure tends to the split planar core.

Probing the interatomic potential of solids with strong-field nonlinear phononics

NASA Astrophysics Data System (ADS)

von Hoegen, A.; Mankowsky, R.; Fechner, M.; Först, M.; Cavalleri, A.

2018-03-01

Nonlinear optical techniques at visible frequencies have long been applied to condensed matter spectroscopy. However, because many important excitations of solids are found at low energies, much can be gained from the extension of nonlinear optics to mid-infrared and terahertz frequencies. For example, the nonlinear excitation of lattice vibrations has enabled the dynamic control of material functions. So far it has only been possible to exploit second-order phonon nonlinearities at terahertz field strengths near one million volts per centimetre. Here we achieve an order-of-magnitude increase in field strength and explore higher-order phonon nonlinearities. We excite up to five harmonics of the A1 (transverse optical) phonon mode in the ferroelectric material lithium niobate. By using ultrashort mid-infrared laser pulses to drive the atoms far from their equilibrium positions, and measuring the large-amplitude atomic trajectories, we can sample the interatomic potential of lithium niobate, providing a benchmark for ab initio calculations for the material. Tomography of the energy surface by high-order nonlinear phononics could benefit many aspects of materials research, including the study of classical and quantum phase transitions.

Dzyaloshinskii-Moriya interaction and magnetic anisotropies in Uranium compounds

NASA Astrophysics Data System (ADS)

Sandratskii, L. M.

2018-05-01

We report on the first-principles study of complex noncollinear magnetic structures in Uranium compounds. We contrast two cases. The first is the periodic magnetic structure of U2Pd2In with exactly orthogonal atomic moments, the second is an incommensurate plane spiral structure of UPtGe where the angle between atomic moments of nearest neighbors is also close to 90°. We demonstrate that the hierarchy of magnetic interactions leading to the formation of the magnetic structure is opposite in the two cases. In U2Pd2In, the magnetic anisotropy plays the leading role, followed by the Dzyaloshinskii-Moriya interaction (DMI) interaction specifying the chirality of the structure. Here, the interatomic exchange interaction does not play important role. In UPtGe the hierarchy of the interactions is opposite. The leading interaction is the interatomic exchange interaction responsible for the formation of the incommensurate spiral structure followed by the DMI responsible for the selected chirality of the helix. The magnetic anisotropy is very weak that is a prerequisite for keeping the distortion of the helical structure weak.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Jingchao, E-mail: zhang@unl.edu, E-mail: yyue@whu.edu.cn; Hong, Yang; Yue, Yanan, E-mail: zhang@unl.edu, E-mail: yyue@whu.edu.cn

As the dimensions of nanocircuits and nanoelectronics shrink, thermal energies are being generated in more confined spaces, making it extremely important and urgent to explore for efficient heat dissipation pathways. In this work, the phonon energy transport across graphene and hexagonal boron-nitride (h-BN) interface is studied using classic molecular dynamics simulations. Effects of temperature, interatomic bond strength, heat flux direction, and functionalization on interfacial thermal transport are investigated. It is found out that by hydrogenating graphene in the hybrid structure, the interfacial thermal resistance (R) between graphene and h-BN can be reduced by 76.3%, indicating an effective approach to manipulatemore » the interfacial thermal transport. Improved in-plane/out-of-plane phonon couplings and broadened phonon channels are observed in the hydrogenated graphene system by analyzing its phonon power spectra. The reported R results monotonically decrease with temperature and interatomic bond strengths. No thermal rectification phenomenon is observed in this interfacial thermal transport. Results reported in this work give the fundamental knowledge on graphene and h-BN thermal transport and provide rational guidelines for next generation thermal interface material designs.« less

On macromolecular refinement at subatomic resolution with interatomic scatterers

DOE Office of Scientific and Technical Information (OSTI.GOV)

Afonine, Pavel V., E-mail: pafonine@lbl.gov; Grosse-Kunstleve, Ralf W.; Adams, Paul D.

2007-11-01

Modelling deformation electron density using interatomic scatters is simpler than multipolar methods, produces comparable results at subatomic resolution and can easily be applied to macromolecules. A study of the accurate electron-density distribution in molecular crystals at subatomic resolution (better than ∼1.0 Å) requires more detailed models than those based on independent spherical atoms. A tool that is conventionally used in small-molecule crystallography is the multipolar model. Even at upper resolution limits of 0.8–1.0 Å, the number of experimental data is insufficient for full multipolar model refinement. As an alternative, a simpler model composed of conventional independent spherical atoms augmented bymore » additional scatterers to model bonding effects has been proposed. Refinement of these mixed models for several benchmark data sets gave results that were comparable in quality with the results of multipolar refinement and superior to those for conventional models. Applications to several data sets of both small molecules and macromolecules are shown. These refinements were performed using the general-purpose macromolecular refinement module phenix.refine of the PHENIX package.« less

Deformation of periodic nanovoid structures in Mg single crystals

NASA Astrophysics Data System (ADS)

Xu, Shuozhi; Su, Yanqing; Zare Chavoshi, Saeed

2018-01-01

Large scale molecular dynamics (MD) simulations in Mg single crystal containing periodic cylindrical voids subject to uniaxial tension along the z direction are carried out. Models with different initial void sizes and crystallographic orientations are explored using two interatomic potentials. It is found that (i) a larger initial void always leads to a lower yield stress, in agreement with an analytic prediction; (ii) in the model with x[\\bar{1}100]-y[0001]-z[11\\bar{2}0] orientations, the two potentials predict different types of tension twins and phase transformation; (iii) in the model with x[0001]-y[11\\bar{2}0]-z[\\bar{1}100] orientations, the two potentials identically predict the nucleation of edge dislocations on the prismatic plane, which then glide away from the void, resulting in extrusions at the void surface; in the case of the smallest initial void, these surface extrusions pinch the void into two voids. Besides bringing new physical understanding of the nanovoid structures, our work highlights the variability and uncertainty in MD simulations arising from the interatomic potential, an issue relatively lightly addressed in the literature to date.

Covalent bonds are created by the drive of electron waves to lower their kinetic energy through expansion

PubMed Central

Schmidt, Michael W.; Ivanic, Joseph; Ruedenberg, Klaus

2014-01-01

An analysis based on the variation principle shows that in the molecules H2+, H2, B2, C2, N2, O2, F2, covalent bonding is driven by the attenuation of the kinetic energy that results from the delocalization of the electronic wave function. For molecular geometries around the equilibrium distance, two features of the wave function contribute to this delocalization: (i) Superposition of atomic orbitals extends the electronic wave function from one atom to two or more atoms; (ii) intra-atomic contraction of the atomic orbitals further increases the inter-atomic delocalization. The inter-atomic kinetic energy lowering that (perhaps counter-intuitively) is a consequence of the intra-atomic contractions drives these contractions (which per se would increase the energy). Since the contractions necessarily encompass both, the intra-atomic kinetic and potential energy changes (which add to a positive total), the fact that the intra-atomic potential energy change renders the total potential binding energy negative does not alter the fact that it is the kinetic delocalization energy that drives the bond formation. PMID:24880263

Covalent bonds are created by the drive of electron waves to lower their kinetic energy through expansion.

PubMed

Schmidt, Michael W; Ivanic, Joseph; Ruedenberg, Klaus

2014-05-28

An analysis based on the variation principle shows that in the molecules H2 (+), H2, B2, C2, N2, O2, F2, covalent bonding is driven by the attenuation of the kinetic energy that results from the delocalization of the electronic wave function. For molecular geometries around the equilibrium distance, two features of the wave function contribute to this delocalization: (i) Superposition of atomic orbitals extends the electronic wave function from one atom to two or more atoms; (ii) intra-atomic contraction of the atomic orbitals further increases the inter-atomic delocalization. The inter-atomic kinetic energy lowering that (perhaps counter-intuitively) is a consequence of the intra-atomic contractions drives these contractions (which per se would increase the energy). Since the contractions necessarily encompass both, the intra-atomic kinetic and potential energy changes (which add to a positive total), the fact that the intra-atomic potential energy change renders the total potential binding energy negative does not alter the fact that it is the kinetic delocalization energy that drives the bond formation.

Universal amorphous-amorphous transition in GexSe100-x glasses under pressure

NASA Astrophysics Data System (ADS)

Yildirim, Can; Micoulaut, Matthieu; Boolchand, Punit; Kantor, Innokenty; Mathon, Olivier; Gaspard, Jean-Pierre; Irifune, Tetsuo; Raty, Jean-Yves

2016-06-01

Pressure induced structural modifications in vitreous GexSe100-x (where 10 ≤ x ≤ 25) are investigated using X-ray absorption spectroscopy (XAS) along with supplementary X-ray diffraction (XRD) experiments and ab initio molecular dynamics (AIMD) simulations. Universal changes in distances and angle distributions are observed when scaled to reduced densities. All compositions are observed to remain amorphous under pressure values up to 42 GPa. The Ge-Se interatomic distances extracted from XAS data show a two-step response to the applied pressure; a gradual decrease followed by an increase at around 15-20 GPa, depending on the composition. This increase is attributed to the metallization event that can be traced with the red shift in Ge K edge energy which is also identified by the principal peak position of the structure factor. The densification mechanisms are studied in details by means of AIMD simulations and compared to the experimental results. The evolution of bond angle distributions, interatomic distances and coordination numbers are examined and lead to similar pressure-induced structural changes for any composition.

Non-patchy strategy for inter-atomic distances from Extended X-ray Absorption Fine Structure

PubMed Central

Xu, Gu; Li, Guifang; LI, Xianya; Liang, Yi; Feng, Zhechuan

2017-01-01

Extended X-ray Absorption Fine Structure (EXAFS) has been one of the few structural probes available for crystalline, non-crystalline and even highly disordered specimens. However, the data analysis involves a patchy and tinkering process, including back-and-forth fitting and filtering, leading to ambiguous answers sometimes. Here we try to resolve this long standing problem, to extract the inter-atomic distances from the experimental data by a single step minimization, in order to replace the tedious and tinkering process. The new strategy is built firmly by the mathematical logic, and made straightforward and undeniable. The finding demonstrates that it is possible to break off from the traditional patchy model fitting, and to remove the logical confusion of a priori prediction of the structure to be matched with experimental data, making it a much more powerful technique than the existing methods. The new method is expected to benefit EXAFS users covering all disciplines. Also, it is anticipated that the current work to be the motivation and inspiration to the further efforts. PMID:28181529

An existence criterion for low-dimensional materials

NASA Astrophysics Data System (ADS)

Chen, Jiapeng; Wang, Biao; Hu, Yangfan

2017-10-01

The discovery of graphene and other two-dimensional (2-D) materials has stimulated a general interest in low-dimensional (low-D) materials. Whereas long time ago, Peierls (1935) and Landau's (1937) theoretical work demonstrated that any one- and two-dimensional materials could not exist in any finite temperature environment. Then, two basic issues became a central concern for many researchers: How can stable low-D materials exist? What kind of low-D materials are stable? Here, we establish an energy stability criterion for low-D materials, which seeks to provide a clear answer to these questions. For a certain kind of element, the stability of its specific low-D structure is determined by several derivatives of its interatomic potential. This atomistic-based approach is then applied to study any straight/planar, low-D, equal-bond-length elemental materials. We found that 1-D monatomic chains, 2-D honeycomb lattices, square lattices, and triangular lattices are the only four permissible structures, and the stability of these structures can only be understood by assuming multi-body interatomic potentials. Using this approach, the stable existence of graphene, silicene and germanene can be explained.

Surface Electrochemistry of Metals

DTIC Science & Technology

1993-04-30

maxima along the 12 directions of open channels .vhich are also the interatomic directions). Elastic scattering angular distributions always contain... scatterer geometric relationships for such samples. Distributions from ordered atomic bilayers reveal that the Auger signal from the underlayer is attenuated...are developing a theoretical model and computational code which include both elastic scattering and inhomogeneous inelastic scattering . We seek

Topological properties of the electron density of solids and molecules. Recent developments in Oviedo.

PubMed

Luaña, Víctor; Costales, Aurora; Mori-Sánchez, Paula; Blanco, Miguel A; Martín Pendás, A

2004-09-01

Some of the latest advances in the analysis of electron density are reviewed, including: (a) topological indices that provide a useful characterization of the global properties of the density; (b) specific results on some prototypical metal and low heteropolarity systems; and (c) calculation of the local curvature of the interatomic surface.

Kinetics of homogeneous nucleation on many-component systems

NASA Technical Reports Server (NTRS)

Hirschfelder, J. O.

1974-01-01

Reiss's (1950) classical treatment of the kinetics of homogeneous nucleation in a system containing two chemical components is extended to many-component systems. The formulation is analogous to the pseudostationary-state theory of chemical reaction rates, with the free energy as a function of the composition of the embryo taking the place of the potential energy as a function of interatomic distances.

Kinetics of homogeneous nucleation in many component systems

NASA Technical Reports Server (NTRS)

Hirschfelder, J. O.

1974-01-01

Reiss's classical treatment of the kinetics of homogeneous nucleation in a system containing two chemical components is extended to many-component systems. The formulation is analogous to the pseudo-stationary state theory of chemical reaction rates with the free energy as a function of the composition of the embryo taking the place of the potential energy as a function of interatomic distances.

Surface Impact Simulations of Helium Nanodroplets

DTIC Science & Technology

2015-06-30

mechanical delocalization of the individual helium atoms in the droplet and the quan- tum statistical effects that accompany the interchange of identical...incorporates the effects of atomic delocaliza- tion by treating individual atoms as smeared-out probability distributions that move along classical...probability density distributions to give effec- tive interatomic potential energy curves that have zero-point averaging effects built into them [25

Interatomic potentials for HeAr, HeKr, and HeXe from multiproperty fits

DOE Office of Scientific and Technical Information (OSTI.GOV)

Danielson, L.J.; Keil, M.

1988-01-15

Crossed molecular beam measurements of differential cross sections (DCS) are reported for elastic scattering of He by Ar, Kr, and Xe at high resolution. Interatomic potentials are determined by simultaneously fitting the DCS's, as well as mixture viscosity and interaction second virial data. Bias due to systematic and potential model errors are examined and are used to estimate the accuracy of the potential energy curves obtained. Attractive well depths are 2.59, 2.67, and 2.64 meV +- 3% for HeAr, HeKr, and HeXe, respectively, agreeing with the best available HeAr potential and a previously proposed HeKr potential, but significantly deeper thanmore » previously reported potentials for HeXe. The HeXe attractive well is also considerably broader than previously reported. Attractive minimum positions are 3.48, 3.70, and 4.00 A ( +- 0.03 A) for HeAr, HeKr, and HeXe, respectively. Including the accurate diffusion data of Dunlop and co-workers (Physica A 95, 561 (1979)) and the absolute integral cross sections of Pirani and Vecchiocattivi (J. Chem. Phys. 66, 372 (1977) and revisions thereto) verify the error bounds for all three potentials.« less

Interatomic interaction effects on second-order momentum correlations and Hong-Ou-Mandel interference of double-well-trapped ultracold fermionic atoms

NASA Astrophysics Data System (ADS)

Brandt, Benedikt B.; Yannouleas, Constantine; Landman, Uzi

2018-05-01

Identification and understanding of the evolution of interference patterns in two-particle momentum correlations as a function of the strength of interatomic interactions are important in explorations of the nature of quantum states of trapped particles. Together with the analysis of two-particle spatial correlations, they offer the prospect of uncovering fundamental symmetries and structure of correlated many-body states, as well as opening vistas into potential control and utilization of correlated quantum states as quantum-information resources. With the use of the second-order density matrix constructed via exact diagonalization of the microscopic Hamiltonian, and an analytic Hubbard-type model, we explore here the systematic evolution of characteristic interference patterns in the two-body momentum and spatial correlation maps of two entangled ultracold fermionic atoms in a double well, for the entire attractive- and repulsive-interaction range. We uncover quantum-statistics-governed bunching and antibunching, as well as interaction-dependent interference patterns, in the ground and excited states, and interpret our results in light of the Hong-Ou-Mandel interference physics, widely exploited in photon indistinguishability testing and quantum-information science.

Size-dependent plastic deformation of twinned nanopillars in body-centered cubic tungsten

NASA Astrophysics Data System (ADS)

Xu, Shuozhi; Startt, Jacob K.; Payne, Thomas G.; Deo, Chaitanya S.; McDowell, David L.

2017-05-01

Compared with face-centered cubic metals, twinned nanopillars in body-centered cubic (BCC) systems are much less explored partly due to the more complicated plastic deformation behavior and a lack of reliable interatomic potentials for the latter. In this paper, the fault energies predicted by two semi-empirical interatomic potentials in BCC tungsten (W) are first benchmarked against density functional theory calculations. Then, the more accurate potential is employed in large scale molecular dynamics simulations of tensile and compressive loading of twinned nanopillars in BCC W with different cross sectional shapes and sizes. A single crystal, a twinned crystal, and single crystalline nanopillars are also studied as references. Analyses of the stress-strain response and defect nucleation reveal a strong tension-compression asymmetry and a weak pillar size dependence in the yield strength. Under both tensile and compressive loading, plastic deformation in the twinned nanopillars is dominated by dislocation slip on {110} planes that are nucleated from the intersections between the twin boundary and the pillar surface. It is also found that the cross sectional shape of nanopillars affects the strength and the initial site of defect nucleation but not the overall stress-strain response and plastic deformation behavior.

Thermal conductivity of silicon using reverse non-equilibrium molecular dynamics

NASA Astrophysics Data System (ADS)

El-Genk, Mohamed S.; Talaat, Khaled; Cowen, Benjamin J.

2018-05-01

Simulations are performed using the reverse non-equilibrium molecular dynamics (rNEMD) method and the Stillinger-Weber (SW) potential to determine the input parameters for achieving ±1% convergence of the calculated thermal conductivity of silicon. These parameters are then used to investigate the effects of the interatomic potentials of SW, Tersoff II, Environment Dependent Interatomic Potential (EDIP), Second Nearest Neighbor, Modified Embedded-Atom Method (MEAM), and Highly Optimized Empirical Potential MEAM on determining the bulk thermal conductivity as a function of temperature (400-1000 K). At temperatures > 400 K, data collection and swap periods of 15 ns and 150 fs, system size ≥6 × 6 UC2 and system lengths ≥192 UC are adequate for ±1% convergence with all potentials, regardless of the time step size (0.1-0.5 fs). This is also true at 400 K, except for the SW potential, which requires a data collection period ≥30 ns. The calculated bulk thermal conductivities using the rNEMD method and the EDIP potential are close to, but lower than experimental values. The 10% difference at 400 K increases gradually to 20% at 1000 K.

Atomistic properties of γ uranium.

PubMed

Beeler, Benjamin; Deo, Chaitanya; Baskes, Michael; Okuniewski, Maria

2012-02-22

The properties of the body-centered cubic γ phase of uranium (U) are calculated using atomistic simulations. First, a modified embedded-atom method interatomic potential is developed for the high temperature body-centered cubic (γ) phase of U. This phase is stable only at high temperatures and is thus relatively inaccessible to first principles calculations and room temperature experiments. Using this potential, equilibrium volume and elastic constants are calculated at 0 K and found to be in close agreement with previous first principles calculations. Further, the melting point, heat capacity, enthalpy of fusion, thermal expansion and volume change upon melting are calculated and found to be in reasonable agreement with experiment. The low temperature mechanical instability of γ U is correctly predicted and investigated as a function of pressure. The mechanical instability is suppressed at pressures greater than 17.2 GPa. The vacancy formation energy is analyzed as a function of pressure and shows a linear trend, allowing for the calculation of the extrapolated zero pressure vacancy formation energy. Finally, the self-defect formation energy is analyzed as a function of temperature. This is the first atomistic calculation of γ U properties above 0 K with interatomic potentials.

Atomistic properties of γ uranium

NASA Astrophysics Data System (ADS)

Beeler, Benjamin; Deo, Chaitanya; Baskes, Michael; Okuniewski, Maria

2012-02-01

The properties of the body-centered cubic γ phase of uranium (U) are calculated using atomistic simulations. First, a modified embedded-atom method interatomic potential is developed for the high temperature body-centered cubic (γ) phase of U. This phase is stable only at high temperatures and is thus relatively inaccessible to first principles calculations and room temperature experiments. Using this potential, equilibrium volume and elastic constants are calculated at 0 K and found to be in close agreement with previous first principles calculations. Further, the melting point, heat capacity, enthalpy of fusion, thermal expansion and volume change upon melting are calculated and found to be in reasonable agreement with experiment. The low temperature mechanical instability of γ U is correctly predicted and investigated as a function of pressure. The mechanical instability is suppressed at pressures greater than 17.2 GPa. The vacancy formation energy is analyzed as a function of pressure and shows a linear trend, allowing for the calculation of the extrapolated zero pressure vacancy formation energy. Finally, the self-defect formation energy is analyzed as a function of temperature. This is the first atomistic calculation of γ U properties above 0 K with interatomic potentials.

Development of a machine learning potential for graphene

NASA Astrophysics Data System (ADS)

Rowe, Patrick; Csányi, Gábor; Alfè, Dario; Michaelides, Angelos

2018-02-01

We present an accurate interatomic potential for graphene, constructed using the Gaussian approximation potential (GAP) machine learning methodology. This GAP model obtains a faithful representation of a density functional theory (DFT) potential energy surface, facilitating highly accurate (approaching the accuracy of ab initio methods) molecular dynamics simulations. This is achieved at a computational cost which is orders of magnitude lower than that of comparable calculations which directly invoke electronic structure methods. We evaluate the accuracy of our machine learning model alongside that of a number of popular empirical and bond-order potentials, using both experimental and ab initio data as references. We find that whilst significant discrepancies exist between the empirical interatomic potentials and the reference data—and amongst the empirical potentials themselves—the machine learning model introduced here provides exemplary performance in all of the tested areas. The calculated properties include: graphene phonon dispersion curves at 0 K (which we predict with sub-meV accuracy), phonon spectra at finite temperature, in-plane thermal expansion up to 2500 K as compared to NPT ab initio molecular dynamics simulations and a comparison of the thermally induced dispersion of graphene Raman bands to experimental observations. We have made our potential freely available online at [http://www.libatoms.org].

Evolution of sp{sup 2} networks with substrate temperature in amorphous carbon films: Experiment and theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gago, R.; Centro de Micro-Analisis de Materiales, Universidad Autonoma de Madrid, Campus de Cantoblanco, 28049 Madrid; Vinnichenko, M.

2005-07-01

The evolution of sp{sup 2} hybrids in amorphous carbon (a-C) films deposited at different substrate temperatures was studied experimentally and theoretically. The bonding structure of a-C films prepared by filtered cathodic vacuum arc was assessed by the combination of visible Raman spectroscopy, x-ray absorption, and spectroscopic ellipsometry, while a-C structures were generated by molecular-dynamics deposition simulations with the Brenner interatomic potential to determine theoretical sp{sup 2} site distributions. The experimental results show a transition from tetrahedral a-C (ta-C) to sp{sup 2}-rich structures at {approx}500 K. The sp{sup 2} hybrids are mainly arranged in chains or pairs whereas graphitic structures aremore » only promoted for sp{sup 2} fractions above 80%. The theoretical analysis confirms the preferred pairing of isolated sp{sup 2} sites in ta-C, the coalescence of sp{sup 2} clusters for medium sp{sup 2} fractions, and the pronounced formation of rings for sp{sup 2} fractions >80%. However, the dominance of sixfold rings is not reproduced theoretically, probably related to the functional form of the interatomic potential used.« less

Molecular dynamics simulations of intergranular fracture in UO2 with nine empirical interatomic potentials

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yongfeng Zhang; Paul C Millett; Michael R Tonks

The intergranular fracture behavior of UO2 was studied using molecular dynamics simulations with a bicrystal model. The anisotropic fracture behavior due to the different grain boundary characters was investigated with the View the MathML source symmetrical tilt S5 and the View the MathML source symmetrical tilt S3 ({1 1 1} twin) grain boundaries. Nine interatomic potentials, seven rigid-ion plus two core–shell ones, were utilized to elucidate possible potential dependence. Initiating from a notch, crack propagation along grain boundaries was observed for most potentials. The S3 boundary was found to be more prone to fracture than the S5 one, indicated bymore » a lower energy release rate associated with the former. However, some potential dependence was identified on the existence of transient plastic deformation at crack tips, and the results were discussed regarding the relevant material properties including the excess energies of metastable phases and the critical energy release rate for intergranular fracture. In general, local plasticity at crack tips was observed in fracture simulations with potentials that predict low excess energies for metastable phases and high critical energy release rates for intergranular fracture.« less

A modified Stillinger-Weber potential for TlBr and its polymorphic extension

DOE PAGES

Zhou, Xiaowang; Foster, Michael E.; Jones, Reese E.; ...

2015-04-30

TlBr is promising for g- and x- radiation detection, but suffers from rapid performance degradation under the operating external electric fields. To enable molecular dynamics (MD) studies of this degradation, we have developed a Stillinger-Weber type of TlBr interatomic potential. During this process, we have also addressed two problems of wider interests. First, the conventional Stillinger-Weber potential format is only applicable for tetrahedral structures (e.g., diamond-cubic, zinc-blende, or wurtzite). Here we have modified the analytical functions of the Stillinger-Weber potential so that it can now be used for other crystal structures. Second, past modifications of interatomic potentials cannot always bemore » applied by a broad community because any new analytical functions of the potential would require corresponding changes in the molecular dynamics codes. Here we have developed a polymorphic potential model that simultaneously incorporates Stillinger-Weber, Tersoff, embedded-atom method, and any variations (i.e., modified functions) of these potentials. As a result, we have implemented this polymorphic model in MD code LAMMPS, and demonstrated that our TlBr potential enables stable MD simulations under external electric fields.« less

On the transferability of electron density in binary vanadium borides VB, V3B4 and VB2.

PubMed

Terlan, Bürgehan; Akselrud, Lev; Baranov, Alexey I; Borrmann, Horst; Grin, Yuri

2015-12-01

Binary vanadium borides are suitable model systems for a systematic analysis of the transferability concept in intermetallic compounds due to chemical intergrowth in their crystal structures. In order to underline this structural relationship, topological properties of the electron density in VB, V3B4 and VB2 reconstructed from high-resolution single-crystal X-ray diffraction data as well as derived from quantum chemical calculations, are analysed in terms of Bader's Quantum Theory of Atoms in Molecules [Bader (1990). Atoms in Molecules: A Quantum Theory, 1st ed. Oxford: Clarendon Press]. The compounds VB, V3B4 and VB2 are characterized by a charge transfer from the metal to boron together with two predominant atomic interactions, the shared covalent B-B interactions and the polar covalent B-M interactions. The resembling features of the crystal structures are well reflected by the respective B-B interatomic distances as well as by ρ(r) values at the B-B bond critical points. The latter decrease with an increase in the corresponding interatomic distances. The B-B bonds show transferable electron density properties at bond critical points depending on the respective bond distances.

A general method for the derivation of the functional forms of the effective energy terms in coarse-grained energy functions of polymers. I. Backbone potentials of coarse-grained polypeptide chains

NASA Astrophysics Data System (ADS)

Sieradzan, Adam K.; Makowski, Mariusz; Augustynowicz, Antoni; Liwo, Adam

2017-03-01

A general and systematic method for the derivation of the functional expressions for the effective energy terms in coarse-grained force fields of polymer chains is proposed. The method is based on the expansion of the potential of mean force of the system studied in the cluster-cumulant series and expanding the all-atom energy in the Taylor series in the squares of interatomic distances about the squares of the distances between coarse-grained centers, to obtain approximate analytical expressions for the cluster cumulants. The primary degrees of freedom to average about are the angles for collective rotation of the atoms contained in the coarse-grained interaction sites about the respective virtual-bond axes. The approach has been applied to the revision of the virtual-bond-angle, virtual-bond-torsional, and backbone-local-and-electrostatic correlation potentials for the UNited RESidue (UNRES) model of polypeptide chains, demonstrating the strong dependence of the torsional and correlation potentials on virtual-bond angles, not considered in the current UNRES. The theoretical considerations are illustrated with the potentials calculated from the ab initio potential-energy surface of terminally blocked alanine by numerical integration and with the statistical potentials derived from known protein structures. The revised torsional potentials correctly indicate that virtual-bond angles close to 90° result in the preference for the turn and helical structures, while large virtual-bond angles result in the preference for polyproline II and extended backbone geometry. The revised correlation potentials correctly reproduce the preference for the formation of β-sheet structures for large values of virtual-bond angles and for the formation of α-helical structures for virtual-bond angles close to 90°.

Electronic band structure effects in the stopping of protons in copper [Electronic band structure non-linear effects in the stopping of protons in copper

DOE PAGES

Quashie, Edwin E.; Saha, Bidhan C.; Correa, Alfredo A.

2016-10-05

Here, we present an ab initio study of the electronic stopping power of protons in copper over a wide range of proton velocities v = 0.02–10a.u. where we take into account nonlinear effects. Time-dependent density functional theory coupled with molecular dynamics is used to study electronic excitations produced by energetic protons. A plane-wave pseudopotential scheme is employed to solve the time-dependent Kohn-Sham equations for a moving ion in a periodic crystal. The electronic excitations and the band structure determine the stopping power of the material and alter the interatomic forces for both channeling and off-channeling trajectories. Our off-channeling results aremore » in quantitative agreement with experiments, and at low velocity they unveil a crossover region of superlinear velocity dependence (with a power of ~1.5) in the velocity range v = 0.07–0.3a.u., which we associate to the copper crystalline electronic band structure. The results are rationalized by simple band models connecting two separate regimes. We find that the limit of electronic stopping v → 0 is not as simple as phenomenological models suggest and it is plagued by band-structure effects.« less

Lattice dynamics and metastability of fcc metals in the hcp structure and the crucial role of spin-orbit coupling in platinum

NASA Astrophysics Data System (ADS)

Schönecker, Stephan; Li, Xiaoqing; Richter, Manuel; Vitos, Levente

2018-06-01

We investigate the lattice dynamical properties of Ni, Cu, Rh, Pd, Ag, Ir, Pt, and Au in the nonequilibrium hcp structure by means of density-functional simulations, wherein spin-orbit coupling (SOC) was considered for Ir, Pt, and Au. The determined dynamical properties reveal that all eight elements possess a metastable hcp phase at zero temperature and pressure. The hcp Ni, Cu, Rh, Pd, and Au previously observed in nanostructures support this finding. We make evident that the inclusion of SOC is mandatory for an accurate description of the phonon dispersion relations and dynamical stability of hcp Pt. The underlying sensitivity of the interatomic force constants is ascribed to a SOC-induced splitting of degenerate band states accompanied by a pronounced reduction of electronic density of states at the Fermi level. To give further insight into the importance of SOC in Pt, we (i) focus on phase stability and examine a lattice transformation related to optical phonons in the hcp phase and (ii) focus on the generalized stacking fault energy (GSFE) of the fcc phase pertinent to crystal plasticity. We show that the intrinsic stable and unstable fault energies of the GSFE scale as in other common fcc metals, provided that the spin-orbit interaction is taken into account.

Molecular dynamics simulations reveal the conformational dynamics of Arabidopsis thaliana BRI1 and BAK1 receptor-like kinases.

PubMed

Moffett, Alexander S; Bender, Kyle W; Huber, Steven C; Shukla, Diwakar

2017-07-28

The structural motifs responsible for activation and regulation of eukaryotic protein kinases in animals have been studied extensively in recent years, and a coherent picture of their activation mechanisms has begun to emerge. In contrast, non-animal eukaryotic protein kinases are not as well understood from a structural perspective, representing a large knowledge gap. To this end, we investigated the conformational dynamics of two key Arabidopsis thaliana receptor-like kinases, brassinosteroid-insensitive 1 (BRI1) and BRI1-associated kinase 1 (BAK1), through extensive molecular dynamics simulations of their fully phosphorylated kinase domains. Molecular dynamics simulations calculate the motion of each atom in a protein based on classical approximations of interatomic forces, giving researchers insight into protein function at unparalleled spatial and temporal resolutions. We found that in an otherwise "active" BAK1 the αC helix is highly disordered, a hallmark of deactivation, whereas the BRI1 αC helix is moderately disordered and displays swinging behavior similar to numerous animal kinases. An analysis of all known sequences in the A. thaliana kinome found that αC helix disorder may be a common feature of plant kinases. © 2017 by The American Society for Biochemistry and Molecular Biology, Inc.

Electron-phonon coupling and thermal transport in the thermoelectric compound Mo 3Sb 7–xTe x

DOE PAGES

Bansal, Dipanshu; Li, Chen W.; Said, Ayman H.; ...

2015-12-07

Phonon properties of Mo 3Sb 7–xTe x (x = 0, 1.5, 1.7), a potential high-temperature thermoelectric material, have been studied with inelastic neutron and x-ray scattering, and with first-principles simulations. The substitution of Te for Sb leads to pronounced changes in the electronic struc- ture, local bonding, phonon density of states (DOS), dispersions, and phonon lifetimes. Alloying with tellurium shifts the Fermi level upward, near the top of the valence band, resulting in a strong suppression of electron-phonon screening, and a large overall stiffening of interatomic force- constants. The suppression in electron-phonon coupling concomitantly increases group velocities and suppresses phononmore » scattering rates, surpassing the effects of alloy-disorder scattering, and re- sulting in a surprising increased lattice thermal conductivity in the alloy. We also identify that the local bonding environment changes non-uniformly around different atoms, leading to variable perturbation strengths for different optical phonon branches. The respective roles of changes in phonon group velocities and phonon lifetimes on the lattice thermal conductivity are quantified. Lastly, our results highlight the importance of the electron-phonon coupling on phonon mean-free-paths in this compound, and also estimates the contributions from boundary scattering, umklapp scattering, and point-defect scattering.« less

Thin films of fullerene-like MoS2 nanoparticles with ultra-low friction and wear

PubMed

Chhowalla; Amaratunga

2000-09-14

The tribological properties of solid lubricants such as graphite and the metal dichalcogenides MX2 (where M is molybdenum or tungsten and X is sulphur or selenium) are of technological interest for reducing wear in circumstances where liquid lubricants are impractical, such as in space technology, ultra-high vacuum or automotive transport. These materials are characterized by weak interatomic interactions (van der Waals forces) between their layered structures, allowing easy, low-strength shearing. Although these materials exhibit excellent friction and wear resistance and extended lifetime in vacuum, their tribological properties remain poor in the presence of humidity or oxygen, thereby limiting their technological applications in the Earth's atmosphere. But using MX2 in the form of isolated inorganic fullerene-like hollow nanoparticles similar to carbon fullerenes and nanotubes can improve its performance. Here we show that thin films of hollow MoS2 nanoparticles, deposited by a localized high-pressure arc discharge method, exhibit ultra-low friction (an order of magnitude lower than for sputtered MoS2 thin films) and wear in nitrogen and 45% humidity. We attribute this 'dry' behaviour in humid environments to the presence of curved S-Mo-S planes that prevent oxidation and preserve the layered structure.

Phonon-defect scattering and thermal transport in semiconductors: developing guiding principles

NASA Astrophysics Data System (ADS)

Polanco, Carlos; Lindsay, Lucas

First principles calculations of thermal conductivity have shown remarkable agreement with measurements for high-quality crystals. Nevertheless, most materials contain defects that provide significant extrinsic resistance and lower the conductivity from that of a perfect sample. This effect is usually accounted for with simplified analytical models that neglect the atomistic details of the defect and the exact dynamical properties of the system, which limits prediction capabilities. Recently, a method based on Greens functions was developed to calculate the phonon-defect scattering rates from first principles. This method has shown the important role of point defects in determining thermal transport in diamond and boron arsenide, two competitors for the highest bulk thermal conductivity. Here, we study the role of point defects on other relatively high thermal conductivity semiconductors, e.g., BN, BeSe, SiC, GaN and Si. We compare their first principles defect-phonon scattering rates and effects on transport properties with those from simplified models and explore common principles that determine these. Efforts will focus on basic vibrational properties that vary from system to system, such as density of states, interatomic force constants and defect deformation. Research supported by the U.S. Department of Energy, Basic Energy Sciences, Materials Sciences and Engineering Division.

Thin films of fullerene-like MoS2 nanoparticles with ultra-low friction and wear

NASA Astrophysics Data System (ADS)

Chhowalla, Manish; Amaratunga, Gehan A. J.

2000-09-01

The tribological properties of solid lubricants such as graphite and the metal dichalcogenides MX2 (where M is molybdenum or tungsten and X is sulphur or selenium) are of technological interest for reducing wear in circumstances where liquid lubricants are impractical, such as in space technology, ultra-high vacuum or automotive transport. These materials are characterized by weak interatomic interactions (van der Waals forces) between their layered structures, allowing easy, low-strength shearing. Although these materials exhibit excellent friction and wear resistance and extended lifetime in vacuum, their tribological properties remain poor in the presence of humidity or oxygen, thereby limiting their technological applications in the Earth's atmosphere. But using MX2 in the form of isolated inorganic fullerene-like hollow nanoparticles similar to carbon fullerenes and nanotubes can improve its performance. Here we show that thin films of hollow MoS2 nanoparticles, deposited by a localized high-pressure arc discharge method, exhibit ultra-low friction (an order of magnitude lower than for sputtered MoS2 thin films) and wear in nitrogen and 45% humidity. We attribute this `dry' behaviour in humid environments to the presence of curved S-Mo-S planes that prevent oxidation and preserve the layered structure.

Heavy-impurity resonance, hybridization, and phonon spectral functions in Fe 1-xM xSi, M=Ir,Os

DOE PAGES

Delaire, O.; Al-Qasir, Iyad I.; May, Andrew F.; ...

2015-03-31

The vibrational behavior of heavy substitutional impurities (M=Ir,Os) in Fe 1-xM xSi (x = 0, 0.02, 0.04, 0.1) was investigated with a combination of inelastic neutron scattering (INS), transport measurements, and first-principles simulations. In this paper, our INS measurements on single-crystals mapped the four-dimensional dynamical structure factor, S(Q;E), for several compositions and temperatures. Our results show that both Ir and Os impurities lead to the formation of a weakly dispersive resonance vibrational mode, in the energy range of the acoustic phonon dispersions of the FeSi host. We also show that Ir doping, which introduces free carriers and increases electron-phonon coupling,more » leads to softened interatomic force-constants compared to doping with Os, which is isoelectronic to Fe. We analyze the phonon S(Q,E) from INS through a Green's function model incorporating the phonon self-energy based on first-principles density functional theory (DFT) simulations. Calculations of the quasiparticle spectral functions in the doped system reveal the hybridization between the resonance and the acoustic phonon modes. Finally, our results demonstrate a strong interaction of the host acoustic dispersions with the resonance mode, likely leading to the large observed suppression in lattice thermal conductivity.« less

Exploring oxidative ageing behaviour of hydrocarbons using ab initio molecular dynamics analysis

NASA Astrophysics Data System (ADS)

Pan, Tongyan; Cheng, Cheng

2016-06-01

With a proper approximate solution to the Schrödinger Equation of a multi-electron system, the method of ab initio molecular dynamics (AIMD) performs first-principles molecular dynamics analysis without pre-defining interatomic potentials as are mandatory in traditional molecular dynamics analyses. The objective of this study is to determine the oxidative-ageing pathway of petroleum asphalt as a typical hydrocarbon system, using the AIMD method. This objective was accomplished in three steps, including (1) identifying a group of representative asphalt molecules to model, (2) determining an atomistic modelling method that can effectively simulate the production of critical functional groups in oxidative ageing of hydrocarbons and (3) evaluating the oxidative-ageing pathway of a hydrocarbon system. The determination of oxidative-ageing pathway of hydrocarbons was done by tracking the generations of critical functional groups in the course of oxidative ageing. The chemical elements of carbon, nitrogen and sulphur all experience oxidative reactions, producing polarised functional groups such as ketones, aldehydes or carboxylic acids, pyrrolic groups and sulphoxides. The electrostatic forces of the polarised groups generated in oxidation are responsible for the behaviour of aged hydrocarbons. The developed AIMD model can be used for modelling the ageing of generic hydrocarbon polymers and developing antioxidants without running expensive experiments.

The NRL Program on Electroactive Polymers.

DTIC Science & Technology

1980-09-15

cell of a point in an aggregate involves selecting the smallest cell formed by planes perpendicularly bisecting all the point to neighbor vectors . Such...plane perpendicular to the interatomic vector is located nearer the smaller atom by bisecting the distance between the sur- faces of spheres whose...density waves (and consequent novel excitations such as solitons (6)). The physical structure as well as the chemical bonding of such polymeric

Quantum chemical calculations of interatomic potentials for computer simulation of solids

NASA Technical Reports Server (NTRS)

1977-01-01

A comprehensive mathematical model by which the collective behavior of a very large number of atoms within a metal or alloy can accurately be simulated was developed. Work was done in order to predict and modify the strength of materials to suit our technological needs. The method developed is useful in studying atomic interactions related to dislocation motion and crack extension.

Mathematical Research in Materials Science: Opportunities and Perspectives. Part 2

DTIC Science & Technology

1993-01-01

spheres and Lennard - Jones potentials , but have not been extended to a general framework that will allow input from more complicated interatomic...focuses on directions for potentially promising collaboration between materials scientists and mathematical scientists, and encourages both communities...interface between the mathematical sciences and other fields. The purpose of this report is not only to focus on directions for potentially promising

Orbital Exponent Optimization in Elementary VB Calculations of the Chemical Bond in the Ground State of Simple Molecular Systems

ERIC Educational Resources Information Center

Magnasco, Valerio

2008-01-01

Orbital exponent optimization in the elementary ab-initio VB calculation of the ground states of H[subscript 2][superscript +], H[subscript 2], He[subscript 2][superscript +], He[subscript 2] gives a fair description of the exchange-overlap component of the interatomic interaction that is important in the bond region. Correct bond lengths and…

An Investigation of the Interatomic Bonding Characteristics of a Ti - 51at.% Al Alloy by X-Ray Diffraction

DTIC Science & Technology

1991-06-01

GROUP SUBGROUP X-ray Diffraction, XRD, TiAI, titanium , aluminum, bonding characteristics, titanium aluminides , Debye-Waller temperature factor...XRD Powder Particles (575X) .............. 47 viii I. INTRODUCTION Titanium aluminides are recognized for their high specific strength, particularly at...bonding characteristics of binary titanium aluminides . Upon the introduction of a third element to the system, a rearrangement of the valence

Development of an interatomic potential for the simulation of defects, plasticity, and phase transformations in titanium

DOE PAGES

Mendelev, M. I.; Underwood, T. L.; Ackland, G. J.

2016-10-17

New interatomic potentials describing defects, plasticity, and high temperature phase transitions for Ti are presented. Fitting the martensitic hcp-bcc phase transformation temperature requires an efficient and accurate method to determine it. We apply a molecular dynamics method based on determination of the melting temperature of competing solid phases, and Gibbs-Helmholtz integration, and a lattice-switch Monte Carlo method: these agree on the hcp-bcc transformation temperatures to within 2 K. We were able to develop embedded atom potentials which give a good fit to either low or high temperature data, but not both. The first developed potential (Ti1) reproduces the hcp-bcc transformationmore » and melting temperatures and is suitable for the simulation of phase transitions and bcc Ti. Two other potentials (Ti2 and Ti3) correctly describe defect properties and can be used to simulate plasticity or radiation damage in hcp Ti. The fact that a single embedded atom method potential cannot describe both low and high temperature phases may be attributed to neglect of electronic degrees of freedom, notably bcc has a much higher electronic entropy. As a result, a temperature-dependent potential obtained from the combination of potentials Ti1 and Ti2 may be used to simulate Ti properties at any temperature.« less

Facet-Dependent Cr(VI) Adsorption of Hematite Nanocrystals.

PubMed

Huang, Xiaopeng; Hou, Xiaojing; Song, Fahui; Zhao, Jincai; Zhang, Lizhi

2016-02-16

In this study, the adsorption process of Cr(VI) on the hematite facets was systematically investigated with synchrotron-based Cr K-edge extended X-ray absorption fine structure (EXAFS) spectroscopy, in situ attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy, density-functional theory calculation, and surface complexation models. Structural model fitting of EXAFS spectroscopy suggested that the interatomic distances of Cr-Fe were, respectively, 3.61 Å for the chromate coordinated hematite nanoplates with exposed {001} facets, 3.60 and 3.30 Å for the chromate coordinated hematite nanorods with exposed {001} and {110} facets, which were characteristic of inner-sphere complexation. In situ ATR-FTIR spectroscopy analysis confirmed the presence of two inner-sphere surface complexes with C3ν and C2ν symmetry, while the C3ν and C2ν species were assigned to monodentate and bidentate inner-sphere surface complexes with average Cr-Fe interatomic distances of 3.60 and 3.30 Å, respectively. On the basis of these experimental and theoretical results, we concluded that HCrO4(-) as dominated Cr(VI) species was adsorbed on {001} and {110} facets in inner-sphere monodentate mononuclear and bidentate binuclear configurations, respectively. Moreover, the Cr(VI) adsorption performance of hematite facets was strongly dependent on the chromate complexes formed on the hematite facets.

Energy of the Isolated Metastable Iron-Nickel FCC Nanocluster with a Carbon Atom in the Tetragonal Interstice.

PubMed

Bondarenko, Natalya V; Nedolya, Anatoliy V

2017-12-01

The energy of the isolated iron-nickel nanocluster was calculated by molecular mechanics method using Lennard-Jones potential. The cluster included a carbon atom that drifted from an inside octahedral interstice to a tetrahedral interstice in [Formula: see text] direction and after that in <222> direction to the surface. In addition, one of 14 iron atoms was replaced by a nickel atom, the position of which was changing during simulation.The energy of the nanocluster was estimated at the different interatomic distances. As a result of simulation, the optimal interatomic distances of Fe-Ni-C nanocluster was chosen for the simulation, in which height of the potential barrier was maximal and face-centered cubic (FCC) nanocluster was the most stable.It is shown that there were three main positions of a nickel atom that significantly affected nanocluster's energy.The calculation results indicated that position of the carbon atom in the octahedral interstice was more energetically favorable than tetrahedral interstice in the case of FCC nanocluster. On the other side, the potential barrier was smaller in the direction [Formula: see text] than in the direction <022>.This indicates that there are two ways for carbon atom to drift to the surface of the nanocluster.

Automated generation of quantum-accurate classical interatomic potentials for metals and semiconductors

NASA Astrophysics Data System (ADS)

Thompson, Aidan; Foiles, Stephen; Schultz, Peter; Swiler, Laura; Trott, Christian; Tucker, Garritt

2013-03-01

Molecular dynamics (MD) is a powerful condensed matter simulation tool for bridging between macroscopic continuum models and quantum models (QM) treating a few hundred atoms, but is limited by the accuracy of available interatomic potentials. Sound physical and chemical understanding of these interactions have resulted in a variety of concise potentials for certain systems, but it is difficult to extend them to new materials and properties. The growing availability of large QM data sets has made it possible to use more automated machine-learning approaches. Bartók et al. demonstrated that the bispectrum of the local neighbor density provides good regression surrogates for QM models. We adopt a similar bispectrum representation within a linear regression scheme. We have produced potentials for silicon and tantalum, and we are currently extending the method to III-V compounds. Results will be presented demonstrating the accuracy of these potentials relative to the training data, as well as their ability to accurately predict material properties not explicitly included in the training data. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Dept. of Energy Nat. Nuclear Security Admin. under Contract DE-AC04-94AL85000.

Molecular dynamics modeling of atomic displacement cascades in 3C-SiC: Comparison of interatomic potentials

DOE Office of Scientific and Technical Information (OSTI.GOV)

Samolyuk, German D.; Osetskiy, Yury N.; Stoller, Roger E.

We used molecular dynamics modeling of atomic displacement cascades to characterize the nature of primary radiation damage in 3C-SiC. We demonstrated that the most commonly used interatomic potentials are inconsistent with ab initio calculations of defect energetics. Both the Tersoff potential used in this work and a modified embedded-atom method potential reveal a barrier to recombination of the carbon interstitial and carbon vacancy which is much higher than the density functional theory (DFT) results. The barrier obtained with a newer potential by Gao and Weber is closer to the DFT result. This difference results in significant differences in the cascademore » production of point defects. We have completed both 10 keV and 50 keV cascade simulations in 3C-SiC at a range of temperatures. In contrast to the Tersoff potential, the Gao-Weber potential produces almost twice as many C vacancies and interstitials at the time of maximum disorder (~0.2 ps) but only about 25% more stable defects at the end of the simulation. Only about 20% of the carbon defects produced with the Tersoff potential recombine during the in-cascade annealing phase, while about 60% recombine with the Gao-Weber potential.« less

Molecular dynamics modeling of atomic displacement cascades in 3C-SiC: Comparison of interatomic potentials

DOE PAGES

Samolyuk, German D.; Osetskiy, Yury N.; Stoller, Roger E.

2015-06-03

We used molecular dynamics modeling of atomic displacement cascades to characterize the nature of primary radiation damage in 3C-SiC. We demonstrated that the most commonly used interatomic potentials are inconsistent with ab initio calculations of defect energetics. Both the Tersoff potential used in this work and a modified embedded-atom method potential reveal a barrier to recombination of the carbon interstitial and carbon vacancy which is much higher than the density functional theory (DFT) results. The barrier obtained with a newer potential by Gao and Weber is closer to the DFT result. This difference results in significant differences in the cascademore » production of point defects. We have completed both 10 keV and 50 keV cascade simulations in 3C-SiC at a range of temperatures. In contrast to the Tersoff potential, the Gao-Weber potential produces almost twice as many C vacancies and interstitials at the time of maximum disorder (~0.2 ps) but only about 25% more stable defects at the end of the simulation. Only about 20% of the carbon defects produced with the Tersoff potential recombine during the in-cascade annealing phase, while about 60% recombine with the Gao-Weber potential.« less

Two Stages of Impact Fracture of Polycrystalline ZnS and ZnSe Compounds

NASA Astrophysics Data System (ADS)

Shcherbakov, I. P.; Dunaev, A. A.; Chmel', A. E.

2018-04-01

Mechanoluminescence (ML) in ductile solids is caused by the motion of charged dislocations in the deformable material. Interatomic bond ruptures followed by electronic structure reconfiguration are the main source of ML in brittle bodies. We studied ML in ceramics composed of mixed ionic/covalent ZnS and ZnSe compounds, which are generated during impact loading higher than the limit deformation. Depending on synthesis method and thermal treatment, the resulting ceramics had different size and geometry of grains and intergrain boundary structure, which presumably may have a significant effect on the dislocation glide. In both materials, the time sweeps of ML pulses have two well-resolved peaks. The position of the peaks along the time axis is substantially dependent on the size of ceramic-forming grains and, to a smaller extent, on the barrier properties of intergrain boundaries. The first peak is associated with plastic deformation preceding disintegration of the crystal structure. The second peak emerges upon crack nucleation as interatomic bonds are ruptured and the material is undergoing local deformation in tips of propagating cracks. The distributions of ML pulse amplitudes (the dependences between the number of pulses and their amplitude) calculated for both peaks individually follow the power law, which demonstrates that the electronic processes having different excitation mechanisms (dislocation motion vs bond rupture) are correlated.

Statistical theory of diffusion in concentrated bcc and fcc alloys and concentration dependencies of diffusion coefficients in bcc alloys FeCu, FeMn, FeNi, and FeCr

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vaks, V. G.; Khromov, K. Yu., E-mail: khromov-ky@nrcki.ru; Pankratov, I. R.

2016-07-15

The statistical theory of diffusion in concentrated bcc and fcc alloys with arbitrary pairwise interatomic interactions based on the master equation approach is developed. Vacancy–atom correlations are described using both the second-shell-jump and the nearest-neighbor-jump approximations which are shown to be usually sufficiently accurate. General expressions for Onsager coefficients in terms of microscopic interatomic interactions and some statistical averages are given. Both the analytical kinetic mean-field and the Monte Carlo methods for finding these averages are described. The theory developed is used to describe sharp concentration dependencies of diffusion coefficients in several iron-based alloy systems. For the bcc alloys FeCu,more » FeMn, and FeNi, we predict the notable increase of the iron self-diffusion coefficient with solute concentration c, up to several times, even though values of c possible for these alloys do not exceed some percent. For the bcc alloys FeCr at high temperatures T ≳ 1400 K, we show that the very strong and peculiar concentration dependencies of both tracer and chemical diffusion coefficients observed in these alloys can be naturally explained by the theory, without invoking exotic models discussed earlier.« less

Intrinsic Resolution of Molecular Electronic Wave Functions and Energies in Terms of Quasi-atoms and Their Interactions.

PubMed

West, Aaron C; Schmidt, Michael W; Gordon, Mark S; Ruedenberg, Klaus

2017-02-09

A general intrinsic energy resolution has been formulated for strongly correlated wave functions in the full molecular valence space and its subspaces. The information regarding the quasi-atomic organization of the molecular electronic structure is extracted from the molecular wave function without introducing any additional postulated model state wave functions. To this end, the molecular wave function is expressed in terms of quasi-atomic molecular orbitals, which maximize the overlap between subspaces of the molecular orbital space and the free-atom orbital spaces. As a result, the molecular wave function becomes the superposition of a wave function representing the juxtaposed nonbonded quasi-atoms and a wave function describing the interatomic electron migrations that create bonds through electron sharing. The juxtaposed nonbonded quasi-atoms are shown to consist of entangled quasi-atomic states from different atoms. The binding energy is resolved as a sum of contributions that are due to quasi-atom formation, quasiclassical electrostatic interactions, and interatomic interferences caused by electron sharing. The contributions are further resolved according to orbital interactions. The various transformations that generate the analysis are determined by criteria that are independent of the working orbital basis used for calculating the molecular wave function. The theoretical formulation of the resolution is quantitatively validated by an application to the C 2 molecule.

Born-Oppenheimer Interatomic Forces from Simple, Local Kinetic Energy Density Functionals

NASA Astrophysics Data System (ADS)

Karasiev, V. V.; Trickey, S. B.; Harris, Frank E.

2006-10-01

Rapid calculation of Born-Oppenheimer (B-O) forces is essential for driving the so-called quantum region of a multi-scale molecular dynamics simulation. The success of density functional theory (DFT) with modern exchange-correlation approximations makes DFT an appealing choice for this role. But conventional Kohn-Sham DFT, even with various linear-scaling implementations, really is not fast enough to meet the challenge of complicated chemo-mechanical phenomena (e.g. stress-induced cracking in the presence of a solvent). Moreover, those schemes involve approximations that are difficult to check practically or to validate formally. A popular alternative, Car-Parrinello dynamics, does not guarantee motion on the B-O surface. Another approach, orbital-free DFT, is appealing but has proven difficult to implement because of the challenge of constructing reliable orbital-free (OF) approximations to the kinetic energy (KE) functional. To be maximally useful for multi-scale simulations, an OF-KE functional must be local (i.e. one-point). This requirement eliminates the two-point functionals designed to have proper linear-response behavior in the weakly inhomogeneous limit. In the face of these difficulties, we demonstrate that there is a way forward. By requiring only that the approximate functional deliver high-quality forces, by exploiting the “conjointness” hypothesis of Lee, Lee, and Parr, by enforcing a basic positivity constraint, and by parameterizing to a carefully selected, small set of molecules we are able to generate a KE functional that does a good job of describing various H q Si m O n clusters as well as CO (providing encouraging evidence of transferability). In addition to that positive result, we discuss several major negative results. First is definitive proof that the conjointness hypothesis is not correct, but nevertheless is useful. The second is the failure of a considerable variety of published KE functionals of the generalized gradient approximation type. Those functionals yield no minimum on the energy surface and give completely incorrect forces. In all cases, the problem can be traced to incorrect behavior of the functionals near the nuclei. Third, the seemingly obvious strategy of direct numerical fitting of OF-KE functional parameters to reproduce the energy surface of selected molecules is unsuccessful. The functionals that result are completely untransferable.

Van-der-Waals interaction of atoms in dipolar Rydberg states

NASA Astrophysics Data System (ADS)

Kamenski, Aleksandr A.; Mokhnenko, Sergey N.; Ovsiannikov, Vitaly D.

2018-02-01

An asymptotic expression for the van-der-Waals constant C 6( n) ≈ -0.03 n 12 K p ( x) is derived for the long-range interaction between two highly excited hydrogen atoms A and B in their extreme Stark states of equal principal quantum numbers n A = n B = n ≫ 1 and parabolic quantum numbers n 1(2) = n - 1, n 2(1) = m = 0 in the case of collinear orientation of the Stark-state dipolar electric moments and the interatomic axis. The cubic polynomial K 3( x) in powers of reciprocal values of the principal quantum number x = 1/ n and quadratic polynomial K 2( y) in powers of reciprocal values of the principal quantum number squared y = 1/ n 2 were determined on the basis of the standard curve fitting polynomial procedure from the calculated data for C 6( n). The transformation of attractive van-der-Waals force ( C 6 > 0) for low-energy states n < 23 into repulsive force ( C 6 < 0) for all higher-energy states of n ≥ 23, is observed from the results of numerical calculations based on the second-order perturbation theory for the operator of the long-range interaction between neutral atoms. This transformation is taken into account in the asymptotic formulas (in both cases of p = 2, 3) by polynomials K p tending to unity at n → ∞ ( K p (0) = 1). The transformation from low- n attractive van-der-Waals force into high- n repulsive force demonstrates the gradual increase of the negative contribution to C 6( n) from the lower-energy two-atomic states, of the A(B)-atom principal quantum numbers n'A(B) = n-Δ n (where Δ n = 1, 2, … is significantly smaller than n for the terms providing major contribution to the second-order series), which together with the states of n″B(A) = n+Δ n make the joint contribution proportional to n 12. So, the hydrogen-like manifold structure of the energy spectrum is responsible for the transformation of the power-11 asymptotic dependence C 6( n) ∝ n 11of the low-angular-momenta Rydberg states in many-electron atoms into the power-12 dependence C 6( n) ∝ n 12 for the dipolar states of the Rydberg manifold.

The evolving quality of frictional contact with graphene.

PubMed

Li, Suzhi; Li, Qunyang; Carpick, Robert W; Gumbsch, Peter; Liu, Xin Z; Ding, Xiangdong; Sun, Jun; Li, Ju

2016-11-24

Graphite and other lamellar materials are used as dry lubricants for macroscale metallic sliding components and high-pressure contacts. It has been shown experimentally that monolayer graphene exhibits higher friction than multilayer graphene and graphite, and that this friction increases with continued sliding, but the mechanism behind this remains subject to debate. It has long been conjectured that the true contact area between two rough bodies controls interfacial friction. The true contact area, defined for example by the number of atoms within the range of interatomic forces, is difficult to visualize directly but characterizes the quantity of contact. However, there is emerging evidence that, for a given pair of materials, the quality of the contact can change, and that this can also strongly affect interfacial friction. Recently, it has been found that the frictional behaviour of two-dimensional materials exhibits traits unlike those of conventional bulk materials. This includes the abovementioned finding that for few-layer two-dimensional materials the static friction force gradually strengthens for a few initial atomic periods before reaching a constant value. Such transient behaviour, and the associated enhancement of steady-state friction, diminishes as the number of two-dimensional layers increases, and was observed only when the two-dimensional material was loosely adhering to a substrate. This layer-dependent transient phenomenon has not been captured by any simulations. Here, using atomistic simulations, we reproduce the experimental observations of layer-dependent friction and transient frictional strengthening on graphene. Atomic force analysis reveals that the evolution of static friction is a manifestation of the natural tendency for thinner and less-constrained graphene to re-adjust its configuration as a direct consequence of its greater flexibility. That is, the tip atoms become more strongly pinned, and show greater synchrony in their stick-slip behaviour. While the quantity of atomic-scale contacts (true contact area) evolves, the quality (in this case, the local pinning state of individual atoms and the overall commensurability) also evolves in frictional sliding on graphene. Moreover, the effects can be tuned by pre-wrinkling. The evolving contact quality is critical for explaining the time-dependent friction of configurationally flexible interfaces.

Enhanced Thermal Diffusion of Li in Graphite by Alternating Vertical Electric Field: A Hybrid Quantum-Classical Simulation Study

NASA Astrophysics Data System (ADS)

Ohba, Nobuko; Ogata, Shuji; Tamura, Tomoyuki; Kobayashi, Ryo; Yamakawa, Shunsuke; Asahi, Ryoji

2012-02-01

Enhancing the diffusivity of the Li ion in a Li-graphite intercalation compound that has been used as a negative electrode in the Li-ion rechargeable battery, is important in improving both the recharging speed and power of the battery. In the compound, the Li ion creates a long-range stress field around itself by expanding the interlayer spacing of graphite. We advance the hybrid quantum-classical simulation code to include the external electric field in addition to the long-range stress field by first-principles simulation. In the hybrid code, the quantum region selected adaptively around the Li ion is treated using the real-space density-functional theory for electrons. The rest of the system is described with an empirical interatomic potential that includes the term relating to the dispersion force between the C atoms in different layers. Hybrid simulation runs for Li dynamics in graphite are performed at 423 K under various settings of the amplitude and frequency of alternating electric fields perpendicular to C-layers. We find that the in-plane diffusivity of the Li ion is enhanced significantly by the electric field if the amplitude is larger than 0.2 V/Å within its order and the frequency is as high as 1.7 THz. The microscopic mechanisms of the enhancement are explained.

Hierarchical Multiscale Modeling of Macromolecules and their Assemblies

PubMed Central

Ortoleva, P.; Singharoy, A.; Pankavich, S.

2013-01-01

Soft materials (e.g., enveloped viruses, liposomes, membranes and supercooled liquids) simultaneously deform or display collective behaviors, while undergoing atomic scale vibrations and collisions. While the multiple space-time character of such systems often makes traditional molecular dynamics simulation impractical, a multiscale approach has been presented that allows for long-time simulation with atomic detail based on the co-evolution of slowly-varying order parameters (OPs) with the quasi-equilibrium probability density of atomic configurations. However, this approach breaks down when the structural change is extreme, or when nearest-neighbor connectivity of atoms is not maintained. In the current study, a self-consistent approach is presented wherein OPs and a reference structure co-evolve slowly to yield long-time simulation for dynamical soft-matter phenomena such as structural transitions and self-assembly. The development begins with the Liouville equation for N classical atoms and an ansatz on the form of the associated N-atom probability density. Multiscale techniques are used to derive Langevin equations for the coupled OP-configurational dynamics. The net result is a set of equations for the coupled stochastic dynamics of the OPs and centers of mass of the subsystems that constitute a soft material body. The theory is based on an all-atom methodology and an interatomic force field, and therefore enables calibration-free simulations of soft matter, such as macromolecular assemblies. PMID:23671457

Allosteric Effect of Adenosine Triphosphate on Peptide Recognition by 3'5'-Cyclic Adenosine Monophosphate Dependent Protein Kinase Catalytic Subunits.

PubMed

Kivi, Rait; Solovjova, Karina; Haljasorg, Tõiv; Arukuusk, Piret; Järv, Jaak

2016-12-01

The allosteric influence of adenosine triphosphate (ATP) on the binding effectiveness of a series of peptide inhibitors with the catalytic subunit of 3'5'-cyclic adenosine monophosphate dependent protein kinase was investigated, and the dependence of this effect on peptide structure was analyzed. The allosteric effect was calculated as ratio of peptide binding effectiveness with the enzyme-ATP complex and with the free enzyme, quantified by the competitive inhibition of the enzyme in the presence of ATP excess, and by the enzyme-peptide complex denaturation assay, respectively It was found that the principle "better binding-stronger allostery" holds for interactions of the studied peptides with the enzyme, indicating that allostery and peptide binding with the free enzyme are governed by the same specificity pattern. This means that the allosteric regulation does not include new ligand-protein interactions, but changes the intensity (strength) of the interatomic forces that govern the complex formation in the case of each individual ligand. We propose that the allosteric regulation can be explained by the alteration of the intrinsic dynamics of the protein by ligand binding, and that this phenomenon, in turn, modulates the ligand off-rate from its binding site as well as the binding affinity. The positive allostery could therefore be induced by a reduction in the enzyme's overall intrinsic dynamics.

Atomistic Design of Favored Compositions for Synthesizing the Al-Ni-Y Metallic Glasses

PubMed Central

Wang, Q.; Li, J. H.; Liu, J. B.; Liu, B. X.

2015-01-01

For a ternary alloy system promising for obtaining the so-called bulk metallic glasses (BMGs), the first priority issue is to predict the favored compositions, which could then serve as guidance for the appropriate alloy design. Taking the Al-Ni-Y system as an example, here we show an atomistic approach, which is developed based on a recently constructed and proven realistic interatomic potential of the system. Applying the Al-Ni-Y potential, series simulations not only clarify the glass formation mechanism, but also predict in the composition triangle, a hexagonal region, in which a disordered state, i.e., the glassy phase, is favored energetically. The predicted region is defined as glass formation region (GFR) for the ternary alloy system. Moreover, the approach is able to calculate an amorphization driving force (ADF) for each possible glassy alloy located within the GFR. The calculations predict an optimized sub-region nearby a stoichiometry of Al80Ni5Y15, implying that the Al-Ni-Y metallic glasses designed in the sub-region could be the most stable. Interestingly, the atomistic predictions are supported by experimental results observed in the Al-Ni-Y system. In addition, structural origin underlying the stability of the Al-Ni-Y metallic glasses is also discussed in terms of a hybrid packing mode in the medium-range scale. PMID:26592568

Ab initio phonon point defect scattering and thermal transport in graphene

NASA Astrophysics Data System (ADS)

Polanco, Carlos A.; Lindsay, Lucas

2018-01-01

We study the scattering of phonons from point defects and their effect on lattice thermal conductivity κ using a parameter-free ab initio Green's function methodology. Specifically, we focus on the scattering of phonons by boron (B), nitrogen (N), and phosphorus substitutions as well as single- and double-carbon vacancies in graphene. We show that changes of the atomic structure and harmonic interatomic force constants locally near defects govern the strength and frequency trends of the scattering of out-of-plane acoustic (ZA) phonons, the dominant heat carriers in graphene. ZA scattering rates due to N substitutions are nearly an order of magnitude smaller than those for B defects despite having similar mass perturbations. Furthermore, ZA phonon scattering rates from N defects decrease with increasing frequency in the lower-frequency spectrum in stark contrast to expected trends from simple models. ZA phonon-vacancy scattering rates are found to have a significantly softer frequency dependence (˜ω0 ) in graphene than typically employed in phenomenological models. The rigorous Green's function calculations demonstrate that typical mass-defect models do not adequately describe ZA phonon-defect scattering rates. Our ab initio calculations capture well the trend of κ vs vacancy density from experiments, though not the magnitudes. This work elucidates important insights into phonon-defect scattering and thermal transport in graphene, and demonstrates the applicability of first-principles methods toward describing these properties in imperfect materials.

Accelerating evaluation of converged lattice thermal conductivity

NASA Astrophysics Data System (ADS)

Qin, Guangzhao; Hu, Ming

2018-01-01

High-throughput computational materials design is an emerging area in materials science, which is based on the fast evaluation of physical-related properties. The lattice thermal conductivity (κ) is a key property of materials for enormous implications. However, the high-throughput evaluation of κ remains a challenge due to the large resources costs and time-consuming procedures. In this paper, we propose a concise strategy to efficiently accelerate the evaluation process of obtaining accurate and converged κ. The strategy is in the framework of phonon Boltzmann transport equation (BTE) coupled with first-principles calculations. Based on the analysis of harmonic interatomic force constants (IFCs), the large enough cutoff radius (rcutoff), a critical parameter involved in calculating the anharmonic IFCs, can be directly determined to get satisfactory results. Moreover, we find a simple way to largely ( 10 times) accelerate the computations by fast reconstructing the anharmonic IFCs in the convergence test of κ with respect to the rcutof, which finally confirms the chosen rcutoff is appropriate. Two-dimensional graphene and phosphorene along with bulk SnSe are presented to validate our approach, and the long-debate divergence problem of thermal conductivity in low-dimensional systems is studied. The quantitative strategy proposed herein can be a good candidate for fast evaluating the reliable κ and thus provides useful tool for high-throughput materials screening and design with targeted thermal transport properties.

Atomistic Design of Favored Compositions for Synthesizing the Al-Ni-Y Metallic Glasses

NASA Astrophysics Data System (ADS)

Wang, Q.; Li, J. H.; Liu, J. B.; Liu, B. X.

2015-11-01

For a ternary alloy system promising for obtaining the so-called bulk metallic glasses (BMGs), the first priority issue is to predict the favored compositions, which could then serve as guidance for the appropriate alloy design. Taking the Al-Ni-Y system as an example, here we show an atomistic approach, which is developed based on a recently constructed and proven realistic interatomic potential of the system. Applying the Al-Ni-Y potential, series simulations not only clarify the glass formation mechanism, but also predict in the composition triangle, a hexagonal region, in which a disordered state, i.e., the glassy phase, is favored energetically. The predicted region is defined as glass formation region (GFR) for the ternary alloy system. Moreover, the approach is able to calculate an amorphization driving force (ADF) for each possible glassy alloy located within the GFR. The calculations predict an optimized sub-region nearby a stoichiometry of Al80Ni5Y15, implying that the Al-Ni-Y metallic glasses designed in the sub-region could be the most stable. Interestingly, the atomistic predictions are supported by experimental results observed in the Al-Ni-Y system. In addition, structural origin underlying the stability of the Al-Ni-Y metallic glasses is also discussed in terms of a hybrid packing mode in the medium-range scale.

Empirical force field-based kinetic Monte Carlo simulation of precipitate evolution and growth in Al-Cu alloys

NASA Astrophysics Data System (ADS)

Joshi, Kaushik; Chaudhuri, Santanu

2016-10-01

Ability to accelerate the morphological evolution of nanoscale precipitates is a fundamental challenge for atomistic simulations. Kinetic Monte Carlo (KMC) methodology is an effective approach for accelerating the evolution of nanoscale systems that are dominated by so-called rare events. The quality and accuracy of energy landscape used in KMC calculations can be significantly improved using DFT-informed interatomic potentials. Using newly developed computational framework that uses molecular simulator LAMMPS as a library function inside KMC solver SPPARKS, we investigated formation and growth of Guiner-Preston (GP) zones in dilute Al-Cu alloys at different temperature and copper concentrations. The KMC simulations with angular dependent potential (ADP) predict formation of coherent disc-shaped monolayers of copper atoms (GPI zones) in early stage. Such monolayers are then gradually transformed into energetically favored GPII phase that has two aluminum layers sandwiched between copper layers. We analyzed the growth kinetics of KMC trajectory using Johnson-Mehl-Avrami (JMA) theory and obtained a phase transformation index close to 1.0. In the presence of grain boundaries, the KMC calculations predict the segregation of copper atoms near the grain boundaries instead of formation of GP zones. The computational framework presented in this work is based on open source potentials and MD simulator and can predict morphological changes during the evolution of the alloys in the bulk and around grain boundaries.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rist, J.; Miteva, T.; Gaire, B.

In this paper we present a comprehensive and detailed study of Interatomic Coulombic Decay (ICD) occurring after irradiating argon dimers with XUV-synchrotron radiation. A manifold of different decay channels is observed and the corresponding initial and final states are assigned. Additionally, the effect of nuclear dynamics on the ICD electron spectrum is examined for one specific decay channel. The internuclear distance-dependent width Γ(R) of the decay is obtained from the measured kinetic energy release distribution of the ions employing a classical nuclear dynamics model.

Lattice Thermal Conductivity of Ultra High Temperature Ceramics ZrB2 and HfB2 from Atomistic Simulations

NASA Technical Reports Server (NTRS)

Lawson, John W.; Murray, Daw S.; Bauschlicher, Charles W., Jr.

2011-01-01

Atomistic Green-Kubo simulations are performed to evaluate the lattice thermal conductivity for single crystals of the ultra high temperature ceramics ZrB2 and HfB2 for a range of temperatures. Recently developed interatomic potentials are used for these simulations. Heat current correlation functions show rapid oscillations which can be identified with mixed metal-Boron optical phonon modes. Agreement with available experimental data is good.

Density and glass forming ability in amorphous atomic alloys: The role of the particle softness

DOE Office of Scientific and Technical Information (OSTI.GOV)

Douglass, Ian; Hudson, Toby; Harrowell, Peter

A key property of glass forming alloys, the anomalously small volume difference with respect to the crystal, is shown to arise as a direct consequence of the soft repulsive potentials between metals. This feature of the inter-atomic potential is demonstrated to be responsible for a significant component of the glass forming ability of alloys due to the decrease in the enthalpy of fusion and the associated depression of the freezing point.

When hydroquinone meets methoxy radical: Hydrogen abstraction reaction from the viewpoint of interacting quantum atoms.

PubMed

Petković, Milena; Nakarada, Đura; Etinski, Mihajlo

2018-05-25

Interacting Quantum Atoms methodology is used for a detailed analysis of hydrogen abstraction reaction from hydroquinone by methoxy radical. Two pathways are analyzed, which differ in the orientation of the reactants at the corresponding transition states. Although the discrepancy between the two barriers amounts to only 2 kJ/mol, which implies that the two pathways are of comparable probability, the extent of intra-atomic and inter-atomic energy changes differs considerably. We thus demonstrated that Interacting Quantum Atoms procedure can be applied to unravel distinct energy transfer routes in seemingly similar mechanisms. Identification of energy components with the greatest contribution to the variation of the overall energy (intra-atomic and inter-atomic terms that involve hydroquinone's oxygen and the carbon atom covalently bound to it, the transferring hydrogen and methoxy radical's oxygen), is performed using the Relative energy gradient method. Additionally, the Interacting Quantum Fragments approach shed light on the nature of dominant interactions among selected fragments: both Coulomb and exchange-correlation contributions are of comparable importance when considering interactions of the transferring hydrogen atom with all other atoms, whereas the exchange-correlation term dominates interaction between methoxy radical's methyl group and hydroquinone's aromatic ring. This study represents one of the first applications of Interacting Quantum Fragments approach on first order saddle points. © 2018 Wiley Periodicals, Inc. © 2018 Wiley Periodicals, Inc.

Intrinsic Resolution of Molecular Electronic Wave Functions and Energies in Terms of Quasi-atoms and Their Interactions

DOE Office of Scientific and Technical Information (OSTI.GOV)

West, Aaron C.; Schmidt, Michael W.; Gordon, Mark S.

A general intrinsic energy resolution has been formulated for strongly correlated wave functions in the full molecular valence space and its subspaces. The information regarding the quasi-atomic organization of the molecular electronic structure is extracted from the molecular wave function without introducing any additional postulated model state wave functions. To this end, the molecular wave function is expressed in terms of quasi-atomic molecular orbitals, which maximize the overlap between subspaces of the molecular orbital space and the free-atom orbital spaces. As a result, the molecular wave function becomes the superposition of a wave function representing the non-bonded juxtaposed quasi-atoms andmore » a wave function describing the interatomic electron migrations that create bonds through electron sharing. The juxtaposed nonbonded quasi-atoms are shown to consist of entangled quasi-atomic states from different atoms. The binding energy is resolved as a sum of contributions that are due to quasi-atom formation, quasiclassical electrostatic interactions and interatomic interferences caused by electron sharing. The contributions are further resolved according to orbital interactions. The various transformations that generate the analysis are determined by criteria that are independent of the working orbital basis used for calculating the molecular wave function. Lastly, the theoretical formulation of the resolution is quantitatively validated by an application to the C 2 molecule.« less

Intrinsic Resolution of Molecular Electronic Wave Functions and Energies in Terms of Quasi-atoms and Their Interactions

DOE PAGES

West, Aaron C.; Schmidt, Michael W.; Gordon, Mark S.; ...

2017-01-30

A general intrinsic energy resolution has been formulated for strongly correlated wave functions in the full molecular valence space and its subspaces. The information regarding the quasi-atomic organization of the molecular electronic structure is extracted from the molecular wave function without introducing any additional postulated model state wave functions. To this end, the molecular wave function is expressed in terms of quasi-atomic molecular orbitals, which maximize the overlap between subspaces of the molecular orbital space and the free-atom orbital spaces. As a result, the molecular wave function becomes the superposition of a wave function representing the non-bonded juxtaposed quasi-atoms andmore » a wave function describing the interatomic electron migrations that create bonds through electron sharing. The juxtaposed nonbonded quasi-atoms are shown to consist of entangled quasi-atomic states from different atoms. The binding energy is resolved as a sum of contributions that are due to quasi-atom formation, quasiclassical electrostatic interactions and interatomic interferences caused by electron sharing. The contributions are further resolved according to orbital interactions. The various transformations that generate the analysis are determined by criteria that are independent of the working orbital basis used for calculating the molecular wave function. Lastly, the theoretical formulation of the resolution is quantitatively validated by an application to the C 2 molecule.« less

Anomalous properties of the acoustic excitations in glasses on the mesoscopic length scale.

PubMed

Monaco, Giulio; Mossa, Stefano

2009-10-06

The low-temperature thermal properties of dielectric crystals are governed by acoustic excitations with large wavelengths that are well described by plane waves. This is the Debye model, which rests on the assumption that the medium is an elastic continuum, holds true for acoustic wavelengths large on the microscopic scale fixed by the interatomic spacing, and gradually breaks down on approaching it. Glasses are characterized as well by universal low-temperature thermal properties that are, however, anomalous with respect to those of the corresponding crystalline phases. Related universal anomalies also appear in the low-frequency vibrational density of states and, despite a longstanding debate, remain poorly understood. By using molecular dynamics simulations of a model monatomic glass of extremely large size, we show that in glasses the structural disorder undermines the Debye model in a subtle way: The elastic continuum approximation for the acoustic excitations breaks down abruptly on the mesoscopic, medium-range-order length scale of approximately 10 interatomic spacings, where it still works well for the corresponding crystalline systems. On this scale, the sound velocity shows a marked reduction with respect to the macroscopic value. This reduction turns out to be closely related to the universal excess over the Debye model prediction found in glasses at frequencies of approximately 1 THz in the vibrational density of states or at temperatures of approximately 10 K in the specific heat.

Resonant interatomic Coulombic decay in HeNe: Electron angular emission distributions

NASA Astrophysics Data System (ADS)

Mhamdi, A.; Trinter, F.; Rauch, C.; Weller, M.; Rist, J.; Waitz, M.; Siebert, J.; Metz, D.; Janke, C.; Kastirke, G.; Wiegandt, F.; Bauer, T.; Tia, M.; Cunha de Miranda, B.; Pitzer, M.; Sann, H.; Schiwietz, G.; Schöffler, M.; Simon, M.; Gokhberg, K.; Dörner, R.; Jahnke, T.; Demekhin, Ph. Â. V.

2018-05-01

We present a joint experimental and theoretical study of resonant interatomic Coulombic decay (RICD) in HeNe employing high resolution cold target recoil ion momentum spectroscopy and ab initio electronic structure and nuclear dynamics calculations. In particular, laboratory- and molecular-frame angular emission distributions of RICD electrons are examined in detail. The exciting-photon energy-dependent anisotropy parameter β (ω ) , measured for decay events that populate bound HeNe+ ions, is in agreement with the calculations performed for the ground ionic state X2Σ1/2 + . A contribution from the a2Π3 /2 final ionic state is found to be negligible. For the He +Ne+ fragmentation channel, the observed laboratory-frame angular distribution of RICD electrons is explained by a slow homogeneous dissociation of bound vibrational levels of the final ionic state A2Π1 /2 into vibrational continua of the lower lying states X2Σ1/2 + and a2Π3 /2 . Our calculations predict that the angular distributions of RICD electrons in the body-fixed dipole plane provide direct access to the electronic character (i.e., symmetry) of intermediate vibronic resonances. However, because of the very slow dissociation of the A2Π1 /2 state, the molecular-frame angular distributions of RICD electrons in the He +Ne+ fragmentation channel are inaccessible to our coincidence experiment.

Inverse scattering transform and soliton solutions for square matrix nonlinear Schrödinger equations with non-zero boundary conditions

NASA Astrophysics Data System (ADS)

Prinari, Barbara; Demontis, Francesco; Li, Sitai; Horikis, Theodoros P.

2018-04-01

The inverse scattering transform (IST) with non-zero boundary conditions at infinity is developed for an m × m matrix nonlinear Schrödinger-type equation which, in the case m = 2, has been proposed as a model to describe hyperfine spin F = 1 spinor Bose-Einstein condensates with either repulsive interatomic interactions and anti-ferromagnetic spin-exchange interactions (self-defocusing case), or attractive interatomic interactions and ferromagnetic spin-exchange interactions (self-focusing case). The IST for this system was first presented by Ieda et al. (2007) , using a different approach. In our formulation, both the direct and the inverse problems are posed in terms of a suitable uniformization variable which allows to develop the IST on the standard complex plane, instead of a two-sheeted Riemann surface or the cut plane with discontinuities along the cuts. Analyticity of the scattering eigenfunctions and scattering data, symmetries, properties of the discrete spectrum, and asymptotics are derived. The inverse problem is posed as a Riemann-Hilbert problem for the eigenfunctions, and the reconstruction formula of the potential in terms of eigenfunctions and scattering data is provided. In addition, the general behavior of the soliton solutions is analyzed in detail in the 2 × 2 self-focusing case, including some special solutions not previously discussed in the literature.

NASA Astrophysics Data System (ADS)

2017-11-01

To deal with these problems investigators usually rely on a calibration method that makes use of a substance with an accurately known set of interatomic distances. The procedure consists of carrying out a diffraction experiment on the chosen calibrating substance, determining the value of the distances with use of the nominal (meter) value of the voltage, and then correcting the nominal voltage by an amount that produces the distances in the calibration substance. Examples of gases that have been used for calibration are carbon dioxide, carbon tetrachloride, carbon disulfide, and benzene; solids such as zinc oxide smoke (powder) deposited on a screen or slit have also been used. The question implied by the use of any standard molecule is, how accurate are the interatomic distance values assigned to the standard? For example, a solid calibrant is subject to heating by the electron beam, possibly producing unknown changes in the lattice constants, and polyatomic gaseous molecules require corrections for vibrational averaging ("shrinkage") effects that are uncertain at best. It has lately been necessary for us to investigate this matter in connection with on-going studies of several molecules in which size is the most important issue. These studies indicated that our usual method for retrieval of data captured on film needed improvement. The following is an account of these two issues - the accuracy of the distances assigned to the chosen standard molecule, and the improvements in our methods of retrieving the scattered intensity data.

First-principles study of point defects at a semicoherent interface

DOE PAGES

Metsanurk, E.; Tamm, A.; Caro, A.; ...

2014-12-19

Most of the atomistic modeling of semicoherent metal-metal interfaces has so far been based on the use of semiempirical interatomic potentials. Here, we show that key conclusions drawn from previous studies are in contradiction with more precise ab-initio calculations. In particular we find that single point defects do not delocalize, but remain compact near the interfacial plane in Cu-Nb multilayers. Lastly, we give a simple qualitative explanation for this difference on the basis of the well known limited transferability of empirical potentials.

Interatomic scattering in energy dependent photoelectron spectra of Ar clusters

DOE Office of Scientific and Technical Information (OSTI.GOV)

Patanen, M.; Benkoula, S.; Nicolas, C.

2015-09-28

Soft X-ray photoelectron spectra of Ar 2p levels of atomic argon and argon clusters are recorded over an extended range of photon energies. The Ar 2p intensity ratios between atomic argon and clusters’ surface and bulk components reveal oscillations similar to photoelectron extended X-ray absorption fine structure signal (PEXAFS). We demonstrate here that this technique allows us to analyze separately the PEXAFS signals from surface and bulk sites of free-standing, neutral clusters, revealing a bond contraction at the surface.

Foraging on the potential energy surface: a swarm intelligence-based optimizer for molecular geometry.

PubMed

Wehmeyer, Christoph; Falk von Rudorff, Guido; Wolf, Sebastian; Kabbe, Gabriel; Schärf, Daniel; Kühne, Thomas D; Sebastiani, Daniel

2012-11-21

We present a stochastic, swarm intelligence-based optimization algorithm for the prediction of global minima on potential energy surfaces of molecular cluster structures. Our optimization approach is a modification of the artificial bee colony (ABC) algorithm which is inspired by the foraging behavior of honey bees. We apply our modified ABC algorithm to the problem of global geometry optimization of molecular cluster structures and show its performance for clusters with 2-57 particles and different interatomic interaction potentials.

Molecular structure, interatomic interactions and vibrational analysis of 1,4-diazabicyclo[3.2.1]octane parent ring system

NASA Astrophysics Data System (ADS)

Britvin, Sergey N.; Rumyantsev, Andrey M.; Zobnina, Anastasia E.; Padkina, Marina V.

2017-02-01

Molecular structure of 1,4-diazabicyclo[3.2.1]octane, a parent ring of TAN1251 family of alkaloids, is herein characterized for the first time in comparison with the structure of nortropane (8-azabicyclo[3.2.1]octane), the parent framework of tropane ring system. The methods of study involve X-ray structural analysis, DFT geometry optimizations with infrared frequency calculations followed by natural bond orbital (NBO) analysis, and vibrational analysis of infrared spectrum.

Foraging on the potential energy surface: A swarm intelligence-based optimizer for molecular geometry

NASA Astrophysics Data System (ADS)

Wehmeyer, Christoph; Falk von Rudorff, Guido; Wolf, Sebastian; Kabbe, Gabriel; Schärf, Daniel; Kühne, Thomas D.; Sebastiani, Daniel

2012-11-01

We present a stochastic, swarm intelligence-based optimization algorithm for the prediction of global minima on potential energy surfaces of molecular cluster structures. Our optimization approach is a modification of the artificial bee colony (ABC) algorithm which is inspired by the foraging behavior of honey bees. We apply our modified ABC algorithm to the problem of global geometry optimization of molecular cluster structures and show its performance for clusters with 2-57 particles and different interatomic interaction potentials.

Tunneling of Two Interacting Fermions

NASA Astrophysics Data System (ADS)

Ishmukhamedov, Ilyas; Ishmukhamedov, Altay

2018-04-01

We consider two interacting atoms subject to a one-dimensional anharmonic trap and magnetic field gradient. This system has been recently investigated by the Heidelberg group in the experiment on two 6Li atoms. In the present paper the tunneling of two cold 6Li atoms, initially prepared in the center-of-mass and relative motion excited state, is explored and full time-dependent simulation of the tunneling dynamics is performed. The dynamics is analyzed for the interatomic coupling strength ranging from strong attraction to strong repulsion.

Automated Algorithms for Quantum-Level Accuracy in Atomistic Simulations: LDRD Final Report.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Thompson, Aidan Patrick; Schultz, Peter Andrew; Crozier, Paul

2014-09-01

This report summarizes the result of LDRD project 12-0395, titled "Automated Algorithms for Quantum-level Accuracy in Atomistic Simulations." During the course of this LDRD, we have developed an interatomic potential for solids and liquids called Spectral Neighbor Analysis Poten- tial (SNAP). The SNAP potential has a very general form and uses machine-learning techniques to reproduce the energies, forces, and stress tensors of a large set of small configurations of atoms, which are obtained using high-accuracy quantum electronic structure (QM) calculations. The local environment of each atom is characterized by a set of bispectrum components of the local neighbor density projectedmore » on to a basis of hyperspherical harmonics in four dimensions. The SNAP coef- ficients are determined using weighted least-squares linear regression against the full QM training set. This allows the SNAP potential to be fit in a robust, automated manner to large QM data sets using many bispectrum components. The calculation of the bispectrum components and the SNAP potential are implemented in the LAMMPS parallel molecular dynamics code. Global optimization methods in the DAKOTA software package are used to seek out good choices of hyperparameters that define the overall structure of the SNAP potential. FitSnap.py, a Python-based software pack- age interfacing to both LAMMPS and DAKOTA is used to formulate the linear regression problem, solve it, and analyze the accuracy of the resultant SNAP potential. We describe a SNAP potential for tantalum that accurately reproduces a variety of solid and liquid properties. Most significantly, in contrast to existing tantalum potentials, SNAP correctly predicts the Peierls barrier for screw dislocation motion. We also present results from SNAP potentials generated for indium phosphide (InP) and silica (SiO 2 ). We describe efficient algorithms for calculating SNAP forces and energies in molecular dynamics simulations using massively parallel computers and advanced processor ar- chitectures. Finally, we briefly describe the MSM method for efficient calculation of electrostatic interactions on massively parallel computers.« less

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

NASA Technical Reports Server (NTRS)

Han, Jie; Globus, Al; Srivastava, Deepak; Chancellor, Marisa K. (Technical Monitor)

1997-01-01

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Molecular dynamics modeling of atomic displacement cascades in 3C-SiC: Comparison of interatomic potentials

NASA Astrophysics Data System (ADS)

Samolyuk, G. D.; Osetsky, Y. N.; Stoller, R. E.

2015-10-01

We used molecular dynamics modeling of atomic displacement cascades to characterize the nature of primary radiation damage in 3C-SiC. We demonstrated that the most commonly used interatomic potentials are inconsistent with ab initio calculations of defect energetics. Both the Tersoff potential used in this work and a modified embedded-atom method potential reveal a barrier to recombination of the carbon interstitial and carbon vacancy which is much higher than the density functional theory (DFT) results. The barrier obtained with a newer potential by Gao and Weber is closer to the DFT result. This difference results in significant differences in the cascade production of point defects. We have completed both 10 keV and 50 keV cascade simulations in 3C-SiC at a range of temperatures. In contrast to the Tersoff potential, the Gao-Weber potential produces almost twice as many C vacancies and interstitials at the time of maximum disorder (∼0.2 ps) but only about 25% more stable defects at the end of the simulation. Only about 20% of the carbon defects produced with the Tersoff potential recombine during the in-cascade annealing phase, while about 60% recombine with the Gao-Weber potential. The Gao-Weber potential appears to give a more realistic description of cascade dynamics in SiC, but still has some shortcomings when the defect migration barriers are compared to the ab initio results.

Probing the interface of doped isotopically mixed helium droplets by the directional anisotropy of interatomic Coulombic decay.

PubMed

Kryzhevoi, Nikolai V; Mateo, David; Pi, Martí; Barranco, Manuel; Cederbaum, Lorenz S

2013-11-07

Interatomic Coulombic decay (ICD) represents an efficient electronic relaxation mechanism of an ionized or an excited system embedded in an environment. The type of this environment and its size have a great impact on the ICD performance. It is stressed that ICD is sensitive to the arrangement of neighboring atoms when the initially created vacancy has a polarization direction. This is demonstrated in the present paper for the case of a 3p-ionized Ca surrounded by He atoms. Useful explicit expressions are derived for the ICD widths which show that the neighbors located along the polarization direction of the ionized orbital have the largest contribution to the ICD rate. By comparison with ab initio results for small clusters, we also show that in a helium environment, the pairwise approximation represents a reliable approach for computing ICD widths. Using this approximation and the density distribution of the helium atoms obtained within density functional theory, we explore ICD in large isotopically mixed helium droplets doped with Ca. A special emphasis is given to the difference between the ICD widths for the Ca3p orbitals directed perpendicular and parallel to the droplet surface. Depending on the size and isotopic composition of the droplet, Ca resides in the interfacial layer between the (4)He core and the (3)He outer shell. Hence, ICD studies in these droplets may provide valuable information on the properties of this interface.

Application of Powder Diffraction Methods to the Analysis of Short- and Long-Range Atomic Order in Nanocrystalline Diamond and SiC: The Concept of the Apparent Lattice Parameter (alp)

NASA Technical Reports Server (NTRS)

Palosz, B.; Grzanka, E.; Gierlotka, S.; Stelmakh, S.; Pielaszek, R.; Bismayer, U.; Weber, H.-P.; Palosz, W.

2003-01-01

Two methods of the analysis of powder diffraction patterns of diamond and SiC nanocrystals are presented: (a) examination of changes of the lattice parameters with diffraction vector Q ('apparent lattice parameter', alp) which refers to Bragg scattering, and (b), examination of changes of inter-atomic distances based on the analysis of the atomic Pair Distribution Function, PDF. Application of these methods was studied based on the theoretical diffraction patterns computed for models of nanocrystals having (i) a perfect crystal lattice, and (ii), a core-shell structure, i.e. constituting a two-phase system. The models are defined by the lattice parameter of the grain core, thickness of the surface shell, and the magnitude and distribution of the strain field in the shell. X-ray and neutron experimental diffraction data of nanocrystalline SiC and diamond powders of the grain diameter from 4 nm up to micrometers were used. The effects of the internal pressure and strain at the grain surface on the structure are discussed based on the experimentally determined dependence of the alp values on the Q-vector, and changes of the interatomic distances with the grain size determined experimentally by the atomic Pair Distribution Function (PDF) analysis. The experimental results lend a strong support to the concept of a two-phase, core and the surface shell structure of nanocrystalline diamond and SiC.

On the role of structure-dynamic relationship in determining the excess entropy of mixing and chemical ordering in binary square-well liquid alloys

NASA Astrophysics Data System (ADS)

Lalneihpuii, R.; Shrivastava, Ruchi; Mishra, Raj Kumar

2018-05-01

Using statistical mechanical model with square-well (SW) interatomic potential within the frame work of mean spherical approximation, we determine the composition dependent microscopic correlation functions, interdiffusion coefficients, surface tension and chemical ordering in Ag-Cu melts. Further Dzugutov universal scaling law of normalized diffusion is verified with SW potential in binary mixtures. We find that the excess entropy scaling law is valid for SW binary melts. The partial and total structure factors in the attractive and repulsive regions of the interacting potential are evaluated and then Fourier transformed to get partial and total radial distribution functions. A good agreement between theoretical and experimental values for total structure factor and the reduced radial distribution function are observed, which consolidates our model calculations. The well-known Bhatia-Thornton correlation functions are also computed for Ag-Cu melts. The concentration-concentration correlations in the long wavelength limit in liquid Ag-Cu alloys have been analytically derived through the long wavelength limit of partial correlation functions and apply it to demonstrate the chemical ordering and interdiffusion coefficients in binary liquid alloys. We also investigate the concentration dependent viscosity coefficients and surface tension using the computed diffusion data in these alloys. Our computed results for structure, transport and surface properties of liquid Ag-Cu alloys obtained with square-well interatomic interaction are fully consistent with their corresponding experimental values.

Heat conduction in diatomic chains with correlated disorder

NASA Astrophysics Data System (ADS)

Savin, Alexander V.; Zolotarevskiy, Vadim; Gendelman, Oleg V.

2017-01-01

The paper considers heat transport in diatomic one-dimensional lattices, containing equal amounts of particles with different masses. Ordering of the particles in the chain is governed by single correlation parameter - the probability for two neighboring particles to have the same mass. As this parameter grows from zero to unity, the structure of the chain varies from regular staggering chain to completely random configuration, and then - to very long clusters of particles with equal masses. Therefore, this correlation parameter allows a control of typical cluster size in the chain. In order to explore different regimes of the heat transport, two interatomic potentials are considered. The first one is an infinite potential wall, corresponding to instantaneous elastic collisions between the neighboring particles. In homogeneous chains such interaction leads to an anomalous heat transport. The other one is classical Lennard-Jones interatomic potential, which leads to a normal heat transport. The simulations demonstrate that the correlated disorder of the particle arrangement does not change the convergence properties of the heat conduction coefficient, but essentially modifies its value. For the collision potential, one observes essential growth of the coefficient for fixed chain length as the limit of large homogeneous clusters is approached. The thermal transport in these models remains superdiffusive. In the Lennard-Jones chain the effect of correlation appears to be not monotonous in the limit of low temperatures. This behavior stems from the competition between formation of long clusters mentioned above, and Anderson localization close to the staggering ordered state.

Ab initio phonon point defect scattering and thermal transport in graphene

DOE Office of Scientific and Technical Information (OSTI.GOV)

Polanco, Carlos A.; Lindsay, Lucas R.

Here, we study the scattering of phonons from point defects and their effect on lattice thermal conductivity κ using a parameter-free ab initio Green's function methodology. Specifically, we focus on the scattering of phonons by boron (B), nitrogen (N), and phosphorus substitutions as well as single- and double-carbon vacancies in graphene. We show that changes of the atomic structure and harmonic interatomic force constants locally near defects govern the strength and frequency trends of the scattering of out-of-plane acoustic (ZA) phonons, the dominant heat carriers in graphene. ZA scattering rates due to N substitutions are nearly an order of magnitudemore » smaller than those for B defects despite having similar mass perturbations. Furthermore, ZA phonon scattering rates from N defects decrease with increasing frequency in the lower-frequency spectrum in stark contrast to expected trends from simple models. ZA phonon-vacancy scattering rates are found to have a significantly softer frequency dependence (~ω 0) in graphene than typically employed in phenomenological models. The rigorous Green's function calculations demonstrate that typical mass-defect models do not adequately describe ZA phonon-defect scattering rates. Our ab initio calculations capture well the trend of κ vs vacancy density from experiments, though not the magnitudes. In conclusion, this work elucidates important insights into phonon-defect scattering and thermal transport in graphene, and demonstrates the applicability of first-principles methods toward describing these properties in imperfect materials.« less

Car and Parrinello meet Green and Kubo: simulating atomic heat transport from equilibrium ab initio molecular dynamics

NASA Astrophysics Data System (ADS)

Baroni, Stefano

Modern simulation methods based on electronic-structure theory have long been deemed unfit to compute heat transport coefficients within the Green-Kubo formalism. This is so because the quantum-mechanical energy density from which the heat flux is derived is inherently ill defined, thus allegedly hampering the use of the Green-Kubo formula. While this objection would actually apply to classical systems as well, I will demonstrate that the thermal conductivity is indeed independent of the specific microscopic expression for the energy density and current from which it is derived. This fact results from a kind of gauge invariance stemming from energy conservation and extensivity, which I will illustrate numerically for a classical Lennard-Jones fluid. I will then introduce an expression for the adiabatic energy flux, derived within density-functional theory, that allows simulating atomic heat transport using equilibrium ab initio molecular dynamics. The resulting methodology is demonstrated by comparing results from ab-initio and classical molecular-dynamics simulations of a model liquid-Argon system, for which accurate inter-atomic potentials are derived by the force-matching method, and applied to compute the thermal conductivity of heavy water at ambient conditions. The problem of evaluating transport coefficients along with their accuracy from relatively short trajectories is finally addressed and discussed with a few representative examples. Partially funded by the European Union through the MaX Centre of Excellence (Grant No. 676598).

Investigating energy-based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1.

PubMed

Zhu, Guanhua; Liu, Wei; Bao, Chenglong; Tong, Dudu; Ji, Hui; Shen, Zuowei; Yang, Daiwen; Lu, Lanyuan

2018-05-01

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure-based and physics-based atomistic force field with an efficient sampling strategy is adopted to simulate a model di-domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low-energy structures and the minimum-size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small-angle X-ray scattering data. It is illustrated that the regularizations of energy and ensemble-size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high-energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure-ensemble optimizations with a topology-based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates. © 2018 Wiley Periodicals, Inc.

POTENTIOLOGY (noun): study focusing on the development of new interatomic pair potential forms; sometimes pursued in an obsessive compulsive manner [The New Yorel Dictionary (2002, unpublished)].} IN SPECTROSCOPY: IT MATTERS

NASA Astrophysics Data System (ADS)

Le Roy, Robert J.

2009-06-01

Spectroscopists have long attempted to summarize what they know about small molecules in terms of a knowledge of potential energy curves or surfaces. For most of the past century, this involved deducing polynomial-expansion force-field coefficients from energy level expressions fitted to experimental data, or for diatomic molecules, by generating tables of many-digit RKR turning points from such expressions. In recent years, however, it has become increasingly common either to use high-level ab initio calculations to compute the desired potentials, or to determine parametrized global analytic potential functions from direct fits to spectroscopic data. In the former case, this invoked a need for robust, flexible, compact, and `portable' analytic potentials for summarizing the information contained in the (sometimes very large numbers of) ab initio points, and making them `user friendly'. In the latter case, the same properties are required for potentials used in the least-squares fitting procedure. In both cases, there is also a cardinal need for potential function forms that extrapolate sensibly, beyond the range of the experimental data or ab initio points. This talk will describe some recent developments in this area, and make a case for what is arguably the `best' general-purpose analytic potential function form now available. Applications to both diatomic molecules and simple polyatomic molecules will be discussed. footnote

Phonon thermal transport in 2H, 4H and 6H silicon carbide from first principles

DOE PAGES

Protik, Nakib Haider; Katre, Ankita; Lindsay, Lucas R.; ...

2017-06-07

Here, silicon carbide (SiC) is a wide band gap semiconductor with a variety of industrial applications. Among its many useful properties is its high thermal conductivity, which makes it advantageous for thermal management applications. In this paper we present ab initio calculations of the in-plane and cross-plane thermal conductivities, κ in and κ out, of three common hexagonal polytypes of SiC: 2H, 4H and 6H. The phonon Boltzmann transport equation is solved iteratively using as input interatomic force constants determined from density functional theory. Both κ in and κ out decrease with increasing n in nH SiC because of additionalmore » low-lying optic phonon branches. These optic branches are characterized by low phonon group velocities, and they increase the phase space for phonon-phonon scattering of acoustic modes. Also, for all n, κ in is found to be larger than κ out in the temperature range considered. At electron concentrations present in experimental samples, scattering of phonons by electrons is shown to be negligible except well below room temperature where it can lead to a significant reduction of the lattice thermal conductivity. This work highlights the power of ab initio approaches in giving quantitative, predictive descriptions of thermal transport in materials. It helps explain the qualitative disagreement that exists among different sets of measured thermal conductivity data and provides information of the relative quality of samples from which measured data was obtained.« less

Ab initio phonon point defect scattering and thermal transport in graphene

DOE PAGES

Polanco, Carlos A.; Lindsay, Lucas R.

2018-01-04

Here, we study the scattering of phonons from point defects and their effect on lattice thermal conductivity κ using a parameter-free ab initio Green's function methodology. Specifically, we focus on the scattering of phonons by boron (B), nitrogen (N), and phosphorus substitutions as well as single- and double-carbon vacancies in graphene. We show that changes of the atomic structure and harmonic interatomic force constants locally near defects govern the strength and frequency trends of the scattering of out-of-plane acoustic (ZA) phonons, the dominant heat carriers in graphene. ZA scattering rates due to N substitutions are nearly an order of magnitudemore » smaller than those for B defects despite having similar mass perturbations. Furthermore, ZA phonon scattering rates from N defects decrease with increasing frequency in the lower-frequency spectrum in stark contrast to expected trends from simple models. ZA phonon-vacancy scattering rates are found to have a significantly softer frequency dependence (~ω 0) in graphene than typically employed in phenomenological models. The rigorous Green's function calculations demonstrate that typical mass-defect models do not adequately describe ZA phonon-defect scattering rates. Our ab initio calculations capture well the trend of κ vs vacancy density from experiments, though not the magnitudes. In conclusion, this work elucidates important insights into phonon-defect scattering and thermal transport in graphene, and demonstrates the applicability of first-principles methods toward describing these properties in imperfect materials.« less

Ab initio computational study on the lattice thermal conductivity of Zintl clathrates [Si19P4] Cl4 and Na4[Al4Si19

NASA Astrophysics Data System (ADS)

Härkönen, Ville J.; Karttunen, Antti J.

2016-08-01

The lattice thermal conductivity of silicon clathrate framework Si23 and two Zintl clathrates, [Si19P4] Cl4 and Na4[Al4Si19] , is investigated by using an iterative solution of the linearized Boltzmann transport equation in conjunction with ab initio lattice dynamical techniques. At 300 K, the lattice thermal conductivities for Si23, [Si19P4] Cl4 , and Na4[Al4Si19] were found to be 43 W/(m K), 25 W/(m K), and 2 W/(m K), respectively. In the case of Na4[Al4Si19] , the order-of-magnitude reduction in the lattice thermal conductivity was found to be mostly due to relaxation times and group velocities differing from Si23 and [Si19P4] Cl4 . The difference in the relaxation times and group velocities arises primarily due to the phonon spectrum at low frequencies, resulting eventually from the differences in the second-order interatomic force constants (IFCs). The obtained third-order IFCs were rather similar for all materials considered here. The present findings are similar to those obtained earlier for some skutterudites. The predicted lattice thermal conductivity of Na4[Al4Si19] is in line with the experimentally measured thermal conductivity of recently synthesized type-I Zintl clathrate Na8[Al8Si38] (polycrystalline samples).

Divide-and-conquer density functional theory on hierarchical real-space grids: Parallel implementation and applications

NASA Astrophysics Data System (ADS)

Shimojo, Fuyuki; Kalia, Rajiv K.; Nakano, Aiichiro; Vashishta, Priya

2008-02-01

A linear-scaling algorithm based on a divide-and-conquer (DC) scheme has been designed to perform large-scale molecular-dynamics (MD) simulations, in which interatomic forces are computed quantum mechanically in the framework of the density functional theory (DFT). Electronic wave functions are represented on a real-space grid, which is augmented with a coarse multigrid to accelerate the convergence of iterative solutions and with adaptive fine grids around atoms to accurately calculate ionic pseudopotentials. Spatial decomposition is employed to implement the hierarchical-grid DC-DFT algorithm on massively parallel computers. The largest benchmark tests include 11.8×106 -atom ( 1.04×1012 electronic degrees of freedom) calculation on 131 072 IBM BlueGene/L processors. The DC-DFT algorithm has well-defined parameters to control the data locality, with which the solutions converge rapidly. Also, the total energy is well conserved during the MD simulation. We perform first-principles MD simulations based on the DC-DFT algorithm, in which large system sizes bring in excellent agreement with x-ray scattering measurements for the pair-distribution function of liquid Rb and allow the description of low-frequency vibrational modes of graphene. The band gap of a CdSe nanorod calculated by the DC-DFT algorithm agrees well with the available conventional DFT results. With the DC-DFT algorithm, the band gap is calculated for larger system sizes until the result reaches the asymptotic value.

Phonon thermal transport in 2H, 4H and 6H silicon carbide from first principles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Protik, Nakib Haider; Katre, Ankita; Lindsay, Lucas R.

Here, silicon carbide (SiC) is a wide band gap semiconductor with a variety of industrial applications. Among its many useful properties is its high thermal conductivity, which makes it advantageous for thermal management applications. In this paper we present ab initio calculations of the in-plane and cross-plane thermal conductivities, κ in and κ out, of three common hexagonal polytypes of SiC: 2H, 4H and 6H. The phonon Boltzmann transport equation is solved iteratively using as input interatomic force constants determined from density functional theory. Both κ in and κ out decrease with increasing n in nH SiC because of additionalmore » low-lying optic phonon branches. These optic branches are characterized by low phonon group velocities, and they increase the phase space for phonon-phonon scattering of acoustic modes. Also, for all n, κ in is found to be larger than κ out in the temperature range considered. At electron concentrations present in experimental samples, scattering of phonons by electrons is shown to be negligible except well below room temperature where it can lead to a significant reduction of the lattice thermal conductivity. This work highlights the power of ab initio approaches in giving quantitative, predictive descriptions of thermal transport in materials. It helps explain the qualitative disagreement that exists among different sets of measured thermal conductivity data and provides information of the relative quality of samples from which measured data was obtained.« less

Infrared spectra of some sulfides and their analogs of binary composition in the long-wave region

NASA Technical Reports Server (NTRS)

Povarennykh, A. S.; Sidorenko, G. A.; Solntseva, L. S.; Solntsev, B. P.

1981-01-01

The far infrared spectra (500-60/cm) of some simple sulfides and their analogs were studied. In all, 22 minerals with different structure types were investigated, out of which 14 are sulfides (galena, alabandite, pyrrhotite, sphalerite, wurtzite, cinnabar, realgar, orpiment, getchelite antimonite, molybdenite, pyrite, marcasite and heazlewoodite) 6 arsenides (niccolite, domeykite, arsenopyrite, lollingite, rammelsbergite and skutterudite), one telluride (tetradymite) and native arsenic. The main bands of infrared absorption spectra of the minerals are compared with the relative strength of the interatomic bonds and their interpretation is given.

Dislocation core structures of tungsten with dilute solute hydrogen

NASA Astrophysics Data System (ADS)

Wang, Yinan; Li, Qiulin; Li, Chengliang; Shu, Guogang; Xu, Ben; Liu, Wei

2017-12-01

In this paper, a combination of quantum mechanical and interatomic potential-based atomistic calculations are used to predict the core structures of screw and edge dislocations in tungsten in the presence of a particular concentration of hydrogen atoms. These configurations of the core structures are the results of two competing energies: the interaction between the partial dislocations and the corresponding generalized stacking fault energy in between the two partial dislocations, which are presented in this work. With this, we can precisely predict the configurations of the hydrogen-doped dislocation core structures.

Multiscale Modeling of Non-crystalline Ceramics (Glass) (FY11)

DTIC Science & Technology

2012-01-01

interatomic potential model for oxides, silicates, and silica-based glasses. J. Phys. Chem. B 2006, 110 (24), 11780–11795. 9. van Beest , B. W. H.; Kramer, G. J...DIRECTOR US ARMY RESEARCH LAB IMNE ALC HRR 2800 POWDER MILL RD ADELPHI MD 20783-1197 1 DIRECTOR US ARMY RESEARCH LAB RDRL CIO LL 2800 POWDER MILL RD...ADELPHI MD 20783-1197 1 DIRECTOR US ARMY RESEARCH LAB RDRL CIO MT 2800 POWDER MILL RD ADELPHI MD 20783-1197 1 DIRECTOR US ARMY RESEARCH LAB RDRL D 2800

Interatomic potential at small internuclear distances. A simple formula for the screening constant

NASA Astrophysics Data System (ADS)

Zinoviev, A. N.

2017-09-01

A simple formula for estimating the screening constant has been proposed. This formula fits well experimental data on the interaction potentials. Quantitative description of the experiment for the effect of electronic screening on the nuclear synthesis reaction cross-section for the D+-D system has been obtained. A conclusion has been made that the differences between the measured cross-sections and their theoretically predicted values, which take place in more complicated cases nuclear synthesis reactions, are not caused by uncertainties in the knowledge of potentials.

Nanopowder synthesis based on electric explosion technology

NASA Astrophysics Data System (ADS)

Kryzhevich, D. S.; Zolnikov, K. P.; Korchuganov, A. V.; Psakhie, S. G.

2017-10-01

A computer simulation of the bicomponent nanoparticle formation during the electric explosion of copper and nickel wires was carried out. The calculations were performed in the framework of the molecular dynamics method using many-body potentials of interatomic interaction. As a result of an electric explosion of dissimilar metal wires, bicomponent nanoparticles having different stoichiometry and a block structure can be formed. It is possible to control the process of destruction and the structure of the formed bicomponent nanoparticles by varying the distance between the wires and the loading parameters.

Relationship between covalence and displacive phase transition temperature in RAO{sub 4} and LiAO{sub 3} (R = rare-earth element and A = Nb and Ta)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tsunekawa, S.; Fukuda, T.; Kimiyama, T.

Crystal structure analyses by TOF neutron powder diffraction are performed for R TaO{sub 4} (R = rare-earth element) and the Ta-O interatomic distances are determined. The relationship between the covalency of A-O bonds (A = Nb and Ta), which show the most shortening upon phase transition, and the transition temperature is discussed for RAO{sub 4} and LiAO{sub 3}, and the parameters of Ta-O covalence are determined.

AB Initi Molecular Dynamics Simulation of the Amorphous Structure of Ca-Mg-Cu and Ca-Mg-Zn Alloys (Postprint)

DTIC Science & Technology

2012-09-26

characteristic coordination poly- hedra present in an amorphous structure.[23,42] A coor- dination polyhedron is defined as an i-centered cluster with...vertices at the first-shell atom positions and edges coinciding with the interatomic bonds in the first shell.[45] Each coordination polyhedron can be...assigned a Voronoi signature (n3, n4, n5, n6), where nm is the number of vertices common to m polyhedron faces (or edges).[46] m is also called the

Voroprot: an interactive tool for the analysis and visualization of complex geometric features of protein structure.

PubMed

Olechnovic, Kliment; Margelevicius, Mindaugas; Venclovas, Ceslovas

2011-03-01

We present Voroprot, an interactive cross-platform software tool that provides a unique set of capabilities for exploring geometric features of protein structure. Voroprot allows the construction and visualization of the Apollonius diagram (also known as the additively weighted Voronoi diagram), the Apollonius graph, protein alpha shapes, interatomic contact surfaces, solvent accessible surfaces, pockets and cavities inside protein structure. Voroprot is available for Windows, Linux and Mac OS X operating systems and can be downloaded from http://www.ibt.lt/bioinformatics/voroprot/.

Enhanced squeezing of a collective spin via control of its qudit subsystems.

PubMed

Norris, Leigh M; Trail, Collin M; Jessen, Poul S; Deutsch, Ivan H

2012-10-26

Unitary control of qudits can improve the collective spin squeezing of an atomic ensemble. Preparing the atoms in a state with large quantum fluctuations in magnetization strengthens the entangling Faraday interaction. The resulting increase in interatomic entanglement can be converted into metrologically useful spin squeezing. Further control can squeeze the internal atomic spin without compromising entanglement, providing an overall multiplicative factor in the collective squeezing. We model the effects of optical pumping and study the tradeoffs between enhanced entanglement and decoherence. For realistic parameters we see improvements of ~10 dB.

Nitroxoline Molecule: Planar or Not? A Story of Battle between π-π Conjugation and Interatomic Repulsion.

PubMed

Tikhonov, Denis S; Sharapa, Dmitry I; Otlyotov, Arseniy A; Solyankin, Peter M; Rykov, Anatolii N; Shkurinov, Alexander P; Grikina, Olga E; Khaikin, Leonid S

2018-02-15

The conformational properties of the nitro group in nitroxoline (8-hydroxy-5-nitroquinoline, NXN) were investigated in the gas phase by means of gas electron diffraction (GED) and quantum chemical calculations, and also with solid-state analysis performed using terahertz time-domain spectroscopy (THz-TDS). The results of the GED refinement show that in the equilibrium structure the NO 2 group is twisted by angle ϕ = 8 ± 3° with respect to the 8-hydroxyoquinoline plane. This is the result of interatomic repulsion of oxygen in the NO 2 group from the closest hydrogen, which overcomes the energy gain from the π-π conjugation of the nitro group and aromatic system of 8-hydroxyoquinoline. The computation of equilibrium geometry using MP2/cc-pVXZ (X = T, Q) shows a large overestimation of the ϕ value, while DFT with the cc-pVTZ basis set performs reasonably well. On the other hand, DFT computations with double-ζ basis sets yield a planar structure of NXN. The refined potential energy surface of the torsion vibration the of nitro group in the condensed phase derived from the THz-TDS data indicates the NXN molecule to be planar. This result stays in good agreement with the previous X-ray structure determination. The strength of the π-system conjugation for the NO 2 group and 8-hydroxyoquinoline is discussed using NBO analysis, being further supported by comparison of the refined semiexperimental gas-phase structure of NXN from GED with other nitrocompounds.

On the continuum mechanics approach for the analysis of single walled carbon nanotubes

NASA Astrophysics Data System (ADS)

Chaudhry, M. S.; Czekanski, A.

2016-04-01

Today carbon nanotubes have found various applications in structural, thermal and almost every field of engineering. Carbon nanotubes provide great strength, stiffness resilience properties. Evaluating the structural behavior of nanoscale materials is an important task. In order to understand the materialistic behavior of nanotubes, atomistic models provide a basis for continuum mechanics modelling. Although the properties of bulk materials are consistent with the size and depends mainly on the material but the properties when we are in Nano-range, continuously change with the size. Such models start from the modelling of interatomic interaction. Modelling and simulation has advantage of cost saving when compared with the experiments. So in this project our aim is to use a continuum mechanics model of carbon nanotubes from atomistic perspective and analyses some structural behaviors of nanotubes. It is generally recognized that mechanical properties of nanotubes are dependent upon their structural details. The properties of nanotubes vary with the varying with the interatomic distance, angular orientation, radius of the tube and many such parameters. Based on such models one can analyses the variation of young's modulus, strength, deformation behavior, vibration behavior and thermal behavior. In this study some of the structural behaviors of the nanotubes are analyzed with the help of continuum mechanics models. Using the properties derived from the molecular mechanics model a Finite Element Analysis of carbon nanotubes is performed and results are verified. This study provides the insight on continuum mechanics modelling of nanotubes and hence the scope to study the effect of various parameters on some structural behavior of nanotubes.

Theoretical study of electron correlation effects in transition metal dimers

NASA Technical Reports Server (NTRS)

Das, G. P.; Jaffe, R. L.

1984-01-01

Introduction of partially localized orbitals is shown to reduce the number of terms needed to describe the bonding in transition metal clusters. Using this formalism, it is possible to compute the various intra- and inter-atomic electron correlation contributions to the bond energy. Calculations demonstrate the relative importance of several kinds of electron correlation terms involving the 3p, 3d, and 4s electrons. Improved interaction potentials are obtained for the dimers V(2) and Cr(2) when additional correlation is added to the CAS SCF results of Walch, Bauschlicher, Roos, and Nelin (1983).

Helium Purge Flow Prevention of Atmospheric Contamination of the Cryogenically Cooled Optics of Orbiting Infrared Telescopes: Calculation of He-O Differential Cross Section,

DTIC Science & Technology

1981-06-05

interactions. Aquilanti and coworkers were able to obtain two analytic forms for the interatomic potential --a Lennard - Jones (12, 6) and an exp(a, 6) function...Sec. UI.D 38 ences between the 3R and 3E- potential functions which described the interac- tions of ground-state oxygen and helium atoms. Instead, for...AO-AIOI 152 AEROSPACE CORP EL SEUMOO CA CHEMISTRY AND PHYSICS LAD r/6 17 HELIUM PLRE FLOW PREVENTION OF ATMOSPHERIC CONTAMINATION OF TAR fTCiO )JN81

Virial Coefficients for the Liquid Argon

NASA Astrophysics Data System (ADS)

Korth, Micheal; Kim, Saesun

2014-03-01

We begin with a geometric model of hard colliding spheres and calculate probability densities in an iterative sequence of calculations that lead to the pair correlation function. The model is based on a kinetic theory approach developed by Shinomoto, to which we added an interatomic potential for argon based on the model from Aziz. From values of the pair correlation function at various values of density, we were able to find viral coefficients of liquid argon. The low order coefficients are in good agreement with theoretical hard sphere coefficients, but appropriate data for argon to which these results might be compared is difficult to find.

Evaporation kinetics of Mg2SiO4 crystals and melts from molecular dynamics simulations

NASA Technical Reports Server (NTRS)

Kubicki, J. D.; Stolper, E. M.

1993-01-01

Computer simulations based on the molecular dynamics (MD) technique were used to study the mechanisms and kinetics of free evaporation from crystalline and molten forsterite (i.e., Mg2SiO4) on an atomic level. The interatomic potential employed for these simulations reproduces the energetics of bonding in forsterite and in gas-phase MgO and SiO2 reasonably accurately. Results of the simulation include predicted evaporation rates, diffusion rates, and reaction mechanisms for Mg2SiO4(s or l) yields 2Mg(g) + 20(g) + SiO2(g).

Dark-soliton dynamics in Bose-Einstein condensates at finite temperature

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jackson, B.; Proukakis, N. P.; Barenghi, C. F.

2007-05-15

The dynamics of a dark soliton in an elongated Bose-Einstein condensate is studied at finite temperatures. In addition to accurately reproducing all stages of the decay of the soliton observed in the experiment of Burger et al. [Phys. Rev. Lett. 83, 5198 (1999)], our numerical simulations reveal the existence of an experimentally accessible parameter regime for which phase-imprinted dark solitons can execute at least one full axial oscillation prior to their decay. The dependence of the decay time scale on temperature and initial soliton depth is analyzed and the role of interatomic collisions quantified.

Formation of 2D nanoparticles with block structure in simultaneous electric explosion of conductors

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kryzhevich, Dmitrij S., E-mail: kryzhev@ispms.ru, E-mail: kost@ispms.ru; Zolnikov, Konstantin P., E-mail: kryzhev@ispms.ru, E-mail: kost@ispms.ru; Abdrashitov, Andrei V.

2014-11-14

A molecular dynamics simulation of nanoparticle formation in simultaneous electric explosion of conductors is performed. Interatomic interaction is described using potentials calculated in the framework of the embedded atom method. High-rate heating results in failure of the conductors with the formation of nanoparticles. The influence of the heating rate, temperature distribution over the specimen cross-section and the distance between simultaneously exploded conductors on the structure of formed nanoparticles is studied. The calculation results show that the electric explosion of conductors allows the formation of nanoparticles with block structure.

Pair-correlation function of a metastable helium Bose-Einstein condensate

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zin, Pawel; Trippenbach, Marek; Gajda, Mariusz

2004-02-01

The pair-correlation function is one of the basic quantities to characterize the coherence properties of a Bose-Einstein condensate. We calculate this function in the experimentally important case of a zero temperature Bose-Einstein condensate in a metastable triplet helium state using the variational method with a pair-excitation ansatz. We compare our result with a pair-correlation function obtained for the hard-sphere potential with the same scattering length. Both functions are practically indistinguishable for distances greater than the scattering length. At smaller distances, due to interatomic interactions, the helium condensate shows strong correlations.

Hydrides of intermetallic compounds with a H/M ratio greater than unity obtained at high hydrogen pressures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Semenenko, K.N.; Klyamkin, S.N.

1993-11-01

Novel hydride phases with H/M > 1 based on Zr{sub 2}Pd, Hf{sub 2}Pd, and Hf{sub 2}Cu (structures of the MoSi{sub 2} type) have been synthesized at high H{sub 2} pressures. The X-ray diffraction investigations of the resulting hydrides have been carried out. Some factors determining the maximum hydrogen content in the hydrides of intermetallic compounds are discussed. A model structure of the hydrides obtained is proposed, which assumes the possibility of direct H-H interactions when the interatomic distances are less than 1 {angstrom}.

Effect of biaxial strain on the magnetism of Fe16N2: Density-functional investigations

NASA Astrophysics Data System (ADS)

Zhou, Wei; Liu, Lijuan; Wu, Ping

2014-02-01

The effect of biaxial strain on the magnetism of α″-Fe16N2 was investigated by the first principles calculations. The GGA, GGA + U and HSE06 calculations give the same result that the magnetic moments increase with the biaxial strain in the ab plane. All non-equivalent Fe atoms contribute to the increase of magnetic moments, although the variations of inter-atomic distances between non-equivalent Fe and N are different. Additionally, the magnetic anisotropy of Fe16N2 could be controlled by the biaxial strain.

Non-Contact Acousto-Thermal Signatures of Plastic Deformation in TI-6AL-4V

NASA Astrophysics Data System (ADS)

Welter, J. T.; Malott, G.; Schehl, N.; Sathish, S.; Jata, K. V.; Blodgett, M. P.

2010-02-01

Plastic deformation introduces changes in a material which include increases in: dislocations, strains, residual stress, and yield stress. However, these changes have a very small impact on the material properties such as elastic modulus, conductivity and ultrasonic wave speed. This is due to the fact that interatomic forces govern these properties, and they are not affected by plastic deformation to any large degree. This is evident from the fact that the changes in electrical resistance and ultrasonic velocity in plastically deformed and virgin samples are very small and can only be determined by highly controlled experiments. Except for X-ray diffraction, there are no direct nondestructive methods for measuring strain and the residual stress. This paper presents an application of the non-contact acousto-thermal signature (NCATS) NDE methodology to detect plastic deformation in flat dog bone Ti-6Al-4V samples. Results of the NCATS measurements on samples subjected to incremental amounts of plastic deformation are presented. The maximum temperature attained by the sample due to acoustic excitation is found to be sensitive to the amount of plastic strain. It is observed that the temperature induced by acoustic excitation increases to a peak followed by a decrease to failure. The maximum temperature peak occurs at plastic strains of 12-14%. It is observed that there is a correlation between the peak in maximum temperature rise and the strain at the experimentally determined ultimate tensile strength. A microstructural based explanation for this will be presented. The results are discussed in reference to utilizing this technique for detection and evaluation of plastic deformation.

Cracks dynamics under tensional stress - a DEM approach

NASA Astrophysics Data System (ADS)

Debski, Wojciech; Klejment, Piotr; Kosmala, Alicja; Foltyn, Natalia; Szpindler, Maciej

2017-04-01

Breaking and fragmentation of solid materials is an extremely complex process involving scales ranging from an atomic scale (breaking inter-atomic bounds) up to thousands of kilometers in case of catastrophic earthquakes (in energy scale it ranges from single eV up to 1024 J). Such a large scale span of breaking processes opens lot of questions like, for example, scaling of breaking processes, existence of factors controlling final size of broken area, existence of precursors, dynamics of fragmentation, to name a few. The classical approach to study breaking process at seismological scales, i.e., physical processes in earthquake foci, is essentially based on two factors: seismic data (mostly) and the continuum mechanics (including the linear fracture mechanics). Such approach has been gratefully successful in developing kinematic (first) and dynamic (recently) models of seismic rupture and explaining many of earthquake features observed all around the globe. However, such approach will sooner or latter face a limitation due to a limited information content of seismic data and inherit limitations of the fracture mechanics principles. A way of avoiding this expected limitation is turning an attention towards a well established in physics method of computational simulations - a powerful branch of contemporary physics. In this presentation we discuss preliminary results of analysis of fracturing dynamics under external tensional forces using the Discrete Element Method approach. We demonstrate that even under a very simplified tensional conditions, the fragmentation dynamics is a very complex process, including multi-fracturing, spontaneous fracture generation and healing, etc. We also emphasis a role of material heterogeneity on the fragmentation process.

Reflection thermal diffuse x-ray scattering for quantitative determination of phonon dispersion relations

DOE PAGES

Mei, A. B.; Hellman, O.; Schlepuetz, C. M.; ...

2015-11-03

Synchrotron reflection x-ray thermal diffuse scattering (TDS) measurements, rather than previously reported transmission TDS, are carried out at room temperature and analyzed using a formalism based upon second-order interatomic force constants and long-range Coulomb interactions to obtain quantitative determinations of MgO phonon dispersion relations (h) over bar omega(j) (q), phonon densities of states g((h) over bar omega), and isochoric temperature-dependent vibrational heat capacities c v (T). We use MgO as a model system for investigating reflection TDS due to its harmonic behavior as well as its mechanical and dynamic stability. Resulting phonon dispersion relations and densities of states are foundmore » to be in good agreement with independent reports from inelastic neutron and x-ray scattering experiments. Temperature-dependent isochoric heat capacities c v (T), computed within the harmonic approximation from (h) over bar omega(j) (q) values, increase with temperature from 0.4 x 10 -4 eV/atom K at 100 K to 1.4 x 10 -4 eV/atom K at 200 K and 1.9 x 10 -4 eV/atom K at 300 K, in excellent agreement with isobaric heat capacity values c p (T) between 4 and 300 K. We anticipate that the experimental approach developed here will be valuable for determining vibrational properties of heteroepitaxial thin films since the use of grazing-incidence (θ ≲ θ c where θ c is the density-dependent critical angle) allows selective tuning of x-ray penetration depths to ≲ 10 nm.« less

Enhanced initial protein adsorption on engineered nanostructured cubic zirconia.

PubMed

Sabirianov, R F; Rubinstein, A; Namavar, F

2011-04-14

Motivated by experimentally-observed biocompatibility enhancement of nanoengineered cubic zirconia (ZrO(2)) coatings to mesenchymal stromal cells, we have carried out computational analysis of the initial immobilization of one known structural fragment of the adhesive protein (fibronectin) on the corresponding surface. We constructed an atomistic model of the ZrO(2) nano-hillock of 3-fold symmetry based on Atom Force Microscopy and Transmission Electron Microscopy images. First principle quantum mechanical calculations show a substantial variation of electrostatic potential at the hillock due to the presence of surface features such as edges and vertexes. Using an implemented Monte Carlo simulated annealing method, we found the orientation of the immobilized protein on the ZrO(2) surface and the contribution of the amino acid residues from the protein sequence to the adsorption energy. Accounting for the variation of the dielectric permittivity at the protein-implant interface, we used a model distance-dependent dielectric function to describe the inter-atom electrostatic interactions in the adsorption potential. We found that the initial immobilization of the rigid protein fragment on the nanostructured pyramidal ZrO(2) surface is achieved with a magnitude of adsorption energy larger than that of the protein on the smooth (atomically flat) surface. The strong attractive electrostatic interactions are a major contributing factor in the enhanced adsorption at the nanostructured surface. In the case of adsorption on the flat, uncharged surface this factor is negligible. We show that the best electrostatic and steric fit of the protein to the inorganic surface corresponds to a minimum of the adsorption energy determined by the non-covalent interactions.

Re-evaluation of low-resolution crystal structures via interactive molecular-dynamics flexible fitting (iMDFF): a case study in complement C4.

PubMed

Croll, Tristan Ian; Andersen, Gregers Rom

2016-09-01

While the rapid proliferation of high-resolution structures in the Protein Data Bank provides a rich set of templates for starting models, it remains the case that a great many structures both past and present are built at least in part by hand-threading through low-resolution and/or weak electron density. With current model-building tools this task can be challenging, and the de facto standard for acceptable error rates (in the form of atomic clashes and unfavourable backbone and side-chain conformations) in structures based on data with dmax not exceeding 3.5 Å reflects this. When combined with other factors such as model bias, these residual errors can conspire to make more serious errors in the protein fold difficult or impossible to detect. The three recently published 3.6-4.2 Å resolution structures of complement C4 (PDB entries 4fxg, 4fxk and 4xam) rank in the top quartile of structures of comparable resolution both in terms of Rfree and MolProbity score, yet, as shown here, contain register errors in six β-strands. By applying a molecular-dynamics force field that explicitly models interatomic forces and hence excludes most physically impossible conformations, the recently developed interactive molecular-dynamics flexible fitting (iMDFF) approach significantly reduces the complexity of the conformational space to be searched during manual rebuilding. This substantially improves the rate of detection and correction of register errors, and allows user-guided model building in maps with a resolution lower than 3.5 Å to converge to solutions with a stereochemical quality comparable to atomic resolution structures. Here, iMDFF has been used to individually correct and re-refine these three structures to MolProbity scores of <1.7, and strategies for working with such challenging data sets are suggested. Notably, the improved model allowed the resolution for complement C4b to be extended from 4.2 to 3.5 Å as demonstrated by paired refinement.

ReaxFF Grand Canonical Monte Carlo simulation of adsorption and dissociation of oxygen on platinum (111)

NASA Astrophysics Data System (ADS)

Valentini, Paolo; Schwartzentruber, Thomas E.; Cozmuta, Ioana

2011-12-01

Atomic-level Grand Canonical Monte Carlo (GCMC) simulations equipped with a reactive force field (ReaxFF) are used to study atomic oxygen adsorption on a Pt(111) surface. The off-lattice GCMC calculations presented here rely solely on the interatomic potential and do not necessitate the pre-computation of surface adlayer structures and their interpolation. As such, they provide a predictive description of adsorbate phases. In this study, validation is obtained with experimental evidence (steric heats of adsorption and isotherms) as well as DFT-based state diagrams available in the literature. The ReaxFF computed steric heats of adsorption agree well with experimental data, and this study clearly shows that indirect dissociative adsorption of O2 on Pt(111) is an activated process at non-zero coverages, with an activation energy that monotonically increases with coverage. At a coverage of 0.25 ML, a highly ordered p(2 × 2) adlayer is found, in agreement with several low-energy electron diffraction observations. Isotherms obtained from the GCMC simulations compare qualitatively and quantitatively well with previous DFT-based state diagrams, but are in disagreement with the experimental data sets available. ReaxFF GCMC simulations at very high coverages show that O atoms prefer to bind in fcc hollow sites, at least up to 0.8 ML considered in the present work. At moderate coverages, little to no disorder appears in the Pt lattice. At high coverages, some Pt atoms markedly protrude out of the surface plane. This observation is in qualitative agreement with recent STM images of an oxygen covered Pt surface. The use of the GCMC technique based on a transferable potential is particularly valuable to produce more realistic systems (adsorbent and adsorbate) to be used in subsequent dynamical simulations (Molecular Dynamics) to address recombination reactions (via either Eley-Rideal or Langmuir-Hinshelwood mechanisms) on variously covered surfaces. By using GCMC and Molecular Dynamics simulations, the ReaxFF force field can be a valuable tool for understanding heterogeneous catalysis on a solid surface. Finally, the use of a reactive potential is a necessary requirement to investigate problems where dissociative adsorption occurs, as typical of many important catalytic processes.

Solar wind ion impacts into ice surfaces: A molecular-dynamics study using the REAX force field

NASA Astrophysics Data System (ADS)

Anders, Christian; Urbassek, Herbert M.

2017-01-01

Molecular dynamics simulation is used to study the effects of solar-wind ion irradiation on an ice target, focusing on the effects of nuclear energy deposition. A reactive force field (REAX) interatomic interaction potential is employed that allows us to model the breaking and formation of molecular bonds and hence to follow the chemistry occurring in the target. As ions we study H and He ions as typical constituents of the solar wind, and Ne as an example of a heavier ion; they impact at the speed of maximum flux in the solar wind, 400 km/s. The ice consists of a mixture of H2O, CO2, CH3OH and NH3. We find that molecular dissociations occur within 0.2 ps after ion impact and new products are formed up to a time of 1 ps; only water has a slower dynamics, due to highly mobile H atoms allowing for late recombinations. The number of dissociations, and hence also of product molecules increases from H over He to Ne ion projectiles and can be quantified by the amount of energy deposited in the target by these ions. The most abundant products formed include CO, OH and NH2. Reaction products are most complex for Ne impact, and include H3O, formaldehyde (H2CO), HO2, and NO. Formaldehyde is important as it is formed relatively frequently and is known as a precursor in the formation of sugars. In addition, molecules containing all CHON elements are formed, among which are CH2NO, CONH, methanolamine (CH5NO), and ethyne (C2H2). Repeated impacts generate novel, and more complex product species; we found CN, CH4, CH3NO, methylamine (CH3NH2), and acetamide (CH3CONH2), among others; the complex species are formed less frequently than the simple fragments. Sputtering occurs for all projectiles, even H. The ejecta are either original molecules - in particular CO2 - or simple fragments; only few product molecules are emitted.

Introducing ab initio based neural networks for transition-rate prediction in kinetic Monte Carlo simulations

NASA Astrophysics Data System (ADS)

Messina, Luca; Castin, Nicolas; Domain, Christophe; Olsson, Pär

2017-02-01

The quality of kinetic Monte Carlo (KMC) simulations of microstructure evolution in alloys relies on the parametrization of point-defect migration rates, which are complex functions of the local chemical composition and can be calculated accurately with ab initio methods. However, constructing reliable models that ensure the best possible transfer of physical information from ab initio to KMC is a challenging task. This work presents an innovative approach, where the transition rates are predicted by artificial neural networks trained on a database of 2000 migration barriers, obtained with density functional theory (DFT) in place of interatomic potentials. The method is tested on copper precipitation in thermally aged iron alloys, by means of a hybrid atomistic-object KMC model. For the object part of the model, the stability and mobility properties of copper-vacancy clusters are analyzed by means of independent atomistic KMC simulations, driven by the same neural networks. The cluster diffusion coefficients and mean free paths are found to increase with size, confirming the dominant role of coarsening of medium- and large-sized clusters in the precipitation kinetics. The evolution under thermal aging is in better agreement with experiments with respect to a previous interatomic-potential model, especially concerning the experiment time scales. However, the model underestimates the solubility of copper in iron due to the excessively high solution energy predicted by the chosen DFT method. Nevertheless, this work proves the capability of neural networks to transfer complex ab initio physical properties to higher-scale models, and facilitates the extension to systems with increasing chemical complexity, setting the ground for reliable microstructure evolution simulations in a wide range of alloys and applications.

EON: software for long time simulations of atomic scale systems

NASA Astrophysics Data System (ADS)

Chill, Samuel T.; Welborn, Matthew; Terrell, Rye; Zhang, Liang; Berthet, Jean-Claude; Pedersen, Andreas; Jónsson, Hannes; Henkelman, Graeme

2014-07-01

The EON software is designed for simulations of the state-to-state evolution of atomic scale systems over timescales greatly exceeding that of direct classical dynamics. States are defined as collections of atomic configurations from which a minimization of the potential energy gives the same inherent structure. The time evolution is assumed to be governed by rare events, where transitions between states are uncorrelated and infrequent compared with the timescale of atomic vibrations. Several methods for calculating the state-to-state evolution have been implemented in EON, including parallel replica dynamics, hyperdynamics and adaptive kinetic Monte Carlo. Global optimization methods, including simulated annealing, basin hopping and minima hopping are also implemented. The software has a client/server architecture where the computationally intensive evaluations of the interatomic interactions are calculated on the client-side and the state-to-state evolution is managed by the server. The client supports optimization for different computer architectures to maximize computational efficiency. The server is written in Python so that developers have access to the high-level functionality without delving into the computationally intensive components. Communication between the server and clients is abstracted so that calculations can be deployed on a single machine, clusters using a queuing system, large parallel computers using a message passing interface, or within a distributed computing environment. A generic interface to the evaluation of the interatomic interactions is defined so that empirical potentials, such as in LAMMPS, and density functional theory as implemented in VASP and GPAW can be used interchangeably. Examples are given to demonstrate the range of systems that can be modeled, including surface diffusion and island ripening of adsorbed atoms on metal surfaces, molecular diffusion on the surface of ice and global structural optimization of nanoparticles.

An automated method for modeling proteins on known templates using distance geometry.

PubMed

Srinivasan, S; March, C J; Sudarsanam, S

1993-02-01

We present an automated method incorporated into a software package, FOLDER, to fold a protein sequence on a given three-dimensional (3D) template. Starting with the sequence alignment of a family of homologous proteins, tertiary structures are modeled using the known 3D structure of one member of the family as a template. Homologous interatomic distances from the template are used as constraints. For nonhomologous regions in the model protein, the lower and the upper bounds for the interatomic distances are imposed by steric constraints and the globular dimensions of the template, respectively. Distance geometry is used to embed an ensemble of structures consistent with these distance bounds. Structures are selected from this ensemble based on minimal distance error criteria, after a penalty function optimization step. These structures are then refined using energy optimization methods. The method is tested by simulating the alpha-chain of horse hemoglobin using the alpha-chain of human hemoglobin as the template and by comparing the generated models with the crystal structure of the alpha-chain of horse hemoglobin. We also test the packing efficiency of this method by reconstructing the atomic positions of the interior side chains beyond C beta atoms of a protein domain from a known 3D structure. In both test cases, models retain the template constraints and any additionally imposed constraints while the packing of the interior residues is optimized with no short contacts or bond deformations. To demonstrate the use of this method in simulating structures of proteins with nonhomologous disulfides, we construct a model of murine interleukin (IL)-4 using the NMR structure of human IL-4 as the template. The resulting geometry of the nonhomologous disulfide in the model structure for murine IL-4 is consistent with standard disulfide geometry.

Intra- and inter-atomic optical transitions of Fe, Co, and Ni ferrocyanides studied using first-principles many-electron calculations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Watanabe, Shinta, E-mail: s-watanabe@nucl.nagoya-u.ac.jp, E-mail: j-onoe@nucl.nagoya-u.ac.jp; Sawada, Yuki; Nakaya, Masato

We have investigated the electronic structures and optical properties of Fe, Co, and Ni ferrocyanide nanoparticles using first-principles relativistic many-electron calculations. The overall features of the theoretical absorption spectra for Fe, Ni, and Co ferrocyanides calculated using a first-principles many-electron method well reproduced the experimental one. The origins of the experimental absorption spectra were clarified by performing a configuration analysis based on the many-electron wave functions. For Fe ferrocyanide, the experimental absorption peaks originated from not only the charge-transfer transitions from Fe{sup 2+} to Fe{sup 3+} but also the 3d-3d intra-transitions of Fe{sup 3+} ions. In addition, the spin crossovermore » transition of Fe{sup 3+} predicted by the many-electron calculations was about 0.24 eV. For Co ferrocyanide, the experimental absorption peaks were mainly attributed to the 3d-3d intra-transitions of Fe{sup 2+} ions. In contrast to the Fe and Co ferrocyanides, Ni ferrocyanide showed that the absorption peaks originated from the 3d-3d intra-transitions of Ni{sup 3+} ions in a low-energy region, while from both the 3d-3d intra-transitions of Fe{sup 2+} ions and the charge-transfer transitions from Fe{sup 2+} to Ni{sup 3+} in a high-energy region. These results were quite different from those of density-functional theory (DFT) calculations. The discrepancy between the results of DFT calculations and those of many-electron calculations suggested that the intra- and inter-atomic transitions of transition metal ions are significantly affected by the many-body effects of strongly correlated 3d electrons.« less

Frenkel pair recombinations in UO2: Importance of explicit description of polarizability in core-shell molecular dynamics simulations

NASA Astrophysics Data System (ADS)

Devynck, Fabien; Iannuzzi, Marcella; Krack, Matthias

2012-05-01

The oxygen and uranium Frenkel pair (FP) recombination mechanisms are studied in UO2 using an empirical interatomic potential accounting for the polarizability of the ions, namely a dynamical core-shell model. The results are compared to a more conventional rigid-ion model. Both model types have been implemented into the cp2k program package and thoroughly validated. The overall picture indicates that the FP recombination mechanism is a complex process involving several phenomena. The FP recombination can happen instantaneously when the distance between the interstitial and the vacancy is small or can be thermally activated at larger separation distances. However, other criteria can prevail over the interstitial-vacancy distance. The surrounding environment of the FP defect, the mechanical stiffness of the matrix, and the orientation of the migration path are shown to be major factors acting on the FP lifetime. The core-shell and rigid-ion models provide a similar qualitative description of the FP recombination mechanism. However, the FP stabilities determined by both models significantly differ in the lower temperature range considered. Indeed, the recombination time of the oxygen and uranium FPs can be up to an order of magnitude lower in the core-shell model at T=600 K and T=1800 K, respectively. These differences highlight the importance of the explicit description of polarizability on some crucial properties such as the resistance to amorphization. This refined description of the interatomic interactions would certainly affect the description of the recrystallization process following a displacement cascade. In turn, the self-healing phase would be better accounted for in the core-shell model and the misestimate inherent to the lack of polarizability in the rigid-ion model corrected.

Stress-induced crystal transition of poly(butylene succinate) studied by terahertz and low-frequency Raman spectroscopy and quantum chemical calculation

NASA Astrophysics Data System (ADS)

Tatsuoka, Seika; Sato, Harumi

2018-05-01

We measured terahertz (THz) and low-frequency Raman spectra of Poly (butylene succinate) (PBS) which shows the crystal transition from α to β by stretching. For the assignment of the absorption peaks in the low-frequency region, we performed quantum chemical calculations with Cartesian-coordinate tensor transfer (CCT) method. Four major peaks appeared in the THz spectra of PBS at around 58, 76, 90, and 100 cm-1, and in the low-frequency Raman spectra a peak was observed at 88 cm-1. The THz peak at 100 cm-1 and the Raman peak at 88 cm-1 show a shift to a lower wavenumber region with increasing temperature. The quantum chemical calculation of β crystal form reveals the new peak appears above 100 cm-1. It was found that two kinds of peaks overlapped at around 100 cm-1 in the THz spectra of PBS. One of them can be assigned to a weak hydrogen bond between the C=O and CH2 groups in the intermolecular chains, which is perpendicular to the molecular chain of the α crystal form. Another one showed a parallel polarization which can be assigned to the intramolecular interaction between O (ether) and H-C groups in the β crystal form. The position of the peak at around 100 cm-1 in the perpendicular polarization changed to a lower wavenumber region with stretching, because of the weakening of the intermolecular hydrogen bonding by increasing the interatomic distances. On the other hand, that of the parallel polarization shifts to a higher wavenumber region because of the shortening of the interatomic distance from α to β crystal form (the strength of the intramolecular hydrogen bonding became stronger) by stretching.

Analysis of oxygen binding-energy variations for BaO on W

NASA Astrophysics Data System (ADS)

Haas, G. A.; Shih, A.; Mueller, D.; Thomas, R. E.

Interatomic Auger analyses have been made of different forms of BaO layers on W substrates. Variations in Auger spectroscopy energies of the Ba4dBa5pO2p interatomic Auger transition were found to be largely governed by the O2p binding energy of the BaO adsorbate. This was illustrated by comparing results of the Auger data values with values derived from O2p binding energies using ultraviolet photoelectron spectroscopy. Very good agreement was observed not only for the W<100> substrate but also for the W<110> substrate which showed two oxygen-induced electronics state. Variations in binding energy were noted for different states of BaO lattice formation and for different amounts of oxidation, ranging from the transition of Ba to BaO and continuing to the BaO 2 stoichiometry and beyond. Effects were also reported for adsorbate alignment and thermal activation (i.e., reduction) of the oxidized state. An empirical relationship was found suggesting that the more tightly bound the O2p states of the BaO adsorbate were, the lower its work function would be. This link between binding energy and work function was observed to be valid not only for cases of poisoning by oxidation, but held as well during reactivation by the subsequent reduction of the oxide. In addition, this relationship also appeared to predict the low work function obtained through the introduction of substances such as Sc to the BaO-W system. Possible qualitative reasons which might contribute to this are discussed in terms of enhanced dipole effects and shifts in band structure.

Elucidation of the atomic-scale mechanism of the anisotropic oxidation rate of 4H-SiC between the (0001) Si-face and ( 000 1 ¯ ) C-face by using a new Si-O-C interatomic potential

NASA Astrophysics Data System (ADS)

Takamoto, So; Yamasaki, Takahiro; Ohno, Takahisa; Kaneta, Chioko; Hatano, Asuka; Izumi, Satoshi

2018-05-01

Silicon carbide (SiC) is an attractive semiconductor material for applications in power electronic devices. However, fabrication of a high-quality SiC/SiO2 interface has been a challenge. It is well-known that there is a great difference in the oxidation rate between the Si-face and the C-face and that the quality of oxide on the Si-face is greater than that on the C-face. However, the atomistic mechanism of the thermal oxidation of SiC remains to be solved. In this paper, a new Si-O-C interatomic potential was developed to reproduce the kinetics of the thermal oxidation of SiC. Using this newly developed potential, large-scale SiC oxidation simulations at various temperatures were performed. The results showed that the activation energy of the Si-face is much larger than that of the C-face. In the case of the Si-face, a flat and aligned interface structure including Si1+ was created. Based on the estimated activation energies of the intermediate oxide states, it is proposed that the stability of the flat interface structure is the origin of the high activation energy of the oxidation of the Si-face. In contrast, in the case of the C-face, it is found that the Si atom at the interface is easily pulled up by the O atoms. This process generates the disordered interface and decreases the activation energy of the oxidation. It is also proposed that many excess C atoms are created in the case of the C-face.

Computational experience with a parallel algorithm for tetrangle inequality bound smoothing.

PubMed

Rajan, K; Deo, N

1999-09-01

Determining molecular structure from interatomic distances is an important and challenging problem. Given a molecule with n atoms, lower and upper bounds on interatomic distances can usually be obtained only for a small subset of the 2(n(n-1)) atom pairs, using NMR. Given the bounds so obtained on the distances between some of the atom pairs, it is often useful to compute tighter bounds on all the 2(n(n-1)) pairwise distances. This process is referred to as bound smoothing. The initial lower and upper bounds for the pairwise distances not measured are usually assumed to be 0 and infinity. One method for bound smoothing is to use the limits imposed by the triangle inequality. The distance bounds so obtained can often be tightened further by applying the tetrangle inequality--the limits imposed on the six pairwise distances among a set of four atoms (instead of three for the triangle inequalities). The tetrangle inequality is expressed by the Cayley-Menger determinants. For every quadruple of atoms, each pass of the tetrangle inequality bound smoothing procedure finds upper and lower limits on each of the six distances in the quadruple. Applying the tetrangle inequalities to each of the (4n) quadruples requires O(n4) time. Here, we propose a parallel algorithm for bound smoothing employing the tetrangle inequality. Each pass of our algorithm requires O(n3 log n) time on a REW PRAM (Concurrent Read Exclusive Write Parallel Random Access Machine) with O(log(n)n) processors. An implementation of this parallel algorithm on the Intel Paragon XP/S and its performance are also discussed.

Abstract ID: 176 Geant4 implementation of inter-atomic interference effect in small-angle coherent X-ray scattering for materials of medical interest.

PubMed

Paternò, Gianfranco; Cardarelli, Paolo; Contillo, Adriano; Gambaccini, Mauro; Taibi, Angelo

2018-01-01

Advanced applications of digital mammography such as dual-energy and tomosynthesis require multiple exposures and thus deliver higher dose compared to standard mammograms. A straightforward manner to reduce patient dose without affecting image quality would be removal of the anti-scatter grid, provided that the involved reconstruction algorithms are able to take the scatter figure into account [1]. Monte Carlo simulations are very well suited for the calculation of X-ray scatter distribution and can be used to integrate such information within the reconstruction software. Geant4 is an open source C++ particle tracking code widely used in several physical fields, including medical physics [2,3]. However, the coherent scattering cross section used by the standard Geant4 code does not take into account the influence of molecular interference. According to the independent atomic scattering approximation (the so-called free-atom model), coherent radiation is indistinguishable from primary radiation because its angular distribution is peaked in the forward direction. Since interference effects occur between x-rays scattered by neighbouring atoms in matter, it was shown experimentally that the scatter distribution is affected by the molecular structure of the target, even in amorphous materials. The most important consequence is that the coherent scatter distribution is not peaked in the forward direction, and the position of the maximum is strongly material-dependent [4]. In this contribution, we present the implementation of a method to take into account inter-atomic interference in small-angle coherent scattering in Geant4, including a dedicated data set of suitable molecular form factor values for several materials of clinical interest. Furthermore, we present scatter images of simple geometric phantoms in which the Rayleigh contribution is rigorously evaluated. Copyright © 2017.

The relation between the Gross Pitaevskii and Bogoliubov descriptions of a dilute Bose gas

NASA Astrophysics Data System (ADS)

Leggett, A. J.

2003-07-01

I formulate a 'pseudo-paradox' in the theory of a dilute Bose gas with repulsive interactions: the standard expression for the ground state energy within the Gross-Pitaevskii (GP) approximation is lower than that in the Bogoliubov approximation, and hence, by the standard variational argument, the former should prima facie be a better approximation than the latter to the true ground state—a conclusion which is of course opposite to the established wisdom concerning this problem. It is shown that the pseudo-paradox is (unsurprisingly) resolved by a correct transcription of the two-body scattering theory to the many-body case; however, contrary to what appears to be a widespread belief, the resolution has nothing to do with any spurious ultraviolet divergences which result from the replacement of the true interatomic potential by a delta-function pseudopotential. Rather, it relates to an infrared divergence which has the consequence that (a) the most obvious form of the GP 'approximation' actually does not correspond to any well-defined ansatz for the many-body wavefunction, and (b) that the 'best shot' at such a wavefunction always produces an energy which exceeds, or at best equals, that calculated in the Bogoliubov approximation. In fact, the necessity of the latter may be seen as a consequence of the need to reduce the Fock term in the energy, which is absent in the two-particle problem but dominant in the many-body case; it does this by increasing the density correlations, at distances less than or approximately equal to the correlation length \\xi , above the value extrapolated from the two-body case. As a by-product I devise an alternative formulation of the Bogoliubov approximation which does not require the explicit replacement of the true interatomic potential by a delta-function pseudopotential.

Time evolution, Lamb shift, and emission spectra of spontaneous emission of two identical atoms

NASA Astrophysics Data System (ADS)

Wang, Da-Wei; Li, Zheng-Hong; Zheng, Hang; Zhu, Shi-Yao

2010-04-01

A unitary transformation method is used to investigate the dynamic evolution of two multilevel atoms, in the basis of symmetric and antisymmetric states, with one atom being initially prepared in the first excited state and the other in the ground state. The unitary transformation guarantees that our calculations are based on the ground state of the atom-field system and the self-energy is subtracted at the beginning. The total Lamb shifts of the symmetric and antisymmetric states are divided into transformed shift and dynamic shift. The transformed shift is due to emitting and reabsorbing of virtual photons, by a single atom (nondynamic single atomic shift) and between the two atoms (quasi-static shift). The dynamic shift is due to the emitting and reabsorbing of real photons, by a single atom (dynamic single atomic shift) and between the two atoms (dynamic interatomic shift). The emitting and reabsorbing of virtual and real photons between the two atoms result in the interatomic shift, which does not exist for the one-atom case. The spectra at the long-time limit are calculated. If the distance between the two atoms is shorter than or comparable to the wavelength, the strong coupling between the two atoms splits the spectrum into two peaks, one from the symmetric state and the other from the antisymmetric state. The origin of the red or blue shifts for the symmetric and antisymmetric states mainly lies in the negative or positive interaction energy between the two atoms. In the investigation of the short time evolution, we find the modification of the effective density of states by the interaction between two atoms can modulate the quantum Zeno and quantum anti-Zeno effects in the decays of the symmetric and antisymmetric states.

Optimal utilization of total elastic scattering cross section data for the determination of interatomic potentials

NASA Technical Reports Server (NTRS)

Bernstein, R. B.; Labudde, R. A.

1972-01-01

The problem of inversion is considered in relation to absolute total cross sections Q(v) for atom-atom collisions and their velocity dependence, and the glory undulations and the transition to high velocity behavior. There is a limit to the amount of information available from Q(v) even when observations of good accuracy (e.g., + or - 0.25%) are in hand over an extended energy range (from thermal energies upward by a factor of greater than 1000 in relative kinetic energy). Methods were developed for data utilization, which take full advantage of the accuracy of the experimental Q(v) measurements.

The Strength of the Metal. Aluminum Oxide Interface

NASA Technical Reports Server (NTRS)

Pepper, S. V.

1984-01-01

The strength of the interface between metals and aluminum oxide is an important factor in the successful operation of devices found throughout modern technology. One finds the interface in machine tools, jet engines, and microelectronic integrated circuits. The strength of the interface, however, should be strong or weak depending on the application. The diverse technological demands have led to some general ideas concerning the origin of the interfacial strength, and have stimulated fundamental research on the problem. Present status of our understanding of the source of the strength of the metal - aluminum oxide interface in terms of interatomic bonds are reviewed. Some future directions for research are suggested.

Atomistic modeling for interfacial properties of Ni-Al-V ternary system

NASA Astrophysics Data System (ADS)

Dong, Wei-ping; Lee, Byeong-Joo; Chen, Zheng

2014-05-01

Interatomic potentials for Ni-Al-V ternary systems have been developed based on the second-nearest-neighbor modified embedded-atom method potential formalism. The potentials can describe various fundamental physical properties of the relevant materials in good agreement with experimental information. The potential is utilized for an atomistic computation of interfacial properties of Ni-Al-V alloys. It is found that vanadium atoms segregate on the γ-fcc/L12 interface and this segregation affects the interfacial properties. The applicability of the atomistic approach to an elaborate alloy design of advanced Ni-based superalloys through the investigation of the effect of alloying elements on interfacial properties is discussed.

A method of solid-solid phase equilibrium calculation by molecular dynamics

NASA Astrophysics Data System (ADS)

Karavaev, A. V.; Dremov, V. V.

2016-12-01

A method for evaluation of solid-solid phase equilibrium curves in molecular dynamics simulation for a given model of interatomic interaction is proposed. The method allows to calculate entropies of crystal phases and provides an accuracy comparable with that of the thermodynamic integration method by Frenkel and Ladd while it is much simpler in realization and less intense computationally. The accuracy of the proposed method was demonstrated in MD calculations of entropies for EAM potential for iron and for MEAM potential for beryllium. The bcc-hcp equilibrium curves for iron calculated for the EAM potential by the thermodynamic integration method and by the proposed one agree quite well.

Final Report for DE-FG02-99ER45795

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wilkins, John Warren

The research supported by this grant focuses on atomistic studies of defects, phase transitions, electronic and magnetic properties, and mechanical behaviors of materials. We have been studying novel properties of various emerging nanoscale materials on multiple levels of length and time scales, and have made accurate predictions on many technologically important properties. A significant part of our research has been devoted to exploring properties of novel nano-scale materials by pushing the limit of quantum mechanical simulations, and development of a rigorous scheme to design accurate classical inter-atomic potentials for larger scale atomistic simulations for many technologically important metals and metalmore » alloys.« less

Coupled-resonator waveguide perfect transport single-photon by interatomic dipole-dipole interaction

NASA Astrophysics Data System (ADS)

Yan, Guo-an; Lu, Hua; Qiao, Hao-xue; Chen, Ai-xi; Wu, Wan-qing

2018-06-01

We theoretically investigate single-photon coherent transport in a one-dimensional coupled-resonator waveguide coupled to two quantum emitters with dipole-dipole interactions. The numerical simulations demonstrate that the transmission spectrum of the photon depends on the two atoms dipole-dipole interactions and the photon-atom couplings. The dipole-dipole interactions may change the dip positions in the spectra and the coupling strength may broaden the frequency band width in the transmission spectrum. We further demonstrate that the typical transmission spectra split into two dips due to the dipole-dipole interactions. This phenomenon may be used to manufacture new quantum waveguide devices.

Origin of Noncubic Scaling Law in Disordered Granular Packing.

PubMed

Xia, Chengjie; Li, Jindong; Kou, Binquan; Cao, Yixin; Li, Zhifeng; Xiao, Xianghui; Fu, Yanan; Xiao, Tiqiao; Hong, Liang; Zhang, Jie; Kob, Walter; Wang, Yujie

2017-06-09

Recent diffraction experiments on metallic glasses have unveiled an unexpected noncubic scaling law between density and average interatomic distance, which led to the speculation of the presence of fractal glass order. Using x-ray tomography we identify here a similar noncubic scaling law in disordered granular packing of spherical particles. We find that the scaling law is directly related to the contact neighbors within the first nearest neighbor shell, and, therefore, is closely connected to the phenomenon of jamming. The seemingly universal scaling exponent around 2.5 arises due to the isostatic condition with a contact number around 6, and we argue that the exponent should not be universal.

Metal-atom Interactions and Clustering in Organic Semiconductor Systems

NASA Astrophysics Data System (ADS)

Tomita, Yoko; Park, Tea-uk; Nakayama, Takashi

2017-07-01

The interatomic interactions and clustering of metal atoms have been studied by first-principles calculations in graphene, pentacene, and polyacetylene as representative organic systems. It is shown that long-range repulsive Coulomb interaction appears between metal atoms with small electronegativity such as Al due to their ionization on host organic molecules, inducing their scattered distribution in organic systems. On the other hand, metal atoms with large electronegativity such as Au are weakly bonded to organic molecules, easily diffuse in molecular solids, and prefer to combine with each other owing to their short-range strong metallic-bonding interaction, promoting metal cluster generation in organic systems.

Consecutive magnetic phase diagram of UCoGe-URhGe-UIrGe system

NASA Astrophysics Data System (ADS)

Pospíšil, Jiří; Haga, Yoshinori; Miyake, Atsushi; Kambe, Shinsaku; Tateiwa, Naoyuki; Tokunaga, Yo; Honda, Fuminori; Nakamura, Ai; Homma, Yoshiya; Tokunaga, Masashi; Aoki, Dai; Yamamoto, Etsuji

2018-05-01

We prepared single crystals in UCo1-xRhxGe and UIr1-xRhxGe systems to establish a complex dU-U-T (dU-U is the shortest interatomic uranium distance and T is temperature) magnetic phase diagram. This recognized a characteristic maximum in magnetic susceptibility at temperature Tmax along the b axis as an important parameter. Three magnetically ordered regions can be distinguished within this scope; first a ferromagnetic region with Curie temperature

Observations of density fluctuations in an elongated Bose gas: ideal gas and quasicondensate regimes.

PubMed

Esteve, J; Trebbia, J-B; Schumm, T; Aspect, A; Westbrook, C I; Bouchoule, I

2006-04-07

We report in situ measurements of density fluctuations in a quasi-one-dimensional 87Rb Bose gas at thermal equilibrium in an elongated harmonic trap. We observe an excess of fluctuations compared to the shot-noise level expected for uncorrelated atoms. At low atomic density, the measured excess is in good agreement with the expected "bunching" for an ideal Bose gas. At high density, the measured fluctuations are strongly reduced compared to the ideal gas case. We attribute this reduction to repulsive interatomic interactions. The data are compared with a calculation for an interacting Bose gas in the quasicondensate regime.

Dynamical Correlation In Some Liquid Alkaline Earth Metals Near Melting

NASA Astrophysics Data System (ADS)

Thakore, B. Y.; Suthar, P. H.; Khambholja, S. G.; Gajjar, P. N.; Jani, A. R.

2010-12-01

The study of dynamical variables: velocity autocorrelation function (VACF) and power spectrum of liquid alkaline earth metals (Ca, Sr, and Ba) have been presented based on the static harmonic well approximation. The effective interatomic potential for liquid metals is computed using our well recognized model potential with the exchange correlation functions due to Hartree, Taylor, Ichimaru and Utsumi, Farid et al. and Sarkar et al. It is observed that the VACF computed using Sarkar et al. gives the good agreement with available molecular dynamics simulation (MD) results [Phys Rev. B 62, 14818 (2000)]. The shoulder of the power spectrum depends upon the type of local field correlation function used.

g_contacts: Fast contact search in bio-molecular ensemble data

NASA Astrophysics Data System (ADS)

Blau, Christian; Grubmuller, Helmut

2013-12-01

Short-range interatomic interactions govern many bio-molecular processes. Therefore, identifying close interaction partners in ensemble data is an essential task in structural biology and computational biophysics. A contact search can be cast as a typical range search problem for which efficient algorithms have been developed. However, none of those has yet been adapted to the context of macromolecular ensembles, particularly in a molecular dynamics (MD) framework. Here a set-decomposition algorithm is implemented which detects all contacting atoms or residues in maximum O(Nlog(N)) run-time, in contrast to the O(N2) complexity of a brute-force approach. Catalogue identifier: AEQA_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEQA_v1_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 8945 No. of bytes in distributed program, including test data, etc.: 981604 Distribution format: tar.gz Programming language: C99. Computer: PC. Operating system: Linux. RAM: ≈Size of input frame Classification: 3, 4.14. External routines: Gromacs 4.6[1] Nature of problem: Finding atoms or residues that are closer to one another than a given cut-off. Solution method: Excluding distant atoms from distance calculations by decomposing the given set of atoms into disjoint subsets. Running time:≤O(Nlog(N)) References: [1] S. Pronk, S. Pall, R. Schulz, P. Larsson, P. Bjelkmar, R. Apostolov, M. R. Shirts, J.C. Smith, P. M. Kasson, D. van der Spoel, B. Hess and Erik Lindahl, Gromacs 4.5: a high-throughput and highly parallel open source molecular simulation toolkit, Bioinformatics 29 (7) (2013).

Mutational analysis of microbial hydroxycinnamoyl-CoA hydratase-lyase (HCHL) towards enhancement of binding affinity: A computational approach.

PubMed

Kumar, Pravin; Ghosh Sachan, Shashwati; Poddar, Raju

2017-10-01

Improving the industrial enzyme for better yield of the product is important and a challenging task. One of such important industrial enzymes is microbial Hydroxycinnamoyl-CoA hydratase-lyase (HCHL). It converts feruloyl-CoA to vanillin. We place our efforts towards the improvement of its catalytic activity with comprehensive computational investigation. Catalytic core of the HCHL was explored with molecular modeling and docking approaches. Site-directed mutations were introduced in the catalytic site of HCHL in a sequential manner to generate different mutants of HCHL. Basis of mutation is to increase the interaction between HCHL and substrate feruloyl-CoA through interatomic forces and hydrogen bond formation. A rigorous molecular dynamics (MD) simulation was performed to check the stability of mutant's structure. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), dynamic cross correlation (DCCM) and principal component analysis (PCA) were also performed to analyze flexibility and stability of structures. Docking studies were carried out between different mutants of HCHL and feruloyl-CoA. Investigation of the different binding sites and the interactions with mutant HCHLs and substrate allowed us to highlight the improved performance of mutants than wild type HCHL. This was further validated with MD simulation of complex consisting of different mutants and substrate. It further confirms all the structures are stable. However, mutant-2 showed better affinity towards substrate by forming hydrogen bond between active site and feruloyl-CoA. We propose that increase in hydrogen bond formation might facilitate in dissociation of vanillin from feruloyl-CoA. The current work may be useful for the future development of 'tailor-made' enzymes for better yield of vanillin. Copyright © 2017 Elsevier Inc. All rights reserved.

Pseudopotential theoretical study of the alkali metals under arbitrary pressure: Density, bulk modulus, and shear moduli

NASA Astrophysics Data System (ADS)

Rasky, Daniel J.; Milstein, Frederick

1986-02-01

Milstein and Hill previously derived formulas for computing the bulk and shear moduli, κ, μ, and μ', at arbitrary pressures, for cubic crystals in which interatomic interaction energies are modeled by pairwise functions, and they carried out the moduli computations using the complete family of Morse functions. The present study extends their work to a pseudopotential description of atomic binding. Specifically: (1) General formulas are derived for determining these moduli under hydrostatic loading within the framework of a pseudopotential model. (2) A two-parameter pseudopotential model is used to describe atomic binding of the alkali metals, and the two parameters are determined from experimental data (the model employs the Heine-Abarenkov potential with the Taylor dielectric function). (3) For each alkali metal (Li, Na, K, Rb, and Cs), the model is used to compute the pressure-versus-volume behavior and, at zero pressure, the binding energy, the density, and the elastic moduli and their pressure derivatives; the theoretical behavior is found to be in excellent agreement with experiment. (4) Calculations are made of κ, μ, and μ' of the bcc alkali metals over wide ranges of hydrostatic compression and expansion. (5) The pseudopotential results are compared with those of arbitrary-central-force models (wherein κ-(2/3)μ=μ'+2P) and with the specific Morse-function results. The pressures, bulk moduli, and zero-pressure shear moduli (as determined for the Morse and pseudopotential models) are in excellent agreement, but important differences appear in the shear moduli under high compressions. The computations in the present paper are for the bcc metals; a subsequent paper will extend this work to include both the bcc and fcc structures, at compressions and expansions where elastic stability or lattice cohesion is, in practice, lost.

A first principles study of the electronic structure, elastic and thermal properties of UB2

NASA Astrophysics Data System (ADS)

Jossou, Ericmoore; Malakkal, Linu; Szpunar, Barbara; Oladimeji, Dotun; Szpunar, Jerzy A.

2017-07-01

Uranium diboride (UB2) has been widely deployed for refractory use and is a proposed material for Accident Tolerant Fuel (ATF) due to its high thermal conductivity. However, the applicability of UB2 towards high temperature usage in a nuclear reactor requires the need to investigate the thermomechanical properties, and recent studies have failed in highlighting applicable properties. In this work, we present an in-depth theoretical outlook of the structural and thermophysical properties of UB2, including but not limited to elastic, electronic and thermal transport properties. These calculations were performed within the framework of Density Functional Theory (DFT) + U approach, using Quantum ESPRESSO (QE) code considering the addition of Coulomb correlations on the uranium atom. The phonon spectra and elastic constant analysis show the dynamic and mechanical stability of UB2 structure respectively. The electronic structure of UB2 was investigated using full potential linear augmented plane waves plus local orbitals method (FP-LAPW+lo) as implemented in WIEN2k code. The absence of a band gap in the total and partial density of states confirms the metallic nature while the valence electron density plot reveals the presence of covalent bond between adjacent B-B atoms. We predicted the lattice thermal conductivity (kL) by solving Boltzmann Transport Equation (BTE) using ShengBTE. The second order harmonic and third-order anharmonic interatomic force constants required as input to ShengBTE was calculated using the Density-functional perturbation theory (DFPT). However, we predicted the electronic thermal conductivity (kel) using Wiedemann-Franz law as implemented in Boltztrap code. We also show that the sound velocity along 'a' and 'c' axes exhibit high anisotropy, which accounts for the anisotropic thermal conductivity of UB2.

Ab initio based study of the ArO- photoelectron spectra: Selectivity of spin-orbit transitions

NASA Astrophysics Data System (ADS)

Buchachenko, A. A.; Jakowski, Jacek; Chałasiński, Grzegorz; Szczȩśniak, M. M.; Cybulski, S. M.

2000-04-01

A combined ab initio atoms-in-molecule approach was implemented to model the photoelectron spectra of the ArO- anion. The lowest adiabatic states of Σ and Π symmetry of ArO and ArO- were investigated using the fourth-order Møller-Plessett perturbation theory including bond functions. The total energies were dissected into electrostatic, exchange, induction, and dispersion components. The complex of Ar with atomic oxygen is only weakly bound, primarily by dispersion interaction. The Π state possesses a deeper minimum (Re=3.4Å,De=380μEh) than the Σ state (Re=3.8Å,De=220μEh). In contrast, the anion complex is fairly strongly bound, primarily by ion-induced dipole induction forces, and the Σ state possesses a deeper minimum at shorter interatomic distances (Re=3.02Å,De=3600μEh) than the Π state (Re=3.35Å,De=2400μEh). The Σ-Π splittings in both systems are mainly due to differences in the exchange repulsion terms. Atoms-in-molecule models were used to account for the spin-orbit interaction, and to generate adiabatic relativistic potentials and wave functions. Collisional properties, diffusion, and mobility coefficients of O and O- in Ar, and absolute total Ar+O scattering cross sections, were calculated and found to agree well with the available experimental data. The photoelectron spectra were simulated within vibronic model, and were found in excellent agreement with the experimental measurements. The bimodal electron kinetic energy distribution was shown to stem from the strong selectivity of spin-orbit transitions, which split into two dense groups, depending on the initial electronic state of the anion. The latter feature cannot be described without explicit consideration of electronic intensity factor.

Students' use of atomic and molecular models in learning chemistry

NASA Astrophysics Data System (ADS)

O'Connor, Eileen Ann

1997-09-01

The objective of this study was to investigate the development of introductory college chemistry students' use of atomic and molecular models to explain physical and chemical phenomena. The study was conducted during the first semester of the course at a University and College II. Public institution (Carnegie Commission of Higher Education, 1973). Students' use of models was observed during one-on-one interviews conducted over the course of the semester. The approach to introductory chemistry emphasized models. Students were exposed to over two-hundred and fifty atomic and molecular models during lectures, were assigned text readings that used over a thousand models, and worked interactively with dozens of models on the computer. These models illustrated various features of the spatial organization of valence electrons and nuclei in atoms and molecules. Despite extensive exposure to models in lectures, in textbook, and in computer-based activities, the students in the study based their explanation in large part on a simple Bohr model (electrons arranged in concentric circles around the nuclei)--a model that had not been introduced in the course. Students used visual information from their models to construct their explanation, while overlooking inter-atomic and intra-molecular forces which are not represented explicitly in the models. In addition, students often explained phenomena by adding separate information about the topic without either integrating or logically relating this information into a cohesive explanation. The results of the study demonstrate that despite the extensive use of models in chemistry instruction, students do not necessarily apply them appropriately in explaining chemical and physical phenomena. The results of this study suggest that for the power of models as aids to learning to be more fully realized, chemistry professors must give more attention to the selection, use, integration, and limitations of models in their instruction.

Origin of Noncubic Scaling Law in Disordered Granular Packing

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xia, Chengjie; Li, Jindong; Kou, Binquan

Recent diffraction experiments on metallic glasses have unveiled an unexpected non-cubic scaling law between density and average interatomic distance, which lead to the speculations on the presence of fractal glass order. Using X-ray tomography we identify here a similar non-cubic scaling law in disordered granular packing of spherical particles. We find that the scaling law is directly related to the contact neighbors within first nearest neighbor shell, and therefore is closely connected to the phenomenon of jamming. The seemingly universal scaling exponent around 2.5 arises due to the isostatic condition with contact number around 6, and we argue that themore » exponent should not be universal.« less

Structural properties and glass transition in Aln clusters

NASA Astrophysics Data System (ADS)

Sun, D. Y.; Gong, X. G.

1998-02-01

We have studied the structural and dynamical properties of several Aln clusters by the molecular-dynamics method combined with simulated annealing. The well-fitted glue potential is used to describe the interatomic interaction. The obtained atomic structures for n=13, 55, and 147 are in agreement with results from ab initio calculations. Our results have demonstrated that the disordered cluster Al43 can be considered as a glass cluster. The obtained thermal properties of glass cluster Al43 are clearly different from the results for high-symmetry clusters, its melting behavior has properties similar to those of a glass solid. The present studies also show that the surface melting behavior does not exist in the studied Aln clusters.

The structural properties of PbF2 by molecular dynamics

NASA Astrophysics Data System (ADS)

Chergui, Y.; Nehaoua, N.; Telghemti, B.; Guemid, S.; Deraddji, N. E.; Belkhir, H.; Mekki, D. E.

2010-08-01

This work presents the use of molecular dynamics (MD) and the code of Dl_Poly, in order to study the structure of fluoride glass after melting and quenching. We are realized the processing phase liquid-phase, simulating rapid quenching at different speeds to see the effect of quenching rate on the operation of the devitrification. This technique of simulation has become a powerful tool for investigating the microscopic behaviour of matter as well as for calculating macroscopic observable quantities. As basic results, we calculated the interatomic distance, angles and statistics, which help us to know the geometric form and the structure of PbF2. These results are in experimental agreement to those reported in literature.

Correlation Functions and Glass Structure

NASA Astrophysics Data System (ADS)

Chergui, Y.; Nehaoua, N.; Telghemti, B.; Guemid, S.; Deraddji, N. E.; Belkhir, H.; Mekki, D. E.

2011-04-01

This work presents the use of molecular dynamics (MD) and the code of Dl Poly, in order to study the structure of fluoride glass after melting and quenching. We are realized the processing phase liquid-phase, simulating rapid quenching at different speeds to see the effect of quenching rate on the operation of the devitrification. This technique of simulation has become a powerful tool for investigating the microscopic behaviour of matter as well as for calculating macroscopic observable quantities. As basic results, we calculated the interatomic distance, angles and statistics, which help us to know the geometric form and the structure of PbF2. These results are in experimental agreement to those reported in literature.

Analytical W-He and H-He interatomic potentials for a W-H-He system

NASA Astrophysics Data System (ADS)

Li, Xiao-Chun; Shu, Xiaolin; Liu, Yi-Nan; Yu, Yi; Gao, F.; Lu, Guang-Hong

2012-07-01

We have constructed W-He and H-He analytical bond-order potentials for a W-H-He system. In combination with the previously self-developed W-H potential [X.-C. Li, X. Shu, Y.-N. Liu, F. Gao, G.-H. Lu, J. Nucl. Mater. 408 (2011) 12] and the Hartree-Fock-dispersion pair potential (Aziz-potential) for He-He interactions, we demonstrate that such potentials behave well for reproducing various properties of the W-H-He system such as defect formation energies, structural properties, and diffusion barriers. Such potentials can be employed to model both the He behaviours and the H-He synergetic effects in the W-H-He system.

Charge ordering in stoichiometric FeTe: Scanning tunneling microscopy and spectroscopy

DOE PAGES

Li, Wei; Yin, Wei -Guo; Wang, Lili; ...

2016-01-04

In this study, we use scanning tunneling microscopy and spectroscopy to reveal a unique stripy charge order in a parent phase of iron-based superconductors in stoichiometric FeTe epitaxy films. The charge order has unusually the same—usually half—period as the spin order. We also found highly anisotropic electron band dispersions being large and little along the ferromagnetic (crystallographic b) and antiferromagnetic (a) directions, respectively. Our data suggest that the microscopic mechanism is likely of the Stoner type driven by interatomic Coulomb repulsion V ij, and that V ij and charge fluctuations, so far much neglected, are important to the understanding ofmore » iron-based superconductors.« less

Thermodynamically accessible titanium clusters TiN, N = 2-32.

PubMed

Lazauskas, Tomas; Sokol, Alexey A; Buckeridge, John; Catlow, C Richard A; Escher, Susanne G E T; Farrow, Matthew R; Mora-Fonz, David; Blum, Volker W; Phaahla, Tshegofatso M; Chauke, Hasani R; Ngoepe, Phuti E; Woodley, Scott M

2018-05-10

We have performed a genetic algorithm search on the tight-binding interatomic potential energy surface (PES) for small TiN (N = 2-32) clusters. The low energy candidate clusters were further refined using density functional theory (DFT) calculations with the PBEsol exchange-correlation functional and evaluated with the PBEsol0 hybrid functional. The resulting clusters were analysed in terms of their structural features, growth mechanism and surface area. The results suggest a growth mechanism that is based on forming coordination centres by interpenetrating icosahedra, icositetrahedra and Frank-Kasper polyhedra. We identify centres of coordination, which act as centres of bulk nucleation in medium sized clusters and determine the morphological features of the cluster.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mishra, Tapan; Das, B. P.; Pai, Ramesh V.

We present a scenario where a supersolid is induced in one of the components of a mixture of two species bosonic atoms where there are no long-range interactions. We study a system of normal and hard-core boson mixture with only the former possessing long-range interactions. We consider three cases: the first where the total density is commensurate and the other two where it is incommensurate to the lattice. By suitable choices of the densities of normal and hard-core bosons and the interaction strengths between them, we predict that the charge density wave and the supersolid orders can be induced inmore » the hard-core species as a result of the competing interatomic interactions.« less

Atomic bonding effects in annular dark field scanning transmission electron microscopy. II. Experiments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Odlyzko, Michael L.; Held, Jacob T.; Mkhoyan, K. Andre, E-mail: mkhoyan@umn.edu

2016-07-15

Quantitatively calibrated annular dark field scanning transmission electron microscopy (ADF-STEM) imaging experiments were compared to frozen phonon multislice simulations adapted to include chemical bonding effects. Having carefully matched simulation parameters to experimental conditions, a depth-dependent bonding effect was observed for high-angle ADF-STEM imaging of aluminum nitride. This result is explained by computational predictions, systematically examined in the preceding portion of this study, showing the propagation of the converged STEM beam to be highly sensitive to net interatomic charge transfer. Thus, although uncertainties in experimental conditions and simulation accuracy remain, the computationally predicted experimental bonding effect withstands the experimental testing reportedmore » here.« less

Strength of the interatomic potential derived from angular scans in LEIS

NASA Astrophysics Data System (ADS)

Primetzhofer, D.; Markin, S. N.; Draxler, M.; Beikler, R.; Taglauer, E.; Bauer, P.

2008-09-01

Angular scans were performed for a Cu(1 0 0) single crystal and He + ions. The results were compared to MARLOWE, KALYPSO and FAN simulations to obtain information on the interaction potential. The influence of the used evaluation procedure on the deduced scattering potential was investigated. The scattering potential is found to be weaker than what is predicted by an uncorrected TFM potential. It was found that the use of a single screening correction factor is applicable in a wide range of impact parameters. It is further shown that selection of single scattering trajectories and a limitation of information depth to the surface layers is possible for neutral and charge integrated spectra.

Monte Carlo simulation of elongating metallic nanowires in the presence of surfactants

NASA Astrophysics Data System (ADS)

Gimenez, M. Cecilia; Reinaudi, Luis; Leiva, Ezequiel P. M.

2015-12-01

Nanowires of different metals undergoing elongation were studied by means of canonical Monte Carlo simulations and the embedded atom method representing the interatomic potentials. The presence of a surfactant medium was emulated by the introduction of an additional stabilization energy, represented by a parameter Q. Several values of the parameter Q and temperatures were analyzed. In general, it was observed for all studied metals that, as Q increases, there is a greater elongation before the nanowire breaks. In the case of silver, linear monatomic chains several atoms long formed at intermediate values of Q and low temperatures. Similar observations were made for the case of silver-gold alloys when the medium interacted selectively with Ag.

Mechanism for transient migration of xenon in UO{sub 2}

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, X.-Y.; Uberuaga, B. P.; Andersson, D. A.

2011-04-11

In this letter, we report recent work on atomistic modeling of diffusion migration events of the fission gas product xenon in UO{sub 2} nuclear fuel. Under nonequilibrium conditions, Xe atoms can occupy the octahedral interstitial site, in contrast to the thermodynamically most stable uranium substitutional site. A transient migration mechanism involving Xe and two oxygen atoms is identified using basin constrained molecular dynamics employing a Buckingham type interatomic potential. This mechanism is then validated using density functional theory calculations using the nudged elastic band method. An overall reduction in the migration barrier of 1.6-2.7 eV is obtained compared to vacancy-mediatedmore » diffusion on the uranium sublattice.« less

Assessment of empirical potential for MOX nuclear fuels and thermomechanical properties

NASA Astrophysics Data System (ADS)

Balboa, Hector; Van Brutzel, Laurent; Chartier, Alain; Le Bouar, Yann

2017-11-01

We assess five empirical interatomic potentials in the approximation of rigid ions and pair interactions for the (U1-y,Puy)O solid solution. The assessment compares available experimental data and Fink's recommendation with simulations on: the structural, thermodynamics, and mechanical properties over the full range of plutonium composition, from pure UO2 to pure PuO2 and for temperatures ranging from 300 K to the melting point. The best results are obtained by potentials referred as Cooper and Potashnikov potentials. The first one reproduces more accurately recommendations for the thermodynamics and mechanical properties exhibiting ductile-like behaviour during crack propagation, while the second one gives brittle behaviour at low temperature.

Quantum Nonlinear Optics without real Photons

NASA Astrophysics Data System (ADS)

Macrí, Vincenzo; Frisk Kockum, Anton; Stassi, Roberto; di Stefano, Omar; Savasta, Salvatore; Nori, Franco

We propose a physical process analogous to spontaneous parametric down-conversion, where one excited atom directly transfers its excitation to a couple of spatially-separated atoms with probability approaching one. The interaction is mediated by the exchange of virtual, rather than real, photons. This nonlinear optical process is coherent and reversible, so that the two excited atoms can transfer back the excitation to the first one: the atomic analogue of sum-frequency generation. The parameters used here correspond to experimentally-demonstrated values in circuit QED. This approach can be extended to consider other nonlinear interatomic processes, e.g. four-qubit mixing, and is an attractive architecture for the realization of quantum devices on a chip. Univ. of Michigan, USA.

Quantum Nonlinear Optics without Photons

NASA Astrophysics Data System (ADS)

Macrı, Vincenzo

Here we propose a physical process analogous to spontaneous parametric down-conversion, where one excited atom directly transfers its excitation to a couple of spatially separated atoms with probability approaching one. The interaction is mediated by the exchange of virtual rather than real photons. This nonlinear optical process is coherent and reversible, so that the couple of excited atoms can transfer back the excitation to the first one: the analogous for atoms of sum-frequency generation. The parameters used here correspond to experimentally-demonstrated values in circuit QED. This approach can be expanded to consider other nonlinear inter-atomic processes as the four-qubit mixing and is an attractive architecture for the realization of quantum devices on a chip.

Molecular dynamics modeling of helium bubbles in austenitic steels

NASA Astrophysics Data System (ADS)

Jelea, A.

2018-06-01

The austenitic steel devices from pressurized water reactors are continuously subjected to neutron irradiation that produces crystalline point defects and helium atoms in the steel matrix. These species evolve into large defects such as dislocation loops and helium filled bubbles. This paper analyzes, through molecular dynamics simulations with recently developed interatomic potentials, the impact of the helium/steel interface on the helium behavior in nanosize bubbles trapped in an austenitic steel matrix. It is shown that the repulsive helium-steel interactions induce higher pressures in the bubble compared to bulk helium at the same temperature and average density. A new equation of state for helium is proposed in order to take into account these interface effects.

Simulation of uniaxial deformation of hexagonal crystals (Mg, Be)

NASA Astrophysics Data System (ADS)

Vlasova, A. M.; Kesarev, A. G.

2017-12-01

Molecular dynamics (MD) simulations were performed for the nanocompression loading of nanocrystalline magnesium and beryllium modeled by an interatomic potential of the embedded atom method (EAM). It is shown that the main deformation modes are prismatic slip and twinning for magnesium, and only prismatic slip for beryllium. The formation of stable configurations of dislocation grids in magnesium and beryllium was observed. Dislocation networks are formed in the habit plane of the twin in a magnesium nanocrystall. Some dislocation reactions are suggested to explain the appearance of such networks. Shockley partial dislocations in a beryllium nanocrystall form grids in the slip plane. A strong anisotropy between slip systems was observed, which is in agreement with experimental data.

Ultrasonic studies of aluminium-substituted Bi(Pb)-2223 superconductors

NASA Astrophysics Data System (ADS)

Solunke, M. B.; Sharma, P. U.; Pandya, M. P.; Lakhani, V. K.; Modi, K. B.; Venugopal Reddy, P.; Shah, S. S.

2005-09-01

The compositional dependence of elastic properties of Al^{3+}-substitu- ted Bi(Pb)-2223 superconducting system with the general formula Bi_{1.7-x}Al_xPb_{0.3}Sr_2Ca_2- Cu_3O_y (x = 0.0, 0.1, 0.2 and 0.3) have been studied by means of ultrasonic pulse transmission (UPT) technique at 1 MHz (300 K). The elastic moduli of the specimens are computed and corrected to zero porosity. The observed variation of elastic constants with aluminium substitution has been explained on the basis of the strength of interatomic bonding. The applicability of heterogeneous metal mixture rule for estimating elastic constants and transition temperature has been tested.

Superfluid Boson-Fermion Mixture: Structure Formation and Collective Periodic Motion

NASA Astrophysics Data System (ADS)

Mitra, A.

2018-01-01

Multiple periodic domain formation due to a modulation instability in a boson-fermion mixture superfluid in the unitary regime has been studied. The periodicity of the structure evolves with time. At the early stage of evolution, bosonic domains show the periodic nature, whereas the periodicity in the fermionic (Cooper pair) domains appears at the late stage of evolution. The nature of interatomic interspecies interactions affects the domain formation. In a harmonic trap, the mixture executes an undamped oscillation. The frequency of the oscillation depends on the relative coupling strength between boson-fermion and fermion-fermion. The repulsive boson-fermion interaction reduces the oscillation frequency, whereas the attractive interaction enhances the frequency significantly.

Navigating at Will on the Water Phase Diagram

NASA Astrophysics Data System (ADS)

Pipolo, S.; Salanne, M.; Ferlat, G.; Klotz, S.; Saitta, A. M.; Pietrucci, F.

2017-12-01

Despite the simplicity of its molecular unit, water is a challenging system because of its uniquely rich polymorphism and predicted but yet unconfirmed features. Introducing a novel space of generalized coordinates that capture changes in the topology of the interatomic network, we are able to systematically track transitions among liquid, amorphous, and crystalline forms throughout the whole phase diagram of water, including the nucleation of crystals above and below the melting point. Our approach, based on molecular dynamics and enhanced sampling or free energy calculation techniques, is not specific to water and could be applied to very different structural phase transitions, paving the way towards the prediction of kinetic routes connecting polymorphic structures in a range of materials.

A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Shimojo, Fuyuki; Hattori, Shinnosuke; Department of Physics, Kumamoto University, Kumamoto 860-8555

We introduce an extension of the divide-and-conquer (DC) algorithmic paradigm called divide-conquer-recombine (DCR) to perform large quantum molecular dynamics (QMD) simulations on massively parallel supercomputers, in which interatomic forces are computed quantum mechanically in the framework of density functional theory (DFT). In DCR, the DC phase constructs globally informed, overlapping local-domain solutions, which in the recombine phase are synthesized into a global solution encompassing large spatiotemporal scales. For the DC phase, we design a lean divide-and-conquer (LDC) DFT algorithm, which significantly reduces the prefactor of the O(N) computational cost for N electrons by applying a density-adaptive boundary condition at themore » peripheries of the DC domains. Our globally scalable and locally efficient solver is based on a hybrid real-reciprocal space approach that combines: (1) a highly scalable real-space multigrid to represent the global charge density; and (2) a numerically efficient plane-wave basis for local electronic wave functions and charge density within each domain. Hybrid space-band decomposition is used to implement the LDC-DFT algorithm on parallel computers. A benchmark test on an IBM Blue Gene/Q computer exhibits an isogranular parallel efficiency of 0.984 on 786 432 cores for a 50.3 × 10{sup 6}-atom SiC system. As a test of production runs, LDC-DFT-based QMD simulation involving 16 661 atoms is performed on the Blue Gene/Q to study on-demand production of hydrogen gas from water using LiAl alloy particles. As an example of the recombine phase, LDC-DFT electronic structures are used as a basis set to describe global photoexcitation dynamics with nonadiabatic QMD (NAQMD) and kinetic Monte Carlo (KMC) methods. The NAQMD simulations are based on the linear response time-dependent density functional theory to describe electronic excited states and a surface-hopping approach to describe transitions between the excited states. A series of techniques are employed for efficiently calculating the long-range exact exchange correction and excited-state forces. The NAQMD trajectories are analyzed to extract the rates of various excitonic processes, which are then used in KMC simulation to study the dynamics of the global exciton flow network. This has allowed the study of large-scale photoexcitation dynamics in 6400-atom amorphous molecular solid, reaching the experimental time scales.« less

State-of-the-art ab initio potential energy curve for the krypton atom pair and thermophysical properties of dilute krypton gas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jäger, Benjamin, E-mail: benjamin.jaeger@uni-rostock.de; Hellmann, Robert, E-mail: robert.hellmann@uni-rostock.de; Bich, Eckard

2016-03-21

A new reference krypton-krypton interatomic potential energy curve was developed by means of quantum-chemical ab initio calculations for 36 interatomic separations. Highly accurate values for the interaction energies at the complete basis set limit were obtained using the coupled-cluster method with single, double, and perturbative triple excitations as well as t-aug-cc-pV5Z and t-aug-cc-pV6Z basis sets including mid-bond functions, with the 6Z basis set being newly constructed for this study. Higher orders of coupled-cluster terms were considered in a successive scheme up to full quadruple excitations. Core-core and core-valence correlation effects were included. Furthermore, relativistic effects were studied not only atmore » a scalar relativistic level using second-order direct perturbation theory, but also utilizing full four-component and Gaunt-effect computations. An analytical pair potential function was fitted to the interaction energies, which is characterized by a depth of 200.88 K with an estimated standard uncertainty of 0.51 K. Thermophysical properties of low-density krypton were calculated for temperatures up to 5000 K. Second and third virial coefficients were obtained from statistical thermodynamics. Viscosity and thermal conductivity as well as the self-diffusion coefficient were computed using the kinetic theory of gases. The theoretical results are compared with experimental data and with results for other pair potential functions from the literature, especially with those calculated from the recently developed ab initio potential of Waldrop et al. [J. Chem. Phys. 142, 204307 (2015)]. Highly accurate experimental viscosity data indicate that both the present ab initio pair potential and the one of Waldrop et al. can be regarded as reference potentials, even though the quantum-chemical methods and basis sets differ. However, the uncertainties of the present potential and of the derived properties are estimated to be considerably lower.« less

Transverse discrete breathers in unstrained graphene

NASA Astrophysics Data System (ADS)

Barani, Elham; Lobzenko, Ivan P.; Korznikova, Elena A.; Soboleva, Elvira G.; Dmitriev, Sergey V.; Zhou, Kun; Marjaneh, Aliakbar Moradi

2017-02-01

Discrete breathers (DB) are spatially localized vibrational modes of large amplitude in defect-free nonlinear lattices. The search for DBs in graphene is of high importance, taking into account that this one atom thick layer of carbon is promising for a number of applications. There exist several reports on successful excitation of DBs in graphene, based on molecular dynamics and ab initio simulations. In a recent work by Hizhnyakov with co-authors the possibility to excite a DB with atoms oscillating normal to the graphene sheet has been reported. In the present study we use a systematic approach for finding initial conditions to excite transverse DBs in graphene. The approach is based on the analysis of the frequency-amplitude dependence for a delocalized, short-wavelength vibrational mode. This mode is a symmetry-dictated exact solution to the dynamic equations of the atomic motion, regardless the mode amplitude and regardless the type of interatomic potentials used in the simulations. It is demonstrated that if the AIREBO potential is used, the mode frequency increases with the amplitude bifurcating from the upper edge of the phonon spectrum for out-of-plane phonons. Then a bell-shaped function is superimposed on this delocalized mode to obtain a spatially localized vibrational mode, i.e., a DB. Placing the center of the bell-shaped function at different positions with respect to the lattice sites, three different DBs are found. Typically, the degree of spatial localization of DBs increases with the DB amplitude, but the transverse DBs in graphene reported here demonstrate the opposite trend. The results are compared to those obtained with the use of the Savin interatomic potential and no transverse DBs are found in this case. The results of this study contribute to a better understanding of the nonlinear dynamics of graphene and they call for the ab initio simulations to verify which of the two potentials used in this study is more precise.

State-of-the-art ab initio potential energy curve for the krypton atom pair and thermophysical properties of dilute krypton gas.

PubMed

Jäger, Benjamin; Hellmann, Robert; Bich, Eckard; Vogel, Eckhard

2016-03-21

A new reference krypton-krypton interatomic potential energy curve was developed by means of quantum-chemical ab initio calculations for 36 interatomic separations. Highly accurate values for the interaction energies at the complete basis set limit were obtained using the coupled-cluster method with single, double, and perturbative triple excitations as well as t-aug-cc-pV5Z and t-aug-cc-pV6Z basis sets including mid-bond functions, with the 6Z basis set being newly constructed for this study. Higher orders of coupled-cluster terms were considered in a successive scheme up to full quadruple excitations. Core-core and core-valence correlation effects were included. Furthermore, relativistic effects were studied not only at a scalar relativistic level using second-order direct perturbation theory, but also utilizing full four-component and Gaunt-effect computations. An analytical pair potential function was fitted to the interaction energies, which is characterized by a depth of 200.88 K with an estimated standard uncertainty of 0.51 K. Thermophysical properties of low-density krypton were calculated for temperatures up to 5000 K. Second and third virial coefficients were obtained from statistical thermodynamics. Viscosity and thermal conductivity as well as the self-diffusion coefficient were computed using the kinetic theory of gases. The theoretical results are compared with experimental data and with results for other pair potential functions from the literature, especially with those calculated from the recently developed ab initio potential of Waldrop et al. [J. Chem. Phys. 142, 204307 (2015)]. Highly accurate experimental viscosity data indicate that both the present ab initio pair potential and the one of Waldrop et al. can be regarded as reference potentials, even though the quantum-chemical methods and basis sets differ. However, the uncertainties of the present potential and of the derived properties are estimated to be considerably lower.

Electronic interaction anisotropy between open-shell lanthanide atoms and helium from cold collision experiment

NASA Astrophysics Data System (ADS)

Krems, R. V.; Buchachenko, A. A.

2005-09-01

Based on measurements of the Zeeman relaxation in a cold gas of He3 [C. I. Hancox, S. C. Doret, M. I. Hummon, L. Luo, and J. M. Doyle, Nature (London) 431, 281 (2004)], we show that the electronic interaction anisotropy between rare-earth atoms with nonzero electronic orbital angular momenta and helium is extremely small. The interaction of the rare-earth atoms with He gives rise to several adiabatic potentials with different electronic symmetries. It is demonstrated that the energy splitting between these potentials does not exceed 0.09cm-1 at interatomic distances larger than the turning point for collisions at 0.8K, including the region of the van der Waals interaction minima.

Theoretical study of the composition pulling effect in InGaN metalorganic vapor-phase epitaxy growth

NASA Astrophysics Data System (ADS)

Inatomi, Yuya; Kangawa, Yoshihiro; Ito, Tomonori; Suski, Tadeusz; Kumagai, Yoshinao; Kakimoto, Koichi; Koukitu, Akinori

2017-07-01

The composition pulling effect in metalorganic vapor-phase InGaN epitaxy was theoretically investigated by thermodynamic analysis. The excess energies of biaxial-strained In x Ga1- x N were numerically calculated using empirical interatomic potentials considering different situations: (i) coherent growth on GaN(0001), (ii) coherent growth on In0.2Ga0.8N(0001), and (iii) bulk growth. Using the excess energies, the excess chemical potentials of InN and GaN alloys were computed. Our results show that compressive strain suppresses In incorporation, whereas tensile strain promotes it. Moreover, assuming chemical equilibrium, the relationship between the solid composition and the growth conditions was predicted. The results successfully reproduced the typical composition pulling effect.

Probing quantum effects in lithium

NASA Astrophysics Data System (ADS)

Deemyad, Shanti; Zhang, Rong

2018-05-01

In periodic table lithium is the first element immediately after helium and the lightest metal. While fascinating quantum nature of condensed helium is suppressed at high densities, lithium is expected to adapt more quantum solid behavior under compression. This is due to the presence of long range interactions in metallic systems for which an increase in the de-Boer parameter (λ/σ, where σ is the minimum interatomic distance and λ is the de-Broglie wavelength) is predicted at higher densities [1,2]. Physics of dense lithium offers a rich playground to look for new emergent quantum phenomena in condensed matter and has been subject of many theoretical and experimental investigations. In this article recent progress in studying the quantum nature of dense lithium will be discussed.

The spectrum of density fluctuations of noble gases probed by THz neutron and x-ray spectroscopy

DOE PAGES

Cunsolo, Alessandro

2016-02-26

Approximately 50 years of inelastic scattering studies of noble gases are reviewed to illustrate the main advances achieved in the understanding of the THz dynamics of simple systems. The gradual departure of the spectral shape from the hydrodynamic regime is discussed with an emphasis on the phenomenology of fast (sub- ps) relaxation processes. This review shows that relaxation phenomena in noble gases have an essentially collisional origin, which is also revealed by the parallelism between their characteristic timescale and the interatomic collision time. In addition, recent THz spectroscopy results on noble gases at extreme thermodynamic conditions are discussed to illustratemore » the need for a revision of our current understanding of the supercritical phase.« less

Size and habit evolution of PETN crystals - a lattice Monte Carlo study

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zepeda-Ruiz, L A; Maiti, A; Gee, R

2006-02-28

Starting from an accurate inter-atomic potential we develop a simple scheme of generating an ''on-lattice'' molecular potential of short range, which is then incorporated into a lattice Monte Carlo code for simulating size and shape evolution of nanocrystallites. As a specific example, we test such a procedure on the morphological evolution of a molecular crystal of interest to us, e.g., Pentaerythritol Tetranitrate, or PETN, and obtain realistic facetted structures in excellent agreement with experimental morphologies. We investigate several interesting effects including, the evolution of the initial shape of a ''seed'' to an equilibrium configuration, and the variation of growth morphologymore » as a function of the rate of particle addition relative to diffusion.« less

The structure of K3C60 and the mechanism of superconductivity.

PubMed Central

Pauling, L

1991-01-01

Analysis of the interatomic distances in the superconducting substance K3C60 indicates that each of the K atoms in tetrahedral interstices between C60 spheres accepts three electrons from C60, thus becoming quadricovalent; its four bonds resonate among the 24 adjacent carbon atoms to give a strong framework in which the negative charges are localized on these K atoms. The electric current is carried by the motion of positive charges (holes) through the network of C60 spheres and the K atoms in octahedral holes. Superconductivity is favored by the localization of the negative charges on the tetrahedral K atoms and their noninvolvement in valence-bond resonance, decreasing the rate of mutual extinction of electrons and holes. PMID:11607222

Crystal structure and chemical bonding in the mixed anion compound BaSF.

PubMed

Driss, D; Cadars, S; Deniard, P; Mevellec, J-Y; Corraze, B; Janod, E; Cario, L

2017-11-28

BaSF was synthesised by a solid state reaction at high temperature and its crystal structure was determined thanks to X-ray diffraction on a single crystal. This transparent yellow fluorochalcogenide has an intergrowth structure built from the stacking of fluorite type layers and sulfur layers. In BaSF sulfur atoms form dimers with interatomic distances as short as 2.1074(10) Å. DFT calculations confirm that this compound is a band insulator with the Fermi level lying in between the antibonding π* and σ* molecular orbitals of the sulfur dimers. Reflectance measurements show that the optical band gap of BaSF is about 2.7 eV in good agreement with the value found from DFT calculations.

Inelastic X-ray Scattering Studies of Plasmons in Carbon Nanotubes

NASA Astrophysics Data System (ADS)

Upton, M. H.; Klie, R. F.; Hill, J. P.; Gog, T.; Casa, D.; Ku, W.; Zhu, Y.; Sfeir, M. Y.; Misewich, J.; Eres, G.; Lowndes, D.

2007-03-01

We investigate the physical parameters controlling the low energy screening in carbon nanotubes via electron energy loss spectroscopy and inelastic x-ray scattering. Two plasmon-like features are observed, one near 9 eV (the so- called π plasmon) and one near 20 eV (the so-called π+σ plasmon). At large nanotube diameters, the π+σ plasmon energies depend exclusively on the number of walls and not on the radius or chiral vector. This shift indicates a change of strength of screening and the effective interaction at inter-atomic distance, and thus suggests an alternative mechanism of tuning the properties of the nanotube in addition to the well-known control provided by chirality and tube diameter.

Saturation and negative temperature coefficient of electrical resistivity in liquid iron-sulfur alloys at high densities from first-principles calculations

NASA Astrophysics Data System (ADS)

Wagle, Fabian; Steinle-Neumann, Gerd; de Koker, Nico

2018-03-01

We report results on electronic transport properties of liquid Fe-S alloys at conditions of planetary cores, computed using first-principle techniques in the Kubo-Greenwood formalism. We describe a combined effect of resistivity saturation due to temperature, compression, and chemistry by comparing the electron mean free path from the Drude response of optical conductivity to the mean interatomic distance. At high compression and high sulfur concentration the Ioffe-Regel condition is satisfied, and the temperature coefficient of resistivity changes sign from positive to negative. We show that this happens due to a decrease in the d density of states at the Fermi level in response to thermal broadening.

Dislocation nucleation from symmetric tilt grain boundaries in body-centered cubic vanadium

NASA Astrophysics Data System (ADS)

Xu, Shuozhi; Su, Yanqing

2018-05-01

We perform molecular dynamics (MD) simulations with two interatomic potentials to study dislocation nucleation from six symmetric tilt grain boundaries (GB) using bicrystal models in body-centered cubic vanadium. The influences of the misorientation angle are explored in the context of activated slip systems, critical resolved shear stress (CRSS), and GB energy. It is found that for four GBs, the activated slip systems are not those with the highest Schmid factor, i.e., the Schmid law breaks down. For all misorientation angles, the bicrystal is associated with a lower CRSS than their single crystalline counterparts. Moreover, the GB energy decreases in compressive loading at the yield point with respect to the undeformed configuration, in contrast to tensile loading.

Analytical bond order potential for simulations of BeO 1D and 2D nanostructures and plasma-surface interactions

NASA Astrophysics Data System (ADS)

Byggmästar, J.; Hodille, E. A.; Ferro, Y.; Nordlund, K.

2018-04-01

An analytical interatomic bond order potential for the Be-O system is presented. The potential is fitted and compared to a large database of bulk BeO and point defect properties obtained using density functional theory. Its main applications include simulations of plasma-surface interactions involving oxygen or oxide layers on beryllium, as well as simulations of BeO nanotubes and nanosheets. We apply the potential in a study of oxygen irradiation of Be surfaces, and observe the early stages of an oxide layer forming on the Be surface. Predicted thermal and elastic properties of BeO nanotubes and nanosheets are simulated and compared with published ab initio data.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Grochola, A.; Kowalczyk, P.; Szczepkowski, J.

Comprehensive spectroscopic studies of hot and ultracold samples of NaCs molecules were combined to complete the investigation of the (3){Omega}=1 <- X {sup 1}{Sigma}{sup +} transition for the NaCs molecule. Polarization labeling, photoassociation, and pulsed laser depletion spectroscopy were used to collect data on rovibrational levels of the (3){Omega}=1 state [here described as the c {sup 3}{Sigma}{sup +} state in Hund's case (a) notation]. The highest observed level was v=72 located {approx}5 GHz below the atomic asymptote Na(3 {sup 2}S{sub 1/2}) + Cs(6 {sup 2}P{sub 3/2}). Approximately 1400 levels were used to construct the potential energy curve of the (3){Omega}=1more » state for the full range of interatomic distances.« less

A molecular dynamics study of melting and dissociation of tungsten nanoparticles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Min; Wang, Jun; Fu, Baoqin

2015-12-15

Molecular dynamics simulations were conducted to study the melting and dissociation of free tungsten nanoparticles. For the various interatomic potentials applied, the melting points of the tungsten nanoparticles increased with increasing nanoparticle diameter. Combining these results with the melting point of bulk tungsten in the experiment, the melting point of nanoparticles with diameters ranging from 4 to 12 nm could be determined. As the temperature increases, free nanoparticles are subject to dissociation phenomena. The dissociation rate was observed to follow Arrhenius behavior, and the Meyer–Neldel rule was obeyed. These results are useful in understanding the behavior of tungsten dust generatedmore » in nuclear fusion devices as well as for the preparation, formation, and application of tungsten powders.« less

Attainment of 40.5 pm spatial resolution using 300 kV scanning transmission electron microscope equipped with fifth-order aberration corrector.

PubMed

Morishita, Shigeyuki; Ishikawa, Ryo; Kohno, Yuji; Sawada, Hidetaka; Shibata, Naoya; Ikuhara, Yuichi

2018-02-01

The achievement of a fine electron probe for high-resolution imaging in scanning transmission electron microscopy requires technological developments, especially in electron optics. For this purpose, we developed a microscope with a fifth-order aberration corrector that operates at 300 kV. The contrast flat region in an experimental Ronchigram, which indicates the aberration-free angle, was expanded to 70 mrad. By using a probe with convergence angle of 40 mrad in the scanning transmission electron microscope at 300 kV, we attained the spatial resolution of 40.5 pm, which is the projected interatomic distance between Ga-Ga atomic columns of GaN observed along [212] direction.

MOLECULAR DYNAMICS OF CASCADES OVERLAP IN TUNGSTEN WITH 20-KEV PRIMARY KNOCK-ON ATOMS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Setyawan, Wahyu; Nandipati, Giridhar; Roche, Kenneth J.

2015-04-16

Molecular dynamics simulations are performed to investigate the mutual influence of two subsequent cascades in tungsten. The influence is studied using 20-keV primary knock-on atoms, to induce one cascade after another separated by 15 ps, in a lattice temperature of 1025 K (i.e. 0.25 of the melting temperature of the interatomic potential). The center of mass of the vacancies at the peak damage during the cascade is taken as the location of the cascade. The distance between this location to that of the next cascade is taken as the overlap parameter. Empirical fits describing the number of surviving vacancies andmore » interstitial atoms as a function of overlap are presented.« less

Kinetics versus thermodynamics in materials modeling: The case of the di-vacancy in iron

NASA Astrophysics Data System (ADS)

Djurabekova, F.; Malerba, L.; Pasianot, R. C.; Olsson, P.; Nordlund, K.

2010-07-01

Monte Carlo models are widely used for the study of microstructural and microchemical evolution of materials under irradiation. However, they often link explicitly the relevant activation energies to the energy difference between local equilibrium states. We provide a simple example (di-vacancy migration in iron) in which a rigorous activation energy calculation, by means of both empirical interatomic potentials and density functional theory methods, clearly shows that such a link is not granted, revealing a migration mechanism that a thermodynamics-linked activation energy model cannot predict. Such a mechanism is, however, fully consistent with thermodynamics. This example emphasizes the importance of basing Monte Carlo methods on models where the activation energies are rigorously calculated, rather than deduced from widespread heuristic equations.

Supramolecular assembly in the epiisopiloturine hydrochloride salt

NASA Astrophysics Data System (ADS)

Mafud, Ana Carolina; Reinheimer, Eric W.; Lima, Filipe Camargo Dalmatti Alves; Batista, Larissa Fernandes; de Paula, Karina; Véras, Leiz Maria Costa; de Souza de Almeida Leite, José Roberto; Venancio, Tiago; Mascarenhas, Yvonne Primerano

2017-05-01

Epiisopiloturine hydrochloride (Epi-HCl) salt was synthetized from epiisopiloturine, an in vivo anthelmintic compound against Schistosoma mansoni worms. Despite there being no acute toxicity in mammalian cells, the compound's water insolubility makes its administration difficult. In this communication, we report the characterization of Epi-HCl its features by spectroscopy, thermal analysis, and PXRD. The single crystals suitable to X-ray diffraction were grown by slow evaporation technique. To better understand the nature of Epi-HCl' solid state, SS-NMR was also used. The salt's intramolecular structure was maintained via cation-pi intramolecular interactions, which in conjunction with hydrogen bonding, gives rise to an extended supramolecular assembly. The interatomic distances within the cations and environment around the chloride anion vary as function of temperature, suggesting a packing relaxation.

Homogenization Models for Carbon Nanotubes

NASA Astrophysics Data System (ADS)

Muc, A.; Jamróz, M.

2004-03-01

Two homogenization models for evaluating Young's modulus of nanocomposites reinforced with single-walled and multi-walled carbon nanotubes are presented. The first model is based on a physical description taking into account the interatomic interaction and nanotube geometry. The elementary cell, here a nanotube with a surrounding resin layer, is treated as a homogeneous body — a material continuum. The second model, similar to a phenomenological engineering one, is obtained by combining the law of mixture with the Cox mechanical model. This model describes the stress distribution along stretched short fibers surrounded by a resin matrix. The similarities between composite materials reinforced with short fibers and nanotubes are elucidated. The results obtained are compared with those for classical microcomposites to demonstrate the advantages and disadvantages of both the composite materials.

Crystal structure of simple metals at high pressures

NASA Astrophysics Data System (ADS)

Degtyareva, Olga

2010-09-01

The effects of pressure on the crystal structure of simple (or sp-) elements are analysed in terms of changes in coordination number, packing density, and interatomic distances, and general rules are established. In the polyvalent elements from groups 14-17, the covalently bonded structures tend to transform to metallic phases with a gradual increase in coordination number and packing density, a behaviour normally expected under pressure. Group 1 and 2 metallic elements, however, show a reverse trend towards structures with low packing density due to intricate changes in their electronic structure. Complex crystal structures such as host-guest and incommensurately modulated structures found in these elements are given special attention in this review in an attempt to determine their role in the observed phase-transition sequences.

Monte Carlo simulation of elongating metallic nanowires in the presence of surfactants

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gimenez, M. Cecilia; Reinaudi, Luis, E-mail: luis.reinaudi@unc.edu.ar; Leiva, Ezequiel P. M.

2015-12-28

Nanowires of different metals undergoing elongation were studied by means of canonical Monte Carlo simulations and the embedded atom method representing the interatomic potentials. The presence of a surfactant medium was emulated by the introduction of an additional stabilization energy, represented by a parameter Q. Several values of the parameter Q and temperatures were analyzed. In general, it was observed for all studied metals that, as Q increases, there is a greater elongation before the nanowire breaks. In the case of silver, linear monatomic chains several atoms long formed at intermediate values of Q and low temperatures. Similar observations weremore » made for the case of silver-gold alloys when the medium interacted selectively with Ag.« less

Phonon interference control of atomic-scale metamirrors, meta-absorbers, and heat transfer through crystal interfaces

NASA Astrophysics Data System (ADS)

Kosevich, Yu. A.; Potyomina, L. G.; Darinskii, A. N.; Strelnikov, I. A.

2018-03-01

The paper theoretically studies the possibility of using the effects of phonon interference between paths through different interatomic bonds for the control of phonon heat transfer through internal crystal interfaces and for the design of phonon metamirrors and meta-absorbers. These metamirrors and meta-absorbers are considered to be defect nanolayers of atomic-scale thicknesses embedded in a crystal. Several analytically solvable three-dimensional lattice-dynamics models of the phonon metamirrors and meta-absorbers at the internal crystal planes are described. It is shown that due to destructive interference in the two or more phonon paths, the internal crystal planes, fully or partially filled with weakly bound or heavy-isotope defect atoms, can completely reflect or completely absorb phonons at the transmission antiresonances, whose wavelengths are larger than the effective thickness of the metamirror or meta-absorber. Due to cooperative superradiant effect, the spectral widths of the two-path interference antiresonances for the plane waves are given by the square of partial filling fraction in the defect crystal plane. Our analysis reveals that the presence of two or more phonon paths plays the dominant role in the emergence of the transmission antiresonances in phonon scattering at the defect crystal planes and in reduction of the thermal interface conductance in comparison with the Fano-resonance concept. We study analytically phonon transmission through internal crystal plane in a model cubic lattice of Si-like atoms, partially filled with Ge-like defect atoms. Such a plane can serve as interference phonon metamirror with the transmission antiresonances in the vicinities of eigenmode frequencies of Ge-like defect atoms in the terahertz frequency range. We predict the extraordinary phonon transmission induced by the two-path constructive interference of the lattice waves in resonance with the vibrations of rare host atoms, periodically distributed in the crystal plane almost completely filled with heavy-isotope defects. We show that the phonon-interference-induced transparency can be produced by the defect nanolayer with the non-nearest-neighbor interactions, filled with two types of isotopes with relatively small difference in masses or binding force constants. In this case, relatively broad transmission antiresonance is accompanied by the narrow transmission peak close to the antiresonance frequency. We describe the softening of the flexural surface acoustic wave, localized at the embedded defect nanolayer, caused by negative surface stress in the layer. The surface wave softening results in spatially periodic static bending deformation of the embedded nanolayer with the definite wave number. The latter effect is estimated for graphene monolayer embedded in a strained matrix of polyethylene. We analyze the effect of nonlinearity in the dynamics of defect atoms on the one- and two-path phonon interference and show that the interference transmission resonances and antiresonances are shifted in frequencies but not completely suppressed by rather strong anharmonicity of interatomic bonds. The reduction of the Kapitza thermal interface conductance caused by the destructive phonon interference in a defect monolayer is described. We show that the additional relatively weak non-nearest-neighbor interactions through the defect crystal plane filled with heavy isotopes substantially reduces the interface thermal conductance, and this effect is stronger in the three-dimensional system than in the quasi-one-dimensional systems studied previously.

Electronic properties of liquid Hg-In alloys : Ab-initio molecular dynamics study

NASA Astrophysics Data System (ADS)

Sharma, Nalini; Thakur, Anil; Ahluwalia, P. K.

2016-05-01

Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-In alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Three liquid Hg-In alloys (Hg10In90, Hg30In70,. Hg50In50, Hg70In30, and Hg90Pb10) at 299 K are considered. The calculated results for liquid Hg (l-Hg) and lead (l-In) are also drawn. Along with the calculated results of considered five liquid alloys of Hg-In alloy. The results obtained from electronic properties namely total density of state and partial density of states help to find the local arrangement of Hg and In atoms and the presence of liquid state in the considered five alloys.

Study of atomic structure of liquid Hg-In alloys using ab-initio molecular dynamics

NASA Astrophysics Data System (ADS)

Sharma, Nalini; Thakur, Anil; Ahluwalia, P. K.

2015-05-01

Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-In alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Five liquid Hg-In mixtures (Hg10In90, Hg30In70, Hg50In50, Hg70In30 and Hg90In10) at 299K are considered. The radial distribution function g(r) and structure factor S(q) of considered alloys are compared with respective experimental results for liquid Hg (l-Hg) and (l-In). The radial distribution function g(r) shows the presence of short range order in the systems considered. Smooth curves of Bhatia-Thornton partial structure factors factor shows the presence of liquid state in the considered alloys.

Equivalent-Continuum Modeling of Nano-Structured Materials

NASA Technical Reports Server (NTRS)

Odegard, Gregory M.; Gates, Thomas S.; Nicholson, Lee M.; Wise, Kristopher E.

2001-01-01

A method has been developed for modeling structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with an equivalent-continuum model. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As an important example with direct application to the development and characterization of single-walled carbon nanotubes, the model has been applied to determine the effective continuum geometry of a graphene sheet. A representative volume element of the equivalent-continuum model has been developed with an effective thickness. This effective thickness has been shown to be similar to, but slightly smaller than, the interatomic spacing of graphite.

Direct photoassociation of halo molecules in ultracold 86 Sr

NASA Astrophysics Data System (ADS)

Aman, J. A.; Hill, Joshua; Killian, T. C.

2017-04-01

We investigate the creation of 1S0 +1S0 halo molecules in strontium 86 through direct photoassociation in an optical dipole trap. We drive two photon Raman transitions near-resonance with a molecular level of the 1S0 +3P1 interatomic potential as the intermediate state. This provides large Frank-Condon factors and allows us to observe resonances for the creation of halo molecules through higher order Raman processes. The halo molecule is bound by EB 85 kHz at low excitation-laser intensity, but experiments show large AC Stark shifts of the molecular binding energy. These conditions suggest that STIRAP should be very effective for improving molecular conversion efficiency. Further experiments in a 3D lattice will explore molecular lifetimes and collision rates. Travel support provided by Shell Corporation.

A study of microkinetic adjustments required to match shock wave experiments and Monte Carlo Direct Simulation for a wide Mach number range

NASA Technical Reports Server (NTRS)

Pham-Van-diep, Gerald C.; Muntz, E. Phillip; Erwin, Daniel A.

1990-01-01

Shock wave thickness predictions from Monte Carlo Direct Simulations, using differential scattering and the Maitland-Smith-Aziz interatomic potential, underpredict experiments as shock Mach numbers increase above about 4. Examination of several sources of data has indicated that at relatively high energies the repulsive portion of accepted potentials such as the Maitland-Smith-Aziz may be too steep. An Exponential-6 potential due to Ross, based on high energy molecular beam scattering data and shock velocity measurements in liquid argon, has been combined with the lower energy portion of the Maitland-Smith-Aziz potential. When this hybrid potential is used in Monte Carlo Direct Simulations, agreement with experiments is improved over the previous predictions using the pure Maitland-Smith-Aziz form.

Manipulating matter rogue waves and breathers in Bose-Einstein condensates.

PubMed

Manikandan, K; Muruganandam, P; Senthilvelan, M; Lakshmanan, M

2014-12-01

We construct higher-order rogue wave solutions and breather profiles for the quasi-one-dimensional Gross-Pitaevskii equation with a time-dependent interatomic interaction and external trap through the similarity transformation technique. We consider three different forms of traps: (i) the time-independent expulsive trap, (ii) time-dependent monotonous trap, and (iii) time-dependent periodic trap. Our results show that when we change a parameter appearing in the time-independent or time-dependent trap the second- and third-order rogue waves transform into the first-order-like rogue waves. We also analyze the density profiles of breather solutions. Here we also show that the shapes of the breathers change when we tune the strength of the trap parameter. Our results may help to manage rogue waves experimentally in a BEC system.

Geometrical eigen-subspace framework based molecular conformation representation for efficient structure recognition and comparison

NASA Astrophysics Data System (ADS)

Li, Xiao-Tian; Yang, Xiao-Bao; Zhao, Yu-Jun

2017-04-01

We have developed an extended distance matrix approach to study the molecular geometric configuration through spectral decomposition. It is shown that the positions of all atoms in the eigen-space can be specified precisely by their eigen-coordinates, while the refined atomic eigen-subspace projection array adopted in our approach is demonstrated to be a competent invariant in structure comparison. Furthermore, a visual eigen-subspace projection function (EPF) is derived to characterize the surrounding configuration of an atom naturally. A complete set of atomic EPFs constitute an intrinsic representation of molecular conformation, based on which the interatomic EPF distance and intermolecular EPF distance can be reasonably defined. Exemplified with a few cases, the intermolecular EPF distance shows exceptional rationality and efficiency in structure recognition and comparison.

Atomic interaction of the MEAM type for the study of intermetallics in the Al-U alloy

NASA Astrophysics Data System (ADS)

Pascuet, M. I.; Fernández, J. R.

2015-12-01

Interaction for both pure Al and Al-U alloys of the MEAM type are developed. The obtained Al interatomic potential assures its compatibility with the details of the framework presently adopted. The Al-U interaction fits various properties of the Al2U, Al3U and Al4U intermetallics. The potential verifies the stability of the intermetallic structures in a temperature range compatible with that observed in the phase diagram, and also takes into account the greater stability of these structures relative to others that are competitive in energy. The intermetallics are characterized by calculating elastic and thermal properties and point defect parameters. Molecular dynamics simulations show a growth of the Al3U intermetallic in the Al/U interface in agreement with experimental evidence.

Simulation of defects in fusion plasma first wall materials

NASA Astrophysics Data System (ADS)

T, Troev; N, Nankov; T, Yoshiie

2014-06-01

Numerical calculations of radiation damages in beryllium, alpha-iron and tungsten irradiated by fusion neutrons were performed using molecular dynamics (MD) simulations. The displacement cascades efficiency has been calculated using the Norgett-Robinson-Torrens (NRT) formula, the universal pair-potential of Ziegler-Biersack-Littmark (ZBL) and the EAM inter-atomic potential. The pair potential overestimates the defects production by a factor of 2. The ZBL pair potential results and the EAM are comparable at higher primary knock-on atom (PKA) energies (E > 100 keV). We found that the most common types of defects are single vacancies, di-vacancies, interstitials and small number of interstitial clusters. On the bases of calculated results, the behavior of vacancies, empty nano-voids and nano-voids with hydrogen and helium were discussed.

The self-consistent calculation of pseudo-molecule energy levels, construction of energy level correlation diagrams and an automated computation system for SCF-X(Alpha)-SW calculations

NASA Technical Reports Server (NTRS)

Schlosser, H.

1981-01-01

The self consistent calculation of the electronic energy levels of noble gas pseudomolecules formed when a metal surface is bombarded by noble gas ions is discussed along with the construction of energy level correlation diagrams as a function of interatomic spacing. The self consistent field x alpha scattered wave (SCF-Xalpha-SW) method is utilized. Preliminary results on the Ne-Mg system are given. An interactive x alpha programming system, implemented on the LeRC IBM 370 computer, is described in detail. This automated system makes use of special PROCDEFS (procedure definitions) to minimize the data to be entered manually at a remote terminal. Listings of the special PROCDEFS and of typical input data are given.

Dynamic stabilities of icosahedral-like clusters and their ability to form quasicrystals

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liang, Xiaogang; Hamid, Ilyar; Duan, Haiming, E-mail: dhm@xju.edu.cn

2016-06-15

The dynamic stabilities of the icosahedral-like clusters containing up to 2200 atoms are investigated for 15 metal elements. The clusters originate from five different initial structures (icosahedron, truncated decahedron, octahedron, closed-shell fragment of an HCP structure, and non-closed-shell fragment of an HCP structure). The obtained order of the dynamic stabilities of the icosahedral-like clusters can be assigned to three groups, from stronger to weaker, according to the size ranges involved: (Zr, Al, Ti) > (Cu, Fe, Co, Ni, Mg, Ag) > (Pb, Au, Pd, Pt, Rh, Ir), which correspond to the predicted formation ability of the quasicrystals. The differences ofmore » the sequences can be explained by analyzing the parameters of the Gupta-type many-body inter-atomic potentials.« less

Molecular-dynamics study of solid-liquid interface migration in fcc metals

NASA Astrophysics Data System (ADS)

Mendelev, M. I.; Rahman, M. J.; Hoyt, J. J.; Asta, M.

2010-10-01

In order to establish a link between various structural and kinetic properties of metals and the crystal-melt interfacial mobility, free-solidification molecular-dynamics simulations have been performed for a total of nine embedded atom method interatomic potentials describing pure Al, Cu and Ni. To fully explore the space of materials properties three new potentials have been developed. The new potentials are based on a previous description of Al, but in each case the liquid structure, the melting point and/or the latent heat are varied considerably. The kinetic coefficient, μ, for all systems has been compared with several theoretical predictions. It is found that at temperatures close to the melting point the magnitude of μ correlates well with the value of the diffusion coefficient in the liquid.

Crystal structure of simple metals at high pressures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Degtyareva, Olga

2010-10-22

The effects of pressure on the crystal structure of simple (or sp-) elements are analysed in terms of changes in coordination number, packing density, and interatomic distances, and general rules are established. In the polyvalent elements from groups 14-17, the covalently bonded structures tend to transform to metallic phases with a gradual increase in coordination number and packing density, a behaviour normally expected under pressure. Group 1 and 2 metallic elements, however, show a reverse trend towards structures with low packing density due to intricate changes in their electronic structure. Complex crystal structures such as host-guest and incommensurately modulated structuresmore » found in these elements are given special attention in this review in an attempt to determine their role in the observed phase-transition sequences.« less

Finite-temperature interatomic exchange and magnon softening in Fe overlayers on Ir(001)

NASA Astrophysics Data System (ADS)

Rodrigues, D. C. M.; Szilva, A.; Klautau, A. B.; Bergman, A.; Eriksson, O.; Etz, C.

2016-07-01

We evaluate how thermal effects soften the magnon dispersion in 6 layers of Fe(001) on top of Ir(001). We perform a systematic study considering noncollinear spin arrangement and calculate configuration-dependent exchange parameters Jij n c following the methodology described by Szilva et al. [Phys. Rev. Lett. 111, 127204 (2013)], 10.1103/PhysRevLett.111.127204. In addition, Monte Carlo simulations were performed in order to estimate the noncollinear spin arrangement as a function of temperature. Hence the Jij n c's related to these configurations were calculated and used in an atomistic spin dynamics approach to evaluate the magnon spectra. Our results show good agreement with recent room-temperature measurements, and highlights how thermal effects produce magnon softening in this, and similar, systems.

Theoretical modeling of zircon's crystal morphology according to data of atomistic calculations

NASA Astrophysics Data System (ADS)

Gromalova, Natalia; Nikishaeva, Nadezhda; Eremin, Nikolay

2017-04-01

Zircon is an essential mineral that is used in the U-Pb dating. Moreover, zircon is highly resistant to radioactive exposure. It is of great interest in solving both fundamental and applied problems associated with the isolation of high-level radioactive waste. There is significant progress in forecasting of the most energetically favorable crystal structures at the present time. Unfortunately, the theoretical forecast of crystal morphology at high technological level is under-explored nowadays, though the estimation of crystal equilibrium habit is extremely important in studying the physical and chemical properties of new materials. For the first time, the thesis about relation of the equilibrium shape of a crystal with its crystal structure was put forward in the works by O.Brave. According to it, the idealized habit is determined in the simplest case by a correspondence with the reticular densities Rhkl of individual faces. This approach, along with all subsequent corrections, does not take into account the nature of atoms and the specific features of the chemical bond in crystals. The atomistic calculations of crystal surfaces are commonly performed using the energetic characteristics of faces, namely, the surface energy (Esurf), which is a measure of the thermodynamic stability of the crystal face. The stable crystal faces are characterized by small positive values of Esurf. As we know from our previous research (Gromalova et al.,2015) one of the constitutive factors affecting the value of the surface energy in calculations is a choice of potentials model. In this regard, we studied several sets of parameters of atomistic interatomic potentials optimized previously. As the first test model («Zircon 1») were used sets of interatomic potentials of interaction Zr-O, Si-O and O-O in the form of Buckingham potentials. To improve playback properties of zircon additionally used Morse potential for a couple of Zr-Si, as well as the three-particle angular harmonic potential. The other sets of interatomic potentials («Zircon 2, Zircon 3») differed from the first in that parameters was found with the help of quantum-chemical calculations of the structure «ab initio».The surface energies for different faces of zircon were calculated using Metadise code (Watson et al., 1996) at P4-3000 personal computer with Windows XP operating system. The computation time for one simple form was from 30 minutes to 12 hours. Calculations have shown that depending on the chosen model the surface energy of zircons faces several changes. For example, Esurf of face (331) obtained using models of potentials «Zircon 2», «Zircon 3» sufficiently similar (2.82 and 3.01 J/mol2 respectively). Meaning of Esurf of this face, calculated on the basis of set «Zircon 1» significantly lower (1,54 J/mol2). With regard to the face (100), it has low surface energies when selecting all three models, with a minimum value (1,14 J/mol2) in the model «Zircon 1». References: Gromalova N.A., Eremin N.N., Urusov V.S. Atomistic computer modeling of the crystal-morpology of corundum group minerals // Zapiski RMO. V. 144. №4. 2015. p. 84-92. Watson G.W., Kelsey E.T., de Leeuw N.H., Harris D.J, Parker S.C. Atomistic simulation of dislocations, surfaces and interfaces in MgO. Journal of the Chemical Society Faraday Transactions. 1996. V.92 P. 433-438.

Atomistically determined phase-field modeling of dislocation dissociation, stacking fault formation, dislocation slip, and reactions in fcc systems

NASA Astrophysics Data System (ADS)

Rezaei Mianroodi, Jaber; Svendsen, Bob

2015-04-01

The purpose of the current work is the development of a phase field model for dislocation dissociation, slip and stacking fault formation in single crystals amenable to determination via atomistic or ab initio methods in the spirit of computational material design. The current approach is based in particular on periodic microelasticity (Wang and Jin, 2001; Bulatov and Cai, 2006; Wang and Li, 2010) to model the strongly non-local elastic interaction of dislocation lines via their (residual) strain fields. These strain fields depend in turn on phase fields which are used to parameterize the energy stored in dislocation lines and stacking faults. This energy storage is modeled here with the help of the "interface" energy concept and model of Cahn and Hilliard (1958) (see also Allen and Cahn, 1979; Wang and Li, 2010). In particular, the "homogeneous" part of this energy is related to the "rigid" (i.e., purely translational) part of the displacement of atoms across the slip plane, while the "gradient" part accounts for energy storage in those regions near the slip plane where atomic displacements deviate from being rigid, e.g., in the dislocation core. Via the attendant global energy scaling, the interface energy model facilitates an atomistic determination of the entire phase field energy as an optimal approximation of the (exact) atomistic energy; no adjustable parameters remain. For simplicity, an interatomic potential and molecular statics are employed for this purpose here; alternatively, ab initio (i.e., DFT-based) methods can be used. To illustrate the current approach, it is applied to determine the phase field free energy for fcc aluminum and copper. The identified models are then applied to modeling of dislocation dissociation, stacking fault formation, glide and dislocation reactions in these materials. As well, the tensile loading of a dislocation loop is considered. In the process, the current thermodynamic picture is compared with the classical mechanical one as based on the Peach-Köhler force.

Modeling material interfaces with hybrid adhesion method

DOE PAGES

Brown, Nicholas Taylor; Qu, Jianmin; Martinez, Enrique

2017-01-27

A molecular dynamics simulation approach is presented to approximate layered material structures using discrete interatomic potentials through classical mechanics and the underlying principles of quantum mechanics. This method isolates the energetic contributions of the system into two pure material layers and an interfacial region used to simulate the adhesive properties of the diffused interface. The strength relationship of the adhesion contribution is calculated through small-scale separation calculations and applied to the molecular surfaces through an inter-layer bond criterion. By segregating the contributions into three regions and accounting for the interfacial excess energies through the adhesive surface bonds, it is possiblemore » to model each material with an independent potential while maintaining an acceptable level of accuracy in the calculation of mechanical properties. This method is intended for the atomistic study of the delamination mechanics, typically observed in thin-film applications. Therefore, the work presented in this paper focuses on mechanical tensile behaviors, with observations in the elastic modulus and the delamination failure mode. To introduce the hybrid adhesion method, we apply the approach to an ideal bulk copper sample, where an interface is created by disassociating the force potential in the middle of the structure. Various mechanical behaviors are compared to a standard EAM control model to demonstrate the adequacy of this approach in a simple setting. In addition, we demonstrate the robustness of this approach by applying it on (1) a Cu-Cu 2O interface with interactions between two atom types, and (2) an Al-Cu interface with two dissimilar FCC lattices. These additional examples are verified against EAM and COMB control models to demonstrate the accurate simulation of failure through delamination, and the formation and propagation of dislocations under loads. Finally, the results conclude that by modeling the energy contributions of an interface using hybrid adhesion bonds, we can provide an accurate approximation method for studies of large-scale mechanical properties, as well as the representation of various delamination phenomena at the atomic scale.« less

``Making the Molecular Movie'': First Frames

NASA Astrophysics Data System (ADS)

Miller, R. J. Dwayne

2011-03-01

Femtosecond Electron Diffraction has enabled atomic resolution to structural changes as they occur, essentially watching atoms move in real time--directly observe transition states. This experiment has been referred to as ``making the molecular movie'' and has been previously discussed in the context of a gedanken experiment. With the recent development of femtosecond electron pulses with sufficient number density to execute single shot structure determinations, this experiment has been finally realized. A new concept in electron pulse generation was developed based on a solution to the N-body electron propagation problem involving up to 10,000 interacting electrons that has led to a new generation of extremely bright electron pulsed sources that minimizes space charge broadening effects. Previously thought intractable problems of determining t=0 and fully characterizing electron pulses on the femtosecond time scale have now been solved through the use of the laser pondermotive potential to provide a time dependent scattering source. Synchronization of electron probe and laser excitation pulses is now possible with an accuracy of 10 femtoseconds to follow even the fastest nuclear motions. The camera for the ``molecular movie'' is well in hand based on high bunch charge electron sources. Several movies depicting atomic motions during passage through structural transitions will be shown. Atomic level views of the simplest possible structural transition, melting, will be presented for a number of systems in which both thermal and purely electronically driven atomic displacements can be correlated to the degree of directional bonding. Optical manipulation of charge distributions and effects on interatomic forces/bonding can be directly observed through the ensuing atomic motions. New phenomena involving strongly correlated electron systems will be presented in which an exceptionally cooperative phase transitions has been observed. The primitive origin of molecular cooperativity has also been discovered in recent studies of molecular crystals. These new developments will be discussed in the context of developing the necessary technology to directly observe the structure-function correlation in biomolecules--the fundamental molecular basis of biological systems.

Short- and medium-range structure of multicomponent bioactive glasses and melts: An assessment of the performances of shell-model and rigid-ion potentials.

PubMed

Tilocca, Antonio

2008-08-28

Classical and ab initio molecular dynamics (MD) simulations have been carried out to investigate the effect of a different treatment of interatomic forces in modeling the structural properties of multicomponent glasses and melts. The simulated system is a soda-lime phosphosilicate composition with bioactive properties. Because the bioactivity of these materials depends on their medium-range structural features, such as the network connectivity and the Q(n) distribution (where Q(n) is a tetrahedral species bonded to n bridging oxygens) of silicon and phosphorus network formers, it is essential to assess whether, and up to what extent, classical potentials can reproduce these properties. The results indicate that the inclusion of the oxide ion polarization through a shell-model (SM) approach provides a more accurate representation of the medium-range structure compared to rigid-ion (RI) potentials. Insight into the causes of these improvements has been obtained by comparing the melt-and-quench transformation of a small sample of the same system, modeled using Car-Parrinello MD (CPMD), to the classical MD runs with SM and RI potentials. Both classical potentials show some limitations in reproducing the highly distorted structure of the melt denoted by the CPMD runs; however, the inclusion of polarization in the SM potential results in a better and qualitatively correct dynamical balance between the interconversion of Q(n) species during the cooling of the melt. This effect seems to reflect the slower decay of the fraction of structural defects during the cooling with the SM potential. Because these transient defects have a central role in mediating the Q(n) transformations, as previously proposed and confirmed by the current simulations, their presence in the melt is essential to produce an accurate final distribution of Q(n) species in the glass.

Improved prediction of heat of mixing and segregation in metallic alloys using tunable mixing rule for embedded atom method

NASA Astrophysics Data System (ADS)

Divi, Srikanth; Agrahari, Gargi; Ranjan Kadulkar, Sanket; Kumar, Sanjeet; Chatterjee, Abhijit

2017-12-01

Capturing segregation behavior in metal alloy nanoparticles accurately using computer simulations is contingent upon the availability of high-fidelity interatomic potentials. The embedded atom method (EAM) potential is a widely trusted interatomic potential form used with pure metals and their alloys. When limited experimental data is available, the A-B EAM cross-interaction potential for metal alloys AxB 1-x are often constructed from pure metal A and B potentials by employing a pre-defined ‘mixing rule’ without any adjustable parameters. While this approach is convenient, we show that for AuPt, NiPt, AgAu, AgPd, AuNi, NiPd, PtPd and AuPd such mixing rules may not even yield the correct alloy properties, e.g., heats of mixing, that are closely related to the segregation behavior. A general theoretical formulation based on scaling invariance arguments is introduced that addresses this issue by tuning the mixing rule to better describe alloy properties. Starting with an existing pure metal EAM potential that is used extensively in literature, we find that the mixing rule fitted to heats of mixing for metal solutions usually provides good estimates of segregation energies, lattice parameters and cohesive energy, as well as equilibrium distribution of metals within a nanoparticle using Monte Carlo simulations. While the tunable mixing rule generally performs better than non-adjustable mixing rules, the use of the tunable mixing rule may still require some caution. For e.g., in Pt-Ni system we find that the segregation behavior can deviate from the experimentally observed one at Ni-rich compositions. Despite this the overall results suggest that the same approach may be useful for developing improved cross-potentials with other existing pure metal EAM potentials as well. As a further test of our approach, mixing rule estimated from binary data is used to calculate heat of mixing in AuPdPt, AuNiPd, AuPtNi, AgAuPd and NiPtPd. Excellent agreement with experiments is observed for AuPdPt.

Size-dependent characterization of embedded Ge nanocrystals: Structural and thermal properties

NASA Astrophysics Data System (ADS)

Araujo, L. L.; Giulian, R.; Sprouster, D. J.; Schnohr, C. S.; Llewellyn, D. J.; Kluth, P.; Cookson, D. J.; Foran, G. J.; Ridgway, M. C.

2008-09-01

A combination of conventional and synchrotron-based techniques has been used to characterize the size-dependent structural and thermal properties of Ge nanocrystals (NCs) embedded in a silica (a-SiO2) matrix. Ge NC size distributions with four different diameters ranging from 4.0 to 9.0 nm were produced by ion implantation and thermal annealing as characterized with small-angle x-ray scattering and transmission electron microscopy. The NCs were well represented by the superposition of bulklike crystalline and amorphous environments, suggesting the formation of an amorphous layer separating the crystalline NC core and the a-SiO2 matrix. The amorphous fraction was quantified with x-ray-absorption near-edge spectroscopy and increased as the NC diameter decreased, consistent with the increase in surface-to-volume ratio. The structural parameters of the first three nearest-neighbor shells were determined with extended x-ray-absorption fine-structure (EXAFS) spectroscopy and evolved linearly with inverse NC diameter. Specifically, increases in total disorder, interatomic distance, and the asymmetry in the distribution of distances were observed as the NC size decreased, demonstrating that finite-size effects govern the structural properties of embedded Ge NCs. Temperature-dependent EXAFS measurements in the range of 15-300 K were employed to probe the mean vibrational frequency and the variation of the interatomic distance distribution (mean value, variance, and asymmetry) with temperature for all NC distributions. A clear trend of increased stiffness (higher vibrational frequency) and decreased thermal expansion with decreasing NC size was evident, confirming the close relationship between the variation of structural and thermal/vibrational properties with size for embedded Ge NCs. The increase in surface-to-volume ratio and the presence of an amorphous Ge layer separating the matrix and crystalline NC core are identified as the main factors responsible for the observed behavior, with the surrounding a-SiO2 matrix also contributing to a lesser extent. Such results are compared to previous reports and discussed in terms of the influence of the surface-to-volume ratio in objects of nanometer dimensions.

Theoretical study of the elasticity, mechanical behavior, electronic structure, interatomic bonding, and dielectric function of an intergranular glassy film model in prismatic β-Si3N4

NASA Astrophysics Data System (ADS)

Ching, W. Y.; Rulis, Paul; Ouyang, Lizhi; Aryal, Sitaram; Misra, Anil.

2010-06-01

Microstructures such as intergranular glassy films (IGFs) are ubiquitous in many structural ceramics. They control many of the important physical properties of polycrystalline ceramics and can be influenced during processing to modify the performance of devices that contain them. In recent years, there has been intense research, both experimentally and computationally, on the structure and properties of IGFs. Unlike grain boundaries or dislocations with well-defined crystalline planes, the atomic scale structure of IGFs, their fundamental electronic interactions, and their bonding characteristics are far more complicated and not well known. In this paper, we present the results of theoretical simulations using ab initio methods on an IGF model in β-Si3N4 with prismatic crystalline planes. The 907-atom model has a dimension of 14.533Å×15.225Å×47.420Å . The IGF layer is perpendicular to the z axis, 16.4Å wide, and contains 72 Si, 32 N, and 124 O atoms. Based on this model, the mechanical and elastic properties, the electronic structure, the interatomic bonding, the localization of defective states, the distribution of electrostatic potential, and the optical dielectric function are evaluated and compared with crystalline β-Si3N4 . We have also performed a theoretical tensile experiment on this model by incrementally extending the structure in the direction perpendicular to the IGF plane until the model fully separated. It is shown that fracture occurs at a strain of 9.42% with a maximum stress of 13.9 GPa. The fractured segments show plastic behavior and the formation of surfacial films on the β-Si3N4 . These results are very different from those of a previously studied basal plane model [J. Chen , Phys. Rev. Lett. 95, 256103 (2005)10.1103/PhysRevLett.95.256103] and add insights to the structure and behavior of IGFs in polycrystalline ceramics. The implications of these results and the need for further investigations are discussed.

Ab-initio molecular dynamics simulations of liquid Hg-Pb alloys

NASA Astrophysics Data System (ADS)

Sharma, Nalini; Thakur, Anil; Ahluwalia, P. K.

2014-04-01

Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-Pb alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Three liquid Hg-Pb mixtures (Hg30Pb70, Hg50Pb50 and Hg90Pb10) at 600K are considered. The radial distribution function g(r) and structure factor S(q) of considered alloys are compared with respective experimental results for liquid Hg (l-Hg) and lead (l-Pb). The radial distribution function g(r) shows the presence of short range order in the systems considered. Smooth curves of Bhatia-Thornton partial structure factors factor shows the presence of liquid state in the considered three alloys. Among the all considered alloys, Hg50Pb50 alloy shows presence of more chemical ordering and presence of hetero-coordination.

Crystal structure prediction supported by incomplete experimental data

NASA Astrophysics Data System (ADS)

Tsujimoto, Naoto; Adachi, Daiki; Akashi, Ryosuke; Todo, Synge; Tsuneyuki, Shinji

2018-05-01

We propose an efficient theoretical scheme for structure prediction on the basis of the idea of combining methods, which optimize theoretical calculation and experimental data simultaneously. In this scheme, we formulate a cost function based on a weighted sum of interatomic potential energies and a penalty function which is defined with partial experimental data totally insufficient for conventional structure analysis. In particular, we define the cost function using "crystallinity" formulated with only peak positions within the small range of the x-ray-diffraction pattern. We apply this method to well-known polymorphs of SiO2 and C with up to 108 atoms in the simulation cell and show that it reproduces the correct structures efficiently with very limited information of diffraction peaks. This scheme opens a new avenue for determining and predicting structures that are difficult to determine by conventional methods.

Gap discrete breathers in strained boron nitride

NASA Astrophysics Data System (ADS)

Barani, Elham; Korznikova, Elena A.; Chetverikov, Alexander P.; Zhou, Kun; Dmitriev, Sergey V.

2017-11-01

Linear and nonlinear dynamics of hexagonal boron nitride (h-BN) lattice is studied by means of molecular dynamics simulations with the use of the Tersoff interatomic potentials. It is found that sufficiently large homogeneous elastic strain along zigzag direction opens a wide gap in the phonon spectrum. Extended vibrational mode with boron and nitrogen sublattices vibrating in-plane as a whole in strained h-BN has frequency within the phonon gap. This fact suggests that a nonlinear spatially localized vibrational mode with frequencies in the phonon gap, called discrete breather (also often termed as intrinsic localized mode), can be excited. Properties of the gap discrete breathers in strained h-BN are contrasted with that for analogous vibrational mode found earlier in strained graphene. It is found that h-BN modeled with the Tersoff potentials does not support transverse discrete breathers.

Detonation Product EOS Studies: Using ISLS to Refine Cheetah

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zaug, J M; Howard, W M; Fried, L E

2001-08-08

Knowledge of an effective interatomic potential function underlies any effort to predict or rationalize the properties of solids and liquids. The experiments we undertake are directed towards determination of equilibrium and dynamic properties of simple fluids at densities sufficiently high that traditional computational methods and semi-empirical forms successful at ambient conditions may require reconsideration. In this paper we present high-pressure and temperature experimental sound speed data on a simple fluid, methanol. Impulsive Stimulated Light Scattering (ISLS) conducted on diamond-anvil cell (DAC) encapsulated samples offers an experimental approach to determine cross-pair potential interactions through equation of state determinations. In addition themore » kinetics of structural relaxation in fluids can be studied. We compare our experimental results with our thermochemical computational model Cheetah. Computational models are systematically improved with each addition of experimental data.« less

Energetic Analysis of Conjugated Hydrocarbons Using the Interacting Quantum Atoms Method.

PubMed

Jara-Cortés, Jesús; Hernández-Trujillo, Jesús

2018-07-05

A number of aromatic, antiaromatic, and nonaromatic organic molecules was analyzed in terms of the contributions to the electronic energy defined in the quantum theory of atoms in molecules and the interacting quantum atoms method. Regularities were found in the exchange and electrostatic interatomic energies showing trends that are closely related to those of the delocalization indices defined in the theory. In particular, the CC interaction energies between bonded atoms allow to rationalize the energetic stabilization associated with the bond length alternation in conjugated polyenes. This approach also provides support to Clar's sextet rules devised for aromatic systems. In addition, the H⋯H bonding found in some of the aromatic molecules studied was of an attractive nature, according to the stabilizing exchange interaction between the bonded H atoms. © 2017 Wiley Periodicals, Inc. © 2017 Wiley Periodicals, Inc.

Influence of Fe/Co ratio on structural and magnetic properties of (Fe100-xCox)84.5Nb5B8.5P2 alloy

NASA Astrophysics Data System (ADS)

Gehlot, K.; Kane, S. N.; Sinha, A. K.; Ghodke, N.; Varga, L. K.

2018-05-01

Structural and magnetic properties of a series of (Fe100-xCox)84.5Nb5B8.5P2 (x = 20, 40, 60) have been investigated respectively by using synchrotron x-ray diffraction and magnetic measurements. Results show that Fe/Co ratio: i) affects stability of the alloy against crystallization, ii) shows evidence for ordering, which has considerable effect on magnetic properties, iii) influences the grain diameter and volume fraction of the formed nano-grains range between 4.8 - 9.5 nm and 1.5 - 9 %, affects magnetic properties considerably. An empirical relation is obtained, which shows linear relationship between interatomic distances for 1st, 2nd co-ordination shell, suggests strong correlation between structural, magnetic properties.

Influence of strain on dislocation core in silicon

NASA Astrophysics Data System (ADS)

Pizzagalli, L.; Godet, J.; Brochard, S.

2018-05-01

First principles, density functional-based tight binding and semi-empirical interatomic potentials calculations are performed to analyse the influence of large strains on the structure and stability of a 60? dislocation in silicon. Such strains typically arise during the mechanical testing of nanostructures like nanopillars or nanoparticles. We focus on bi-axial strains in the plane normal to the dislocation line. Our calculations surprisingly reveal that the dislocation core structure largely depends on the applied strain, for strain levels of about 5%. In the particular case of bi-axial compression, the transformation of the dislocation to a locally disordered configuration occurs for similar strain magnitudes. The formation of an opening, however, requires larger strains, of about 7.5%. Furthermore, our results suggest that electronic structure methods should be favoured to model dislocation cores in case of large strains whenever possible.

Atomic quantum corrals for Bose-Einstein condensates

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xiong Hongwei; Kavli Institute for Theoretical Physics China, Chinese Academy of Sciences, Beijing 100190; Wu Biao

2010-11-15

We consider the dynamics of Bose-Einstein condensates in a corral-like potential. Compared to the electronic quantum corrals, the atomic quantum corrals have the advantages of allowing direct and convenient observation of the wave dynamics, together with adjustable interaction strength. Our numerical study shows that these advantages not only allow exploration of the rich dynamical structures in the density distribution but also make the corrals useful in many other aspects. In particular, the corrals for atoms can be arranged into a stadium shape for the experimental visualization of quantum chaos, which has been elusive with electronic quantum corrals. The density correlationmore » is used to describe quantitatively the dynamical quantum chaos. Furthermore, we find that the interatomic interaction can greatly enhance the dynamical quantum chaos, for example, inducing a chaotic behavior even in circle-shaped corrals.« less

Adiabatic passage of radio-frequency-assisted Förster resonances in Rydberg atoms for two-qubit gates and the generation of Bell states

NASA Astrophysics Data System (ADS)

Beterov, I. I.; Hamzina, G. N.; Yakshina, E. A.; Tretyakov, D. B.; Entin, V. M.; Ryabtsev, I. I.

2018-03-01

High-fidelity entangled Bell states are of great interest in quantum physics. Entanglement of ultracold neutral atoms in two spatially separated optical dipole traps is promising for implementation of quantum computing and quantum simulation and for investigation of Bell states of material objects. We propose a method to entangle two atoms via long-range Rydberg-Rydberg interaction. Alternative to previous approaches, based on Rydberg blockade, we consider radio-frequency-assisted Stark-tuned Förster resonances in Rb Rydberg atoms. To reduce the sensitivity of the fidelity of Bell states to the fluctuations of interatomic distance, we propose to use the double adiabatic passage across the radio-frequency-assisted Stark-tuned Förster resonances, which results in a deterministic phase shift of the collective two-atom state.

A database of new zeolite-like materials.

PubMed

Pophale, Ramdas; Cheeseman, Phillip A; Deem, Michael W

2011-07-21

We here describe a database of computationally predicted zeolite-like materials. These crystals were discovered by a Monte Carlo search for zeolite-like materials. Positions of Si atoms as well as unit cell, space group, density, and number of crystallographically unique atoms were explored in the construction of this database. The database contains over 2.6 M unique structures. Roughly 15% of these are within +30 kJ mol(-1) Si of α-quartz, the band in which most of the known zeolites lie. These structures have topological, geometrical, and diffraction characteristics that are similar to those of known zeolites. The database is the result of refinement by two interatomic potentials that both satisfy the Pauli exclusion principle. The database has been deposited in the publicly available PCOD database and in www.hypotheticalzeolites.net/database/deem/. This journal is © the Owner Societies 2011

Study of the high-pressure helium phase diagram using molecular dynamics

NASA Astrophysics Data System (ADS)

Koci, L.; Ahuja, R.; Belonoshko, A. B.; Johansson, B.

2007-01-01

The rich occurrence of helium and hydrogen in space makes their properties highly interesting. By means of molecular dynamics (MD), we have examined two interatomic potentials for 4He. Both potentials are demonstrated to reproduce high-pressure solid and liquid equation of state (EOS) data. The EOS, solid-solid transitions and melting at high pressures (P) were studied using a two-phase method. The Buckingham potential shows a good agreement with theoretical and experimental EOS, but does not reproduce experimental melting data. The Aziz potential shows a perfect match with theoretical melting data. We conclude that there is a stable body-centred-cubic (bcc) phase for 4He at temperatures (T) above 340 K and pressures above 22 GPa for the Buckingham potential, whereas no bcc phase is found for the Aziz potential in the applied PT range.

Entropic vs. elastic models of fragility of glass-forming liquids: Two sides of the same coin?

NASA Astrophysics Data System (ADS)

Sen, Sabyasachi

2012-10-01

The two most influential atomistic models that have been proposed in the literature to explain the temperature dependent activation energy of viscous flow of a glass-forming liquid, i.e., its fragility, are the configurational entropy model of Adam and Gibbs [J. Chem. Phys. 43, 139 (1965), 10.1063/1.1696442] and the elastic "shoving" model of Dyre et al. [J. Non-Cryst. Solids 352, 4635 (2006), 10.1016/j.jnoncrysol.2006.02.173]. Here we demonstrate a qualitative equivalence between these two models starting from the well-established general relationships between the interatomic potentials, elastic constants, structural rearrangement, and entropy in amorphous materials. The unification of these two models provides important predictions that are consistent with experimental observations and shed new light into the problem of glass transition.

A revised set of values of single-bond radii derived from the observed interatomic distances in metals by correction for bond number and resonance energy

PubMed Central

Pauling, Linus; Kamb, Barclay

1986-01-01

An earlier discussion [Pauling, L. (1947) J. Am. Chem. Soc. 69, 542] of observed bond lengths in elemental metals with correction for bond number and resonance energy led to a set of single-bond metallic radii with values usually somewhat less than the corresponding values obtained from molecules and complex ions. A theory of resonating covalent bonds has now been developed that permits calculation of the number of resonance structures per atom and of the effective resonance energy per bond. With this refined method of correcting the observed bond lengths for the effect of resonance energy, a new set of single-bond covalent radii, in better agreement with values from molecules and complex ions, has been constructed. PMID:16593698

Feedback control of an interacting Bose-Einstein condensate using phase-contrast imaging

NASA Astrophysics Data System (ADS)

Szigeti, S. S.; Hush, M. R.; Carvalho, A. R. R.; Hope, J. J.

2010-10-01

The linewidth of an atom laser is limited by density fluctuations in the Bose-Einstein condensate (BEC) from which the atom laser beam is outcoupled. In this paper we show that a stable spatial mode for an interacting BEC can be generated using a realistic control scheme that includes the effects of the measurement backaction. This model extends the feedback theory, based on a phase-contrast imaging setup, presented by Szigeti, Hush, Carvalho, and Hope [Phys. Rev. APLRAAN1050-294710.1103/PhysRevA.80.013614 80, 013614 (2009)]. In particular, it is applicable to a BEC with large interatomic interactions and solves the problem of inadequacy of the mean-field (coherent state) approximation by utilizing a fixed number state approximation. Our numerical analysis shows the control to be more effective for a condensate with a large nonlinearity.

Feedback control of an interacting Bose-Einstein condensate using phase-contrast imaging

DOE Office of Scientific and Technical Information (OSTI.GOV)

Szigeti, S. S.; Hush, M. R.; Carvalho, A. R. R.

2010-10-15

The linewidth of an atom laser is limited by density fluctuations in the Bose-Einstein condensate (BEC) from which the atom laser beam is outcoupled. In this paper we show that a stable spatial mode for an interacting BEC can be generated using a realistic control scheme that includes the effects of the measurement backaction. This model extends the feedback theory, based on a phase-contrast imaging setup, presented by Szigeti, Hush, Carvalho, and Hope [Phys. Rev. A 80, 013614 (2009)]. In particular, it is applicable to a BEC with large interatomic interactions and solves the problem of inadequacy of the mean-fieldmore » (coherent state) approximation by utilizing a fixed number state approximation. Our numerical analysis shows the control to be more effective for a condensate with a large nonlinearity.« less

Molecular dynamics of shock loading of metals with defects

DOE Office of Scientific and Technical Information (OSTI.GOV)

Belak, J.F.

1997-12-31

The finite rise time of shock waves in metals is commonly attributed to dissipative or viscous behavior of the metal. This viscous or plastic behavior is commonly attributed to the motion of defects such as dislocations. Despite this intuitive understanding, the experimental observation of defect motion or nucleation during shock loading has not been possible due to the short time scales involved. Molecular dynamics modeling with realistic interatomic potentials can provide some insight into defect motion during shock loading. However, until quite recently, the length scale required to accurately represent a metal with defects has been beyond the scope ofmore » even the most powerful supercomputers. Here, the author presents simulations of the shock response of single defects and indicate how simulation might provide some insight into the shock loading of metals.« less

Structural and electronic properties of chiral single-wall copper nanotubes

NASA Astrophysics Data System (ADS)

Duan, YingNi; Zhang, JianMin; Xu, KeWei

2014-04-01

The structural, energetic and electronic properties of chiral ( n, m) (3⩽ n⩽6, n/2⩽ m⩽ n) single-wall copper nanotubes (CuNTs) have been investigated by using projector-augmented wave method based on density-functional theory. The (4, 3) CuNT is energetically stable and should be observed experimentally in both free-standing and tip-suspended conditions, whereas the (5, 5) and (6, 4) CuNTs should be observed in free-standing and tip-suspended conditions, respectively. The number of conductance channels in the CuNTs does not always correspond to the number of atomic strands comprising the nanotube. Charge density contours show that there is an enhanced interatomic interaction in CuNTs compared with Cu bulk. Current transporting states display different periods and chirality, the combined effects of which lead to weaker chiral currents on CuNTs.

Decomposition of Proteins into Dynamic Units from Atomic Cross-Correlation Functions.

PubMed

Calligari, Paolo; Gerolin, Marco; Abergel, Daniel; Polimeno, Antonino

2017-01-10

In this article, we present a clustering method of atoms in proteins based on the analysis of the correlation times of interatomic distance correlation functions computed from MD simulations. The goal is to provide a coarse-grained description of the protein in terms of fewer elements that can be treated as dynamically independent subunits. Importantly, this domain decomposition method does not take into account structural properties of the protein. Instead, the clustering of protein residues in terms of networks of dynamically correlated domains is defined on the basis of the effective correlation times of the pair distance correlation functions. For these properties, our method stands as a complementary analysis to the customary protein decomposition in terms of quasi-rigid, structure-based domains. Results obtained for a prototypal protein structure illustrate the approach proposed.

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

NASA Technical Reports Server (NTRS)

Lawson, JOhn W.; Daw, Murray S.; Bauschlicher, Charles W.

2011-01-01

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are candidate materials for applications in extreme environments because of their high melting point, good mechanical properties and reasonable oxidation resistance. Unlike many ceramics, these materials have high thermal conductivity which can be advantageous, for example, to reduce thermal shock. Recently, we developed Tersoff style interatomic potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current. Results at room temperature and at elevated temperatures will be reported.

Development of interatomic potentials appropriate for simulation of devitrification of Al 90Sm 10 alloy

DOE PAGES

Mendelev, M. I.; Zhang, F.; Ye, Z.; ...

2015-04-23

In this study, a semi-empirical potential for the Al 90Sm 10 alloy is presented. The potential provides satisfactory reproduction of pure Al properties, the formation energies of a set of Al–Sm crystal phases with Sm content about 10%, and the structure of the liquid Al 90Sm 10 alloy. During molecular dynamics simulation in which the liquid alloy is cooled at a rate of 10 10 K/s, the developed potential produces a glass structure with lower ab initio energy than that produced by ab initio molecular dynamics (AIMD) itself using a typical AIMD cooling rate of 8 ∙10 13 K/s. Basedmore » on these facts the developed potential should be suitable for simulations of phase transformations in the Al 90Sm 10 alloy.« less

Hetero-diffusion of Au epitaxy on stepped Ag(110) surface: Study of the jump rate and diffusion coefficient

NASA Astrophysics Data System (ADS)

Benlattar, M.; El koraychy, E.; Kotri, A.; Mazroui, M.

2017-12-01

We have used molecular dynamics simulations combined with an interatomic potential derived from the embedded atom method, to investigate the hetero-diffusion of Au adatom near a stepped Ag(110) surface with the height of one monoatomic layer. The activation energies for different diffusion processes, which occur on the terrace and near the step edge, are calculated both by molecular statics and molecular dynamics simulations. Static energies are found by the drag method, whereas the dynamic barriers are computed at high temperature from the Arrhenius plots. Our numerical results reveal that the jump process requires very high activation energy compared to the exchange process either on the terrace or near the step edge. In this work, other processes, such as upward and downward diffusion at step edges, have also been discussed.

cDF Theory Software for mesoscopic modeling of equilibrium and transport phenomena

DOE Office of Scientific and Technical Information (OSTI.GOV)

2015-12-01

The approach is based on classical Density Functional Theory ((cDFT) coupled with the Poisson-Nernst-Planck (PNP) transport kinetics model and quantum mechanical description of short-range interaction and elementary transport processes. The model we proposed and implemented is fully atomistic, taking into account pairwise short-range and manybody long-range interactions. But in contrast to standard molecular dynamics (MD) simulations, where long-range manybody interactions are evaluated as a sum of pair-wise atom-atom contributions, we include them analytically based on wellestablished theories of electrostatic and excluded volume interactions in multicomponent systems. This feature of the PNP/cDFT approach allows us to reach well beyond the length-scalesmore » accessible to MD simulations, while retaining the essential physics of interatomic interactions from first principles and in a parameter-free fashion.« less

Potential energy curves of the Na2+ molecular ion from all-electron ab initio relativistic calculations

NASA Astrophysics Data System (ADS)

Bewicz, Anna; Musiał, Monika; Kucharski, Stanisław A.

2017-11-01

The equation-of-motion coupled-cluster method for electron affinity calculations has been used to study potential energy curves (PECs) for the Na+2 molecular ion. Although the studied molecule represents the open shell system the applied approach employs the closed shell Na+ 22 ion as the reference. In addition the Na+ 22 system dissociates into the closed shell fragments; hence, the restricted Hartree-Fock scheme can be used within the whole range of interatomic distances, from 2 to 45 Å. We used large basis set engaging 268 basis functions with all 21 electrons correlated. The relativistic effects are included via second-order Douglas-Kroll method. The computed PECs, spectroscopic molecular constants and vibrational energy levels agree well with experimental values if the latter are available or with other theoretical data.

Evaluation of melting point of UO 2 by molecular dynamics simulation

NASA Astrophysics Data System (ADS)

Arima, Tatsumi; Idemitsu, Kazuya; Inagaki, Yaohiro; Tsujita, Yuichi; Kinoshita, Motoyasu; Yakub, Eugene

2009-06-01

The melting point of UO 2 has been evaluated by molecular dynamics simulation (MD) in terms of interatomic potential, pressure and Schottky defect concentration. The Born-Mayer-Huggins potentials with or without a Morse potential were explored in the present study. Two-phase simulation whose supercell at the initial state consisted of solid and liquid phases gave the melting point comparable to the experimental data using the potential proposed by Yakub. The heat of fusion was determined by the difference in enthalpy at the melting point. In addition, MD calculations showed that the melting point increased with pressure applied to the system. Thus, the Clausius-Clapeyron equation was verified. Furthermore, MD calculations clarified that an addition of Schottky defects, which generated the local disorder in the UO 2 crystal, lowered the melting point.

Physics of Shock Compression and Release: NEMD Simulations of Tantalum and Silicon

NASA Astrophysics Data System (ADS)

Hahn, Eric; Meyers, Marc; Zhao, Shiteng; Remington, Bruce; Bringa, Eduardo; Germann, Tim; Ravelo, Ramon; Hammerberg, James

2015-06-01

Shock compression and release allow us to evaluate physical deformation and damage mechanisms occurring in extreme environments. SPaSM and LAMMPS molecular dynamics codes were employed to simulate single and polycrystalline tantalum and silicon at strain rates above 108 s-1. Visualization and analysis was accomplished using OVITO, Crystal Analysis Tool, and a redesigned orientation imaging function implemented into SPaSM. A comparison between interatomic potentials for both Si and Ta (as pertaining to shock conditions) is conducted and the influence on phase transformation and plastic relaxation is discussed. Partial dislocations, shear induced disordering, and metastable phase changes are observed in compressed silicon. For tantalum, the role of grain boundary and twin intersections are evaluated for their role in ductile spallation. Finally, the temperature dependent response of both Ta and Si is investigated.

Pairing mechanism in Bi-O superconductors: A finite-size chain calculation

NASA Astrophysics Data System (ADS)

Aligia, A. A.; Nuez Regueiro, M. D.; Gagliano, E. R.

1989-09-01

We have studied the pairing mechanism in BiO3 systems by calculating the binding energy of a pair of holes in finite Bi-O chains, for parameters that simulate three-dimensional behavior. In agreement with previous results using perturbation theory in the hopping t, for covalent Bi-O binding and parameters for which the parent compound has a disproportionate ground state, pairing induced by the presence of biexcitons is obtained for sufficiently large interatomic Coulomb repulsion. The analysis of appropriate correlation functions shows a rapid metallization of the system as t and the number of holes increase. This fact shrinks the region of parameters for which the finite-size calculations can be trusted without further study. The same model for other parameters yields pairing in two other regimes: bipolaronic and magnetic excitonic.

An empirical potential for simulating vacancy clusters in tungsten.

PubMed

Mason, D R; Nguyen-Manh, D; Becquart, C S

2017-12-20

We present an empirical interatomic potential for tungsten, particularly well suited for simulations of vacancy-type defects. We compare energies and structures of vacancy clusters generated with the empirical potential with an extensive new database of values computed using density functional theory, and show that the new potential predicts low-energy defect structures and formation energies with high accuracy. A significant difference to other popular embedded-atom empirical potentials for tungsten is the correct prediction of surface energies. Interstitial properties and short-range pairwise behaviour remain similar to the Ackford-Thetford potential on which it is based, making this potential well-suited to simulations of microstructural evolution following irradiation damage cascades. Using atomistic kinetic Monte Carlo simulations, we predict vacancy cluster dissociation in the range 1100-1300 K, the temperature range generally associated with stage IV recovery.

Specific features of nonvalent interactions in orthorhombic perovskites

NASA Astrophysics Data System (ADS)

Serezhkin, V. N.; Pushkin, D. V.; Serezhkina, L. B.

2014-07-01

It is established that isostructural orthorhombic perovskites ABO3 (sp. gr. Pnma in different systems, no. 62, Z = 4), depending on the specificity of nonvalent interactions (which determine the combinatorial-topological type of the Voronoi-Dirichlet polyhedra (VDPs) of four basis atoms), are divided into ten different stereotypes. It is shown by the example of 259 perovskites belonging to the DyCrO3 stereotype that VDP characteristics can be used to quantitatively estimate the distortion of BO6 octahedra, including that caused by the Jahn-Teller effect. It is found that one of the causes of the distortion of the coordination polyhedra of atoms in the structure of orthorhombic perovskites is heteroatomic metal-metal interactions, for which the interatomic distances are much shorter than the sum of the Slater radii of A and B atoms.

Study of atomic structure of liquid Hg-In alloys using ab-initio molecular dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sharma, Nalini; Ahluwalia, P. K.; Thakur, Anil

2015-05-15

Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-In alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Five liquid Hg-In mixtures (Hg{sub 10}In{sub 90}, Hg{sub 30}In{sub 70}, Hg{sub 50}In{sub 50}, Hg{sub 70}In{sub 30} and Hg{sub 90}In{sub 10}) at 299K are considered. The radial distribution function g(r) and structure factor S(q) of considered alloys are compared with respective experimental results for liquid Hg (l-Hg) and (l-In). The radial distribution function g(r) shows the presence of short range order in the systems considered. Smooth curves of Bhatia-Thornton partial structure factors factormore » shows the presence of liquid state in the considered alloys.« less

Crucial role of nuclear dynamics for electron injection in a dye–semiconductor complex

DOE PAGES

Monti, Adriano; Negre, Christian F. A.; Batista, Victor S.; ...

2015-06-05

In this study, we investigate the electron injection from a terrylene-based chromophore to the TiO 2 semiconductor bridged by a recently proposed phenyl-amide-phenyl molecular rectifier. The mechanism of electron transfer is studied by means of quantum dynamics simulations using an extended Hückel Hamiltonian. It is found that the inclusion of the nuclear motion is necessary to observe the photoinduced electron transfer. In particular, the fluctuations of the dihedral angle between the terrylene and the phenyl ring modulate the localization and thus the electronic coupling between the donor and acceptor states involved in the injection process. The electron propagation shows characteristicmore » oscillatory features that correlate with interatomic distance fluctuations in the bridge, which are associated with the vibrational modes driving the process. The understanding of such effects is important for the design of functional dyes with optimal injection and rectification properties.« less

DOE Office of Scientific and Technical Information (OSTI.GOV)

Smedskjaer, Morten M., E-mail: mos@bio.aau.dk; Bauchy, Mathieu; Mauro, John C.

The properties of glass are determined not only by temperature, pressure, and composition, but also by their complete thermal and pressure histories. Here, we show that glasses of identical composition produced through thermal annealing and through quenching from elevated pressure can result in samples with identical density and mean interatomic distances, yet different bond angle distributions, medium-range structures, and, thus, macroscopic properties. We demonstrate that hardness is higher when the density increase is obtained through thermal annealing rather than through pressure-quenching. Molecular dynamics simulations reveal that this arises because pressure-quenching has a larger effect on medium-range order, while annealing hasmore » a larger effect on short-range structures (sharper bond angle distribution), which ultimately determine hardness according to bond constraint theory. Our work could open a new avenue towards industrially useful glasses that are identical in terms of composition and density, but with differences in thermodynamic, mechanical, and rheological properties due to unique structural characteristics.« less

Hydrogen bonds in concreto and in computro

NASA Astrophysics Data System (ADS)

Stouten, Pieter F. W.; Kroon, Jan

1988-07-01

Molecular dynamics simulations of liquid water and liquid methanol have been carried out. For both liquids an effective pair potential was used. The models were fitted to the heat of vaporization, pressure and various radial distribution functions resulting from diffraction experiments on liquids. In both simulations 216 molecules were put in a cubic periodical ☐. The system was loosely coupled to a temperature bath and to a pressure bath. Following an initial equilibration period relevant data were sampled during 15 ps. The distributions of oxygen—oxygen distances in hydrogen bonds obtained from the two simulations are essentially the same. The distribution obtained from crystal data is somewhat different: the maximum has about the same position, but the curve is much narrower, which can be expected merely from the fact that diffraction experiments only supply average atomic positions and hence average interatomic distances. When thermal motion is taken into account a closer likeness is observed.

Ab initio study of the electrostatic multipole nature of torsional potentials in CH3SSCH3, CH3SSH, and HOOH

NASA Technical Reports Server (NTRS)

Sokalski, W. A.; Lai, J.; Luo, N.; Sun, S.; Shibata, M.; Ornstein, R.; Rein, R.

1991-01-01

The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.