Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dispersion function”

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 91 records · Page 5

Crustal Structure of the Mesopotamian Plain, East of Iraq

The crustal structure of Iraq was investigated through analyzing teleseismic data from 12 new seismic stations. Three seismic stations are located within the Zagros Fold-Thrust Belt, and nine are within the Mesopotamian Plain. Joint inversion of P wave receiver function and Rayleigh wave dispersion data were employed to resolve S wave velocity structure models beneath each station. Combining these models with available Moho depths within and around the study area reveals that the Moho depth is smoothly increasing from the Arabian platform toward the Zagros Mountains. Further, an exceptional deeper root is observed near the eastern edge of the Mesopotamian Plain where sedimentary pile is the thickest. This root was interpreted as a structure inherited from the successive Mesozoic rifting of the NE Arabian platform enhanced by sedimentary loading and progressive subsidence. A low-velocity uppermost mantle was observed beneath the Arabian Foredeep attesting the existence of a low-strength lithospheric mantle beneath the region southwest from the Zagros deformation front. The weak rheology of the uppermost mantle may have allowed local sinking of the crust and deepening of the Moho boundary due to vertical loading and the Late Miocene lateral contraction.

58 GEOSCIENCES↗

Machine-learning-informed scattering correlation analysis of sheared colloids

We have carried out theoretical analysis, Monte Carlo simulations and machine-learning analysis to quantify microscopic rearrangements of dilute dispersions of spherical colloidal particles from coherent scattering intensity. Both monodisperse and polydisperse dispersions of colloids were created and underwent a rearrangement consisting of an affine simple shear and non-affine rearrangement using the Monte Carlo method. We calculated the coherent scattering intensity of the dispersions and the correlation function of intensity before and after the rearrangement and generated a large data set of angular correlation functions for varying system parameters, including number density, polydispersity, shear strain and non-affine rearrangement. Singular value decomposition of the data set shows the feasibility of machine-learning inversion from the correlation function for the polydispersity, shear strain and non-affine rearrangement using only three parameters. A Gaussian process regressor is then trained on the data set and can retrieve the affine shear strain, non-affine rearrangement and polydispersity with relative errors of 3%, 1% and 6%, respectively. Altogether, our model provides a framework for quantitative studies of both steady and non-steady microscopic dynamics of colloidal dispersions using coherent scattering methods.

Gaussian process regression↗

Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost

Hybrid or “extended” symmetry-adapted perturbation theory (XSAPT) replaces traditional SAPT’s treatment of dispersion with better performing alternatives while at the same time extending two-body (dimer) SAPT to a many-body treatment of polarization using a self-consistent charge embedding procedure. The present work presents a systematic study of how XSAPT interaction energies and energy components converge with respect to the choice of Gaussian basis set. Errors can be reduced in a systematic way using correlation-consistent basis sets, with aug-cc-pVTZ results converged within <0.1 kcal/mol. Similar (if slightly less systematic) behavior is obtained using Karlsruhe basis sets at much lower cost, and we introduce new versions with limited augmentation that are even more efficient. Pople-style basis sets, which are more efficient still, often afford good results if a large number of polarization functions are included. The dispersion models used in XSAPT afford much faster basis-set convergence as compared to the perturbative description of dispersion in conventional SAPT, meaning that “compromise” basis sets (such as jun-cc-pVDZ) are no longer required and benchmark-quality results can be obtained using triple-ζ basis sets. The use of diffuse functions proves to be essential, especially for the description of hydrogen bonds. As a result, the “δ(Hartree–Fock)” correction for high-order induction can be performed in double-ζ basis sets without significant loss of accuracy, leading to a mixed-basis approach that offers 4× speedup over the existing (cubic scaling) XSAPT approach.

74 ATOMIC AND MOLECULAR PHYSICS↗

The HectoMAP Cluster Survey: Spectroscopically Identified Clusters and their Brightest Cluster Galaxies (BCGs)

