Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stabilization method”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Observational evidence against strongly stabilizing tropical cloud feedbacks

In this paper, we present a method to attribute cloud radiative feedbacks to convective processes, using subcloud layer buoyancy as a diagnostic of stable and deep convective regimes. Applying this approach to tropical remote sensing measurements over years 2000–2016 shows that an inferred negative short-term cloud feedback from deep convection was nearly offset by a positive cloud feedback from stable regimes. The net cloud feedback was within statistical uncertainty of the National Center for Atmospheric Research Community Atmosphere Model (CAM5) with historical forcings, with discrepancies in the partitioning of the cloud feedback into convective regimes. Compensation between high-cloud responses to tropics-wide warming in stable and unstable regimes resulted in smaller net changes in high-cloud fraction with warming. In addition, deep convection and associated high clouds set in at warmer temperatures in response to warming, as a consequence of nearly invariant subcloud buoyancy. Lastly, this invariance further constrained the magnitude of cloud radiative feedbacks and is consistent with climate model projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Self‐Assembly of 0D/3D Perovskite Bi‐Layer from a Micro‐Emulsion Ink

Abstract 2D/3D bilayer perovskite synthesized using sequential deposition methods has shown effectiveness in enhancing the stability of perovskite solar devices. However, these approaches present several limitations such as uncontrolled chemical processes, disordered interfacial states, and microscale heterogeneities that can chemically, structurally, and electronically compromise the performance of solar modules. Here, this work demonstrates an emulsion‐based self‐assembly approach using natural lipid biomolecules in a nonionic solution system to form a 0D/3D bilayer structure. The new capping layer is composed of 0D‐entity nanoparticles of perovskite encapsulated by a hydrophobic lipid membrane, analogous to a cell structure, formed through a molecular self‐assembly process. This 0D layer provides a strong water repellent characteristics, optimum interface microstructure, and excellent homogeneity that drives significant enhancement in stability. Solar modules with a large active area of 70 cm 2 fabricated using films comprising of 0D/3D bilayer structure are found to show consistent efficiency of >19% for 2800 h of continuous illumination in the air (60% relative humidity). This emulsion‐based self‐assembly approach is expected to have a transformative impact on the design and development of stable perovskite‐based devices.

0D/3D bilayers↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Boundary Corrections for Kernel Approximation to Differential Operators

The kernel-based approach to operator approximation for partial differential equations has been shown to be unconditionally stable for linear PDEs and numerically exhibit unconditional stability for non-linear PDEs. These methods have the same computational cost as an explicit finite difference scheme but can exhibit order reduction at boundaries. In previous work on periodic domains, order reduction was addressed, yielding high-order accuracy. The issue addressed in this work is the elimination of order reduction of the kernel-based approach for a more general set of boundary conditions. Further, we consider the case of both first and second order operators. To demonstrate the theory, we provide not only the mathematical proofs but also experimental results by applying various boundary conditions to different types of equations. The results agree with the theory, demonstrating a systematic path to high order for kernel-based methods on bounded domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Crack Path Selection in Thermally Loaded Borosilicate/Steel Bibeam Specimen

Here, we have developed a novel specimen for studying crack paths in glass. Under certain conditions, the specimen reaches a state where the crack must select between multiple paths satisfying the K II = 0 condition. This path selection is a simple but challenging benchmark case for both analytical and numerical methods of predicting crack propagation. We document the development of the specimen, using an uncracked and instrumented test case to study the effect of adhesive choice and validate the accuracy of both a simple beam theory model and a finite element model. In addition, we present preliminary fracture test results and provide a comparison to the path predicted by two numerical methods (mesh restructuring and XFEM). The directional stability of the crack path and differences in kink angle predicted by various crack kinking criteria is analyzed with a finite element model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Control of High-Temperature Static and Transient Thermomechanical Behavior of SiMo Ductile Iron by Al Alloying

