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.
Following publication of the original article, the author identified an error in Graphical Abstract and updated Acknowledgment section. In Graphical abstract, there is a type on the unit E B (eV) which has been updated with this correction.
Recent work on quaternary semiconductors Cu 2 BaSn(S,Se) 4 and Ag 2 BaSnSe 4 for photovoltaic and thermoelectric applications, respectively, has shown the promise of exploring the broader family of defect-resistant I 2 -II-IV-X 4 materials (where I, II, and IV refer to the formal oxidation state of the metal cations and X is a chalcogen anion) with tetrahedrally coordinated I/IV cations and larger II cations (i.e., Sr, Ba, Pb, and Eu) for optoelectronic and energy-related applications. Chemical dissimilarity among the II and I/IV atoms represents an important design motivation because it presents a barrier to antisite formation, which otherwise may act as electronically harmful defects. We herein show how all 31 experimentally reported I 2 -II-IV-X 4 examples (with large II cations and tetrahedrally coordinated smaller I/IV cations), which form within five crystal structure types, are structurally linked. Based on these structural similarities, we derive a set of tolerance factors that serve as descriptors for phase stability within this family. Despite common usage in the well-studied perovskite system, Shannon ionic radii are found to be insufficient for predicting metal–chalcogen bond lengths, pointing to the need for experimentally derived correction factors as part of an empirically driven learning approach to structure prediction. Here, we use the tolerance factors as a predictive tool and demonstrate that four new I 2 -II-IV-X 4 compounds, Ag 2 BaSiS 4 , Ag 2 PbSiS 4 , Cu 2 PbGeS 4 , and Cu 2 SrSiS 4 , can be synthesized in correctly predicted phases. One of these compounds, Ag 2 PbSiS 4 , shows potentially promising optoelectronic properties for photovoltaic applications.
Not Available
With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.
We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.
The performance of quantum error correcting (QEC) codes is often studied under the assumption of spatiotemporally uniform error rates. On the other hand, experimental implementations almost always produce heterogeneous error rates, in either space or time, as a result of effects such as imperfect fabrication and/or cosmic rays. It is therefore important to understand if and how their presence can affect the performance of QEC in qualitative ways. Here, in this work, we study the effects of nonuniform error rates in the representative examples of the 1D repetition code and the 2D toric code, focusing on when they have extended spatiotemporal correlations; these may arise, for instance, from rare events (such as cosmic rays) that temporarily elevate error rates over the entire code patch. These effects can be described in the corresponding statistical mechanics models for decoding, where long-range correlations in the error rates lead to extended rare regions of weaker coupling. For the 1D repetition code where the rare regions are linear, we find two distinct decodable phases: a conventional ordered phase in which logical failure rates decay exponentially with the code distance, and a rare-region dominated Griffiths phase in which failure rates are parametrically larger and decay as a stretched exponential. In particular, the latter phase is present when the error rates in the rare regions are above the bulk threshold. For the 2D toric code where the rare regions are planar, we find no decodable Griffiths phase: rare events which boost error rates above the bulk threshold lead to an asymptotic loss of threshold and failure to decode. Unpacking the failure mechanism implies that techniques for suppressing extended sequences of repeated rare events (which, without intervention, will be statistically present with high probability) will be crucial for QEC with the toric code.
In high-energy and astroparticle physics, event generators play an essential role, even in the simplest data analyses. As analysis techniques become more sophisticated, e.g. based on deep neural networks, their correct description of the observed event characteristics becomes even more important. Physical processes occurring in hadronic collisions are simulated within a Monte Carlo framework. A major challenge is the modeling of hadron dynamics at low momentum transfer, which includes the initial and final phases of every hadronic collision. QCD-inspired phenomenological models used for these phases cannot guarantee completeness or correctness over the full phase space. These models usually include parameters which must be tuned to suitable experimental data. Until now, event generators have been developed and tuned mainly on the basis of data from high-energy physics experiments at accelerators. The wealth of data available from the latest generation of astroparticle experiments has not yet been fully exploited, and in many cases is not satisfactorily described. Both kinds of data sets are complementary as astroparticle experiments provide sensitivity especially to hadrons produced nearly parallel to the collision axis and cover center-of-mass energies up to several hundred TeV, well beyond those reached at colliders so far. In this report, we provide an overview of state-of-the-art event generators and their tuning, including the most relevant inputs from high-energy accelerator and astroparticle experiments. We present a road map that shows, for the first time, how the unified tuning of event generators with accelerator-based and astroparticle data can be performed.
A distributed control system uses a central controller in Internet communication with a local controller to manage grid tie attachment with a battery to form an integrated battery energy storage system (BESS). The BESS is capable of charging or discharging the battery, as well as correcting grid phase with volt amp reactive (VAR) leading or lagging operation modes. Examples shown include simple BESS charging and discharging, BESS integrated with renewable energy sources (here photovoltaic), and direct current fast charge (DCFC) connections with an electric vehicle.
The Fermilab Linac delivers 400 MeV H- beam to the rest of the accelerator chain. Providing stable intensity, energy, and emittance is key since it directly affects downstream machines. To operate high current beam, accelerators must minimize uncontrolled particle loss; this can be accomplished by minimizing beam longitudinal emittance via RF parameter optimization. However, RF tuning is required daily since the resonance frequency of the accelerating cavities is affected by ambient temperature and humidity variations and thus drifts with time. In addition, the energy and phase space distribution of particles emerging from the ion source are subject to fluctuations. Such drift is not unique to Fermilab, but rather affects most laboratories. We are exploring machine learning (ML) algorithms for automated RF tuning for 2 objectives: optimization of Linac output energy and phase oscillation correction, with an emphasis on time-drift aware modeling that can account for conditions changing over time.
Semiempirical methods like density functional tight-binding (DFTB) allow extensive phase space sampling, making it possible to generate free energy surfaces of complex reactions in condensed-phase environments. Such a high efficiency often comes at the cost of reduced accuracy, which may be improved by developing a specific reaction parametrization (SRP) for the particular molecular system. Thiol–disulfide exchange is a nucleophilic substitution reaction that occurs in a large class of proteins. Its proper description requires a high-level ab initio method, while DFT-GAA and hybrid functionals were shown to be inadequate, and so is DFTB due to its DFT-GGA descent. We develop an SRP for thiol–disulfide exchange based on an artificial neural network (ANN) implementation in the DFTB+ software and compare its performance to that of a standard SRP approach applied to DFTB. Furthermore, as an application, we use both new DFTB-SRP as components of a QM/MM scheme to investigate thiol–disulfide exchange in two molecular complexes: a solvated model system and a blood protein. Demonstrating the strengths of the methodology, highly accurate free energy surfaces are generated at a low cost, as the augmentation of DFTB with an ANN only adds a small computational overhead.
We introduce a new theoretical framework based on Feynman diagrams to compute phase shifts in matter wave interferometry. The method allows for analytic computation of higher order quantum corrections, beyond the traditional semi-classical approximation. These additional terms depend on the finite size of the initial matter wavefunction and/or have higher order dependence on ℏ. We apply the method to compute the response of matter wave interferometers to power law potentials and potentials with an arbitrary spatial dependence. The analytic expressions are validated by comparing to numerical simulations, and estimates are provided for the scale of the quantum corrections to the phase shift response to the gravitational field of the earth, anharmonic trapping potentials, and gravitational fields from local proof masses. We also find that for certain experimentally feasible parameters, these corrections are large enough to be measured and could lead to systematic errors if they are not mitigated. We find that to first order in a spatially dependent potential, quantum corrections vanish when the initial matter wavepacket has spherical symmetry and the potential satisfies Laplace's equation. We anticipate these quantum corrections will be especially important for trapped matter wave interferometers and for free-space matter wave interferometers in the presence of proof masses. These interferometers are becoming increasingly sensitive tools for mobile inertial sensing, gravity surveying, tests of gravity and its interplay with quantum mechanics, and searches for dark energy.
Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.
Ptychography is a powerful imaging technique that is used in a variety of fields, including materials science, biology, and nanotechnology. However, the accuracy of the reconstructed ptychography image is highly dependent on the accuracy of the recorded probe positions which often contain errors. These errors are typically corrected jointly with phase retrieval through numerical optimization approaches. When the error accumulates along the scan path or when the error magnitude is large, these approaches may not converge with satisfactory result. We propose a fundamentally new approach for ptychography probe position prediction for data with large position errors, where a neural network is used to make single-shot phase retrieval on individual diffraction patterns, yielding the object image at each scan point. The pairwise offsets among these images are then found using a robust image registration method, and the results are combined to yield the complete scan path by constructing and solving a linear equation. We show that our method can achieve good position prediction accuracy for data with large and accumulating errors on the order of 10 2 pixels, a magnitude that often makes optimization-based algorithms fail to converge. For ptychography instruments without sophisticated position control equipment such as interferometers, our method is of significant practical potential.
Piezoresponse force microscopy (PFM) is a powerful characterization technique to readily image and manipulate the ferroelectric domains. PFM gives an insight into the strength of local piezoelectric coupling and polarization direction through PFM amplitude and phase, respectively. Converting measured arbitrary units into units of effective piezoelectric constant remains a challenge, and insufficient methods are often used. While most quantification efforts have been spent on quantifying the PFM amplitude signal, little attention has been given to the PFM phase, which is often arbitrarily adjusted to fit expectations. This is problematic when investigating materials with unknown or negative sign of the probed effective electrostrictive coefficient or strong frequency dispersion of electromechanical responses, because assumptions about the PFM phase cannot be reliably made. The PFM phase can, however, provide important information on the polarization orientation and the sign of the effective electrostrictive coefficient probed by PFM. Most notably, the orientation of the PFM hysteresis loop is determined by the PFM phase. Moreover, when presenting PFM data as a combined signal, the resulting response can be artificially lowered or asymmetric if the phase data have not been correctly processed. Here, we explain the PFM amplitude quantification process and demonstrate a path to identify the phase offset required to extract correct meaning from the PFM phase data. We explore different sources of phase offsets including the experimental setup, instrumental contributions, and data analysis. We discuss the physical working principles of PFM and develop a strategy to extract physical meaning from the PFM amplitude and phase.
Explore the source record for details and available documents.
Alternate stacking of a highly conducting metallic layer with a magnetic triangular layer found in delafossite PdCrO 2 provides an excellent platform for discovering intriguing correlated quantum phenomena. Thin film growth of delafossites may enable not only the tuning of the basic physical properties beyond what bulk materials can exhibit, but also the development of novel hybrid materials by interfacing with dissimilar materials, yet this has proven to be extremely challenging. Here, we report the epitaxial growth of metallic delafossite PdCrO 2 films by pulsed laser epitaxy (PLE). The fundamental role of the PLE growth conditions, epitaxial strain, and chemical and structural characteristics of the substrate is investigated by growing under various growth conditions and on various types of substrates. While strain plays a large role in improving the crystallinity, the direct growth of epitaxial PdCrO 2 films without impurity phases was not successful. We attribute this difficulty to both the chemical and structural dissimilarities with the substrate and volatile nature of the PdO sublayer, which make nucleation of the right phase difficult. This difficulty was overcome by growing CuCrO 2 buffer layers before PdCrO 2 films were grown. Unlike PdCrO 2 , CuCrO 2 films were readily grown with a relatively wide growth window. Only a monolayer thick buffer layer was sufficient to grow the correct PdCrO 2 phase. This result indicates that the epitaxy of Pd-based delafossites is extremely sensitive to the chemistry and structure of the interface, necessitating near perfect substrate materials. The resulting films are commensurately strained and show an antiferromagnetic transition at 40 K that persists down to as thin as 3.6 nm in thickness. This work provides key insights into advancing the epitaxial growth of the broader class of metallic delafossites for both studying the basic physical properties and developing new spintronic and computing devices.
The long-held conclusion that magnetic order in delta phase Pu makes the structure mechanically unstable in density functional theory calculations is incorrect if the magnetic order is allowed to be non-collinear. The non-collinear 3Q spin structure is intrinsically cubic and makes all structurally equivalent bonds equivalent in their bonding character. Applied to density functional theory calculations on delta phase Pu, the 3Q spin structure results in elastic constants and phonons with the correct symmetry. Finally, the calculated phonon dispersion agrees better with experiment than previous calculations, and the calculated thermal expansion shows the unique behavior seen experimentally: it is negative.