Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix”

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 937 records · Page 52

Coal-derived graphene foam and micron-sized silicon composite anodes for lithium-ion batteries

Silicon-based materials demonstrate significant potential as lithium-ion batteries (LIBs) anode, but their expansion and degradation present engineering design challenges for commercial application. In this study, a porous three-dimensional (3D) graphene and micron-sized silicon composite anode (Si@G foam) was synthesized using humic acid (HA) derived from coal as a graphene precursor. In-situ formation of graphene structure through reducing HA was confirmed by Raman spectra. The reduced HA (rHA) shows similar electrical conductivity compared to the commercial conductive carbon. SEM images depict a 3D skeleton of coal-derived graphene, with silicon particles distributed on the 3D graphene foam's internal surface. The Si@G composite-anode displays a good reversible capacity of ~656 mAh/g at a current density of 50 mA/g, as well as a high-rate capability of ~433 mAh/g at a current density of 800 mA/g, and outstanding cycling stability –89.8% capacity retention after 300 cycles, which is significantly higher than that of other foam structures. During lithiation and de-lithiation, the graphene foam serves as a matrix of electrical conductors and a volume expansion support for silicon. This 3D graphene network will be beneficial for developing advanced silicon-based anodes for high-performance LIBs.

25 ENERGY STORAGE↗

Flow annealed importance sampling bootstrap meets differentiable particle physics

High-energy physics requires the generation of large numbers of simulated data samples from complex but analytically tractable distributions called matrix elements. Surrogate models, such as normalizing flows, are gaining popularity for this task due to their computational efficiency. We adopt an approach based on flow annealed importance sampling bootstrap (FAB) that evaluates the differentiable target density during training and helps avoid the costly generation of training data in advance. We show that FAB reaches higher sampling efficiency with fewer target evaluations in high dimensions in comparison to other methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mechanical Properties of Permanent Foaming Fixatives for Deactivation and Decommissioning Activities

DOE-EM have identified an operational requirement for a fixative that can immobilize and/or encapsulate residual contamination in 3D void volumes (pipes, gloveboxes, waste containers, etc.) during D and D activities. Failure to safely and effectively immobilize residual contamination can: Put workers at risk, Contaminate the public and environment, Drive up operational costs. Commercial-off-the-shelf (COTS) polyurethane (PU) foams is one possible solution and is currently being investigated in collaboration with Savannah River National Lab (SRNL). Mechanical property testing will dictate how well PU foams can withstand physical stressors outlined in safety documents while immobilizing residual contamination. Results can provide a performance criteria for Section 5 of ASTM E3191: Standard Specification for Permanent Foaming Fixatives Used to Mitigate Spread of Radioactive Contamination. Polyurethanes consists of repeating units of urethanes which are produced from an exothermic reaction between a polyol (an alcohol with two or more hydroxyl groups in a molecule) and an isocyanate [2]. PUs consists of 3D cells that are formed by a blowing agent which is normally water. Water interacting with the isocyanate creates carbon dioxide and urea which contributes to cells' growth and expansion. The cross-linking between the polyol and isocyanate solidify the cell faces and walls. Micro scale properties that can influence mechanical properties include: Relative density (ratio between density of cellular material and density of the solid of which the cell walls are made of). Cells connectivity between edges and faces, and number of cells contact neighbors [3]. Whether a foam is open-celled or closed-celled. Closed-celled foams have a thin membrane that encloses each cell which can contribute to the stiffness of the foam. Incorporating fillers in the PU matrix like carbon-based nanoparticles in the form of expandable graphite (EG) can increase the thermal and mechanical performance [4]. When exposed to heat, EG produces a char layer that prevents heat and oxygen transfer which is known as intumescing. Mechanical performance is dependent on amount of EG loading since the cell sizes are linearly dependent on it. Overall, the rigid foams (I-R2) proved superior. SEM Analysis: I-R2 had the most consistent spherical sizes (average cell diameter: 489 μm), which enhance the relative density the most [5]. Tensile Testing Results: I-R2 experienced the largest breaking load (103.65 lbf) and stress (289.42 psi). I-R2 also had a higher Young's Modulus compared to R1 (2.99 ksi to 2.09 ksi). Compression Testing Results: I-R2 had the highest stress at 80% deflection (360.67 psi), load (1387.48 lbf) and modulus (3.92 ksi)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Charge Transport in Zirconium-Based Metal–Organic Frameworks

