Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

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,621 records · Page 90

A driven similarity renormalization group approach to quantum many-body problems

Applications of the similarity renormalization group (SRG) approach [F. Wegner, Ann. Phys. 506, 77 (1994) and S. D. Głazek and K. G. Wilson, Phys. Rev. D 49, 4214 (1994)] to the formulation of useful many-body theories of electron correlation are considered. In addition to presenting a production-level implementation of the SRG based on a single-reference formalism, a novel integral version of the SRG is reported, in which the flow of the Hamiltonian is driven by a source operator. It is shown that this driven SRG (DSRG) produces a Hamiltonian flow that is analogous to that of the SRG. Compared to the SRG, which requires propagating a set of ordinary differential equations, the DSRG is computationally advantageous since it consists of a set of polynomial equations. The equilibrium distances, harmonic vibrational frequencies, and vibrational anharmonicities of a series of diatomic molecules computed with the SRG and DSRG approximated with one- and two-body normal ordered operators are in good agreement with benchmark values from coupled cluster with singles, doubles, and perturbative triples. Particularly surprising results are found when the SRG and DSRG methods are applied to C{sub 2} and F{sub 2}. In the former case, both methods fail to converge, while in the latter case an unbound potential energy curve is obtained. A modified commutator approximation is shown to correct these problems in the case of the DSRG method.

74 ATOMIC AND MOLECULAR PHYSICS↗

Optimal Water-Power Flow Problem: Formulation and Distributed Optimal Solution

This paper formalizes an optimal water-power flow (OWPF) problem to optimize the use of controllable assets across power and water systems while accounting for the couplings between the two infrastructures. Tanks and pumps are optimally managed to satisfy water demand while improving power grid operations; {for the power network, an AC optimal power flow formulation is augmented to accommodate the controllability of water pumps.} Unfortunately, the physics governing the operation of the two infrastructures and coupling constraints lead to a nonconvex (and, in fact, NP-hard) problem; however, after reformulating OWPF as a nonconvex, quadratically-constrained quadratic problem, a feasible point pursuit-successive convex approximation approach is used to identify feasible and optimal solutions. In addition, a distributed solver based on the alternating direction method of multipliers enables water and power operators to pursue individual objectives while respecting the couplings between the two networks. The merits of the proposed approach are demonstrated for the case of a distribution feeder coupled with a municipal water distribution network.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Exact relativistic expressions for polarization of incoherent Thomson scattering

We present a derivation of the degree of polarization for incoherent Thomson scattering (TS) using Mueller matrix formalism. An exact analytic solution is obtained for spectrum-integrated matrix elements. The solution is valid for the full range of incident polarizations, scattering angles, and electron thermal motion from non-relativistic to ultra-relativistic. It is based on a newly developed theoretical model, a finite transit time (FTT) correction to previous theoretical work on TS polarization. The Mueller matrix elements are substantially different from previous calculations without the FTT correction, even to the lowest linear order in T e =m e c 2 <<1. Mathematically, the derivation is a unique example of fully analytical integration of the 3D scattering operator over a relativistic Maxwellian distribution function; experimentally, the results have application to the use of the polarization properties of Thomson scattered light as a method of electron temperature measurement. The results can also be used as a reliable tool for benchmarking and verification of numerical codes for frequency resolved properties of TS polarization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generalization of the Kohn-Sham system that can represent arbitrary one-electron density matrices

Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of noninteracting particles, is the workhorse of the theory. The particular form of the Kohn-Sham wave function admits only idempotent one-electron density matrices whereas wave functions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept and introducing a suitable dot product as well as a probability density, a noninteracting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans' theorem are available. Lastly, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single-determinant wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

97 MATHEMATICS AND COMPUTING↗

Orbits for 18 Visual Binaries and Two Double-line Spectroscopic Binaries Observed with HRCAM on the CTIO SOAR 4 m Telescope, Using a New Bayesian Orbit Code Based on Markov Chain Monte Carlo

