Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster characterization”

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

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

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DFT Analysis of the Binding of Rare Earth Nitrates at Internal and External Surfaces of MCM-22

MCM-22-type zeolites constitute a well-characterized tunable class of aluminosilicates suitable for elucidating the fundamental aspects of the binding of rare earth elements (REEs) in layered materials. Starting from the pure silica version ITQ-1 with a unit cell of Si 72 O 144 , a model for periodic bulk crystalline MCM-22 with a finite Al concentration is provided by replacing a Si atom with an Al atom in the unit cell at suitable tetrahedral sites near an internal pore surface. Then, a H atom is added to an O atom bridging Si and Al atoms to create a Brønsted acid site (BAS). There are no internal silanol groups in this bulk model. To generate a model for an external surface, we adopt the fully hydroxylated surface structure of a layer within the ITQ-1 precursor with two silanols per lateral unit cell. A BAS on the external surface can be generated by replacing a near-surface Si atom with an Al atom and adding a H atom, as above. The strength of binding at a BAS of REE, X, taken to be present in the solution phase as nitrates, is determined from the energy change in the reaction X(NO 3 ) 3 + ≡Si–{OH}–Al≡ → ≡Si–{OX(NO 3 ) 2 }–Al≡+ HNO 3 . The strength of binding at the silanols is determined similarly. Binding energies are determined from two approaches. The first performs periodic plane-wave density functional theory (DFT) total energy analysis for an entire unit cell of MCM-22. The second utilizes cluster models capturing the local environment of REE binding sites and performs DFT analysis with localized basis sets. The two approaches yield consistent results for Nd, revealing similarly strong binding at either an external or internal BAS, but much weaker binding at a silanol site. This is consistent with the picture deduced from recent experiments. Here, we also comment on binding at Al-bridged siloxane sites, which have been suggested as alternative binding sites to BAS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
DOE OSTI2023

A full year of aerosol size distribution data from the central Arctic under an extreme positive Arctic Oscillation: insights from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition

Abstract. The Arctic environment is rapidly changing due to accelerated warming in the region. The warming trend is driving a decline in sea ice extent, which thereby enhances feedback loops in the surface energy budget in the Arctic. Arctic aerosols play an important role in the radiative balance and hence the climate response in the region, yet direct observations of aerosols over the Arctic Ocean are limited. In this study, we investigate the annual cycle in the aerosol particle number size distribution (PNSD), particle number concentration (PNC), and black carbon (BC) mass concentration in the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. This is the first continuous, year-long data set of aerosol PNSD ever collected over the sea ice in the central Arctic Ocean. We use a k-means cluster analysis, FLEXPART simulations, and inverse modeling to evaluate seasonal patterns and the influence of different source regions on the Arctic aerosol population. Furthermore, we compare the aerosol observations to land-based sites across the Arctic, using both long-term measurements and observations during the year of the MOSAiC expedition (2019–2020), to investigate interannual variability and to give context to the aerosol characteristics from within the central Arctic. Our analysis identifies that, overall, the central Arctic exhibits typical seasonal patterns of aerosols, including anthropogenic influence from Arctic haze in winter and secondary aerosol processes in summer. The seasonal pattern corresponds to the global radiation, surface air temperature, and timing of sea ice melting/freezing, which drive changes in transport patterns and secondary aerosol processes. In winter, the Norilsk region in Russia/Siberia was the dominant source of Arctic haze signals in the PNSD and BC observations, which contributed to higher accumulation-mode PNC and BC mass concentrations in the central Arctic than at land-based observatories. We also show that the wintertime Arctic Oscillation (AO) phenomenon, which was reported to achieve a record-breaking positive phase during January–March 2020, explains the unusual timing and magnitude of Arctic haze across the Arctic region compared to longer-term observations. In summer, the aerosol PNCs of the nucleation and Aitken modes are enhanced; however, concentrations were notably lower in the central Arctic over the ice pack than at land-based sites further south. The analysis presented herein provides a current snapshot of Arctic aerosol processes in an environment that is characterized by rapid changes, which will be crucial for improving climate model predictions, understanding linkages between different environmental processes, and investigating the impacts of climate change in future Arctic aerosol studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Predicting 3D moisture sorption behavior of materials from 1D investigations

