SEARCH · Engineering Papers

Results for “energy use behavior”

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.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2024

NUMERICAL MODELING OF A SOLID OXIDE FUEL CELL FOR USE IN REAL-TIME SIMULATION AND CYBER-PHYSICAL SYSTEMS

Cyber-physical systems provide a mechanism with which to investigate the physical phenomena and behavior of traditionally cost-prohibitive or otherwise fragile equipment. For the National Energy Technology Laboratory (NETL), this approach resulted in the Hybrid Performance (Hyper) facility which features a gas turbine-SOFC hybrid cycle utilizing real turbomachinery and a simulated SOFC stack. This allows for the investigation of combined cycle performance and control strategies, in an exhaustive manner, both without fear of destroying delicate state-of-the-art fuel cells, and with the full accuracy of real-world turbomachinery. Issues arose between the transient response of the SOFC model being limited to a sample time of 80 milliseconds, due to the calculation time of the SOFC model taking on average 40 milliseconds to calculate for a given timestep with spikes in calculation time reaching the 80 millisecond threshold. In order to be able to match the speed of transients from the turbomachinery and likewise better discern transient behavior, it was determined that the SOFC model must be optimized to operate at a sample time of 5 milliseconds. Therefore, it is necessary to optimize the SOFC model in order to decrease the calculation time from around 40 milliseconds, down to at the most 5 milliseconds. To do this, both the electrochemical algorithm and the thermal algorithm used to simulate the physical behavior of the SOFC are investigated to determine where improvements can be made. To this end the rootfinding numerical recipes of the electrochemical algorithm are investigated as the complex electrochemistry requires a highly iterative nested dual convergence loop to resolve the voltage-current relationship, and likewise the temporal discretization of the thermal algorithm is modified for the sake of higher accuracy and stability. Ultimately the new electrochemical algorithm featuring higher order rootfinding schemes proves to be efficient enough to reach the sub 5 millisecond target, signifying an order of magnitude reduction in calculation time, and when coupled with the new temporal discretization similar calculation time characteristics show that a fully implicit, higher order temporal discretization can also successfully be used if desired. Ultimately this result means that the cyber-physical simulation system can operate at higher sample rates, and resolve transient events at significantly higher resolution and fidelity.

Arias, Jesus↗

DOE OSTI2020

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1991

Relaxation processes in a low-order three-dimensional magnetohydrodynamics model

The time asymptotic behavior of a Galerkin model of 3D magnetohydrodynamics (MHD) has been interpreted using the selective decay and dynamic alignment relaxation theories. A large number of simulations has been performed that scan a parameter space defined by the rugged ideal invariants, including energy, cross helicity, and magnetic helicity. It is concluded that time asymptotic state can be interpreted as a relaxation to minimum energy. A simple decay model, based on absolute equilibrium theory, is found to predict a mapping of initial onto time asymptotic states, and to accurately describe the long time behavior of the runs when magnetic helicity is present. Attention is also given to two processes, operating on time scales shorter than selective decay and dynamic alignment, in which the ratio of kinetic to magnetic energy relaxes to values 0(1). The faster of the two processes takes states initially dominant in magnetic energy to a state of near-equipartition between kinetic and magnetic energy through power law growth of kinetic energy. The other process takes states initially dominant in kinetic energy to the near-equipartitioned state through exponential growth of magnetic energy.

Stribling, Troy↗

NASA NTRS2018

Evaluation of a Thermodynamic Ice Crystal Icing Model Using Experimental Ice Accretion Data