Metal–organic frameworks (MOFs) are a class of crystalline porous materials characterized by inorganic nodes and multitopic organic linkers. Because of their molecular-scale porosity and periodic intraframework chemical functionality, MOFs are attractive scaffolds for supporting and/or organizing catalysts, photocatalysts, chemical-sensing elements, small enzymes, and numerous other functional-property-imparting, nanometer-scale objects. Notably, these objects can be installed after the synthesis of the MOF, eliminating the need for chemical and thermal compatibility of the objects with the synthesis milieu. Thus, postsynthetically functionalized MOFs can present three-dimensional arrays of high-density, yet well-separated, active sites. Depending on the application and corresponding morphological requirements, MOF materials can be prepared in thin-film form, pelletized form, isolated single-crystal form, polycrystalline powder form, mixed-matrix membrane form, or other forms. For certain applications, most obviously catalytic hydrolysis and electro- or photocatalytic water splitting, but also many others, an additional requirement is water stability. MOFs featuring hexa-zirconium(IV)-oxy nodes satisfy this requirement. For applications involving electrocatalysis, charge storage, photoelectrochemical energy conversion, and chemiresistive sensing, a further requirement is electrical conductivity, as embodied in electron or hole transport. As most MOFs, under most conditions, are electrically insulating, imparting controllable charge-transport behavior is both a chemically intriguing and chemically compelling challenge. Herein, we describe three strategies to render zirconium-based metal–organic frameworks (MOFs) tunably electrically conductive and, therefore, capable of transporting charge on the few nanometers (i.e., several molecular units) to few micrometers (i.e., typical dimensions for MOF microcrystallites) scale. The first strategy centers on redox-hopping between periodically arranged, chemically equivalent sites, essentially repetitive electron (or hole) self-exchange. Zirconium nodes are electrically insulating, but they can function as grafting sites for (a) redox-active inorganic clusters or (b) molecular redox couples. Alternatively, charge hopping based on linker redox properties can be exploited. Marcus’s theory of electron transfer has proven useful for understanding/predicting trends in redox-hopping based conductivity, most notably, in accounting for variations as great as 3000-fold depending on the direction of charge propagation through structurally anisotropic MOFs. In MOF environments, propagation of electronic charge via redox hopping is necessarily accompanied by movement of charge-compensating ions. Consequently, rates of redox hopping can depend on both the identity and concentration of ions permeating the MOF. In the context of electrocatalysis, an important goal is to transport electronic charge fast enough to match or exceed the inherent activity of MOF-based or MOF-immobilized catalysts. Bandlike electronic conductivity is the focus of an alternative strategy: one based on the introduction of molecular guests capable of forming donor–acceptor charge transfer complexes with the host framework. Theory again can be applied predictively to alter conductivity. A third strategy similarly emphasizes electronic conductivity, but it makes use of added bridges in the form of molecular oligomers or inorganic clusters that can then be linked to span the length of a MOF crystallite. For all strategies, retention of molecular-scale porosity is emphasized, as this property is key to many applications. Lastly, while our focus is on Zr-MOFs, the described approaches clearly are extendable to other MOF compositions, as has already been demonstrated, in part, in studies by others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗

Layered Heterostructure Ionogel Electrolytes for High‐Performance Solid‐State Lithium‐Ion Batteries