We present orbital elements and mass sums for 18 visual binary stars of spectral types B to K (five of which are new orbits) with periods ranging from 20 to more than 500 yr. For two double-line spectroscopic binaries with no previous orbits, the individual component masses, using combined astrometric and radial velocity data, have a formal uncertainty of ∼0.1 M{sub ⊙}. Adopting published photometry and trigonometric parallaxes, plus our own measurements, we place these objects on an H-R diagram and discuss their evolutionary status. These objects are part of a survey to characterize the binary population of stars in the Southern Hemisphere using the SOAR 4 m telescope+HRCAM at CTIO. Orbital elements are computed using a newly developed Markov chain Monte Carlo (MCMC) algorithm that delivers maximum-likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate the uncertainty of our derived parameters in a robust way. For spectroscopic binaries, using our approach, it is possible to derive a self-consistent parallax for the system from the combined astrometric and radial velocity data (“orbital parallax”), which compares well with the trigonometric parallaxes. We also present a mathematical formalism that allows a dimensionality reduction of the feature space from seven to three search parameters (or from 10 to seven dimensions—including parallax—in the case of spectroscopic binaries with astrometric data), which makes it possible to explore a smaller number of parameters in each case, improving the computational efficiency of our MCMC code.

79 ASTRONOMY AND ASTROPHYSICS↗

Synthesis and magnetic structure of the layered manganese oxide selenide Sr{sub 2}MnO{sub 2}Ag{sub 1.5}Se{sub 2}

The synthesis of a high-purity sample of the layered oxide selenide Sr{sub 2}MnO{sub 2}Ag{sub 1.5}Se{sub 2} is reported. At ambient temperature it crystallises in the space group I4/mmm with two formula units in the unit cell and lattice parameters a=4.08771(1) Å, c=19.13087(8) Å. The compound displays mixed-valent manganese in a formal oxidation state close to +2.5 and powder neutron diffraction measurements reveal that below the Néel temperature of 63(1) K this results in an antiferromagnetic structure which may be described as A-type, modelled in the magnetic space group P{sub I}4/mnc (128.410 in the Belov, Neronova and Smirnova (BNS) scheme) in which localised Mn moments of 3.99(2) μ{sub B} are arranged in ferromagnetic layers which are coupled antiferromagnetically. In contrast to the isostructural compound Sr{sub 2}MnO{sub 2}Cu{sub 1.5}S{sub 2}, Sr{sub 2}MnO{sub 2}Ag{sub 1.5}Se{sub 2} does not display long range ordering of coinage metal ions and vacancies, nor may significant amounts of the coinage metal readily be deintercalated using soft chemical methods. - Graphical abstract: Sr{sub 2}MnO{sub 2}Ag{sub 1.5}Se{sub 2} containing mixed valent Mn ions undergoes magnetic ordering with ferromagnetic coupling within MnO{sub 2} sheets and antiferromagnetic coupling between MnO{sub 2} sheets. - Highlights: • High purity sample of Sr{sub 2}MnO{sub 2}Ag{sub 1.5}Se{sub 2} obtained. • Magnetic structure determined. • Compared with related mixed-valent manganite oxide chalcogenides.

36 MATERIALS SCIENCE↗

Ten-Year ALARA Review of Dosimetry Results 1 January 2008 through 31 December 2017

Dosimetry results from 1 January 2008 through 31 December 2017 were reviewed to demonstrate that radiation protection methods used at Sandia National Laboratories are compliant with regulatory limits and consistent with the philosophy to keep exposures to radiation As Low As is Reasonably Achievable (ALARA). Personnel dosimetry (external and internal) and environmental thermoluminescent dosimeter results were reviewed for the Sandia National Laboratories in New Mexico, California, and Nevada. ALARA is a philosophical approach to radiation protection by managing and controlling radiation exposures (individual and collective) to the work force and to the public to levels that are As Low As is Reasonably Achievable taking social, technical, economic, practical, and public policy considerations into account. ALARA is not a dose limit but a planning tool with the objective to keep doses below applicable dose limits As Low As is Reasonably Achievable. In the case of Sandia National Laboratories, formal ALARA goals are not needed since Collective and individual doses are well below applicable limits through operational ALARA practices implemented during Work Planning and Control activities.

61 RADIATION PROTECTION AND DOSIMETRY↗

