Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex geometry”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Bis(imidazole)-based cysteine labeling tool for metalloprotein assembly

Precise metal-protein coordination by design remains a considerable challenge. Polydentate, high-metal-affinity protein modifications, both chemical and recombinant, can enable metal localization. However, these constructs are often bulky, conformationally and stereochemically ill-defined, or coordinately saturated. Here, we expand the biomolecular metal-coordination toolbox with the irreversible attachment to cysteine of bis(1-methylimidazol-2-yl)ethene (“BMIE”), which generates a compact imidazole-based metal-coordinating ligand. Conjugate additions of small-molecule thiols (thiocresol and N-Boc-Cys) with BMIE confirm general thiol reactivity. The BMIE adducts are shown to complex the divalent metal ions Cu ++ and Zn ++ in bidentate (N 2 ) and tridentate (N 2 S*) coordination geometries. Cysteine-targeted BMIE modification (>90% yield at pH 8.0) of a model protein, the S203C variant of carboxypeptidase G2 (CPG2), measured with ESI-MS, confirms its utility as a site-selective bioconjugation method. ICP-MS analysis confirms mono-metallation of the BMIE-modified CPG2 protein with Zn ++ , Cu ++ , and Co ++ . EPR characterization of the BMIE-modified CPG2 protein reveals the structural details of the site selective 1:1 BMIE-Cu ++ coordination and symmetric tetragonal geometry under physiological conditions and in the presence of various competing and exchangeable ligands (H 2 O/HO – , tris, and phenanthroline). An X-ray protein crystal structure of BMIE-modified CPG2-S203C demonstrates that the BMIE modification is minimally disruptive to the overall protein structure, including the carboxypeptidase active sites, although Zn ++ metalation could not be conclusively discerned at the resolution obtained. The carboxypeptidase catalytic activity of BMIE-modified CPG2-S203C was also assayed and found to be minimally affected. Finally, these features, combined with ease of attachment, define the new BMIE-based ligation as a versatile metalloprotein design tool, and enable future catalytic and structural applications.

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

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Synthesis, Properties, and Electrochemical Proton Reduction of a Homoleptic Tetrathiolato Ni-Site Model of [NiFe]-Hydrogenase

[NiFe]-hydrogenase enzymes process H 2 at a nonplanar tetracysteinato-Ni site, the sole participator in proton binding/redox chemistry during turnover. With the objective of assessing whether a simple tetrahedral/ tetrathiolato-Ni 2+ core could promote H 2 evolution reaction (HER), we synthesized (Et 4 N) 2 [Ni(S-p-CF 3 −Ph) 4 ] (1) employing para-trifluoromethylbenzenethiolate ( − S-p-CF 3 −Ph) as a Ni-site analog of [NiFe]-hydrogenase. Spectroscopic measurements and X-ray crystallography confirm the distorted tetrahedral geometry of 1. Dissolution of 1 results in partial thiolate dissociation and formation of S,S-bridged complexes such as (Et 4 N) 2 [Ni 2 (S-p-CF 3 −Ph) 6 ] (3) among other ill-defined species. Dissociation is further accelerated in the presence of Brønsted acids, complicating the assessment of 1 for proton reduction. However, this dissociation/proton instability is suppressed in the presence of additional thiolate ligand to ensure tetrahedral/tetrathiolato 1 persists in solution. Electrochemical HER activity was evaluated by monitoring the current response of an MeCN solution of 1/ excess thiolate after sequential titration with a weak Brønsted acid (acetic acid). The results suggest that 1 is a modest electrocatalyst for the HER with a turnover frequency of 14.5 ± 3.6 s −1 and an overpotential of 0.72 ± 0.02 V. Control experiments and supplementary DFT computations indicate that 1, or a species derived from 1, is responsible for the HER and suggest an ECCE-type mechanism.

Evolution reactions↗
DOE OSTI2022

Efficient Intermolecular Charge Transport in π-Stacked Pyridinium Dimers Using Cucurbit[8]uril Supramolecular Complexes

Intermolecular charge transport through π-conjugated molecules plays an essential role in biochemical redox processes and energy storage applications. In this work, we observe highly efficient intermolecular charge transport upon dimerization of pyridinium molecules in the cavity of a synthetic host (cucurbit[8]uril, CB[8]). Stable, homoternary complexes are formed between pyridinium molecules and CB[8] with high binding affinity, resulting in an offset stacked geometry of two pyridiniums inside the host cavity. The charge transport properties of free and dimerized pyridiniums are characterized using a scanning tunneling microscope-break junction (STM-BJ) technique. Our results show that π-stacked pyridinium dimers exhibit comparable molecular conductance to isolated, single pyridinium molecules, despite a longer transport pathway and a switch from intra- to intermolecular charge transport. Control experiments using a CB[8] homologue (cucurbit[7]uril, CB[7]) show that the synthetic host primarily serves to facilitate dimer formation and plays a minimal role on molecular conductance. Molecular modeling using density functional theory (DFT) reveals that pyridinium molecules are planarized upon dimerization inside the host cavity, which facilitates charge transport. In addition, the π-stacked pyridinium dimers possess large intermolecular LUMO–LUMO couplings, leading to enhanced intermolecular charge transport. Altogether, this work demonstrates that supramolecular assembly can be used to control intermolecular charge transport in π-stacked molecules.

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