This paper presents the evaluation of a thermodynamic ice crystal icing model, previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. It has been proposed that there are two types of distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. In the former, ice accretion occurs where a freeze fraction (0 to 1) of melted ice crystals freezes on a surface, along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized by having strong adhesion to the surface. In the latter, icing occurs from accumulated unmelted ice on a surface, where a melt fraction (0 to 1) dictates the amount of unmelted impinged ice. This melt-dominated icing is characterized by weakly bonded surface adhesion. The experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. This mass loss parameter must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface. This loss parameter, however, along with the freeze and melt fraction, are the only experimental parameters that are currently not measured directly. Using reported icing growth rates from published ice crystal icing experiments, a methodology is proposed to determine these unknown parameters. This work takes reported ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) in 2016 and tests NASA collaborated on with the National Research Council (NRC) of Canada in 2012 that examined the fundamental physics of ice crystal icing. Those research efforts sought to generate icing conditions representative of those that occur inside a jet engine when ingesting ice crystals. This paper presents the fundamental equations of the thermodynamic model, the methodology used to determine the aforementioned unknown icing parameters, and results from model evaluation using experimental data. In addition, this paper builds on the previously proposed model by adding a transient conduction term to explain ice growth behavior at the onset of experimental tests that was observed to be different from steady-state ice growth that occurred later in the test run.With the addition of this energy term, this becomes a quasi-steady model. A key finding from this work suggests that mass loss fractions can exceed 0.90 for steady ice growth periods. In addition, due to conductive heat fluxes when using a warmer-than-freezing airfoil, lower mass loss fraction values were calculated during the initial transient period.

Ice-Crystal Ice Accretion↗

DOE OSTI2020

Asymptotic Behavior of Cosmologies with $$\Lambda >0$$ in $$2+1$$ Dimensions

We study, using mean curvature flow methods, 2+1 dimensional cosmologies with a positive cosmological constant and matter satisfying the dominant and the strong energy conditions. If the spatial slices are compact with non-positive Euler characteristic and are initially expanding everywhere, then we prove that the spatial slices reach infinite volume, asymptotically converge on average to de Sitter and they become, almost everywhere, physically indistinguishable from de Sitter. This holds true notwithstanding the presence of initial arbitrarily-large density fluctuations and the formation of black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Topography-Induced TKE Budget Behavior Over an Amazon Forest

Data from three different heights (35, 50 and 81 m) of one of the ATTO Project towers in the Amazon forest were used to calculate the TKE (turbulence kinetic energy) budget and some other statistics within the RSL (roughness sublayer). The statistical analyses were carried out for unstable and stable cases. The vertical transport and the vertical advection terms do not explain the imbalances found in the TKE budget, highlighting the impor- tance of horizontal terms in complex terrain. Indeed, patterns independent of the time of day were found for vertical transport, horizontal and vertical velocity skewness, the mean vertical component velocity, and momentum and vertical TKE flux. Here, the highest values of mean vertical velocity were observed from the direction where there is a valley near the tower, which may indicate a topographical effect. Monin-Obukhov Similarity Theory (MOST) predictions for the mean velocity gradient, strictly applicable only in the inertial sublayer (ISL), gave reasonable predictions at 35 and 50 m, but dimensionless dissipation rates displayed large scatter. At 81 m, there were clear deviations from MOST for the two functions, disclosing the effect of topography. In addition, a neutral LES (large-eddy simu- lation) study showed good agreement with tower data for the same general wind direction.

ATTO project↗

DOE OSTI2026

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗

DOE OSTI2020

Evaluation of a V 8 C 7 Anode for Oxygen Evolution in Alkaline Media: Unusual Morphological Behavior

Metallic vanadium carbide (V 8 C 7 ) with cubic symmetry is examined as an oxygen evolution reaction (OER) precatalyst in alkaline media. Herein, we used quasi in situ scanning electron microscopy and energy-dispersive X-ray spectrometry to investigate the structural transformation of the precatalyst V 8 C 7 microparticles during extended cyclic voltammetric (CV) OER testing. Interestingly, an anisotropic morphological transformation (from a distorted sphere to a cuboid) of V 8 C 7 was observed. Our theoretical and experimental results strongly suggest that this morphological change happens due to the selective self-oxidation and dissolution of the V 8 C 7 (110) and (111) surfaces and the subsequent exposure of the relatively stable (100), (010), and (001) surfaces. Lastly, these results also suggest that these stable facets are preferable for the OER, and a current density of 10 mA·cm –2 was delivered at an overpotential of 503 mV after the extended CV OER testing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Cutting out the middleman: calibrating and validating a dynamic vegetation model (ED2-PROSPECT5) using remotely sensed surface reflectance

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2021