Abstract Ionogel electrolytes based on ionic liquids and gelling matrices offer several advantages for solid‐state lithium‐ion batteries, including nonflammability, wide processing compatibility, and favorable electrochemical and thermal properties. However, the absence of ionic liquids that are concurrently stable at low and high potentials constrains the electrochemical windows of ionogel electrolytes and thus their high‐energy‐density applications. Here, ionogel electrolytes with a layered heterostructure are introduced, combining high‐potential (anodic stability: >5 V vs Li/Li + ) and low‐potential (cathodic stability: <0 V vs Li/Li + ) imidazolium ionic liquids in a hexagonal boron nitride nanoplatelet matrix. These layered heterostructure ionogel electrolytes lead to extended electrochemical windows, while preserving high ionic conductivity (>1 mS cm −1 at room temperature). Using the layered heterostructure ionogel electrolytes, full‐cell solid‐state lithium‐ion batteries with a nickel manganese cobalt oxide cathode and a graphite anode are demonstrated, exhibiting voltages that are unachievable with either the high‐potential or low‐potential ionic liquid alone. Compared to ionogel electrolytes based on mixed ionic liquids, the layered heterostructure ionogel electrolytes enable higher stability operation of full‐cell lithium‐ion batteries, resulting in significantly enhanced cycling performance.

Hyun, Woo Jin↗

Scalable, Ultrathin, and High-Temperature-Resistant Solid Polymer Electrolytes for Energy-Dense Lithium Metal Batteries

We report all-solid-state batteries (ASSBs) demonstrate great promise, offering high energy density, good thermal stability, and safe operation compared with traditional Li-ion batteries. Among various solid-state electrolytes (SSEs), solid polymer electrolytes (SPEs) offer an attractive choice due to their thinness, low density, and good manufacturability. However, ultrathin SPEs that work with practical current densities or at high temperatures remain challenging, limiting applicable conditions of SPE-based batteries. Here, the authors report a novel scalable, ultrathin, and high-temperature-resistant SPE for ASSBs. This design includes an electrospun polyacrylonitrile (PAN) matrix and polyethylene oxide (PEO)/Li salt ionic conductor, which offers a stable LiF and Li 3 N containing SSE/Li interface. The unique interface—as well as the good mechanical strength—inhibits lithium dendrites and prevents short circuiting. As a result, symmetrical Li-Li cells deliver more than 300 h cyclability at 0.5 mA cm -2 . ASSBs fabricated with only 5 µm-thickness PAN-PEO/lithium bis(trifluoromethanesulfonyl)imide reach 300 cycles at 0.3 C rate at 60 °C. The excellent thermal stability of PAN also results in safer SPEs at high temperatures. The design extends battery operation up to temperatures of 120 and 150 °C, where it achieves 500 cycles at C/2 rate and 100 cycles at 2C rate, respectively.

25 ENERGY STORAGE↗

Medium-Range Ordering in the Ionic Glass Electrolytes LiPON and LiSiPON

In this work, we provide an in-depth structural characterization of the amorphous ionic glasses LiPON and LiSiPON with high Li content. Based on ab initio molecular dynamics simulations, the structure of these materials is an inverted structure with either isolated polyanion tetrahedra or polyanion dimers in a Li+ matrix. Based on neutron scattering data, this type of inverted structure leads to a significant amount of medium-range ordering in the structure, as demonstrated by two sharp diffraction peaks and a periodic structural oscillation in the density function G(r). While this medium-range ordering is commonly observed in liquids and metallic glasses, it has not previously been observed in oxides. On a local scale, adding N and Si increases the number of anion bridges and polyanion dimer structures, leading to higher ionic conductivity. In the medium-range ordering, the addition of Si leads to more disorder in the polyanion substructure but a significant increase in the ordering of the O substructure. Finally, we demonstrate that this inverted structure with medium-range ordering results in a glassy material that is both mechanically stiff and ductile on the nanoscale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Speeding Up Hartree–Fock in JuliaChem with Density Fitting