We apply a friends-of-friends (FoF) algorithm to identify galaxy clusters and we use the catalog to explore the evolutionary synergy between brightest cluster galaxies (BCGs) and their host clusters. We base the cluster catalog on the dense HectoMAP redshift survey (2000 redshifts deg -2 ). The HectoMAP FoF catalog includes 346 clusters with 10 or more spectroscopic members within the range 0.05 < z < 0.55 and with a median z = 0.29. We list these clusters and their members. We also include central velocity dispersions (σ*, BCG ) for the FoF cluster BCGs, a distinctive feature of the HectoMAP FoF catalog. HectoMAP clusters with higher galaxy number density (80 systems) are all genuine clusters with a strong concentration and a prominent BCG in Subaru/Hyper Suprime-Cam images. The phase-space diagrams show the expected elongation along the line of sight. Lower-density systems include some low reliability systems. We establish a connection between BCGs and their host clusters by demonstrating that σ*, BCG /σ cl decreases as a function of cluster velocity dispersion (σ cl ), in contrast, numerical simulations predict a constant σ*, BCG /σ cl . Sets of clusters at two different redshifts show that BCG evolution in massive systems is slow over the redshift range z < 0.4. The data strongly suggest that minor mergers may play an important role in BCG evolution in clusters with σ cl ≳ 300 km s -1 . For lower mass systems (σ cl < 300 km s -1 ), major mergers may play a significant role. The coordinated evolution of BCGs and their host clusters provides an interesting test of simulations in high-density regions of the universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Tailoring the Near-Surface Environment of Rh Single-Atom Catalysts for Selective CO 2 Hydrogenation

For this study, we used a combination of experimental spectroscopies, density functional theory calculations, and CO 2 hydrogenation studies to investigate the effects of modifying single-atom Rh 1 /TiO 2 catalysts with functionalized phosphonic acid monolayers. We found that the deposition of specific amine-functionalized ligands resulted in an ~8× increase in site-specific CO 2 reduction turnover frequency at 150 °C and a ~2× increase at 250 °C. On-stream stability also improved following ligand deposition. The effect of the modifier on reactivity was highly sensitive to the proximity of the amine functional group to the surface, which was controlled by adjusting the length of the phosphonic acid tail. Furthermore, deposition of alkyl phosphonic acids without an amine functional group resulted in blocked CO 2 adsorption and a near-complete loss of catalytic activity. Infrared spectroscopy studies suggested that the amine group provided binding sites for CO 2 that enabled hydrogenation when the amine was positioned near a Rh1 site. Phosphonic acid-modified catalysts also exhibited high selectivity to CO over the series product methane; the selectivity effect was traced to modification of the Rh1 sites to favor CO desorption. Phosphonic acid deposition resulted in 80–90% loss of accessible Rh1 sites, likely due to blocking by tail groups. However, even with the loss of sites, under low-temperature reaction conditions, the rates of CO 2 hydrogenation were improved with the coatings, indicating that the remaining sites are highly efficient. Organic functionalization of the supports for atomically dispersed catalysts offers the opportunity to precisely control the positioning of functional groups in the vicinity of a well-defined active site, potentially enabling an additional level of control over active site design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Amorphous dispersions of flubendazole in hydroxypropyl methylcellulose: Formulation stability assisted by pair distribution function analysis

Here, we use X-ray pair distribution function (PDF) analysis applied to high-energy synchrotron X-ray powder diffraction data to evaluate the amorphous solid dispersions interactions and their aging stability. The obtained systems are based on hydroxypropyl methylcellulose (hypromellose) derivatives and flubendazole (FBZ) drug dispersions prepared using a spray-dryer technique. We carry out stability studies under aging parameters (40 °C/75% relative humidity) to tune the systems’ recrystallization. The results reveal that ion-base interactions between the drug-polymer matrix are responsible for reducing clustering processes yielding slower recrystallization and different ordering in the hypromellose phthalate (HPMCP/FBZ) and hypromellose acetate succinate (HPMC-AS/FBZ) systems and complete drug clustering in hypromellose (HPMC-E3/FBZ). The structural ordering was accessed using differential X-ray PDFs that revealed the region between 3.5 Å and 5.0 Å could be related to FBZ intermolecular interactions and is more ordered for the least stable system (HPMC-E3/FBZ) and less ordered for the most stable system (HPMCP/FBZ). These results show that the ion-base interactions between drug and matrix occur at these intermolecular distances.

60 APPLIED LIFE SCIENCES↗

Efficient Pd on carbon catalyst for ammonium formate dehydrogenation: Effect of surface oxygen functional groups

Formate solution is one of a kind promising liquid organic hydrogen carrier but suffers from kinetic challenges. This study investigates the kinetics of the surfacefunctionalized Pd on carbon catalysts for formate dehydrogenation and the impact of O-functional groups. Here, the fraction of the distinguished O-functional groups was modulated by the different concentrated HNO 3 solutions treatment or by H 2 reduction. This study shows that the O-functional groups play critical roles in dispersing Pd nanoparticles and decreasing the activation energy for dehydrogenation. Density functional theory (DFT) calculation reveals that most O-functional groups enhance formate adsorption on the Pd active site. However, the existence of C=O groups consumes reducing agents and hinders the formation of metallic Pd. The electron transfer from Pd to oxygen functional groups is unfavorable to dehydrogenation. The as-prepared Pd5/re-ACA (reduced activated carbon washed by acid) exhibited significant activity with a higher turnover frequency of 13,511 h -1 than commercial Pd/C.