Cutting out the Middleman: Calibrating and Validating a Dynamic Vegetation Model (ED2-PROSPECT5) Using Remotely Sensed Surface Reflectance

Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant processes as well as understanding land–atmosphere feedbacks. In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty. Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified. Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state. We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types. Using a single common set of parameters across all sites, the calibrated model was able to accurately reproduce surface reflectance for sites with highly varied forest composition and structure. However, the calibrated model's predictions of leaf area index (LAI) were less robust, capturing only 46 % of the variability in the observations. Comparing the ED2 radiative transfer model with another two-stream soil–leaf–canopy radiative transfer model commonly used in remote sensing studies (PRO4SAIL) illustrated structural errors in the ED2 representation of direct radiation backscatter that resulted in systematic underestimation of reflectance. In addition, we also highlight that, to directly compare with a two-stream radiative transfer model like EDR, we had to perform an additional processing step to convert the directional reflectance estimates of AVIRIS to hemispherical reflectance (also known as “albedo”). In future work, we recommend that vegetation models add the capability to predict directional reflectance, to allow them to more directly assimilate a wide range of airborne and satellite reflectance products. We ultimately conclude that despite these challenges, using dynamic vegetation models to predict surface reflectance is a promising avenue for model calibration and validation using remote sensing data.

Alexey N Shiklomanov↗

DOE OSTI2021

Low mode implosion symmetry sensitivity in low gas-fill NIF cylindrical hohlraums

Achieving an efficient capsule implosion in National Ignition Facility indirect-drive target experiments requires symmetric hohlraum x-ray drive for the duration of the laser pulse. This is commonly achieved using two-sided two-cone laser irradiation of cylindrical hohlraums that, in principle, can zero the time average of all spherical harmonic asymmetry modes <6 as well as the time dependence of the usually dominant mode 2. In practice, experimental evidence indicates that maintaining symmetric drive becomes limited late in the pulse due to the inward expansion of the hohlraum wall and outward expansion of the capsule ablator plasmas impairing the propagation of the inner-cone laser beams. This effect is enhanced in hohlraums employing low gas-fill, now used almost exclusively as these provide the highest performing implosions and reduce Stimulated Brillouin and Raman backscatter losses, since the gas plasma provides less back pressure to limit blow-in of the hohlraum wall and capsule ablator plasmas. In order to understand this dynamic behavior, we combined multi-keV X-ray imaging of the wall and imploded fuel plasmas as we changed a single parameter at a time: hohlraum gas-fill, laser outer cone picket energy, radius of high density carbon capsules used, and laser beam polar and azimuthal pointing geometry. In this work, we developed a physics-based multi-parameter experimental scaling to explain the results that extend prior scalings and compare those to radiation hydrodynamic simulations to develop a more complete picture of how hohlraum, capsule, and laser parameters affect pole vs equator drive symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Dynamic Modeling of Sequential Service Restoration in Islanded Single Master Microgrids

This paper presents a methodology for generating restoration sequences in distribution systems and microgrids with emphasis on deriving optimal switching times for remote control switches (RCS). Microgrid architecture considered are assumed to be operating in islanded single master mode. Control methods for handling the unbalance inherent in distribution systems using single-phase controlled PQ inverters were developed and incorporated in the proposed method. The problem of considering simultaneously both energy optimization and dynamic stability of the system, which is lacking in existing methods, was addressed by studying the transient behavior of microgrids based on extended EMTP simulations and then approximate and incorporate these behaviors into the optimization formulation for restoration. The approach is based on multi-objective optimization of final restored load and total restoration time. The methodology was studied on a modified IEEE 123 node test feeder. An EMTP simulation in PSCAD of the formed microgrids was used to validate the accuracy of the dynamic and steady-state characteristics of the resulting systems. By considering the transient behavior of DGs during restoration, a better restoration sequence that does not trigger protection relays could be realized.

