Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effect of hydride vapor phase epitaxy growth conditions on the degree of atomic ordering in GaInP

We report a systematic study of CuPt-type ordering in hydride vapor phase epitaxy (HVPE)-grown Ga 0.5 In 0.5 P. Selected-area electron diffraction reveals ordering in samples grown on (001) GaAs substrates offcut toward (1$\overline{1}$1)B. The ordering is single-variant, occurring only on (1$\overline{1}$1)B planes and not on ($\overline{1}$11)B. Quantitative analysis of the order parameter by high-resolution x-ray diffraction (HRXRD) indicates that ordering increases with deposition temperature in samples grown at 600–700 °C with a constant gas-phase V/III ratio ~3. Ordering increases with V/III ratio in the range of 1.3–6.7 at a constant deposition temperature of 650°C. Photoluminescence measurements correlate the order parameter with Ga0.5In0.5P bandgap contraction, though the contraction is larger than expected based on the magnitude of order parameters measured by HRXRD. A possible reason for this discrepancy is that the photoluminescence emission occurs in the lower bandgap ordered domains, which are small and evenly dispersed throughout the material. We also show that the degree of ordering decreases with growth rate, disappearing at ~60 µm/h. The HVPE-grown material exhibits a generally weaker ordering than organometallic vapor phase epitaxy (OMVPE)-grown material, likely due to the moderate V/III ratios employed, in contrast to the V/III ratios in the 100s typical of OMVPE. However, the tendency for ordering to increase with V/III ratio suggests that the same dimer-induced stress mechanism used to explain the occurrence of ordering in OMVPE-grown material also applies to HVPE. The tendencies for ordering to increase with deposition temperature and decrease with growth rate show that kinetics limit the degree of ordering, also in agreement with OMVPE trends.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Machine learning prediction of self-diffusion in Lennard-Jones fluids

In this work, different machine learning (ML) methods were explored for the prediction of self-diffusion in Lennard-Jones (LJ) fluids. Using a database of diffusion constants obtained from the molecular dynamics simulation literature, multiple Random Forest (RF) and Artificial Neural Net (ANN) regression models were developed and characterized. The role and improved performance of feature engineering coupled to the RF model development was also addressed. The performance of these different ML models was evaluated by comparing the prediction error to an existing empirical relationship used to describe LJ fluid diffusion. It was found that the ANN regression models provided superior prediction of diffusion in comparison to the existing empirical relationships.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Isotope effects on the high pressure viscosity of liquid water measured by differential dynamic microscopy

In this paper, differential dynamic microscopy is performed in diamond anvil cells to measure the viscosity of water along the 24 °C isotherm to high-pressure by determination of the tracer diffusion coefficient of monodisperse silica spheres of known diameter and application of the Stokes-Einstein-Sutherland equation. This technique allows liquid samples to be compressed to greater pressure prior to freezing than with other viscometry methods. The highest-pressure measurement was made at 1.67 GPa, considerably deeper into the supercompressed regime than previously reported. The effect of isotopic composition is investigated with samples of normal water, heavy water, and partially deuterated water. When data below 0.25 GPa are excluded a free volume model fits the observed viscosities well yielding a theoretical glass transition density close to that observed in very-high-density amorphous ice. The improved fit above 0.25 GPa coincides with the loss of other anomalous behaviors in liquid water caused by hydrogen bonding and represents a transition to properties closer to those of a simple liquid.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

An Entropy-Maximization Approach to Automated Training Set Generation for Interatomic Potentials

Machine learning-based interatomic potentials are currently garnering a lot of attention as they strive to achieve the accuracy of electronic structure methods at the computational cost of empirical potentials. Given their generic functional forms, the transferability of these potentials is highly dependent on the quality of the training set, the generation of which can be highly labor-intensive. Good training sets should at once contain a very diverse set of configurations while avoiding redundancies that incur cost without providing benefits. We formalize these requirements in a local entropy-maximization framework and propose an automated sampling scheme to sample from this objective function. We show that this approach generates much more diverse training sets than unbiased sampling and is competitive with hand-crafted training sets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Displacement of carbon atoms in few-layer graphene

Molecular dynamics simulations were performed to study the susceptibility of carbon atom displacement under electron irradiation. The mapping of threshold displacement energies at different recoiling directions showed that the energies are very sensitive to the layer configurations and positions of neighboring atoms. Carbon atoms on the top and the bottom layers of few-layer graphene are most vulnerable to irradiation damage due to lack of constraints from the neighboring graphene layers. As indirect experiment evidence, transmission electron microscopy was performed on the edge of folded few-layer graphene, which made it possible to reveal “the inside” and compare irradiation tolerance of atoms at different layers, by using an electron analysis beam for both displacement creation and in situ characterization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Preparation of high orbital angular momentum Rydberg states by optical-millimeter-wave STIRAP