54 ENVIRONMENTAL SCIENCES↗

Cost function for low-dimensional manifold topology assessment

Abstract In reduced-order modeling, complex systems that exhibit high state-space dimensionality are described and evolved using a small number of parameters. These parameters can be obtained in a data-driven way, where a high-dimensional dataset is projected onto a lower-dimensional basis. A complex system is then restricted to states on a low-dimensional manifold where it can be efficiently modeled. While this approach brings computational benefits, obtaining a good quality of the manifold topology becomes a crucial aspect when models, such as nonlinear regression, are built on top of the manifold. Here, we present a quantitative metric for characterizing manifold topologies. Our metric pays attention to non-uniqueness and spatial gradients in physical quantities of interest, and can be applied to manifolds of arbitrary dimensionality. Using the metric as a cost function in optimization algorithms, we show that optimized low-dimensional projections can be found. We delineate a few applications of the cost function to datasets representing argon plasma, reacting flows and atmospheric pollutant dispersion. We demonstrate how the cost function can assess various dimensionality reduction and manifold learning techniques as well as data preprocessing strategies in their capacity to yield quality low-dimensional projections. We show that improved manifold topologies can facilitate building nonlinear regression models.

42 ENGINEERING↗

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dispersion formulas in QFTs, CFTs and holography

We study momentum space dispersion formulas in general QFTs and their applications for CFT correlation functions. We show, using two independent methods, that QFT dispersion formulas can be written in terms of causal commutators. The first derivation uses analyticity properties of retarded correlators in momentum space. The second derivation uses the largest time equation and the defining properties of the time-ordered product. At four points we show that the momentum space QFT dispersion formula depends on the same causal double-commutators as the CFT dispersion formula. At n -points, the QFT dispersion formula depends on a sum of nested advanced commutators. For CFT four-point functions, we show that the momentum space dispersion formula is equivalent to the CFT dispersion formula, up to possible semi-local terms. We also show that the Polyakov-Regge expansions associated to the momentum space and CFT dispersion formulas are related by a Fourier transform. In the process, we prove that the momentum space conformal blocks of the causal double-commutator are equal to cut Witten diagrams. Finally, by combining the momentum space dispersion formulas with the AdS Cutkosky rules, we find a complete, bulk unitarity method for AdS/CFT correlators in momentum space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Theory meets experiment: Insights into structure and magnetic properties of Fe 1-x Ni x B alloy

Here, we studied the structural and magnetic properties of the solid solution Fe 1-x Ni x B through theoretical and experimental approaches. Powder X-ray diffraction, X-ray Pair Distribution Function analysis, and energy dispersive X-ray spectroscopy reveal that the Fe 1-x Ni x B solid solution crystallizes in the β-FeB structure type up to x = 0.6–0.7 and exhibits anisotropic unit cell volume contraction with increasing Ni concentration. Magnetic measurements showed a transition from ferromagnetism to paramagnetism around x = 0.7. For x = 0.5, the low (< 0.3 μ B ) magnetic moments suggest itinerant magnetism despite the relatively high Curie temperature (up to 225 K). Theoretical calculations indicated different types of magnetic orderings depending on the Fe/Ni atomic order, with the antiferromagnetic state being stable for ordered FeNiB 2 , whereas the ground state is ferromagnetic for the disordered alloy. Calculations also predicted the coexistence of low- and high-spin states in Fe atoms around the composition with x = 0.5, in line with the experimental evidence from 57 Fe Mössbauer spectroscopy. The two magnetically distinct Fe sites for x = 0.3, 0.4, and 0.5 observed by 57 Fe Mössbauer spectroscopy can also be interpreted as two magnetically different regions or clusters. The formation of these clusters could affect the critical behavior near a quantum magnetic transition based on a potential ferromagnetic quantum critical point identified computationally and experimentally near x = 0.64. This work highlights the complex interplay between structure and magnetism in Fe 1-x Ni x B alloys, suggesting areas for future research on quantum critical behavior.