Surface deformation coupled with self-organized pattern on a liquid anode of 1 atm DC glow discharge

Intricate, self-organized plasma structures observed above the surface of a liquid anode of atmospheric DC glow discharge were found to give rise to coherent, organized surface deformation and mechanical wave formation at the plasma–liquid interface. This new phenomenon indicates that the liquid is closely coupled to the plasma by the anode sheath’s electrohydrodynamic (EHD) force. A scientific question then arises: Do surface perturbations, coupled with the nonuniform surface charge distribution, enhance the electric field and induce self-organization? Using the reflective background-oriented schlieren technique, the liquid surface profile under the plasma pattern was measured for the first time. The results show that surface distortions are driven by the repulsive Coulomb force of nonuniform net-negative surface charge acted by the anode sheath field. The impacts of various operating parameters on the patterns and surface waves were examined, revealing the significance of gas heating and liquid charge relaxation time in the pattern formation mechanism. Time-resolved dynamics of a pulsed DC discharge indicated that the surface deformation only appeared after the establishment of plasma patterns. Statistically, the surface wave under the plasma has high wave numbers (8000–16000 m −1 ) and small amplitudes ($<$ 10 µm), generally found in the capillary wave regime. Yet the motion of surface deformations is in tandem with the plasma pattern and exhibits a nondispersive nature of constant phase velocity (0.1–0.4 m s −1 ), suggesting the dominant role of EHD force over the surface tension in the observed surface wave. These results indicate that the deformed liquid surface is driven by the EHD force of the plasma sheath. Although they share a similar geometry, the deformation and wave dynamics of the liquid surface do not stimulate a plasma pattern. Importantly, the complex EHD coupling in the plasma–liquid system raises awareness and new challenges for plasma control engineering, and the quantitative characteristics of the nondispersive surface wave are informative for advancing relevant theory and modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Controlling quantum many-body dynamics in driven Rydberg atom arrays

Dynamic stabilization of an array Large-scale systems comprising one-dimensional chains and two-dimensional arrays of excited atoms held in a programmable optical lattice are a powerful platform with which to simulate emergent phenomena. Bluvstein et al. built an array of up to 200 Rydberg atoms and subjected the system to periodic excitation. Under such driven excitation, they found that the array of atoms stabilized, freezing periodically into what looked like time crystals. Understanding and controlling the dynamic interactions in quantum many-body systems lies at the heart of contemporary condensed matter physics and the exotic phenomena that can occur. Science , this issue p. 1355

Science & Technology - Other Topics↗
DOE OSTI2020

Understanding Bifacial PV Modeling: Raytracing and View Factor Models

As part of the PV Magazine Webinar 'View-factor vs. ray tracing - which bifacial modelling techniques should you use?,' the second part lead by NREL will explore the different rear-irradiance calculation software, which fall into two categories: view-factor and ray-tracing models. View factor models assume isotropic scattering of reflected rays, allowing for calculation of irradiance by integration. Due-diligence software such as PVSyst or SAM use the view-factor model. Ray-tracing models simulate multipath reflection and absorption of individual rays entering a scene. Raytracing software such as bifacial_radiance, which is the only available open-source toolkit, offers the possibility of reproducing complex scenes, including shading or finite-system edge effects. This can be useful to evaluate specific system geometries, shading and obstructions, areas of different albedo and electrical mismatch, all topics of interest to the prediction of bifacial performance. Model agreement for view factors and bifacial_radiance software is improving and has been shown to be between 2% (absolute) when compared with measured results.

14 SOLAR ENERGY↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2023

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗
DOE OSTI2020

Connecting particle interactions to agglomerate morphology and rheology of boehmite nanocrystal suspensions