Rydberg states of molecules are intrinsically challenging to study due to the presence of fast non-radiative decay pathways, such as predissociation. However, selectively exciting Rydberg states with values of the orbital angular momentum (ℓ) ℓ ≳ 3 is a productive strategy to minimize this rapid decay and to populate molecular Rydberg states with lifetimes that approach those of atoms. Here, in this proof-of-principle demonstration, we transfer population to an nf Rydberg state of the calcium atom by stimulated Raman adiabatic passage, in which an optical and a millimeter-wave field couple the initial and final states via an intermediate nd Rydberg state. Numerical simulations reproduce the observed time and frequency dependences of the population transfer and suggest the utility of this scheme to populate high-ℓ Rydberg states of molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Convergence acceleration of Monte Carlo many-body perturbation methods by using many control variates

Herein, the use of many control variates is proposed as a method to accelerate the second- and third-order Monte Carlo (MC) many-body perturbation (MC-MP2 and MC-MP3) calculations. A control variate is an exactly integrable function that is strongly correlated or anti-correlated with the target function to be integrated by the MC method. Evaluating both integrals as well as their covariances in the same MC run, one can effect a mutual cancellation of the statistical uncertainties and biases in the MC integrations, thereby accelerating its convergence considerably. Six and thirty-six control variates, whose integrals are known a priori, are generated for MC-MP2 and MC-MP3, respectively, by systematically replacing one or more two-electron-integral vertexes of certain configurations by zero-valued overlap-integral vertexes in their Goldstone diagrams. The variances and co- variances of these control variates are computed at a marginal cost, enhancing the overall efficiency of the MC-MP2 and MC-MP3 calculations by a factor of up to 14 and 20, respectively.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Correlated atomic dynamics in liquid seen in real space and time

In liquids, the timescales for structure, diffusion, and phonon are all similar, of the order of a pico-second. This not only makes characterization of liquid dynamics difficult but also renders it highly questionable to describe liquids in these terms. In particular, the current definition of the structure of liquids by the instantaneous structure may need to be expanded because the liquid structure is inherently dynamic. Here, we advocate describing the liquid structure through the distinct-part of the Van Hove function, which can be determined by inelastic neutron and x-ray scattering measurements as well as by simulation. It depicts the dynamic correlation between atoms in space and time, starting with the instantaneous correlation function at t = 0. The observed Van Hove functions show that the atomic dynamics is strongly correlated in some liquids, such as water. The effect of atomic correlation on various transport properties of fluid, including viscosity and diffusivity, is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Consequences of sp 2 –sp 3 boron isomerization in supercooled liquid borates

In this study, time-resolved high-energy synchrotron x-ray total scattering measurements on supercooled molten lithium metaborate (LiBO 2 ) reveal an isomerization reaction involving conversion of trigonal sp 2 boron to tetrahedral sp 3 boron during quenching and glass formation. Van't Hoff analysis yields an accurate enthalpy change, ΔH = 21(1) kJ mol –1 boron, from which we develop an analytical model for the sp 3 isomer fraction and its contribution to configurational heat capacity (C p conf ) and entropy as a function of temperature and composition. Isomerization constitutes 40% of the total calorimetric C p conf at the glass transition for LiBO 2 and directly contributes to the observed rise in liquid fragility with the lithium content.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Accuracy of XAS theory for unraveling structural changes of adsorbates: CO on Ni(100)

Studying surface reactions using ultrafast optical pump and X-ray probe experiments relies on accurate calculations of X-ray spectra of adsorbates for correct identification of the spectral signatures and their dynamical evolution. We show that experimental XAS can be well reproduced for the different binding sites in a static prototype system CO/Ni(100) at a standard DFT (GGA) level of theory, using a plane-wave basis and pseudopotentials. This validates its utility in analyzing ultrafast X-ray probe experiments. The accuracy of computed relative core level binding energies is about 0.2 eV, representing a lower limit for which spectral features can be resolved with this method. We also show that the commonly used Z +1 approximation gives very good core binding energy shifts overall. However, we find a discrepancy for CO adsorbed in the hollow site, which we assign to the significantly stronger hybridization in hollow bonding than in on-top.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Prospects for sub-nanometer scale imaging of optical phenomena using electron microscopy

Imaging of optical phenomena at the sub-nanometer scale can offer fundamental insights into the electronic or vibrational states in atomic-scale defects, molecules, and nanoparticles, which are important in quantum information, heterogeneous catalysis, optoelectronics, and structural biology. Several techniques have surpassed the traditional Abbe diffraction limit and attained spatial resolutions down to a few nanometers, but sub-nanometer scale optics has remained elusive. Here, we propose an approach that combines spectrally specific photoabsorption with sub-nanometer scale resolution transmission electron microscopy (TEM) of photoexcited electrons. We first estimate the signal level and conditions required for imaging nanoscale optical phenomena in core-shell quantum dots (QDs) like CdS/CdTe. Furthermore, we show the possibility of imaging photoexcited states of atomic-scale defects in a monolayer hexagonal boron nitride (h-BN) using ab initio and high resolution (HR)TEM simulations. The ability to directly visualize photoexcited states at the sub-nanometer scale opens opportunities to study properties of individual quantum dots and atomic defects.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Tuning spin–orbit coupling in (6,5) single-walled carbon nanotube doped with sp 3 defects