42 ENGINEERING↗

DOE OSTI2021

Hadronic resonance production and interaction in p-Pb collisions at LHC energies in EPOS3

Using the EPOS3 model with the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) model to describe the hadronic phase, we study the production of short-lived hadronic resonances and the modification of their yields and p T spectra in p-Pb collisions at √sNN= 5.02 TeV. High-multiplicity p-Pb collisions exhibit similar behavior to midperipheral Pb-Pb collisions at LHC energies, and we find indications of a short-lived hadronic phase in p-Pb collisions that can modify resonance yields and p T spectra through scattering processes. The evolution of resonance production is investigated as a function of the system size, which is related to the lifetime of the hadronic phase in order to study the onset and evolution of scattering effects in p-Pb collisions. We also study hadron production separately in the core and corona parts of these collisions and explore how this division affects the total particle yields as the system size increases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Decoding the Desorption Mechanism of 2LiH:1Mg(NH2)2 Using Metal Borohydrides

The complex metal hydride 2LiH:1Mg(NH2)2 has emerged as a promising material for stationary hydrogen storage applications, such as seasonal storage or energy backup systems, due to its high volumetric and gravimetric capacities and robust reversibility. However, its widespread adoption is hindered by sluggish reaction rates, performance degradation upon cycling, and improper end-use cases. To address these problems and better understand the desorption pathway, we used metal borohydrides (MBH4; M = Li–Cs) as chemical probes. A thorough analysis of the bulk behavior of all six materials, including hydrogen cycling experiments, X-ray absorption spectroscopy, FTIR, pXRD, solid-state NMR, and ab initio DFT simulations, shows that the borohydride additives decrease the activation energy of hydrogen release by about 20 kJ/mol for MBH4@2:1 materials versus pristine. Furthermore, more surface-sensitive studies show that the amide-to-imide desorption pathway in these materials, while essentially complete in the bulk, is incomplete in the near-surface region, suggesting an “inverse core–shell” desorption mechanism for amide dehydrogenation to imide. The kinetic enhancements produced by MBH4 additives (M = K, Rb, and Cs) are attributed to the destabilization of the amide N–H bond and interaction with the LiH/Mg(NH2)2 interface to promote H–H bond formation. An inverse core-shell mechanism is also operative in the hydrogen desorption for the 2LiH:1LiNH2 system, suggesting this may be a general feature of amides. Given the fast dehydrogenation rate and large gravimetric capacity, these materials satisfy these requirements for telecom backups and seasonal microgrid storage applications.

Absorption↗

DOE OSTI2021

Effect of water frustration on water oxidation catalysis in the nanoconfined interlayers of layered manganese oxides birnessite and buserite

The role of geometric frustration of water molecules on the rate of water oxidation in the nanoconfined interlayer of manganese-oxide layered materials (birnessite, buserite) is examined in a well-controlled experiment. Calcium buserite is prepared, and used in a split-batch synthetic protocol to prepare calcium birnessite, sodium buserite, and sodium birnessite, and partially dehydrated sodium birnessite. Thus, four samples are prepared in which features effecting catalytic efficiency (defect density, average manganese oxidation state) are controlled, and the main difference is the degree of hydration of the interlayer (two layers of water in buserites vs one layer of water in birnessite). Here, molecular dynamics simulations predict birnessite samples to exhibit geometric water frustration, which facilitates redox catalysis by lowering the Marcus reorganizaiton energy of electron transfer, while buserite samples exhibit traditional intermolecular hydrogen bonding among the two-layer aqeuous region, leading to slower catalytic behavior akin to redox reactions in bulk water. Water oxdiation activity is investigated using chemical and electrochemical techniques, demonstrating and quantifying the role of water frustration in enhancing catalysis. Calculation and experiment demonstrate dehydrated sodium birnessite to be most effective, and calcium buserite the least effective, with a difference in electrocatlytic overpotential of ~750 mV and a ~20-fold difference in turnover number.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Energy Services Interface: Architecture, Requirements, and Agreements