Rheology imposes significant challenges on processing of complex suspensions such as nuclear waste slurries at the Hanford and Savannah River sites. Understanding rheology connecting to microstructures and underlying particle interactions in complex slurries is therefore important for both fundamental knowledge and practical applications. Here, we use suspensions of aluminum oxyhydroxide minerals in the form of boehmite as an analog of the radioactive waste slurry to gain physical insights on the correlation between particle interactions, microstructures, and slurry rheology. Specifically, we use a combination of Couette rheometry and small-angle scattering techniques (independently and simultaneously) to understand how the slurry microstructure changes under flow and how these structural changes manifest themselves in the bulk rheology of the suspensions. Our experiments show that the boehmite slurries are thixotropic, with the rheology and structure of the suspensions changing with increasing exposure to flow. In the slurries, particle aggregates begin as loose, system-spanning clusters, but exposure to moderate shear rates causes the aggregates to irreversibly consolidate into denser clusters of finite size. The microstructural changes directly influence the rheological properties of the slurries such as viscosity and viscoelasticity. More importantly, our study shows that solution pH affects the amount of structural rearrangement and the kinetics of the rearrangement process, with an increase in pH leading to faster and more dramatic changes in the bulk rheology. Such dynamic microstructural changes and resultant rheology are understood via correlations between particle interactions and strength of particle network, coupled with surface chemistry and anisotropic nature of particle interactions. Nearly identical structural changes are also observed in Poiseuille flow geometries, implying that the observed changes are relevant in the pipe flow conditions present during waste processing.

Weston, Javen S.↗
DOE OSTI2020

Enhancing nuclear data validation analysis by using machine learning

In this work, we demonstrate how machine learning (ML) techniques can be used as an auxiliary tool for nuclear data validation analysis. The output of the ML analysis can inform evaluators and validators on the quality, or lack thereof, of specific nuclear data and benchmark experiments with respect to simulating these experimental benchmark values. To this end, measured and simulated effective neutron multiplication factors, k eff , of 875 selected ICSBEP critical assemblies and the latter's sensitivities with respect to nuclear data as well as benchmarks' features (e.g., material nuclide, core geometry) are used as input for a random forest (RF) regression model. The RF is used to encode the complex inter-dependencies between thousands of nuclear data features (i.e., sensitivity profiles and aspects describing the measurements) and 875 simulated and experimental k eff values in order to understand what nuclear data features are most informative for predicting bias. The complexity of relationships and high-dimensional space is difficult-to-impossible to search fully using simply expert judgment. As a first proof-of-concept—step, it is demonstrated that this technique is able to correctly trace large discrepancies between simulated and experimental k eff back to fabricated shortcomings in nuclear data that were perturbed to simulated k eff values. In a second, real-case scenario, step, the RF algorithm is used to validate the ENDF/B-VIII.0 library in comparison to ENDF/B-VII.1 nuclear data. One case is showcased where the chosen ML algorithms highlighted nuclear data (the 19 F(n,inl) cross section from 0.4–0.9 MeV) that are shown to be problematic by comparing them to associated differential experimental data and nuclear data from other libraries. In addition to that, it is shown that the RF results point towards poor benchmark experiments and associated underestimated uncertainties (e.g., the PU-SOL-THERM-028 series). However, using the RF algorithm for validating nuclear data with respect to k eff is currently limited to pinpointing groups of questionable nuclear data due to the inherent correlations between features introduced by the nuclear data themselves and how k eff is simulated. Due to this, we recommend that the ML methods presented be used to augment—rather than replace—the expert knowledge of evaluators and validators.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

In situ imaging of three dimensional freeze printing process using rapid x-ray synchrotron radiography

Three dimensional freeze printing (3DFP) combines the advantages of freeze casting and additive manufacturing to fabricate multifunctional aerogels. Freeze casting is a cost-effective, efficient, and versatile method capable of fabricating micro-scale porous structures inside the aerogels for many different applications. The 3DFP provided the capability of fabricating highly customized geometries with controlled microporous structures as well. However, there are still many unexplained phenomena and features because of the complexity of post-processes and indirect observation methods. This study demonstrates the design and construction of the in situ imaging systems, which use the x-ray synchrotron radiography to observe freeze casting and 3DFP processes. With the advantages provided by the in situ x-ray imaging techniques, the ice crystal growth with its unique lamellar structures can be observed during the freeze casting process. Here, the movement of freeze front, material deposition, and growth of ice crystals can also be visualized during the inkjet-based 3DFP process.

3D printing↗
DOE OSTI2022

Effect of insulator length and fill pressure on filamentation and neutron production in a 4.6 kJ dense plasma focus

Optimization of neutron yields from dense plasma focus devices is a complex multi-faceted challenge that necessitates the prudent selection of mechanical constraints such as the electrode and insulator geometries. Here, the neutron yield is found to significantly depend on the insulator length. As the length of the insulator increases, the exposed anode length traveled by the sheath during the run-down phase decreases. This suggests an increase in the optimal fill pressure with increasing insulator length to maintain the pinch time near peak current. However, in the present study, the opposite trend is observed—the optimal fill pressure for neutron production decreases with increasing insulator length. Optical probing of the sheath from run-down to the pinch reveals significant plasma filamentation with increasing pressure and a dependence of insulator length on filamentation onset. A direct consequence of increased filamentation is a reduction in mass sweeping efficiency, directly quantified as a function of fill pressure for the first time.

