Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complexity theory”

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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Near Zone: Basic scattering code user's manual with space station applications

The Electromagnetic Code - Basic Scattering Code, Version 3, is a user oriented computer code to analyze near and far zone patterns of antennas in the presence of scattering structures, to provide coupling between antennas in a complex environment, and to determine radiation hazard calculations at UHF and above. The analysis is based on uniform asymptotic techniques formulated in terms of the Uniform Geometrical Theory of Diffraction (UTD). Complicated structures can be simulated by arbitrarily oriented flat plates and an infinite ground plane that can be perfectly conducting or dielectric. Also, perfectly conducting finite elliptic cylinder, elliptic cone frustum sections, and finite composite ellipsoids can be used to model the superstructure of a ship, the body of a truck, and airplane, a satellite, etc. This manual gives special consideration to space station modeling applications. This is a user manual designed to give an overall view of the operation of the computer code, to instruct a user in how to model structures, and to show the validity of the code by comparing various computed results against measured and alternative calculations such as method of moments whenever available.

Marhefka, R. J.↗
NASA NTRS1996

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Britcher, Colin P.↗
NASA NTRS2005

Planetary Gravity Fields and Their Impact on a Spacecraft Trajectory

The present work touches an interdisciplinary aspect of space exploration: the improvement of spacecraft navigation by means of enhanced planetary interior model derivation. The better the bodies in our solar system are known and modelled, the more accurately (and safely) a spacecraft can be navigated. In addition, the information about the internal structure of a planet, moon or any other planetary body can be used in arguments for different theories of solar system evolution. The focus of the work lies in a new approach for modelling the gravity field of small planetary bodies: the implementation of complex ellipsoidal coordinates (figure 1, [4]) for irregularly shaped bodies that cannot be represented well by a straightforward spheroidal approach. In order to carry out the required calculations the computer programme GRASP (Gravity Field of a Planetary Body and its Influence on a Spacecraft Trajectory) has been developed [5]. The programme furthermore allows deriving the impact of the body s gravity field on a spacecraft trajectory and thus permits predictions for future space mission flybys.

Weinwurm, G.↗
DOE OSTI2020

C–H Bond Activation via U(II) in the Reduction of Heteroleptic Bis(trimethylsilyl)amide U(III) Complexes