This report continues and extends the development of the Energy Services Interface (ESI) that recently has been led by the U.S. Department of Energy’s Grid Modernization Laboratory Consortium. The ESI is intended to facilitate the coordination of multiple, flexible energy resources that work in tandem to satisfy grid objectives by using performance attributes encoded within an ESI Service Template. The ESI relies on a pair of interfaces, representing the service requestor (the consumer of a service) and a service provider (who manages its resources per the agreed performance attributes) to achieve what is commonly known as a “grid service”. Due to its performance-driven approach, the ESI is hypothesized to be able to satisfy all common grid needs via only six common ESI service types. This report first reviews the status of ESI development, including its fundamental tenets. The report makes three important contributions to ESI development: First, it recognizes the similarity between service level agreements and the contract-like agreements that would be needed between energy service requestors and providers and recommends that ESI service agreements be modeled after the web services agreement specification. Second, whereas prior development efforts had focused on a flat, five-stage lifecycle, this report asserts that the ESI should have three behavioral layers in its architecture—the discovery, agreement, and service layers. Finally, this report offers concrete requirements that should be useful toward the development of the service requestors’ and service providers’ respective communications interfaces.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

A Variable Pressure Multi-Pressure Rail System Design for Agricultural Applications

This paper presents a solution for reducing energy loss in the hydraulic control system of agricultural tractors and their implements. The solution is referred to as a multi-pressure rail (MPR) and provides power to the hydraulic functions following a pressure control logic, as opposed to the traditional flow control logic typical of hydraulic systems used in off-road vehicles. The proposed hydraulic control system allows for elimination of redundant flow control valves in the state-of-the-art system, which cause excessive throttling losses leading to poor overall energy efficiency. Related work on MPR technology targets construction vehicles, where the MPR solution can allow energy recovery during overrunning loads and better engine management. This paper alternatively addresses the case of agricultural applications where functions mostly operate under resistive load conditions with slow dynamics, which offers an opportunity to target throttle losses. For this purpose, the paper introduces a variable pressure control strategy to handle the instantaneous pressure at each rail. To develop both the controller and the hydraulic system architecture, a stationary test rig is conceived and used to validate a numerical simulation model of the MPR system and its control strategy. Particular focus is given to the dynamic behavior of the system during the switches of a function between different pressure rails, which needs to ensure reduced oscillations of the flow provided to each hydraulic function. Once validated, the simulation model is used to predict the energy savings of the MPR solution in an actual application: a 435 hp hydraulic tractor powering a 16-row planter, for which operating features during typical drive cycles were available to the authors. The results show up to 59% total power reduction at the pump shaft, corresponding to 89.8% system efficiency gain.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Modulating the Extent of Anisotropic Cuprophilicity via High Pressure with Piezochromic Luminescence Sensitization

We report metallophilicity has been widely studied as a fundamental supramolecular interaction. However, the extent and directionality thereof remain controversial. A major obstacle lies in the difficulty to separately control the geometry and chemical composition. Herein, we address this challenge by modulating metallophilicity with mechanical pressure. Using a multinuclear Cu(I) complex as model system, we report anomalous anisotropies of (supra)molecular structures, vibrations, and interaction energies upon isotropic compression as well as concomitant (essentially turn-on) piezochromic luminescence enhancement with ~10 3 modulation. The in situ characterizations indicate opposite behaviors of contact distances and cuprophilic interactions for intermolecular vs intramolecular Cu-Cu pairs under pressure. Theoretical calculations break down the attractive and repulsive forces associated with cuprophilicity, its spontaneous 4p-3d hybridization origin, and direction-dependent interaction strength. The use of isotropic mechanical force reveals the intrinsic anisotropy of metallophilicity in multinuclear systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