In this work, the density fitting (DF) approximation is added to the restricted Hartree–Fock (RHF) implementation in the JuliaChem computational chemistry code. Utilizing a DF algorithm that uses symmetry and integral screening, a significant reduction in time to compute the Fock matrix is achieved. The symmetry and screening DF-RHF techniques were adapted to be performed on graphics processing units (GPUs), which are well suited to perform the matrix multiplications that comprise the bulk of the Fock build time in DF-RHF. The JuliaChem DF-RHF GPU algorithm employs a novel approach that automatically switches between two DF-RHF algorithms depending on the number of basis functions in the calculation. The JuliaChem GPU DF-RHF implementation demonstrates up to 2× speedup for Fock build times compared to the existing best-in-class GPU DF-RHF implementation by operating directly on screened intermediate matrices. Due to the high portability of the Julia language code, the JuliaChem CPU and GPU DF-RHF implementations could be benchmarked on a variety of CPU and GPU architectures from multiple hardware vendors.

Hayes, John J. [Ames Laboratory, and Iowa State Un↗

Dynamics of Polymerization and Gelation in Epoxy Nanocomposites via X-ray Photon Correlation Spectroscopy

The details of the curing process of epoxy resins are notoriously difficult to ascertain, hampering polymer matrix composite manufacturing. Here, the curing of a series of six epoxy resins containing dilute spherical nanoparticles was investigated via X-ray photon correlation spectroscopy (XPCS). The resin formulation was varied to achieve a range of crosslink densities and topologies. Prior to gelation, the viscosities obtained from XPCS agreed with bulk rheology. A clear dynamic transition was detected in all samples, where the q-scaling changes from approximately q2 .2 to q 1 , indicating a transition from sub-diffusive to ballistic-like nanoparticle motion. For the samples with a high crosslink density, this change was abrupt and coincided with gelation, and there was no evidence of crosslink heterogeneity. Chain entanglement played a role in the transition for samples with lower crosslink density. Finally, this work reveals the potential of XPCS to deliver new quantitative insights into the curing process of thermosetting resins.

36 MATERIALS SCIENCE↗

van der Waals driven anharmonic melting of the 3D charge density wave in VSe 2

Understanding of charge-density wave (CDW) phases is a main challenge in condensed matter due to their presence in high-T c superconductors or transition metal dichalcogenides (TMDs). Among TMDs, the origin of the CDW in VSe 2 remains highly debated. Here, by means of inelastic x-ray scattering and first-principles calculations, we show that the CDW transition is driven by the collapse at 110 K of an acoustic mode at q CDW = (2.25 0 0.7) r.l.u. The softening starts below 225 K and expands over a wide region of the Brillouin zone, identifying the electron-phonon interaction as the driving force of the CDW. This is supported by our calculations that determine a large momentum-dependence of the electron-phonon matrix-elements that peak at the CDW wave vector. Our first-principles anharmonic calculations reproduce the temperature dependence of the soft mode and the T CDW onset only when considering the out-of-plane van der Waals interactions, which reveal crucial for the melting of the CDW phase. The nature of the charge density wave transition in VSe 2 is still debated. Here, the authors demonstrate that the transition is mainly driven by electron-phonon interactions, despite the presence of the Fermi-surface nesting, and that Wan-der-Waals forces are responsible for melting of the charge density wave order.

74 ATOMIC AND MOLECULAR PHYSICS↗

The effective number of parameters in kernel density estimation

We devise a new formula for measuring the effective degrees of freedom (EDoF) in kernel density estimation (KDE). Starting from the orthogonal polynomial sequence (OPS) expansion for the ratio of the empirical to the oracle density, we show how convolution with the kernel leads to a new OPS with respect to which one may express the resulting KDE. The expansion coefficients of the two OPS systems can then be related via a kernel sensitivity matrix, which leads to a natural oracle definition of EDoF through the trace operator. Asymptotic properties of the (empirical) plug-in EDoF are worked out through influence functions, and connections with other empirical EDoFs are established. Minimization of Kullback-Leibler divergence is investigated as an alternative to integrated squared error based bandwidth selection rules, yielding a new normal scale rule. The methodology, which arises from a proper oracle formulation and is not restricted to convolution kernels, suggests the possibility of a new bandwidth selection rule based on an information criterion such as AIC.

bandwidth selection↗

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extracting the electronic structure of light elements in bulk materials through a Compton scattering method in the readily accessible hard x-ray regime

Here, our Compton profile measurements of Ti and TiH 2 using readily available hard X-ray radiation at 27.5 keV, detected by both a Hitachi Vortex silicon-drift detector and a high-resolution superconducting transition-edge sensor array, are found to be in excellent accord with state-of-the-art density functional theory based calculations. The spherically averaged difference between the Compton profiles of TiH 2 and Ti is well described by an inverted parabola, supporting an itinerant behavior of the electron gas screening the protons in the Ti matrix. Our experimental approach, validated by two different detectors, extends the applicability of Compton scattering technique to the readily accessible hard x-ray regime (below 30 keV). Our study suggests possibilities for experiments at low-flux bending magnet synchrotron beamlines and paves the way for the development of tabletop Compton experiments with x-ray tubes.

36 MATERIALS SCIENCE↗

A differentiable perturbation-based weak lensing shear estimator

Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-per cent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (fpfs), using jax, which is publicly available as part of a new suite of analytic shear algorithms called anacal. This code can analytically calibrate the shear response of any non-linear observable constructed with the fpfs shapelets and detection modes utilizing autodifferentiation (ad), generalizing the formalism to include a family of shear estimators with corrections for detection and selection biases. Using the ad capability of jax, it calculates the full Hessian matrix of the non-linear observables, which improves the previously presented second-order noise bias correction in the shear estimation. As an illustration of the power of the new anacal framework, we optimize the effective galaxy number density in the space of the generalized shear estimators using an LSST-like galaxy image simulation for the 10 yr LSST. For the generic shear estimator, the magnitude of the multiplicative bias |m| is below 3 × 10 –3 (99.7 per cent confidence interval), and the effective galaxy number density is improved by 5 per cent. We also discuss some planned future additions to the anacal software suite to extend its applicability beyond the fpfs measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

First-principles study of electron transport in ScN

We investigate the conduction-band structure and electron mobility in rocksalt ScN based on density functional theory. The first-principles band structure allows us to obtain band velocities and effective masses as a function of energy. Electron-phonon scattering is assessed by explicitly computing the q-dependent electron-phonon matrix elements, with the inclusion of the long-range electrostatic interaction. The influence of free-carrier screening on the electron transport is assessed using the random phase approximation. We find a notable enhancement of electron mobility when the carrier concentration exceeds 10 20 cm -3 . We calculate the room-temperature electron mobility in ScN to be 587 cm2 /Vs at low carrier concentrations. When the carrier concentration is increased, the electron mobility starts to decrease significantly around n = 10 19 cm -3 , and drops to 240 cm 2 /Vs at n = 10 21 cm -3 . We also explore the influence of strain in (111)- and (100)-oriented ScN films. For (111) films, we find that a 1.0% compressive epitaxial strain increases the in-plane mobility by 72 cm 2 /Vs and the out-of-plane mobility by 50 cm 2 /Vs. For (100) films, a 1.0% compressive epitaxial strain increases the out-of-plane mobility by as much as 172 cm2 /Vs, but has a weak impact on the in-plane mobility. Our study sheds light on electron transport in ScN at different electron concentrations and shows how strain engineering could increase the electron mobility.

36 MATERIALS SCIENCE↗

Enhanced Hanford High-Aluminum Waste Glass Property Data Development

This study investigated the effects of aluminum concentration on simulated high level waste glass properties to eventually establish an aluminum limit (as single-component or multiple-component constraints) for glass formulations for high-aluminum Hanford wastes. A test matrix of 25 high-aluminum glasses (20 = Al 2 O 3 = 28.57 wt%) was generated, and the chemical compositions were measured. The following properties were measured: crystal formation after centerline canister cooling, crystallinity as a function of temperature, density, viscosity, electrical conductivity, toxic leaching characteristics using the Toxicity Characteristic Leaching Procedure, product consistency using the Product Consistency Test, and SO 3 solubility. These results are reported here.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