Abstract Silicon and molybdenum ( SiMo ) ductile iron is commonly used for exhaust manifolds because these components experience thermal cycling in oxidizing environment, which requires resistance to fatigue during transient thermomechanical loads. Previous studies have demonstrated that alloying elements, such as Al , to SiMo ductile iron reduces the amount of surface degradation during static high-temperature exposure. However, deterioration of sphericity of the graphite nodules and a decrease in ductility could affect the tendency of cracking during thermal cycling. In this article, the effect of Al alloying on static and transient thermomechanical behavior of SiMo ductile iron was investigated to optimize the amount of Al alloying. A thermodynamic approach was used to confirm the effect of the Al alloying on the phase transformations in two SiMo cast irons, alloyed by 1.8% Al and 3% Al . These two alloys were cast in a laboratory along with the baseline SiMo ductile iron. Several experimental methods were used to evaluate the dimensional stability, physical properties, static oxidation, and failure resistance during constrained thermal cycling testing to compare their high-temperature capability. Experimental results verified that Al alloying increases the temperature range and decreases volume change during eutectoid transformation, which together with enhancement of oxidation protection improved the dimensional stability. Thermocycling tests showed that the number of cycles to failure depends on the amount of Al alloying and the applied high-temperature exposure during each cycle. SEM/EDX, high-resolution TEM and µCT analysis were used to verify the mechanism resulting from the Al alloying protection. It was shown that an optimal level of Al alloying for balancing oxidation and thermal cracking resistance depends on thermomechanical conditions of application.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Non-equilibrium molecular dynamics simulation to evaluate the effect of confinement on fluid flow in silica nanopores

The flow of fluids in nano-confinement has applications in separations, water purification, medical systems and in the recovery of fluids in petroleum systems. It is believed that fluids do not obey continuum laws when flowing in nanopores (or in confinement). This paper attempts to shed light on various nano confinement effects such as fluid-wall interactions, pore size and molecular geometry, using molecular dynamic simulation. Here, water, hexane, and methanol flow behaviors were simulated for pore diameters ranging from 1 to 8 nm. In addition to the density analysis of the confined fluids, the fluid flows in saturated nanopores were simulated by the sectional flow method. Water and methanol molecules were completely stabilized in the 1 nm pore. Hexane molecules could not enter the 1 nm pore, due to geometric considerations. For the 2 – 8 nm pores, all fluids showed reduced flow rate compared to the Hagen–Poiseuille flow, due to an interfacial molecular layer stabilized at the pore surface. Flow reduction observed in these studies is contrary to significant flow enhancements observed for fluid flow in carbon nanotubes of similar dimensions. Water flow showed almost constant stick length for all the pore sizes. Methanol flow had the largest stick length. Hexane flow was reduced because of overcrowding of molecules at the pore surface. The effect of confinement diminished with an increase in pore diameter.

42 ENGINEERING↗
DOE OSTI2023

Accuracy and precision of ESI-Orbitrap-IRMS observations of hours to tens of hours via reservoir injection

Here, orbitrap isotope ratio mass spectrometry (Orbitrap-IRMS) has recently been applied to high-precision, natural-abundance isotope ratio measurements of a diverse range of compounds, including amino acids, oxyanions, fatty acids, and metals. These measurements can characterize many isotope ratios simultaneously at high (≈1.0‰) precision. In a successful experiment, observed precision will track the shot-noise limit and be limited by experimental time. Some isotope ratios, for example those involving 17 O in organic compounds or multiply-substituted (‘clumped’) isotopologues, require experimental times of hours to tens of hours to achieve desired precision, while current sample introduction techniques focus on observations on the order of seconds to tens of minutes. In this study, we characterize Orbitrap-IRMS performance for three long duration measurements (individual acquisitions ≥1 h and as long as 24 h) using an automated reservoir injection system coupled to a Q Exactive HF Orbitrap with an electrospray ionization (ESI) source. First, we characterize long-term intra-measurement stability through a 24-h long measurement of acetone. We report the following isotope ratios and precisions (as acquisition errors, errors on the observed ratio within this measurement (σ AE )): 13 C/ 12 C (σ AE = 0.07‰), 17 O/ 16 O (σ AE = 1.1‰), 18 O/ 16 O (σ AE = 0.3‰), and 13 C 13 C/ 12 C (σ AE = 0.65‰). The σ AE of each tracks the shot noise limit throughout and is limited by the challenging conditions (high resolution and low numbers of ions per scan) required for 17 O/ 16 O measurement in the presence of 13 C via Orbitrap. Second, we characterize inter-measurement stability via a sequence of seven 75-min analyses of perchlorate. We observe the following ratios and acquisition errors: 37 Cl/ 35 Cl (σ AE = 0.09‰); 17 O/ 16 O (σ AE = 1.6‰); 18 O/ 16 O (σ AE = 0.7‰), 37 Cl 17 O/ 35 Cl 16 O (σ AE = 2.7‰), and 37 Cl 18 O/ 35 Cl 16 O (σ AE = 1.2‰). However, we find that inter-measurement drift between acquisitions limits our accuracy and precision for standardized measurements (i.e., error on reported δ values) to ≈1‰ for the 37 Cl/ 35 Cl measurement. Hence, the benefits of low σ AE may not be fully realized. Third, we demonstrate accuracy via sample/standard comparisons of a methionine sample with 13 C enrichment of ≈20‰ relative to a known standard. Using a sequence of seven 60-min analyses, we recover the following isotope ratios and standardized precisions (i.e., error on reported δ values, denoted propagated acquisition errors, σ PAE ): 33 S/ 32 S (σ PAE = 1.0‰), 34 S/ 32 S (σ PAE = 0.7‰), 15 N/ 14 N (σ PAE = 2.1‰), 2 H/ 1 H (σ PAE = 3.2‰), 13 C/ 12 C (σ PAE = 0.4‰), 18 O/ 16 O (σ PAE = 1.6‰), & 13 C 13 C/ 12 C (σ PAE = 2.8‰) with confirmation of accurate results for the known 13 C/ 12 C and 13 C 13 C/ 12 C enrichments. Together, our results demonstrate the viability of Orbitrap-IRMS for long duration measurements of diverse sample types via an automated reservoir injection system. Inter-measurement stability remains a challenge; we expect our methods to be most applicable to extended measurements of hard-to-observe properties, such as 17 O in organics and clumped isotopologues.