Bhaskar, Gourab [Iowa State Univ., Ames, IA (Unite↗

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase Stability, Band Gap Tuning, and Rashba Splitting in Selenium-Alloyed Bournonite: CuPbSb(S 1– x Se x ) 3

Recently, bournonite (CuPbSbS 3 ) has been identified as a potential ferroelectric photovoltaic (PV) material with a 1.3 eV band gap, which falls in an appropriate range for single junction PV devices. Progress in applying bournonite has yielded several studies reporting successful thin-film processing, the most recent of which demonstrated a 2.65% power conversion efficiency for a bournonite-based PV device. In an effort to explore bournonite band gap engineering for PV and other prospective applications (e.g., thermoelectric, spintronic), we here consider the solid-state synthesis of selenium-alloyed bournonite, CuPbSb(S 1-x Se x ) 3 , across the full range of x (0.0 ≤ x ≤ 1.0), and report phase purity for 0.0 ≤ x ≤ 0.5. We characterize the crystal structure and band gap of the samples using X-ray diffraction and diffuse reflectance spectroscopy, determining a band gap decrease for single-phase samples from 1.25 eV at x = 0.0 to 1.06 eV at x = 0.5. Formation energies determined using dispersion corrected hybrid density functional theory (DFT) support experimental findings and are consistent with the stability (instability) of the x = 0.5 (x = 1.0) structure relative to starting materials. The computed band structures show a decreasing trend in the band gap with increasing x, again consistent with experiment, with the states at the conduction (valence) band minima (maxima) being principally derived from Pb (S/Se) states. Spin texture analysis and detailed comparison of spin splitting parameters show that Se-alloying modifies the band gap while maintaining Rashba spin splitting character within the electronic structure. Rashba splitting is most pronounced for the conduction band minimum along the Γ-X path in the Brillouin zone. Energy separation between spin states is maintained at ΔE~0.2 eV for the various Se contents, with σz as the spin direction. Furthermore, a smaller spin splitting occurs along Γ-Z (decreasing value with increasing x), with σ x as the spin direction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ALD with Alternative Coreactants: Which Work, Which Do Not, and Why

In this paper, we present a combined experimental and theoretical study that systematically evaluates a series of alcohols (i.e., primary, secondary, and tertiary) as coreactants for trimethylaluminum (TMA)-based atomic layer deposition (ALD) of Al 2 O 3 . We employed in situ quartz crystal microbalance techniques and ex situ X-ray photoelectron spectroscopy to probe growth across a range of temperatures from T = 120 to 285 °C. Dispersion-inclusive hybrid density functional theory was employed to identify potential reaction pathways and compute corresponding Gibbs free-energy barriers and rate coefficients. Experimentally on an Al 2 O 3 surface, sustained thin-film growth, where the rate is constant over many cycles, was only observed with tertiary alcohols (tert-butanol and 2-methyl-2-butanol) at elevated temperatures (e.g., T = 285 °C); for tert-butanol specifically, sustained growth also occurred at T = 240 °C, with no sustained growth at or below T = 210 °C. Growth using primary and secondary alcohols decayed to a negligible value after only 2–3 cycles on an Al 2 O 3 surface at all tested temperatures. Intentional trace additions of H 2 O to anhydrous alcohols restored sustained growth but at reduced deposition rates relative to pure H 2 O. Theoretical analysis supports a bimolecular alkoxy β-H elimination mechanism as the primary reaction pathway: initiated by the formation of a bound alkoxy species, the subsequent β-H elimination step is rate-determining, with a barrier that systematically decreases as the degree of the alcohol coreactant is increased from primary to secondary to tertiary. Computed rate coefficients (k) also indicate a strong temperature dependence (with k increasing by ∼4–5 orders of magnitude from T = 120 to 285 °C) and are consistent with the experimental observation that only tertiary alcohols lead to steady thin-film growth at elevated temperatures. This integrated approach establishes why only certain alcohols can sustain Al 2 O 3 ALD and delivers a predictive framework for identifying effective alcohol coreactants.

36 MATERIALS SCIENCE↗

Mechanochemical Solid Form Screening of Zeolitic Imidazolate Frameworks Using Structure-Directing Liquid Additives

We demonstrate a systematic application of the mechanochemical liquid-assisted grinding (LAG) methodology to screen for forms of zinc imidazolate (ZnIm 2 ), of fundamental importance as the simplest member of the zeolitic imidazolate framework materials family. The exploration of 45 different liquid additives, selected based on their molecular structure and physicochemical properties has resulted in eight different ZnIm 2 topological forms, appearing in 13 crystallographically distinct solid forms (including two previously unknown forms of the crb (BCT) topology), amorphous phases, and the interrupted moc-Zn 4 Im 8 HIm. All prepared topological forms were also explored computationally, using dispersion-corrected periodic density functional theory (DFT) calculations, enabling the rationalization of screening outcomes, and setting the stage for future prediction of additive-directed metal–organic framework (MOF) synthesis. This first systematic exploration of LAG in screening for three-dimensional MOFs demonstrates the potential of the liquid additive to not only accelerate materials synthesis, but also to direct it toward topologically different MOFs. The discovery of novel forms of a material that already exhibits at least 21 crystallographically and functionally different forms provides a strong testimony on the power of mechanochemistry in metal–organic materials discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