Hahn, E. N. (ORCID:0000000223050532)↗
DOE OSTI2020

Charge density waves on a half-filled decorated honeycomb lattice

Tight binding models like the Hubbard Hamiltonian are most often explored in the context of uniform intersite hopping $\textit{t}$. The electron-electron interactions, if sufficiently large compared to this translationally invariant $\textit{t}$, can give rise to ordered magnetic phases and Mott insulator transitions, especially at commensurate filling. The more complex situation of nonuniform $\textit{t}$ has been studied within a number of situations, perhaps most prominently in multiband geometries where there is a natural distinction of hopping between orbitals of different degree of overlap. In this paper we explore related questions arising from the interplay of multiple kinetic energy scales and electron-phonon interactions. Specifically, we use determinant quantum Monte Carlo (DQMC) to solve the half-filled Holstein Hamiltonian on a “decorated honeycomb lattice,” consisting of hexagons with internal hopping $\textit{t}$ coupled together by $\textit{t'}$. This modulation of the hopping introduces a gap in the Dirac spectrum and affects the nature of the topological phases. Here, we determine the range of $\textit{t/t'}$ values which support a charge density wave phase about the Dirac point of uniform hopping $\textit{t = t'}$, as well as the critical transition temperature $T_c$. The QMC simulations are compared with the results of mean field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Micro-architected material design for mechanical response

Rapid advances in additive manufacturing (AM) have enabled the creation of micro-architected materials—also known as mechanical metamaterials—with unprecedented control over fine-scale geometries and arrangements of multiple material constituents. These “materials” can achieve unique and extraordinary effective mechanical properties through their complex architectures rather than composition alone. A key challenge is to design for these bespoke effective mechanical responses within the constraints of available AM techniques (i.e., given a set of desired effective properties), identify a (often nonunique) micro-architecture and selection of material constituents that achieves them. Two main strategies have emerged. Gradient-based methods use sensitivity analysis to iteratively refine candidate designs, while data-driven methods learn micro-architecture-constituent relationships from existing examples to propose new designs. This article reviews these design approaches for micro-architected materials with tailored mechanical responses that can be fabricated by AM as well as their applications.

Spadaccini, Christopher M [Lawrence Livermore Nati↗
DOE OSTI2020

Integrable deformation of $\mathbb{C}$P n and generalised Kähler geometry

We build on the results of [1] for generalised frame fields on generalised quotient spaces and study integrable deformations for CP n . In particular we show how, when the target space of the Principal Chiral Model is a complex projective space, a two-parameter deformation can be introduced in principle. The second parameter can however be removed via a diffeomorphism, which we construct explicitly, in accordance with the results stemming from a thorough integrability analysis we carry out. We also elucidate how the deformed target space can be seen as an instance of generalised Kähler, or equivalently bi-Hermitian, geometry. In this respect, we find the generic form of the pure spinors for CP n and the explicit expression for the generalised Kähler potential for n = 1, 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Zinc-based cyclens containing pyridine and cross-bridges: X-ray and DFT structures, Lewis acidity, gas-phase acidity, and p K a values

Herein it is well known that catalytic centers containing the zinc(II) ion can act as both Lewis and Bronsted-Lowry acids. In addition to coordination number, the Zn coordination geometry can also strongly impact the acidity of the active site, no matter what measure of acidity is considered. Herein, we report the first pentacoordinate, zinc-ammonia complex containing a pyridine based tetraazamacrocycle, [(pyclen)Zn(NH 3 )](PF 6 ) 2 , that has applications in Lewis acid catalysis. From this structure, we obtain binding energies and acidities for a series of related pyclen and cross-bridged cyclen type tetraazamacrocycles comprising the pentacoordinate N 4 Zn(II)–OH 2 entity, collectively referred to as [(R-cyclen)Zn–OH 2 ] 2+ . Results from gas- and aqueous-phase density functional theory (M05-2X) and ab initio (MP2) calculations reported herein demonstrate that molecular geometry has a substantial influence on both the Zn–OH x binding strengths and the deprotonation energy of coordinated H 2 O, but not necessarily in a predicable way. While generally more constrained N–Zn–N coordination leads to greater Zn–OH 2 binding energies (Lewis acidities), the corresponding Lewis acidities of the complexes don’t always correlate with the Zn–OH 2 bond lengths nor the (Bronsted) acidity of the coordinated H 2 O. Additionally, the order of the Lewis acid strength of the [(R-cyclen)Zn(II)] complexes changes as the basic –OH 2 ligand is replaced with its –OH counterpart.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