58 GEOSCIENCES↗
DOE OSTI2017

Room temperature biosynthesis of greatly stable fluorescent ZnO quantum dots for the selective detection of Cr{sup 3+} ions

Graphical abstract: Biosynthesized fluorescent ZnO quantum dots are extremely stable and are able to sense up to 2 nM of Cr3+ ions in aqueous solution. - Highlights: • ZnO quantum dotes (QDs) are rapidly synthesized from the bio waste at room temperature. • The capping and reducing agent is from natural source hence biodegradable. • The water soluble fluorescent ZnO quantum dots are highly stable. • The QDs are able to sense 2 nM of Cr{sup 3+} ions in aqueous solution. - Abstract: Synthesis of extremely stable and water soluble quantum dots (QDs) would be smart if the method is rapid and eco-friendly. Here in, we have employed a simple strategy for the preparation of well dispersed, luminescent ZnO quantum dots by using pomegranate peel extract at room temperature. This method excludes the use of strong commercial stabilizing agents, different pH conditions, higher temperature and inert atmosphere during the synthesis. The microscopic study has revealed the ZnO QDs to be of a narrow size with an average size distribution of about 4 nm. The water soluble QDs exhibited a strong emission peak at 490 nm and were stable for twelve months. The emission from ZnO QDs are selectively quenched upon addition of Cr{sup 3+} ions. The interaction between the Cr{sup 3+} and the biomolecules on the exterior of ZnO QDs, enables the system to detect up to 2 nM Cr{sup 3+} ions in aqueous solution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Risk Assessment in a Chemical Laboratory Following an Explosive Incident Involving a Novel Diazonium Compound: Retrospective Analysis and Lessons Learned

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. Furthermore, the steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Computer simulations↗
DOE OSTI2017

Nonadiabatic Ab Initio Molecular Dynamics with the Floating Occupation Molecular Orbital-Complete Active Space Configuration Interaction Method [Non-Adiabatic Ab Initio Molecular Dynamics with Floating Occupation Molecular Orbitals CASCI Method]

Here, we show that the floating occupation molecular orbital complete active space configuration interaction (FOMO-CASCI) method is a promising alternative to the widely used complete active space self-consistent field (CASSCF) method in direct nonadiabatic dynamics simulations. We have simulated photodynamics of three archetypal molecules in photodynamics: ethylene, methaniminium cation, and malonaldehyde. We compared the time evolution of electronic populations and reaction mechanisms as revealed by the FOMO-CASCI and CASSCF approaches. Generally, the two approaches provide similar results. Some dynamical differences are observed, but these can be traced back to energetically minor differences in the potential energy surfaces. We suggest that the FOMO-CASCI method represents, due to its efficiency and stability, a promising approach for direct ab initio dynamics in the excited state.

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

Morphology Control of Carbon-Free Spinel NiCo 2 O 4 Catalysts for Enhanced Bifunctional Oxygen Reduction and Evolution in Alkaline Media