Abstract Moisture in materials can be a source of future outgassing and exacerbate unwanted changes in physical and chemical properties. Here, we investigate the effect of sample size and shape on the moisture transport phenomena through a combined experimental and modeling approach. Several different materials varying in size and shape were investigated over a wide range of relative humidities (0–90%) and temperatures ( $$30{-}70^{\,\circ } \hbox {C}$$ 30 - 70 ∘ C ) using gravimetric type dynamic vapor sorption (DVS). A dynamic triple-mode sorption model, developed previously, was employed to describe the experimental results with good success; the model includes absorption, adsorption, pooling (clustering) of species, and molecular diffusion. Here we show that the full triple-mode sorption model is robust enough to predict the dynamic uptake and outgassing of 3-dimensional (3D) samples using parameters derived from quasi-1D samples. This successful demonstration on three different materials (filled polydimethylsiloxane (PDMS), unfilled PDMS, and ceramic inorganic composite) illustrates that the model is robust at describing the scale-independent physics and chemistry of moisture sorption and diffusion materials. This work demonstrates that while sorption mechanisms manifest in testing of all sample sizes, some of these mechanisms were so subtle that they were overlooked in our initial modeling and assessment, illustrating the importance of multi-scale experiments in the development of robust predictive capabilities. Our study also outlines the challenges and viable solutions for global optimization of a multi-parameter model. The ability to quantify moisture sorption and diffusion, independent of scale, using 1D lab-scale experiments enables prediction of long-term bulk materials behavior in real applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating correlated defects in metal organic frameworks using theory-guided inelastic neutron scattering spectroscopy

Metal organic frameworks (MOFs) that incorporate metal oxide cluster nodes, exemplified by UiO-66, have been widely studied, especially in terms of their deviations from the ideal, defect-free crystalline structures. Although defects such as missing linkers, missing nodes, and the presence of adventitious synthesis-derived node ligands (such as acetates and formates) have been proposed, their exact structures remain unknown. Previously, it was demonstrated that defects are correlated and span multiple unit cells. The highly specialized techniques used in these studies are not easily applicable to other MOFs. Thus, there is a need to develop new experimental and computational approaches to understand the structure and properties of defects in a wider variety of MOFs. Here, we show how low-frequency phonon modes measured by inelastic neutron scattering (INS) spectroscopy can be combined with density functional theory (DFT) simulations to provide unprecedented insights into the defect structure of UiO-66. We are able to identify and assign peaks in the fingerprint region (<100 cm –1 ) which correspond to phonon modes only present in certain defective topologies. Specifically, this analysis suggests that our sample of UiO-66 consists of predominantly defect-free fcu regions with smaller domains corresponding to a defective bcu topology with 4 and 2 acetate ligands bound to the Zr 6 O 8 nodes. Importantly, the INS/DFT approach provides detailed structural insights (e.g., relative positions and numbers of acetate ligands) that are not accessible with microscopy-based techniques. Furthermore, the quantitative agreement between DFT simulations and the experimental INS spectrum combined with the relative simplicity of sample preparation, suggests that this methodology may become part of the standard and preferred protocol for the characterization of MOFs, and, in particular, for elucidating the structure defects in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hydrogen–Deuterium Exchange Mass Spectrometry Identifies Local and Long-Distance Interactions within the Multicomponent Radical SAM Enzyme, PqqE

