pre-postprocess elastic constants 2021
One scripts generates input and submission script files for third party DFT software VASP to compute elastic constants. The other script postprocesses the results from VASP.
Blaschke, Daniel↗
SEARCH · Engineering Papers
Results for “Postprocessing”
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.
Z-Target Radiography Postprocessing With A Deep Convolution Neural Network
Analyzing X-ray radiographs is crucial for understanding target behavior in Inertial Confinement Fusion (ICF) and High Energy Density (HED) platforms. However, the density of Magneto Raleigh Taylor (MRT) bands and limitations of target materials often obscure relevant spike growth and density information. To address this issue, machine learning postprocessing techniques can be applied to remove darkened regions in radiography images. In this study, a novel method is presented for removing MRT darkened regions from z-target radiographs using a convolutional neural network (CNN). The CNN, consisting of six layers, treats the darkened regions as noise and employs a mixed loss function and end-to-end frameworks to suppress them while preserving sharpness. The six-layer architecture is designed to effectively learn features when provided with a larger volume of learning space. Each layer is optimized using a mixed loss function that combines a standard loss pixel approach with a multi-scaled structural similarity index loss, which considers luminance, contrast, and structure in local neighborhoods. This approach is particularly beneficial for capturing the stochastic structure of MRT limbs. Due to the limited availability of experimental data, training is conducted using synthetic target radiography from 3D Alegra simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ultrathin VO 2 Films on Functional Substrates
The metal–insulator transition (MIT) in vanadium dioxide (VO 2 ) thin films is strongly affected by grain size, thickness, and interfacial properties. Typically, a minimum thickness around 50 nm is required for VO 2 to exhibit a significant MIT when functional substrates like sapphire and silicon are used. Several works have shown that thin films below 20 nm, with up to 2–3 decades of change in the resistance across the MIT, can be achieved but require complex pre- or postprocessing of the samples. We show that predeposition substrate condition control facilitates the direct growth of VO 2 ultrathin 15 nm films, exhibiting a resistance change between 3 and 4 decades across the MIT. Our findings indicate that the interface between the film and the substrate is crucial in determining the initial growth layers and the structural evolution. With appropriate substrate surface treatment, the desired VO 2 MIT can be enhanced regardless of the substrate crystallographic orientation. Moreover, we propose a novel approach to obtain large resistance changes across the MIT in ultrathin VO 2 films by incorporating a predeposited 1.5 nm vanadium oxide buffer layer, thereby eliminating the need to use different materials or complex pre- or postprocessing of the samples. Here we also demonstrate that this method improves the transition of 25–50 nm VO 2 thin films on silicon substrates. Our study reveals a simple approach for direct growth of ultrathin VO 2 films exhibiting a significant MIT, which is commonly accepted unattainable over substrates of technological importance, such as sapphire and silicon.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Evaluating a quantum-classical quantum Monte Carlo algorithm with Matchgate shadows
Solving the electronic structure problem of molecules and solids to high accuracy is a major challenge in quantum chemistry and condensed matter physics. The rapid emergence and development of quantum computers offer a promising route to systematically tackle this problem. Recent work by [Huggins et al ., Nature (London) 603 , 416 (2022)] proposed a hybrid quantum-classical quantum Monte Carlo (QC-QMC) algorithm using Clifford shadows to determine the ground state of a Fermionic Hamiltonian. This approach displayed inherent noise resilience and the potential for improved accuracy compared to its purely classical counterpart. Nevertheless, the use of Clifford shadows introduces an exponentially scaling postprocessing cost. In this work, we investigate an improved QC-QMC scheme utilizing the recently developed Matchgate shadows technique [Commun. Math. Phys. 404 , 629 (2023)], which removes the aforementioned exponential bottleneck. We observe from experiments on quantum hardware that the use of Matchgate shadows in QC-QMC is inherently noise robust. We show that this noise resilience has a more subtle origin than in the case of Clifford shadows. Nevertheless, we find that classical postprocessing, while asymptotically efficient, requires hours of runtime on thousands of classical CPUs for even the smallest chemical systems, presenting a major challenge to the scalability of the algorithm.
Monte Carlo methods↗
$Z$-Pinch Interferometry Analysis With the Fourier-Based TNT Code
Here, we present the analysis of interferometry diagnostics with the user-friendly Talbot Numerical Tool (TNT), a Fourier-based postprocessing code that enables real-time assessment of plasma systems. TNT performance was explored with visible and infrared interferometry in pulsed-power-driven Z -pinch configurations to expand its capabilities beyond Talbot X-ray interferometry in the high-intensity laser environment. TNT enabled accurate electron density characterization of magnetically driven plasma flows and shocks through phase-retrieval methods that did not require data modification or masking. TNT demonstrated enhanced resolution, detecting below 4 % fringe shift, which corresponds to 8.7 × 10 15 cm –2 within 28 μ m, approaching the laser probing system limit. TNT was tested against a well-known interferometry analysis software, delivering an average resolving power nearly ten times better (~28 μ m versus ~ 210 μ m) when resolving plasma ablation features. TNT demonstrated higher sensitivity when probing sharp electron density gradients in supersonic shocks. A maximum electron areal density of 4.1 × 10 17 cm –2 was measured in the shocked plasma region, and a minimum electron density detection of ~ 1.0 × 10 15 cm –2 was achieved. When probing colliding plasma flows, the calculations of the effective adiabatic index and the associated errors were improved from γ* = 2.6 ± 1.6 – 1.4 ± 0.2 with TNT postprocessing, contributing valuable data for the interpretation of radiative transport. Additional applications of TNT in the characterization of pulsed-power plasmas and beyond are discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Robust Design Under Uncertainty in Quantum Error Mitigation
Error mitigation techniques are crucial to achieving near-term quantum advantage. Classical postprocessing of quantum computation outcomes is a popular approach for error mitigation, which includes methods, such as zero noise extrapolation, virtual distillation, and learning-based error mitigation. However, these techniques have limitations due to the propagation of uncertainty resulting from the finite shot number of a quantum measurement. In this work, we introduce general and unbiased methods for quantifying the uncertainty and error of error-mitigated observables based on the strategic sampling of error mitigation outcomes. We then extend our approach to demonstrate the optimization of performance and robustness of error mitigation under uncertainty. To illustrate our methods, we apply them to zero noise extrapolation and Clifford date regression in the ground state of the XY model simulated using depolarizing and International Business Machines Corporation (IBM) Toronto noise models, respectively. In particular, we optimize the choice of noise levels and the allocation of shots for zero noise extrapolation and the distribution of the training circuits for Clifford data regression. While our methods are readily applicable to any postprocessing-based error mitigation approach, in practice they must not be prohibitively expensive—even though they perform optimizations of the error mitigation hyperparameters requiring sampling of a statistical distribution of error mitigation outcomes. By leveraging surrogate-based optimization, we show that our methods can efficiently perform optimal design for a zero noise extrapolation implementation. We then further demonstrate the transferability of learned zero noise extrapolation hyperparameters to other similar circuits.
97 MATHEMATICS AND COMPUTING↗
PFLOTRAN modeling data and scripts associated with “Refining the Hydrogeologic Framework of a Large River Corridor Model Using Waterborne Transient Electromagnetics”
NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the publication “Refining the Hydrogeologic Framework of a Large River Corridor Model Using Waterborne Transient Electromagnetics” submitted to Water Resources Research (Terry et al. 2025). The data package contains the groundwater modeling dataset from PFLOTRAN software. It includes the python script for mesh generation, boundary condition setting, PFLOTRAN input deck formation and postprocessing. It couples groundwater flow and species transport for Hanford Reach river corridor and pipelines the model generation and processing. This model can be used to easily generate the model and analysis for Hanford site. It can also be adjusted to other hydrologic area with ease. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package consists of 6 folders: (1) “data” contains all necessary data as input and intermediate data for processing; (2) “mesh” contains all mesh related files to generate mesh in Hanford Reach river corridor; (3) “model_run” contains the generated script for PFLOTRAN modeling; (4) “notebooks” contains all the Python script to generate the model; (5) “output” contains all the output from the computation; (6) “postprocessing” contains the Python script to generate scientific figure for manuscript. All files are .csv (comma-separated values), .h5 (HDF5 format), .in (input files), .ipynb (Jupyter notebooks), .p (Python pickle), .png (images), .PNG (images), .py (Python scripts), .pyc (Python bytecode), .r (R scripts), .sh (shell scripts), .txt (text files), .vtu (3D mesh/visualization format), .xz (compressed archive), or .zip (compressed archive).
54 ENVIRONMENTAL SCIENCES↗
Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization
This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.
97 MATHEMATICS AND COMPUTING↗
Demonstrating the impact of ro-vibrationally excited H 2 on divertor detachment via population modelling
EDGE2D-EIRENE profiles characteristic for low recycling, high recycling and detached plasmas are postprocessed applying a vibrationally resolved collisional radiative (CR) model based on the Yacora solver to predict excited state densities and thus the molecular Lyman-, Werner- and Fulcher band emission intensities. These emission intensities are compared to predictions from the AMJUEL database standardly used in EIRENE. Lower Fulcher band and higher Lyman- and Werner band emission is predicted by the Yacora model in comparison to AMJUEL suggesting that inconsistent results from model based evaluation of measured spectra are expected depending on which model and molecular emission band is considered. The implications on detachment rate predictions are discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A Unified Workflow for Sensitivity-Based Kinetic Analysis in Microkinetic Models
Degrees of rate control (DRC), apparent activation energies, and apparent reaction orders are established local sensitivity diagnostics for interpreting microkinetic models, but applying them routinely to large mechanisms often requires substantial reaction-specific bookkeeping, perturbation design, and postprocessing. Here, in this study, we present a unified derivative-based workflow that evaluates these quantities from a single compiled reaction-network model and target-rate definition. For any user-provided microkinetic model, the workflow compiles the mechanism into stoichiometrically consistent mass-action rate equations, solves the surface dynamics, and uses automatic differentiation to compute sensitivities with respect to rate constants, temperature, and gas partial pressures. By combining their calculations in the same framework, the workflow clearly demonstrates the relationships between different DRCs and the apparent activation energy. Using existing examples of propylene partial oxidation and methane oxidation on Pd(100), we verify expected transient redistribution of rate control, distinguish net Campbell DRCs from one-sided directional sensitivities, and show how apparent activation energy can be reconstructed either from one-sided DRCs or from state-based DRCs while critical mechanistic insights are obtained consistently. In the methane oxidation case, a pathway-subset test further illustrates how a simplified mechanism preserves key kinetic signatures of a full model, showing the potential of our user-friendly tool for model construction beyond kinetic analysis.
36 MATERIALS SCIENCE↗
Self-Sensing Composites via an Embedded 3D-Printed PVDF-MoS 2 Nanosensor for Structural Health Monitoring
Carbon fiber (CF)-reinforced epoxy composites are widely used in vehicle applications, where early damage detection is crucial for reliability and safety. To address this need, we developed a self-sensing epoxy/CF composite by embedding a PVDF-MoS 2 nanosensor via an embedded 3D printing method. By harnessing the intrinsic curing kinetics of epoxy, we tailored its rheological properties to optimize the embedded printing process, enabling precise and reliable support for sensor filaments without compromising the composite’s structural and functional integrity. Through comprehensive rheological and kinetic analysis, we established a quantitative relationship among curing temperature, conversion rate, and resulting yield modulus─defining a narrow processing window essential for successful sensor integration. Specifically, we identified that an epoxy yield modulus range of 180–294 Pa and a conversion rate below 10% are critical to support the PVDF-MoS 2 filament architecture. Here, this embedded 3D printing method produces complex and multimaterial PVDF-MoS 2 sensors within an epoxy matrix with minimal deformation and reduced postprocessing, which is scalable and adaptable for industrial applications. Under cyclic loading, the embedded sensors exhibited stable signals under constant loads and increased voltage signals in response to crack formation (17–35% higher) and catastrophic failure (1 order of magnitude higher), effectively capturing structural changes in real time. This study demonstrates the potential of PVDF-MoS 2 nanocomposite sensor materials for real-time structural health monitoring in epoxy–CF composite systems, enabling early detection of defects and stress anomalies, significantly reducing the risk of unexpected failures, and enhancing structural reliability.
PVDF-MoS2 sensor↗
Tailoring Piezoelectricity of 3D Printing PVDF-MoS 2 Nanocomposite via In Situ Induced Shear Stress
3D printing of unique structures with tunable properties offers significant advantages in the fabrication of complex and customized electronic devices. Here, this study introduces a process-microstructure–property-guided manufacturing route to fabricate PVDF-2D MoS 2 piezoelectric nanocomposites with tunable piezoelectric properties without having a postprocess. We control PVDF’s microstructure through direct ink writing (DIW) 3D printing while tuning PVDF-MoS 2 interfacial strain by controlling rheology and 3D printing parameters, such as nozzle size and printing speed. Our approach demonstrates tunable piezoelectricity in PVDF-MoS 2 , achieving a 15-fold increase in the piezoelectric coefficient (d 33 ) at a printing-induced shear stress of 6685 Pa. This enhancement arises from the electrostatic interactions between PVDF and MoS 2 and the filler distribution and alignment caused by the in situ shear stress in 3D printing, as confirmed by XPS and Raman mapping analyses. Our findings advance the understanding of piezoelectric mechanisms in PVDF-based nanocomposites, laying the foundation for 3D printing of piezoelectric sensors in wearable device applications with enhanced performance and customization capabilities.
2D MoS2↗
3D Printing of Functional Hydrogel Devices for Screenings of Membrane Permeability and Selectivity
Developing a fundamental understanding of the effects of varying ligand chemistries on mass transport rates is key to designing membranes with solute-specific selectivity. While permeation cells offer a robust method to characterize membrane performance, they are limited to assessing a single membrane chemistry or salt solution per test. As a result, investigating the effects of varying ligand chemistries on membrane performance can be a tedious process, involving both the preparation of multiple samples and numerous, time-consuming permeation tests. This study uses digital light processing (DLP) 3D printing to fabricate a millifluidic flow-based permeation device made from a hydrogel active ester network that can be easily functionalized with ion-selective ligands. Without the need for bonding or assembly steps, ligands can be introduced and tested in the permeation device by simply injecting a small volume of a ligand solution. Various salt concentrations and molecular species can be cycled through a single device by switching the solution feeding into the salt reservoir, thereby reducing the number of samples needed for permeability and selectivity screenings. This research sets the groundwork for formulation development and postprocessing methods to 3D-print functional millifluidic devices capable of assessing solute selectivity in membranes and polymer adsorbents for aqueous separations. In this work, comparable salt permeability trends were observed with both 3D-printed devices and traditional assays. Devices were functionalized with an imidazole ligand to investigate salt permeability and selectivity of monovalent and divalent salts. Measurements showed increasing permeability for monovalent salts (NaCl) relative to divalent salts (MgCl 2 , CuCl 2 ) in functionalized membranes, with higher monovalent/divalent selectivity at increasing imidazole grafting densities. Here, the methods and findings described here represent a step toward developing higher-throughput methods with 3D-printed devices for screening the effects of ligand chemistry on mass transport rates in membrane materials.
36 MATERIALS SCIENCE↗
Dual-Wavelength Simultaneous Patterning of Degradable Thermoset Supports for One-Pot Embedded 3D Printing
Vat photopolymerization (VP) techniques have enabled the fabrication of complex geometries while balancing high precision and fast processing times. 3D printed objects are traditionally built layer-by-layer with newly cured layers being structurally supported by previous ones. Fabricating unsupported features such as overhangs and arches risks misalignment and sagging, limiting the range of accessible designs. To overcome this issue, support structures are fabricated along with the primary object as temporary scaffolds that provide stability and conserve print fidelity. For VP specifically, patterning dissolvable sacrificial supports is attractive to avoid manual removal after printing. In this study, we demonstrate a base-degradable thermoset to pattern print supports in a one-pot formulation along with the primary structural material. Efficient printing is enabled using a dual-wavelength negative imaging (DWNI) DLP printer that patterns the degradable thermoset with visible light and the permanent network with UV light, which are simultaneously projected using a single digital micromirror device (DMD). Printed objects undergo thermal postprocessing to enhance the final conversion of the primary material, after which thermoset supports are degraded in a basic, aqueous solution. This approach provides a robust method for the dual-wavelength patterning of sacrificial thermoset supports, broadening the range of accessible 3D printable materials and geometries.
3D printing↗
Learning earthquake ground motions via conditional generative modeling
Predicting high-fidelity ground motions for future earthquakes is crucial for seismic hazard assessment and infrastructure resilience. Conventional empirical simulations suffer from sparse sensor distribution and geographically localized earthquake locations, while physics-based methods are computationally intensive and require accurate representations of Earth structures and earthquake sources. We propose an artificial intelligence (AI) spectrogram generator, Conditional Generative Modeling for Ground Motion (CGM-GM). CGM-GM leverages earthquake magnitudes and geographic coordinates of earthquakes and sensors as inputs, when postprocessed with phase information, capturing spatially continuous Fourier amplitude spectra (FAS) as well as properties such as P and S arrivals, and waveform durations, without explicit physics constraints. This is achieved through a probabilistic autoencoder that extracts latent distributions in the time-frequency domain and variational sequential models for prior and posterior distributions. We evaluate the performance of CGM-GM using small-magnitude earthquake records from the San Francisco Bay Area, a region with high seismic risks. Here, we report that CGM-GM demonstrates potential for complementing physics-based simulations and non-ergodic empirical ground motion models, as well as shows promise in seismology and beyond.
geophysics↗
Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN
In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.
dielectric properties↗
Error-mitigated nonorthogonal quantum eigensolver via shadow tomography
We present a shadow-tomography-enhanced nonorthogonal quantum eigensolver (NOQE) for more efficient and accurate electronic structure calculations on near-term quantum devices. By integrating shadow tomography into the NOQE, the measurement cost scales linearly rather than quadratically with the number of reference states, while also reducing the required qubits and circuit depth by half. This approach enables extraction of all matrix elements via randomized measurements and classical postprocessing. We analyze its sample complexity and show that, for small systems, it remains constant in the high-precision regime, while for larger systems, it scales linearly with the system size. We further apply shadow-based error mitigation to suppress noise-induced bias without increasing quantum resources. Demonstrations on the hydrogen molecule in the strongly correlated regime achieve chemical accuracy under realistic noise, showing that our method is both resource-efficient and noise-resilient for practical quantum chemistry simulations in the near term.
quantum algorithms & computation↗
Revealing the Structure and Dynamics of Self-Generated Electric and Magnetic Fields Near Plasma Stagnation in Laser-Driven Hohlraums
By coupling newly developed triparticle charged particle radiography with radiography reconstruction algorithms and novel reconstruction postprocessing techniques, the spatial structure and time evolution of self-generated electric and magnetic fields in laser-driven vacuum hohlraums have been quantitatively revealed. Through high-fidelity data from a series of experiments, it is shown that late in the hohlraum evolution (after the end of laser drive) these fields are strongly correlated in both space and time, providing evidence that their evolution is primarily dominated by advection with the plasma flow. At these late times, plasma flow velocities inferred from both gross radiography analysis and field reconstructions (and corroborated with Thomson scattering measurements) indicate that the plasma is approaching stagnation near the hohlraum axis. Finally, these experiments provide not only new physical insight into spontaneously generated hohlraum fields, but also provide important spatially and temporally resolved information for future benchmarking of numerical codes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