Spinel NiCo 2 O 4 is considered a promising precious metal-free catalyst that is also carbon-free for oxygen electrocatalysis. Current efforts mainly focus on optimal chemical doping and substituent to tune its electronic structures for enhanced activity. Here, we study its morphology control and elucidate the morphology-dependent catalyst performance for bifunctional oxygen reduction reaction (ORR) and oxygen evolution reaction (OER). Three types of NiCo 2 O 4 catalysts with significantly distinct morphologies were prepared using temple-free, Pluronic-123 (P-123) soft, and SiO 2 hard templates, respectively, via hydrothermal methods following by a calcination. While the hard-template yields sphere-like dense structures, soft-template assists the formation of a unique nano-needle cluster assembly containing abundant meso- and macro pores. Furthermore, the effect of morphology of NiCo 2 O 4 on their corresponding bifunctional catalytic performance was systematically investigated. The flower-like nano-needle assembly NiCo 2 O 4 catalyst via the soft template method exhibited the highest catalytic activity and stability for both ORR and OER. In particular, it exhibited an onset and half-wave potentials of 0.94 and 0.82 V vs. RHE, respectively, for the ORR in alkaline media. Although it is still inferior to Pt, the NiCo 2 O 4 represents one of the best ORR catalyst compared to other reported carbon-free oxides. Meanwhile, remarkable OER activity and stability were achieved with an onset potential of 1.48 V and a current density of 15 mA/cm 2 at 1.6 V, showing no activity loss after 20,000 potential cycles (0 to 1.9 V). The demonstrated stability is even superior to Ir for the OER. The morphology-controlled approach provides an effective solution to create a robust 3D architecture with increased surface areas and enhanced mass transfer. More importantly, the soft template can yield high degree of spinel crystallinity with ideal stoichiometric ratios between Ni and Co, thus promoting structural integrity with enhanced electrical conductivity and catalytic properties.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2017

Chemical Potential of Metal Atoms in Supported Nanoparticles: Dependence upon Particle Size and Support

Metal nanoparticles supported on oxide surfaces form the basis for many industrial catalysts and promise to play an ever-increasing role in future energy and environmental technologies. The chemical potential (μ) of the metal atoms in these particles depends strongly on particle size and support and is an important factor that determines their catalytic properties, including their binding strengths to adsorbed reaction intermediates and their long-term stability against sintering. We present here a method for estimating this chemical potential as a function of particle size for different metal/oxide combinations. We show that this chemical potential for late transition metals is well approximated for a particle of diameter D by μ( D ) = [(3γ m – E adh )(1 + D o / D )](2 V m / D ), where γ m is the surface energy of the bulk metal, E adh is the adhesion energy at the bulk metal/oxide interface, and D o is ~1.5 nm, and V m is the molar volume of the bulk metal. We further show that E adh increases with (1) increasing heat of formation of the most stable oxide of the metal from metal gas atoms plus O 2 (gas) per mole of metal atoms, (2) decreasing enthalpy of reduction of the oxide to its next lower oxidation state plus O 2 (gas), per mole of oxygen atoms, and (3) increasing density of surface oxygen atoms on the oxide surface. The linear scaling of E adh with these properties allows estimations of E adh for a variety of metal/oxide combinations. Using this E adh estimate in the above equation with known values for γ m allows estimates of metal chemical potential versus particle size for late transition metals on various oxide supports. WE conclude this will improve our ability to understand structure–property relations in catalysis and design better catalysts.

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

All 2D Material Printed Diodes and Circuits on Paper for Sustainable Electronics

Sustainable electronics aim to reduce environmental impact by using ecofriendly materials, energy-efficient manufacturing, and recyclable components. However, existing approaches rely on complex, resource-intensive methods, rare metals, or nanomaterials with limited stability as well as plastic substrates, raising sustainability issues. Solution-processed two-dimensional (2D) materials offer a promising alternative: water-based and biocompatible conductive, semiconductive, and insulating 2D material inks can be produced with scalable techniques and are suitable for the fabrication of fully printed devices on low-cost and biodegradable paper substrates. However, 2D material only and fully printed diodes on paper have not yet been reported. Here, we demonstrate fully inkjet-printed 2D material-based diodes on paper using metal-insulator-semiconductor and metal-insulator-metal-semiconductor architectures. Water-based graphene and MoS2 inks, prepared by liquid-phase exfoliation, are used for the metallic and insulating films, while electrochemical exfoliation is used to produce the semiconducting MoS2 ink. The highest forward-to-reverse current ratio obtained is 330 (at ±2 V), while the forward current density is 1 mA/cm2 (at 1 V), making the diode performance comparable to the best solution-processed diodes reported so far. However, in contrast to previous works, fabrication occurs entirely at room temperature and ambient pressure, without using high-pressure sputtering, thermal evaporation, and any precious metal ink. The devices maintain stable performance under bending up to strain of 4% over 10,000 cycles. Finally, the diodes are successfully integrated with other 2D-material based electrical components to realize fully printed RC circuits, differentiators, integrators, and AC-to-DC converters onto paper, hence demonstrating the suitability of our approach for sustainable and disposable integrated circuits.