Controlled calculation of the thermal conductivity for a spinon Fermi surface coupled to a U(1) gauge field

Motivated by recent transport measurements on the candidate spin-liquid phase of the organic triangular lattice insulator EtMe{sub 3}Sb[Pd(dmit){sub 2}]{sub 2}, we perform a controlled calculation of the thermal conductivity at intermediate temperatures in a spin liquid system where a spinon Fermi surface is coupled to a U(1) gauge field. The present computation builds upon the double expansion approach developed by Mross et al. (2010) for small ϵ=z{sub b}−2 (where z{sub b} is the dynamical critical exponent of the gauge field) and large number of fermionic species N. Using the so-called memory matrix formalism that most crucially does not assume the existence of well-defined quasiparticles at low energies in the system, we calculate the temperature dependence of the thermal conductivity κ of this model due to non-critical Umklapp scattering of the spinons for a finite N and small ϵ. Then we discuss the physical implications of such theoretical result in connection with the experimental data available in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electronic Conduction through Atomic Chains, Quantum Well and Quantum Wire

Charge transport is dynamically and strongly linked with atomic structure, in nanostructures. We report our ab-initio calculations on electronic transport through atomic chains and the model calculations on electron-electron and electron-phonon scattering rates in presence of random impurity potential in a quantum well and in a quantum wire. We computed synthesis and ballistic transport through; (a) C and Si based atomic chains attached to metallic electrodes, (b) armchair (AC), zigzag (ZZ), mixed, rotated-AC and rotated-ZZ geometries of small molecules made of 2S, 6C and 4H atoms attaching to metallic electrodes, and (c) carbon atomic chain attached to graphene electrodes. Computed results show that synthesis of various atomic chains are practically possible and their transmission coefficients are nonzero for a wide energy range. The ab-initio calculations on electronic transport have been performed with the use of Landauer-type scattering formalism formulated in terms of Grben's functions in combination with ground-state DFT. The electron-electron and electron-phonon scattering rates have been calculated as function of excitation energy both at zero and finite temperatures for disordered 2D and 1D systems. Our model calculations suggest that electron scattering rates in a disordered system are mainly governed by effective dimensionality of a system, carrier concentration and dynamical screening effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hybrid framework for the simulation of stochastic chemical kinetics

Stochasticity plays a fundamental role in various biochemical processes, such as cell regulatory networks and enzyme cascades. Isothermal, well-mixed systems can be modelled as Markov processes, typically simulated using the Gillespie Stochastic Simulation Algorithm (SSA) [25]. While easy to implement and exact, the computational cost of using the Gillespie SSA to simulate such systems can become prohibitive as the frequency of reaction events increases. This has motivated numerous coarse-grained schemes, where the “fast” reactions are approximated either using Langevin dynamics or deterministically. While such approaches provide a good approximation when all reactants are abundant, the approximation breaks down when one or more species exist only in small concentrations and the fluctuations arising from the discrete nature of the reactions become significant. This is particularly problematic when using such methods to compute statistics of extinction times for chemical species, as well as simulating non-equilibrium systems such as cell-cycle models in which a single species can cycle between abundance and scarcity. In this paper, a hybrid jump-diffusion model for simulating well-mixed stochastic kinetics is derived. It acts as a bridge between the Gillespie SSA and the chemical Langevin equation. For low reactant reactions the underlying behaviour is purely discrete, while purely diffusive when the concentrations of all species are large, with the two different behaviours coexisting in the intermediate region. A bound on the weak error in the classical large volume scaling limit is obtained, and three different numerical discretisations of the jump-diffusion model are described. The benefits of such a formalism are illustrated using computational examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Global reverse supply chain design for solid waste recycling under uncertainties and carbon emission constraint