Interactions among proteins and peptides are essential for many biological activities including the tailoring of peptide substrates to produce natural products. The first step in the production of the bacterial redox cofactor pyrroloquinoline quinone (PQQ) from its peptide precursor is catalyzed by a radical SAM (rSAM) enzyme, PqqE. We describe the use of hydrogen–deuterium exchange mass spectrometry (HDX-MS) to characterize the structure and conformational dynamics in the protein–protein and protein–peptide complexes necessary for PqqE function. HDX-MS-identified hotspots can be discerned in binary and ternary complex structures composed of the peptide PqqA, the peptide-binding chaperone PqqD, and PqqE. Structural conclusions are supported by size-exclusion chromatography coupled to small-angle X-ray scattering (SEC-SAXS). HDX-MS further identifies reciprocal changes upon the binding of substrate peptide and S-adenosylmethionine (SAM) to the PqqE/PqqD complex: long-range conformational alterations have been detected upon the formation of a quaternary complex composed of PqqA/PqqD/PqqE and SAM, spanning nearly 40 Å, from the PqqA binding site in PqqD to the PqqE active site Fe 4 S 4 . Interactions among the various regions are concluded to arise from both direct contact and distal communication. The described experimental approach can be readily applied to the investigation of protein conformational communication among a large family of peptide-modifying rSAM enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Measurement of Kinetic-Scale Current Filamentation Dynamics and Associated Magnetic Fields in Interpenetrating Plasmas

Here, we present the first local, quantitative measurements of ion current filamentation and magnetic field amplification in interpenetrating plasmas, characterizing the dynamics of the ion Weibel instability. The interaction of a pair of laser-generated, counterpropagating, collisionless, supersonic plasma flows is probed using optical Thomson scattering (TS). Analysis of the TS ion-feature revealed anticorrelated modulations in the density of the two ion streams at the spatial scale of the ion skin depth c/ω pi = 120 μm, and a correlated modulation in the plasma current. The inferred current profile implies a magnetic field amplitude ~30 ± 6 T , corresponding to ~1% of the flow kinetic energy, indicating that magnetic trapping is the dominant saturation mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Parametrizing the detector response with neural networks

In high energy physics, characterizing the response of a detector to radiation is one of the most important and fundamental experimental tasks. In many cases, this task is accomplished by parameterizing summary statistics of the full detector response probability density. The parameterized detector response can then be used for calibration as well as for directly improving physics analysis sensitivity. This paper reports on how to parameterize summary statistics of the detector response using neural networks. Specifically, neural networks are powerful tools for incorporating multidimensional data and the loss function used during training determines which summary statistic is learned. One common summary statistic that has not been combined with deep learning (as far as the authors are aware) is the mode. A neural network-based approach to mode learning is proposed and empirically demonstrated in the context of high energy jet calibrations. Altogether, the neural network-based toolkit for detector response parameterization can enhance the utility of data collected at high energy physics experiments and beyond.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy-particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). In this work, we apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Young Stellar Objects, Accretion Disks, and Their Variability with Rubin Observatory LSST

Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multiband filters option available in Rubin LSST, in particular the u, g, r, i filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g., field crowdedness, massive neighbors, metallicity) by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Optimization Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide–Fast–Deep observing mode that will be used to survey the sky over the full duration of the main survey (≈10 yr). The combination of these two modes will be vital to investigate the connection between the inner-disk dynamics and longer-term eruptive variability behaviors, such as those observed on EX Lupi–type objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mapping causal patterns in crystalline solids

The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO 3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy. Localized properties, including polarization, lattice parameter, and chemical composition, are parameterized from atomic-scale imaging, and their causal relationships are reconstructed using a linear non-Gaussian acyclic model. This approach is further extended to explore the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged at both the expected locations, such as domain walls and interfaces, and at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase-separated regions in the morphotropic phase boundaries. Overall, we posit that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only the generative processes that can lead to observed microscopic states but also the causal links between multiple interacting subsystems.

Causal inference↗
DOE OSTI2020

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↗
DOE OSTI2021

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Intense windstorms in the northeastern United States