Reduction of (C 5 Me 5 ) 2 U III (NR 2 ) and (C 5 Me 5 )U III (NR 2 ) 2 (R = SiMe 3 ) with potassium graphite in the presence of 2.2.2-cryptand (crypt) generates dark solutions that have UV-visible spectra consistent with time-dependent density functional theory (TDDFT) calculations on the U(II) products, [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - . However, the solutions quickly change color and form the U(III) C–H bond activation products [K(crypt)][(C 5 Me 5 ) 2 U III (CH 2 SiMe 2 NSiMe 3 –κC,κN)], 1, and [K(crypt)][(C 5 Me 5 )U III (NR 2 )(CH 2 SiMe 2 NSiMe 3 –κC,κN)], 2, that were identified by X-ray crystallography. DFT calculations on the putative [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - complexes revealed 5f 3 6d 1 ground state electron configurations as previously found in isolable [(C 5 H 4 SiMe 3 ) 3 U II ] - , which indicated that these low symmetry heteroleptic complexes are reasonable precursors for new U(II) complexes.

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

Ferroelectric switching pathways and domain structure of SrBi 2 ⁢(Ta,Nb) 2 ⁢O 9 from first principles

Several families of layered perovskite oxide ferroelectrics exhibit a coupling between polarization and structural order parameters, such as octahedral rotation distortions. This coupling provides opportunities for novel electric field-based manipulation of material properties, and it also stabilizes complex domain patterns and domain-wall vortices. Among layered perovskites with such coupled orders, the Aurivillius-phase oxides SrBi 2 ⁢B 2 ⁢O 9 (B = Ta,Nb) are well-known for their excellent room-temperature ferroelectric performance. This work combines group theoretic analysis with density functional theory calculations to examine the ferroelectric switching processes of SrBi2⁢B2⁢O9. Low-energy two-step ferroelectric switching paths are identified, with polarization reversal facilitated by structural order-parameter rotations. Analysis of the domain structure reveals how the relative energetics of the coupled order parameters translates into a network of several distinct domain-wall types linked by domain-wall vortex structures. Comparisons are made between the ferroelectric switching and domain structure of SrBi 2 ⁢B 2 ⁢O 9 and those of the layered n = 2 Ruddlesden-Popper hybrid improper ferroelectrics. The results provide insight into how ferroelectric properties may be optimized by engineering the complex crystal structures of Aurivillius-phase oxides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Resolving local ordering and structure in Mn x Ge 1- x Te alloys through thermodynamic ensembles of pair distribution functions

Characterizing local bonding environments in complex materials is essential for understanding and optimizing their properties. Equally as important is the ability to predict local motifs as a function of synthesis conditions, enhancing chemists’ ability to design properties into materials. In this study, we present an approach to leverage statistical mechanics to generate temperature- and energy-informed ensemble averaged pair distribution functions (PDFs). This method, which we have named Thermodynamic Ensemble Averages of PDFs for Ordering and Transformations (TEAPOT), utilizes density functional theory (DFT) to relax supercells while incorporating energetic penalties for local order, enabling accurate and computationally efficient analysis of local structure. We apply this method to the neutron PDF measurements of the pseudobinary MnTe–GeTe (MGT) alloy, demonstrating its capability to resolve complex local distortions and chemical ordering. Our results reveal detailed insights into phase transformations and local distortions driven by Mn substitution. For compositions that globally present as rock salt, our analysis reveals that Ge coordination geometry is heavily impacted by synthesis temperature. We propose that high temperature synthesis conditions promote a lowered Ge polyhedra distortion, promoting high charge carrier mobility due to the alignment of local and global structure. Incorporating statistical mechanics and computation into experimental analysis thus guides synthesis of tailored local structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Variational design principles for nonequilibrium colloidal assembly

Using large deviation theory and principles of stochastic optimal control, we show that rare molecular dynamics trajectories conditioned on assembling a specific target structure encode a set of interactions and external forces that lead to enhanced stability of that structure. Such a relationship can be formulated into a variational principle, for which we have developed an associated optimization algorithm and have used it to determine optimal forces for targeted self-assembly within nonequilibrium steady-states. We illustrate this perspective on inverse design in a model of colloidal cluster assembly within linear shear flow. We find that colloidal clusters can be assembled with high yield using specific short-range interactions of tunable complexity. Shear decreases the yields of rigid clusters, while small values of shear increase the yields of nonrigid clusters. The enhancement or suppression of the yield due to shear is rationalized with a generalized linear response theory. Furthermore, by studying 21 unique clusters made of six, seven, or eight particles, we uncover basic design principles for targeted assembly out of equilibrium.

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

Strategically Modified Ligand Incorporating Mixed Phosphonate and Carboxylate Groups to Enhance Performance in All-Iron Redox Flow Batteries

Iron redox flow batteries (Fe-RFBs) hold significant promise for achieving cost-effectiveness and utilizing abundant materials for stationary energy storage applications. Here, a design of a novel Fe complex utilizing a nitrogenous phosphonate/carboxylate mixed ligand, N,N-Bis(phosphonomethyl)glycine (BMPG), is presented to achieve high performance Fe anolyte. Compared to its all-phosphonate form, nitrilotri(methylphosphonic acid) (NTMPA), the new complex Fe(BPMG) 2 demonstrates a negatively shifted redox potential, resulting in ≈0.07 V (≈10%) increase in battery output voltage. Full battery testing paired with ferrocyanide catholyte demonstrates stable cycling (capacity degradation <0.0001%/cycle) over 730 consecutive charge/discharge cycles with Coulombic Efficiency of 100% at a current density of 20 mA cm -2 under near neutral pH (≈8). Of particular interest, density functional theory (DFT) studies and operando Raman measurements provide strong evidence supporting a molecular structure in BPMG, which reveals the mixed phosphonate/carboxylate groups in BPMG maintain the octahedral coordination of the Fe ion center with phosphonates exclusively, while leaving the carboxylate unbound for both Fe(II) and Fe(III) complexes. This structural similarity between BPMG-based Fe(II) and Fe(III) complexes effectively mitigates the slow redox reaction kinetics observed in Fe(NTMPA) 2 anolyte, where significant ligand reorientation occurs between Fe(II) and Fe(III) complexes.

25 ENERGY STORAGE↗
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

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

Physically Informed Machine Learning Prediction of Electronic Density of States

The electronic structure of a material, such as its density of states (DOS), provides key insights into its physical and functional properties and serves as a valuable source of high-quality features for many materials screening and discovery workflows. Still, the computational cost of calculating the DOS, most commonly with density functional theory (DFT), becomes prohibitive for meeting high-fidelity or high-throughput requirements, necessitating a cheaper but sufficiently accurate surrogate. To fulfill this demand, we develop a general machine learning method based on graph neural networks for predicting the DOS purely from atomic positions, six orders of magnitude faster than DFT. This approach can effectively use large materials databases and be applied generally across the entire periodic table to materials classes of arbitrary compositional and structural diversity. We furthermore devise a highly adaptable scheme for physically informed learning which encourages the DOS prediction to favor physically reasonable solutions defined by any set of desired constraints. This functionality provides a means for ensuring that the predicted DOS is reliable enough to be used as an input to downstream materials screening workflows to predict more complex functional properties, which rely on accurate physical features.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct Transformation of SiH 4 to a Molecular L(H) 2 Co=Si=Co(H) 2 L Silicide Complex

The synthesis of bimetallic molecular silicide complexes is reported, based on the use of multiple Si–H bond activations in SiH 4 at the metal centers of 14-electron LCo I fragments (L = Tp", HB(3,5-diisopropylpyrazolyl) 3 – ; [BP 2 tBu Pz], PhB(CH 2 P t Bu 2 ) 2 (pyrazolyl)). Upon exposure of (Tp"Co) 2 (μ-N 2 ) (1) to SiH 4 , a mixture of (Tp"Co) 2 (μ-H) (2) and (Tp"Co) 2 (μ-H) 2 (3) was formed and no evidence for Si–H oxidative addition products was observed. In contrast, [BP 2 tBu Pz]-supported Co complexes led to Si–H oxidative additions with the generation of silylene and silicide complexes as products. Notably, the reaction of ([BP 2 tBu Pz]Co) 2 (μ-N 2 ) (5) with SiH 4 gave the dicobalt silicide complex [BP 2 tBu Pz](H) 2 Co=Si=Co(H) 2 [BP 2 tBu Pz] (8) in high yield, representing the first direct route to a symmetrical bimetallic silicide. Here, the effect of the [BP 2 tBu Pz] ligand on Co–Si bonding in 7 and 8 was explored by analysis of solid-state molecular structures and density functional theory (DFT) investigations. Upon exposure to CO or DMAP (DMAP = 4-dimethylaminopyridine), 8 converted to the corresponding [BP 2 tBu Pz]Co(L) x adducts (L = CO, x = 2; L = DMAP, x = 1) with concomitant loss of SiH 4 , despite the lack of significant Si–H interactions in the starting complex. On heating to 60 °C, 8 underwent reaction with MeCl to produce small quantities of Me x SiH 4–x (x = 1–3), demonstrating functionalization of the μ-silicon atom in a molecular silicide to form organosilanes.

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

Importance of Dispersion in the Molecular Geometries of Mn(III) Spin-Crossover Complexes

The computational investigation of the molecular geometries of a pair of manganese(III) spin-crossover complexes is reported. For the geometry of the quintet high-spin state, density functionals significantly overestimate Mn–Namine bond distances, although the geometry for the triplet intermediate-spin state is well described. Here, comparisons with several wave function-based methods demonstrate that this error is due to the limited ability of commonly used density functionals to recover dispersion beyond a certain extent. Among the methods employed for geometry optimization, restricted open-shell Møller–Plesset perturbation theory (MP2) appropriately describes the high-spin geometry but results in a slightly shorter Mn–O distance in both spin states. On the other hand, extended multistate complete active space second-order perturbation theory (XMS-CASPT2) provides a good description of the geometry for the intermediate-spin state but also sufficiently recovers dispersion, performing well for the high-spin state. Despite the fact that the electronic structure of both spin states is dominated by one-electron configuration, XMS-CASPT2 offers a balanced approach, leading to molecular geometries with much better agreement with experiment than MP2 and DFT. A scan along the Mn–N amine bond demonstrates that for these complexes coupled cluster methods (i.e., DLPNO-CCSD(T)) also yield bond distances in agreement with experiment while multiconfiguration pair density functional theory (MC-PDFT) is unable to recover dispersion well enough, analogous to single-reference DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Chemical differentiation in regions of massive star formation

We have reexamined the origin of the apparent differentiation between nitrogen-bearing molecules and complex oxygen-bearing molecules that is observed in hot molecular cores associated with massive protostars. Observations show that methanol is an ubiquitous and abundant component of protostellar ices. Recent observations suggest that ammonia may constitute an appreciable fraction of the ices toward some sources. In contrast to previous theories that suggested that N/O differentiation was caused by an anticorrelation between methanol and ammonia in the precursor grain mantles, we show that the presence of ammonia in mantles and the core temperature are key quantities in determining N/O differentiation. Calculations are presented which show that when large amounts of ammonia are evaporated alkyl cation transfer reactions are suppressed and the abundances of complex O-bearing organic molecules greatly reduced. Cooler cores (100 K) eventually evolve to an oxygen-rich chemical state similar to that attained when no ammonia was injected, but on a timescale that is an order of magnitude longer (~10(5) yr). Hotter cores (300 K) never evolve an O-rich chemistry unless ammonia is almost absent from the mantles. In this latter case, a complex O-rich chemistry develops on a timescale of ~10(4) yr, as in previous models, but disappears in about 2 x 10(5) yr, after which time the core is rich in NH3, HCN, and other N-bearing molecules. There are thus two ways in which N-rich cores can occur. We briefly discuss the implications for the determination of hot-core ages and for explaining N/O differentiation in several well-studied sources.

NASA Center ARC↗
DOE OSTI2025

Probing Effects of Electron Addition to Polycyclic Conjugated Hydrocarbons Containing Four‐Membered Rings

The rectangular cyclobutadiene (CBD, C 4 H 4 ) is a unique moiety for building nonbenzenoid polycyclic conjugated hydrocarbons with interesting electron‐accepting properties. Herein, the investigation on chemical reduction of several CBD‐containing polycyclic hydrocarbons with increasing conjugation length is reported: biphenylene (C 12 H 8 ), dimethyl[2]naphthalene (C 22 H 16 ), and tetramethyl‐dibenzo‐[3]phenylene (C 30 H 22 ). The two‐step sequential reduction is first demonstrated by in situ spectroscopic investigation and then confirmed by the isolation of single crystals of the reduced products. The X‐ray crystallographic analysis reveals the formation of several mono‐ and doubly reduced products in solvent‐separated and complexed forms. The crystal structures for both neutral parents and corresponding reduced products unravel the changes in bond alternation in each ring of the fused systems. Density functional theory (DFT) and nucleus‐independent chemical shift (NICS) scan calculations reveal that the two‐electron addition reduces the aromatic character in the benzenoid rings but has minor influence on the antiaromatic CBD rings.

Zhu, Yikun [Department of Chemistry University at ↗
DOE OSTI2026

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2022

Analysis of Intermediates and Products from the Dehydrogenation of Mg(BH 4 ) 2

The thermodynamic properties of key compounds, Mg(B 3 H 8 ) 2 , MgB 2 H 6 , MgB 10 H 10 , Mg(B 11 H 14 ) 2 , Mg 3 (B 3 H 6 ) 2 , and MgB 12 H 12 , proposed to be formed in the release of hydrogen from magnesium borohydride Mg(BH 4 ) 2 and uptake of hydrogen by MgB 2 , have been investigated using solid–state density functional theory (DFT) calculations. Further, more accurate treatment of cell–size effects to the entropies were also investigated, in order to improve the accuracy of the thermodynamic properties of complex borohydrides. We find that the zero–point energy corrections can lower the electronic energies of reaction by 20–30 kJ/(mol H 2 ) for these intermediates, while adding the thermal and entropy contributions results in the total decrease up to ~50 kJ/(mol H 2 ). Although our treatment lowers the calculated formation energy of Mg(B 3 H 8 ) 2 , it is still too high to explain the experimental observation of B 3 H 8 – . We discuss possible reasons for this disparity and propose that the formation of B 3 H 8 – and H – in a disordered amorphous phase has a large energy difference compared to the phase–separated Mg(B 3 H 8 ) 2 and MgH 2 considered in calculations. Comparison of the experimental and NMR chemical shifts calculated within a DFT approach for known species Mg(BH 4 ) 2 , Mg(B 3 H 8 ) 2 , Mg(B 11 H 14 ) 2 , MgB 10 H 10 and MgB 12 H 12 provides validation for predicting the chemical shifts of the other compounds which are yet to be confirmed experimentally. These include MgB 2 H 6 and the proposed tri–anion species Mg 3 (B 3 H 6 ) 2 that both have favorable thermodynamics for reversible hydrogen storage in Mg(BH 4 ) 2 without the formation of MgH 2 as a co–product which could phase–separate and inhibit rehydrogenation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