Highlights: • We propose a robust model for global reverse supply chain design. • The uncertain waste collection rate, currency exchange rate, and the maritime cost are incorporated. • The total emission and cost are affected by both the uncertainties and the global carbon cap. • The waste products flows between countries vary with the levels of uncertainties and the global carbon cap. • There is an opportunity to reduce the maritime emission. - Abstract: The emergence of concerns over environmental protection, resource conservation as well as the development of logistics operations and manufacturing technology has led several countries to implement formal collection and recycling systems of solid waste. Such recycling system has the benefits of reducing environmental pollution, boosting the economy by creating new jobs, and generating income from trading the recyclable materials. This leads to the formation of a global reverse supply chain (GRSC) of solid waste. In this paper, we investigate the design of such a GRSC with a special emphasis on three aspects; (1) uncertainty of waste collection levels, (2) associated carbon emissions, and (3) challenges posed by the supply chain’s global aspect, particularly the maritime transportation costs and currency exchange rates. To the best of our knowledge, this paper is the first attempt to integrate the three above-mentioned important aspects in the design of a GRSC. We have used mixed integer-linear programming method along with robust optimization to develop the model which is validated using a sample case study of e-waste management. Our results show that using a robust model by taking the complex interactions characterizing global reverse supply chain networks into account, we can create a better GRSC. The effect of uncertainties and carbon constraints on decisions to reduce costs and emissions are also shown.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SU-F-E-17: A Dedicated Teaching and Research Linac as a Stepping Stone to Introduce Medical Physics to Students

Purpose: This work describes how a non-clinical, research and teaching Linac is used as an extremely motivating and exciting way to introduce students to medical physics. Methods: The dedicated facility was inaugurated in 2014. The facility is composed of a fully equipped and functional state-of-the-art Varian TrueBeam Linac and a complete set of physics instruments and QA phantoms for the Linac and onboard imaging. The Linac bunker and treatment console are oversized such that a class of 12–15 can comfortably fit, seated if needed for longer sessions. A 3cr undergraduate laboratory course that includes medical imaging, x-ray source characterization (mAs, kVp, and filtering) and many others including an introductory Linac laboratory was created. The latter is composed of one general 4-hours session and a weekly 4-hours session for teams of two students. The general session includes a hands-on presentation of the Linac, its environment and a formal safety and radiation protection course (with an exam). Results: Since the winter of 2015, senior undergraduate (total of 15) pursuing either the medical physics or the biomedical engineering tracks can register. At the Linac, the students are allowed full control of the experiments, including set-up and irradiation. Supervisor intervention is limited to safety concerns for students or equipment. Measurements of output factors using two chambers (regular and small field) for various field sizes (1×1 to 30×30 cm{sup 2}) and of detailed depth-dose curves for 6 MV, 6 and 12 MeV beams are to be performed and discussed in a formal report. Conclusion: Full access to, and control of, a Linac is the high point of this course. It provides a glimpse of medical physics and generates an experimental background for those continuing to CAMPEP programs. This dedicated, non-clinical facility further enable enhance CAMPEP graduate teaching and research activities not possible with a clinical device.

60 APPLIED LIFE SCIENCES↗

Electric response of a metal-molecule-metal junction to laser pulse by solving hierarchical equations of motion

We have combined the quantum dissipative theory and the time dependent density functional theory to perform the first principle calculation of laser induced quantum dynamical electron transport through a molecule weak bridged to two electrodes. The formalism of hierarchical equations of motion based on non-equilibrium Green’s function theory has been taken in this work. Numerical simulations of optical absorption spectra of benzene, laser induced transient current without and with bias, charge pumping effect, as well as the spectrum analysis from the current in Au-benzene-Au molecular junction are presented and discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomistic modeling of the metallic-to-semiconducting phase boundaries in monolayer MoS{sub 2}

Recent experimental demonstration on the coexistence of metallic and semiconducting phases in the same monolayer MoS{sub 2} crystal has attracted much attention for its use in ultra-low contact resistance-MoS{sub 2} transistors. However, the electronic structures of the metallic-to-semiconducting phase boundaries, which appear to dictate the carrier injection in such transistors, are not yet well understood. In this letter, interfacing the 2H and 1T′ polytypes appropriately, we first model the “beta” (β) and the “gamma” (γ) phase boundaries, and demonstrate good agreement with experiential results. We then apply first-principles based density functional theory to calculate the electronic structures for those optimized geometries. We further employ non equilibrium Green's function formalism to evaluate the transmission spectra and the local density of states (LDOS) in order to assess the Schottky barrier nature of the phase boundaries. Our study reveals that while the γ boundary yields p-type Schottky barrier, the β boundary leads to the distinct symmetric Schottky barrier with an atomically sharp transition region. This understanding could be useful for designing high performance transistors using phase-engineered MoS{sub 2} crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-E-P-22: AAPM Task Group 263 Tackling Standardization of Nomenclature for Radiation Therapy