2D materials↗
DOE OSTI2021

Computation-guided optimization of split protein systems

Splitting bioactive proteins into conditionally reconstituting fragments is a powerful strategy for building tools to study and control biological systems. However, split proteins often exhibit a high propensity to reconstitute, even without the conditional trigger, limiting their utility. Current approaches for tuning reconstitution propensity are laborious, context-specific or often ineffective. Here, we report a computational design strategy grounded in fundamental protein biophysics to guide experimental evaluation of a sparse set of mutants to identify an optimal functional window. We hypothesized that testing a limited set of mutants would direct subsequent mutagenesis efforts by predicting desirable mutant combinations from a vast mutational landscape. This strategy varies the degree of interfacial destabilization while preserving stability and catalytic activity. We validate our method by solving two distinct split protein design challenges, generating both design and mechanistic insights. Furthermore, this new technology will streamline the generation and use of split protein systems for diverse applications.

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

Insights into the structure–property–activity relationship in molybdenum-doped octahedral molecular sieve manganese oxides for catalytic oxidation

Hexavalent molybdenum ions substituted into the mixed-valent framework of octahedral molecular sieve manganese oxides (1–10 mol% Mo–K-OMS-2) were systematically prepared via a single-step refluxing method. Here, the structure, composition, morphology, thermal stability, and textural and redox properties of the Mo–K-OMS-2 materials were characterized by various experimental techniques. Density functional theory (DFT) calculations were performed to study the electronic properties of K-OMS-2 materials, including the influence of molybdenum on such properties. The catalytic activity of K-OMS-2 towards CO oxidation dramatically increased with Mo dopant concentration, until an optimum incorporation of 5 mol% was reached, giving full CO conversion at 120 °C (1 vol% CO, 1 vol% O2). More importantly, the activity of the 5% Mo–K-OMS-2 catalyst was (reversibly) inhibited by water vapour (3%) in the feed only at low temperatures and exhibited better tolerance and stability during long-run experiments (28 h) as compared to undoped K-OMS-2. The superior performance of Mo–K-OMS-2 catalyst for CO oxidation can be attributed to the following factors: (1) morphological evolution from nanofibers to nanospheres, as well as the slightly distorted structure, generated a sharp increase in the specific surface area up to 258 m2 g–1, hence, increased the number of surface active sites; (2) improved mobility of reactive oxygen species at the surface and enhanced redox properties; and (3) tailored electronic properties promoted more active sites for oxidation reactions, which is supported by DFT calculations.

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

The many-body expansion for metals. I. The alkaline earth metals Be, Mg, and Ca

We examine the many-body expansion (MBE) for alkaline earth metal clusters, Be n , Mg n , Ca n ( n = 4, 5, 6), at the Møller–Plesset second order perturbation theory, coupled-cluster singles and doubles with perturbative triples, multi-reference perturbation theory, and multi-reference configuration interaction levels of theory. The magnitude of each term in the MBE is evaluated for several geometrical configurations. We find that the behavior of the MBE for these clusters depends strongly on the geometrical arrangement and, to a lesser extent, on the level of theory used. Another factor that affects the MBE is the in situ (ground or excited) electronic state of the individual atoms in the cluster. For most geometries, the three-body term is the largest, followed by a steady decrease in absolute energy for subsequent terms. Though these systems exhibit non-negligible multi-reference effects, there was little qualitative difference in the MBE when employing single vs multi-reference methods. Useful insights into the connectivity and stability of these clusters have been drawn from the respective potential energy surfaces and quasi-atomic orbitals for the various dimers, trimers, and tetramers. Through these analyses, we investigate the similarities and differences in the binding energies of different-sized clusters for these metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