Single-walled carbon nanotubes (SWCNTs) containing sp 3 defects are a promising class of optoelectronic materials with bright photoluminescence and demonstrated single-photon emission. Using density functional theory simulations, complemented by measurements, we investigate the electronic structure of a series of quantum defects attached to (6,5) SWCNT with the goal of tuning the spin–orbit coupling by introduction of a heavy atom in the defect structure. We characterize the ground state electronic and spin properties of four synthesized and three potential defects on the tube and find that all of the synthesized defects considered introduce a localized midgap defect-centered state containing a single electron, ≈0.2–0.3 eV above the valence band. The spin density is located at the sp 3 defect site with negligible spin–orbit coupling even with the presence of a Pd atom. Three additional functional groups were tested via computation to increase spin localization near the metal, thereby increasing spin–orbit coupling. We predict that only the chlorodiphosphanepalladium(II)– [Cl(PH 3 ) 2 Pd(II)–] defect results in increased spin–orbit splitting of the defect state and the conduction band associated with the pristine-like SWCNT, a measure of the spin–orbit coupling of excited state transitions. This study suggests that for unpassivated sp 3 defects in (6,5) SWCNT, forming a direct bond between a heavy atom and the sp 3 carbon allows for tuning of spin–orbit coupling.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Phase noise limited frequency shift impulsive Raman spectroscopy

We introduce a method to enable optical amplification of a coherent Raman spectroscopy signal, which we call radio frequency (RF) Doppler Raman spectroscopy. In this article, we consider the perturbation of a probe pulse in a sample due to an excited Raman vibrational coherence as a generalized Doppler shift, which connects a time-varying optical path length (the product of the propagation length and refractive index, OPL = n ℓ) with an optical frequency shift. Amplification of a Raman signal outside of the focused interaction is enabled by converting the Doppler frequency shift experienced by a laser probe pulse into a periodic timing jitter. This transit time perturbation is detected through the phase of a RF electronic signal measured at a harmonic of the probe pulse train with a method adapted from precision metrology techniques used to measure laser pulse train timing jitter. Measurement of a timing jitter allows access to much lower noise floors than other coherent Raman techniques, and by exploiting the new capability to scale the signal of a coherent Raman spectroscopic signal, this method opens the potential to detect very weak Raman signals that are currently not observable due to limits of illumination intensity imposed by laser damage to the specimen and noise.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

A comparison of Raman and pyrometry dynamic temperature measurements of shocked cyclohexane

Optical pyrometry is a well-developed, broadly applicable method to determine temperature in dynamic compression experiments. However, it measures radiation from only the skin depth of the radiant material, is highly susceptible to spurious backgrounds, and requires an often-unmeasurable free parameter—the sample dynamic emissivity. Raman spectroscopy offers a method to determine temperature directly from fundamental theory, but it is practically limited in applications to a small subset of shock experiments. In this work, we describe an experiment allowing simultaneous application of Raman and pyrometric temperature techniques to benchmark a specific instance of pyrometry. We attempted both measurements on multiple experiments and compare the successful temperature results between the separate but highly repeatable shots. Here, we find that dynamic shock temperatures inferred from Raman spectroscopy and pyrometry generally agree within the uncertainties of each measurement, although those uncertainties are larger than typical for a specially designed pyrometry experiment on a shocked metal sample.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Static dipole polarizabilities of polyacenes using self-interaction-corrected density functional approximations

Density functional approximations are known to significantly overestimate the polarizabilities of long chain-like molecules. We study the static electric dipole polarizabilities and the vertical ionization potentials of polyacenes from benzene to pentacene using the Fermi-L\"owdin orbital based self-interaction corrected (FLOSIC) density functional method. The orbital-by-orbital self-interaction correction corrects for the overestimation tendency of density functional approximations. The polarizabilities calculated with FLOSIC-DFA are however overly corrected. We also tested the recently developed locally-scaled self-interaction correction (LSIC) method on the polyacenes. The local-scaling method applies full SIC in the one-electron regions and restores the proper behavior of the SIC exchange-correlation functionals in the uniform density limit. The results show that LSIC removes the overcorrection tendency of the FLOSIC-DFA and produces results that are in excellent agreement with reference CCSD values. The vertical ionization potentials with LSIC also show good agreement with available experimental values.

74 ATOMIC AND MOLECULAR PHYSICS↗

Refine your search to explore more results.