Purpose: There is growing recognition of need for increased clarity and consistency in the nomenclatures used for body and organ structures, DVH metrics, toxicity, dose and volume units, etc. Standardization has multiple benefits; e.g. facilitating data collection for clinical trials, enabling the pooling of data between institutions, making transfers (i.e. hand-offs) between centers safer, and enabling vendors to define “default” settings. Towards this goal, the American Association of Physicists in Medicine (AAPM) formed a task group (TG263) in July of 2014, operating under the Work Group on Clinical Trials to develop consensus statements. Guiding principles derived from the investigation and example nomenclatures will be presented for public feedback. Methods: We formed a multi-institutional and multi-vendor collaborative group of 39 physicists, physicians and others involved in clinical use and electronic transfer of information. Members include individuals from IROC, NRG, IHE-RO, DICOM WG-7, ASTRO and EORTC groups with overlapping interests to maximize the quality of the consensus and increase the likelihood of adoption. Surveys of group and NRG members were used to define current nomenclatures and requirements. Technical requirements of vendor systems and the proposed DICOM standards were examined. Results: There is a marked degree of inter and intra institutional variation in current approaches, resulting from inter-vendor differences in capabilities, clinic specific conceptualizations and inconsistencies. Using a consensus approach, the group defined optimal formats for the naming of targets and normal structures. A formal objective assessment of 13 existing clinically-used software packages show that all had capabilities to accommodate these recommended nomenclatures. Conclusions: A multi-stakeholder effort is making significant steps forward in developing a standard nomenclature that will work across platforms. Our current working list includes > 550 structures. Outreach efforts are ongoing to ensure broader participation in evaluating and testing the principles as they are developed by TG263.

60 APPLIED LIFE SCIENCES↗

Novel Carbon (C)-Boron (B)-Nitrogen (N)-Containing H2 Storage Materials