Abstract. Windstorms are a major natural hazard in many countries. The objective of this study is to identify and characterize intense windstorms during the last 4 decades in the US Northeast and determine both the sources of cyclones responsible for these events and the manner in which those cyclones differ from the cyclone climatology. The windstorm detection is based on the spatial extent of locally extreme wind speeds at 100 m height from the ERA5 reanalysis database. During the top 10 windstorms, wind speeds exceed their local 99.9th percentile over at least one-third of land-based ERA5 grid cells in this high-population-density region of the USA. Maximum sustained wind speeds at 100 m during these windstorms range from 26 to over 43 ms−1, with wind speed return periods exceeding 6.5 to 106 years (considering the top 5 % of grid cells during each storm). Property damage associated with these storms, with inflation adjusted to January 2020, ranges from USD 24 million to over USD 29 billion. Two of these windstorms are linked to decaying tropical cyclones, three are Alberta clippers, and the remaining storms are Colorado lows. Two of the 10 re-intensified off the east coast, leading to development of nor'easters. These windstorms followed frequently observed cyclone tracks but exhibit maximum intensities as measured using 700 hPa relative vorticity and mean sea level pressure that is 5–10 times the mean values for cyclones that followed similar tracks over this 40-year period. The time evolution of wind speeds and concurrent precipitation for those windstorms that occurred after the year 2000 exhibit good agreement with in situ ground-based and remote sensing observations, plus storm damage reports, indicating that the ERA5 reanalysis data have a high degree of fidelity for large, damaging windstorms such as these. A larger pool of the top 50 largest windstorms exhibit evidence of only weak serial clustering, which is in contrast to the relatively strong serial clustering of windstorms in Europe.

58 GEOSCIENCES↗
DOE OSTI2021

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗
DOE OSTI2025

Thermal Weight Determination and Interstate Coupling in State-Averaged ADAPT-VQE

Characterizing electronic thermal states at low temperatures is an important but challenging task in quantum chemistry and condensed matter physics, making it a prime candidate for a useful application in quantum computing. One of the most successful methods for state preparation on quantum computers is the Adaptive, Problem-Tailored (ADAPT) Variational Quantum Eigensolver (VQE), which has recently been generalized to treat excited states within a state-averaged framework as well as Gibbs states. In this work, we introduce Helmholtz-Optimized Thermal (HOT) ADAPT-VQE, an ancilla-free strategy for preparing Gibbs states that directly minimizes the Helmholtz free energy by targeting the dominant eigenstates of the thermal ensemble. We demonstrate the usefulness of HOT-ADAPT-VQE by predicting the free energy of two model systems with strongly correlated ground states: (1) the Fe 2+ cation in a magnetic field and (2) a [Cu 2 O 7 ] 10– fragment of the Mott insulator La 2 CuO 4 . Our results demonstrate that HOT-ADAPT-VQE significantly improves upon Gibbs-state estimates from multistate variants of ADAPT-VQE, often with substantially shallower quantum circuits, making it a promising candidate for thermal-state calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Folding Coarse-Grained Oligomer Models with PyRosetta

Non-biological foldamers are a promising class of macromolecules that share similarities to classical biopolymers such as proteins and nucleic acids. Currently, designing novel foldamers is a non-trivial process, often involving many iterations of trial synthesis and characterization until folded structures are observed. In this work, we aim to tackle these foldamer design challenges using computational modeling techniques. We developed CG PyRosetta, an extension to the popular protein folding python package, PyRosetta, which introduces coarse-grained (CG) residues into PyRosetta, enabling the folding of toy CG foldamer models. Although these models are simplified, they can help explore overarching physical hypotheses about how oligomers can form. Through systematic variation of CG parameters in these models, we can investigate various folding hypotheses at the CG scale to inform the design process of new foldamer chemistries. In this study, we demonstrate CG PyRosetta’s ability to identify minimum energy structures with a diverse structural search over a range of simple models, as well as two hypothesis-driven parameter scans investigating the effects of side-chain size and internal backbone angle on secondary structures. We are able to identify several types of secondary structures from single- and double-helices to sheet-like and knot-like structures. Here, we show how side-chain size and backbone bond angle both play an important role in the structure and energetics of these toy models. Optimal side-chain sizes promote favorable packing of side chains, while specific backbone bond angles influence the specific helix type found in folded structures.

36 MATERIALS SCIENCE↗