SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Theoretical calculations predict the formation of phonon vortices at heavy impurities in two-dimensional materials. The vortex morphology depends on the impurity mass and the local symmetry of the defect.
In an one-pot in situ trans/esterification of low-grade feedstocks, the simultaneous presence of triglycerides (TAGs) and free fatty acids (FFAs) presents a dual catalytic challenge. While TAGs undergo transesterification to form biodiesel and glycerol, FFAs present in the feedstocks are esterified simultaneously, producing water as a by-product. This water can deactivate acid catalysts by interacting with their active sites, reducing catalytic efficiency and reusability. To overcome this, we developed a hydrophobic sulfonic acid-functionalized biochar catalyst derived from Citrus sinensis (orange peel) via sulfonation with H 2 SO 4 and subsequent silylation using hexamethyldisilazane (HMDS). The optimized catalyst exhibited a high surface area (355.09 m 2 g −1 ), sulfur content (5.33 wt%), and strong hydrophobicity (water contact angle: 154°). Compared to non-hydrophobic analogs, it showed enhanced activity and reusability (up to 10 cycles). Using response surface methodology with central composite design, a biodiesel yield of 99.1 ± 0.4% was achieved under optimal conditions. Life cycle assessment was performed to evaluate the environmental impacts of biodiesel production utilising the synthesized catalyst, considering 1000 kg of biodiesel produced as 1 functional unit. In conclusion, the recorded results showed the cumulative abiotic depletion of fossil resources over the entire biodiesel production process as 87 243.423 MJ and global warming potential as 4103.494 kg CO 2 equivalent.
In this work a 300 keV transmission electron microscope was modified to produce broadband pulsed beams that can be, in principle, between 40 MHz and 12 GHz, corresponding to temporal resolution in the nanosecond to picosecond range without an excitation laser. The key enabling technology is a pair of phase-matched modulating and de-modulating traveling wave metallic comb striplines (pulsers). An initial temporal resolution of 30 ps was achieved with a strobe frequency of 6.0 GHz. The placement of the pulsers, mounted immediately below the gun, allows for preservation of all optical configurations, otherwise available to the unmodified instrument, and therefore makes such a post-modified instrument for dual-use, i.e., both pulsed-beam mode (i.e., stroboscopic time-resolved) and conventional continuous waveform mode. In this article, we describe the elements inserted into the beam path, challenges encountered during integration with an in-service microscope, and early results from an electric-field-driven pump–probe experiment. We conclude with ideas for making this class of instruments broadly applicable for examining cyclical and repeatable phenomena.
Sandia’s Z Pulsed Power Facility is able to dynamically compress matter to extreme states with exceptional uniformity, duration, and size, which are ideal for investigating fundamental material properties of high energy density conditions. X-ray diffraction (XRD) is a key atomic scale probe since it provides direct observation of the compression and strain of the crystal lattice and is used to detect, identify, and quantify phase transitions. Because of the destructive nature of Z-Dynamic Material Property (DMP) experiments and low signal vs background emission levels of XRD, it is very challenging to detect a diffraction signal close to the Z-DMP load and to recover the data. We have developed a new Spherical Crystal Diffraction Imager (SCDI) diagnostic to relay and image the diffracted x-ray pattern away from the load debris field. The SCDI diagnostic utilizes the Z-Beamlet laser to generate 6.2-keV Mn–He α x rays to probe a shock-compressed material on the Z-DMP load. Finally, a spherically bent crystal composed of highly oriented pyrolytic graphite is used to collect and focus the diffracted x rays into a 1-in. thick tungsten housing, where an image plate is used to record the data.
Here, we describe recent efforts to improve our predictive modeling of rate-dependent behavior at, or near, a phase transition using molecular dynamics simulations. Cadmium sulfide (CdS) is a well-studied material that undergoes a solid-solid phase transition from wurtzite to rock salt structures between 3 and 9 GPa. Atomistic simulations are used to investigate the dominant transition mechanisms as a function of orientation, size and rate. We found that the final rock salt orientations were determined relative to the initial wurtzite orientation, and that these orientations were different for the two orientations and two pressure regimes studied. The CdS solid-solid phase transition is studied, for both a bulk single crystal and for polymer-encapsulated spherical nanoparticles of various sizes.
Pulsed power accelerators compress electrical energy in space and time to provide versatile experimental platforms for high energy density and inertial confinement fusion science. The 80-TW “Z” pulsed power facility at Sandia National Laboratories is the largest pulsed power device in the world today. Z discharges up to 22 MJ of energy stored in its capacitor banks into a current pulse that rises in 100 ns and peaks at a current as high as 30 MA in low-inductance cylindrical targets. Considerable progress has been made over the past 15 years in the use of pulsed power as a precision scientific tool. This paper reviews developments at Sandia in inertial confinement fusion, dynamic materials science, x-ray radiation science, and pulsed power engineering, with an emphasis on progress since a previous review of research on Z in Physics of Plasmas in 2005.
Pulsed-power generators can produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies using the magnetic loading technique. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength–strain rate dependence, structural phase transitions, and density of crystal defects, such as dislocations. Here, we present a cost-effective, compact, pulsed x-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically driven ramp compression of materials with a single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials’ community to investigate in situ dynamic phase transitions critical to equation of states. Finally, we present results using this new diagnostic to evaluate lattice compression in Zr and Al and to capture signatures of phase transitions in CdS.
Alkali antimonide photocathodes have wide applications in free-electron lasers and electron cooling. The short lifetime of alkali antimonide photocathodes necessitates frequent replacement of the photocathodes during a beam operation. Furthermore, exposure to mediocre vacuum causes loss of photocathode quantum efficiency due to the chemical reaction with residual gas molecules. Theoretical analyses have shown that covering an alkali antimonide photocathode with a monolayer graphene or hexagonal boron nitride protects it in a coarse vacuum environment due to the inhibition of chemical reactions with residual gas molecules. Alkali antimonide photocathodes require an ultra-high vacuum environment, and depositing a monolayer 2D material on it poses a serious challenge. In the present work, we have incorporated a novel method known as intercalation, in which alkali atoms pass through the defects of a graphene thin film to create a photocathode material underneath. Initially, Sb was deposited on a Si substrate, and a monolayer graphene was transferred on top of the Sb film. Heat cleaning around 550–600 °C effectively removed the Sb oxides, leaving metallic Sb underneath the graphene layer. Depositing Cs on top of a monolayer graphene enabled the intercalation process. Atomic force microscopy, Raman spectroscopy, x-ray photoelectron spectroscopy, low energy electron microscopy, and x-ray diffraction measurements were performed to evaluate photocathode formation underneath the monolayer graphene. Our analysis shows that Cs penetrated the graphene and reacted with Sb and formed Cs 3 Sb.
Hydrocarbon polymers are used in a wide variety of practical applications. In the field of dynamic compression at extreme pressures, these polymers are used at several high energy density (HED) experimental facilities. One of the most common polymers is poly(methyl methacrylate) or PMMA, also called Plexiglass® or Lucite®. Here, we present high-fidelity, hundreds of GPa range experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for PMMA to more than threefold compression up to 650 GPa and re-shock Hugoniot states up to 1020 GPa in an off-Hugoniot regime, where experimental data are even sparser. These data can be used to put additional constraints on tabular equation of state (EOS) models. The present results provide clear evidence for the need to re-examine the existing tabular EOS models for PMMA above ~120 GPa as well as perhaps revisit EOSs of similar hydrocarbon polymers commonly used in HED experiments investigating dynamic compression, hydrodynamics, or inertial confinement fusion.
Photonic Doppler velocimetry (PDV) is an established technique for measuring the velocities of fast-moving surfaces in high-energy-density experiments. In the standard approach to PDV analysis, the short-time Fourier transform (STFT) is used to generate a spectrogram from which the velocity history of the target is inferred. The user chooses the form, duration, and separation of the window function. Here, in this study, we present a Bayesian approach to infer the velocity directly from the PDV oscilloscope trace, without using the spectrogram for analysis. This is clearly a difficult inference problem due to the highly periodic nature of the data, but we find that with carefully chosen prior distributions for the model parameters, we can accurately recover the injected velocity from synthetic data. We validate this method using PDV data collected at the STAR two-stage light gas gun at Sandia National Laboratories, recovering shock-front velocities in quartz that are consistent with those inferred using the STFT-based approach and are interpolated across regions of low signal-to-noise data. Although this method does not rely on the same user choices as the STFT, we caution that it can be prone to misspecification if the chosen model is not sufficient to capture the velocity behavior. Analysis using posterior predictive checks can be used to establish whether a better model is required, although more complex models come with additional computational cost, often taking more than several hours to converge when sampling the Bayesian posterior. We, therefore, recommend it be viewed as a complementary method to that of the STFT-based approach.
An important component of dynamic compression x-ray diffraction (XRD) experiment analysis is geometry calibration: proper data interpretation requires knowledge of the precise detector position and orientation and, if the experiment involves a single-crystal sample, knowledge of the lattice orientation. The determination of these parameters in the arbitrary three-dimensional (3D) scattering geometries often present in dynamic compression facilities is challenging, as the associated optimization problem can be highly nonlinear, nonsmooth, and discontinuous. We present a genetic algorithm-based approach for performing dynamic compression XRD calibrations that overcomes these obstacles. We provide details regarding the image processing, algorithm implementation, and open-source software deployment and demonstrate the capability of the approach to calibrate the detector and crystal parameters in 3D geometries. Notably, we demonstrate the solver’s capacity to find the crystal orientation without a priori rotation constraints.
Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.
We present pressure amplifiers that have been developed to enable the study of material behavior at terapascal (TPa) pressures using two-stage light-gas gun (2SLGG) facilities. When impacted by a hyper-velocity projectile, Mach reflections within the pressure amplifier eventually merge to form a Mach stem, generating a planar output shock that greatly exceeds the initial drive. Validation experiments demonstrated that peak pressures of 0.6 and 1.073 TPa were achieved in quartz samples for two different amplifier designs. These pressures represent an approximately fivefold increase over the sample pressures typically attainable at 2SLGG facilities. The temporal and spatial uniformity of the pressure drive is suitable for high-precision equation-of-state measurements of materials relative to a shock standard. This work greatly extends the capabilities of 2SLGG facilities and opens new avenues of study in high-energy-density physics using these drivers.
An independent set (IS) is a set of vertices in a graph such that no edge connects any two vertices. In adiabatic quantum computation [E. Farhi, et al ., Science 292, 472–475 (2001); A. Das, B. K. Chakrabarti, Rev. Mod. Phys. 80, 1061–1081 (2008)], a given graph G ( V , E ) can be naturally mapped onto a many-body Hamiltonian H IS G ( V , E ) , with edges E being the two-body interactions between adjacent vertices V . Thus, solving the IS problem is equivalent to finding all the computational basis ground states of H IS G ( V , E ) . Very recently, non-Abelian adiabatic mixing (NAAM) has been proposed to address this task, exploiting an emergent non-Abelian gauge symmetry of H IS G ( V , E ) [B. Wu, H. Yu, F. Wilczek, Phys. Rev. A 101, 012318 (2020)]. Here, we solve a representative IS problem G ( 8 , 7 ) by simulating the NAAM digitally using a linear optical quantum network, consisting of three C-Phase gates, four deterministic two-qubit gate arrays (DGA), and ten single rotation gates. The maximum IS has been successfully identified with sufficient Trotterization steps and a carefully chosen evolution path. Remarkably, we find IS with a total probability of 0.875(16), among which the nontrivial ones have a considerable weight of about 31.4%. Our experiment demonstrates the potential advantage of NAAM for solving IS-equivalent problems.
The primary goal of protein science is to understand how proteins function, which requires understanding the functional dynamics responsible for transitions between different functional structures of a protein. A central concept is the exact reaction coordinates that can determine the value of committor for any protein configuration, which provide the optimal description of functional dynamics. Despite intensive efforts, identifying the exact reaction coordinates (RCs) in complex molecules remains a formidable challenge. Using the recently developed generalized work functional, we report the discovery of the exact RCs for an important functional process—the flap opening of HIV-1 protease. Our results show that this process has six RCs, each one is a linear combination of ~240 backbone dihedrals, providing the precise definition of collectivity and cooperativity in the functional dynamics of a protein. Applying bias potentials along each RC can accelerate flap opening by 10 3 to 10 4 folds. The success in identifying the RCs of a protein with 198 residues represents a significant progress beyond that of the alanine dipeptide, currently the only other complex molecule for which the exact RCs for its conformational changes are known. Our results suggest that the generalized work functional (GWF) might be the fundamental operator of mechanics that controls protein dynamics.