The following summarizes the research conducted for DOE project DE-EE0005658 “Novel Carbon(C)-Boron(B)-Nitrogen(N)-Containing H2 Storage Materials”. This work focused in part on the continued study of two materials identified from the preceding project DE-FG360GO18143 (“Hydrogen Storage by Novel CBN Heterocycle Materials”) as lead candidates to meet the DOE technical targets for either vehicular or non-automotive hydrogen storage applications. Specifically, a room-temperature liquid, 3-methyl-1,2-cyclopentane (B), and a high H2 capacity solid, 1,2-BN-cyclohexane (J), were selected for further characterization and performance optimization. In addition to these compounds, the current project also aimed to prepare several new materials predicted to be disposed towards direct reversibility of H2 release and uptake, a feature deemed critical to achieving efficient recycling of spent fuel end products. To assist in the rational design of these and other next-generation materials, this project undertook to investigate the mechanism of hydrogen release from established compounds (mainly B and J) using a combined experimental/computational approach. Among this project’s signature accomplishments, the preliminary synthetic route to B was optimized for production on decagram scale. With such quantities of material available, its performance in powering an actual 30 W proton exchange membrane (PEM) fuel cell stack was tested and found to be identical to that of facility H2. Despite this positive proof-of-concept achievement, however, further consideration of neat B as a potential hydrogen storage material was abandoned due to evidence of thermal instability. Specifically, mass spectrometry-coupled thermogravimetric analysis (TGA-MS) revealed significant H2 release from B to initiate at 50 °C, well below the 60 °C minimum threshold set by the DOE. This result prompted a more extensive investigation in the decomposition mechanism of B vis-à-vis that of J, which exhibited in neat form a substantially higher onset temperature for spontaneous H2 release (70 °C). Solution-phase kinetic experiments using ReactIR established a second-order dependence for the initial loss of H2 from both B and J; Arrhenius analysis, however, revealed the activation barrier for this reaction was lower for B than for J, which presumably contributes to the diminished thermal stability of the former. On the basis of these and other experimental results, extensive computational efforts yielded a reasonable mechanistic model for the dehydrogenation of 1,2-BN-cycloalkane materials. While the prospect of neat B as a suitable hydrogen storage material was discarded, it was proposed that the combination of B with more thermally stable amine-borane-based materials might afford mixtures with improved properties. Indeed, when B was combined with ammonia borane (AB) in a 2:1 molar ratio, the two materials formed a liquid. More significantly, this mixture remained liquid even after complete dehydrogenation, thus establishing the potential for a single-phase fuel cycle. (In contrast, the dehydrogenation product of neat B is a low melting solid (mp = 28-30 °C).) Another advantage conferred by the blend formulation was a dramatic reduction in the amount of borazine produced by AB. Borazine is a well-known contaminant of H2 produced by the thermal decomposition of neat AB, and exerts deleterious effects on fuel cell performance. Residual gas analysis (RGA) of the gas stream generated from the B-AB blend, however, detected just 0.01% borazine content when a Pt-Ni nanoparticle dehydrogenation catalyst was used. In all ii then, the 2:1 B-AB blend marks a major achievement in the effort to develop a suitable liquid amine-borane hydrogen storage material, and merits further investigation into the optimization for practical adoption. Similar realization of the potential of J as a high % wt. H2 material required a method to dehydrogenate the carbonaceous components of the molecule without the use of a sacrificial hydrogen acceptor, as had been reported in the previous project. Ultimately, this reaction was achieved for a B,N-disubstituted BN-cyclohexene model substrate using a gas flow system with a fixed Pd/C catalyst bed. Considerable work remains, however, to translate these initial results into a general protocol for complete dehydrogenation of fully saturated BN-cycloalkane materials such as J. With concrete confirmation of the possibility to perform both BN and CC dehydrogenation on a single theoretical substrate, COMSOL modeling was used to evaluate the effects of thermodynamically coupling the two reactions. It was hypothesized that the heat generated from exothermic BN dehydrogenation would partially drive the endothermic CC dehydrogenation reaction; this additional heat consumption was expected to in turn confer the benefit of lowering the maximum reactor temperature. A two-dimensional model of an axisymmetric reactor including experimental kinetic and calculated thermodynamic parameters for both reactions did indeed predict these outcomes. The extent to which the effects of thermodynamic coupling actually manifested, however, were also revealed to depend strongly on the relative rates of the two reactions, as well as the magnitude of the equilibrium constant governing the progress of the endothermic process. Given the evident complexity of attaining high effective % wt. H2 capacity with J, alternative systems were investigated for greater facility of extensive H2 release. Among those studied, 1,2,4,5-bis-BN-cyclohexane (H) demonstrated the most favorable properties, particularly with respect to thermal stability: rather than decompose, a neat sample instead sublimed when heated above 150 °C. Nevertheless, two commercially available catalytic systems were identified to effect release of two H2 equivalents from H. Release of further equivalents were apparently impeded by the formation of either polymeric material or one of two dimeric cage compounds depending on the catalyst used. Notably, a method to regenerate H from these product mixtures remains to be developed. Thus, while H may prove useful for certain long-term energy storage needs, it is currently less suited applications involving frequent fuel consumption. Similar difficulties were also encountered in attempts to realize the complete fuel cycle of 1,3-BN-cyclohexane (E) and B,N-substituted derivatives thereof. It had been initially proposed that E would provide for readily reversible BN dehydrogenation through a measure of frustrated Lewis pair-type character. Indeed, computations predicted this reaction would be essentially thermoneutral in solution. In the course of attempts to fully hydrogenate the spent fuel, however, dimeric species formed and proved resistant to further BN reduction. While a number of monomeric cyclic compounds were also successfully synthesized as formal boron-nitrogen frustrated Lewis pairs, none demonstrated any capacity to split H2 across the BN unit. The challenge of developing a practical amine-borane-based material for readily reversible hydrogen storage thus remains unresolved at this time. As such, it deserves consideration as a major objective of any future work.

08 HYDROGEN↗