SEARCH · Search NASA

Results for “deformation function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2017

The analytic structure of conformal blocks and the generalized Wilson-Fisher fixed points

Here, we describe in detail the method used in our previous work arXiv:1611.10344 to study the Wilson-Fisher critical points nearby generalized free CFTs, exploiting the analytic structure of conformal blocks as functions of the conformal dimension of the exchanged operator. Our method is equivalent to the mechanism of conformal multiplet recombination set up by null states. We also compute, to the first non-trivial order in the ε-expansion, the anomalous dimensions and the OPE coefficients of infinite classes of scalar local operators using just CFT data. We study single-scalar and O(N)-invariant theories, as well as theories with multiple deformations. When available we agree with older results, but we also produce a wealth of new ones. Furthermore, unitarity and crossing symmetry are not used in our approach and we are able to apply our method to non-unitary theories as well. Some implications of our results for the study of the non-unitary theories containing partially conserved higher-spin currents are briefly mentioned.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Structural study of the effects of mutations in proteins to identify the molecular basis of the loss of local structural fluidity leading to the onset of autoimmune diseases

Protein-Protein Interactions (PPIs) are crucial in most of the biological processes and PPI dysfunctions are known to be associated with the onsets of various diseases. One of such diseases is the auto-immune disease. Auto-immune diseases are one among the less studied group of diseases with very high mortality rates. Thus, we tried to correlate the appearances of mutations with their probable biochemical basis of the molecular mechanisms leading to the onset of the disease phenotypes. We compared the effects of the Single Amino Acid Variants (SAVs) in the wild type and mutated proteins to identify any structural deformities that might lead to altered PPIs leading ultimately to disease onset. For this we used Relative Solvent Accessibility (RSA) as a spatial parameter to compare the structural perturbation in mutated and wild type proteins. We observed that the mutations were capable to increase intra-chain PPIs whereas inter-chain PPIs would remain mostly unaltered. This might lead to more intra-molecular friction causing a deleterious alteration of protein's normal function. A Lyapunov exponent analysis, using the altered RSA values due to polymorphic and disease causing mutations, revealed polymorphic mutations have a positive mean value for the Lyapunov exponent while disease causing mutations have a negative mean value. Thus, local spatial stochasticity has been lost due to disease causing mutations, indicating a loss of structural fluidity. The amino acid conversion plot also showed a clear tendency of altered surface patch residue conversion propensity than polymorphic conversions. So far, this is the first report that compares the effects of different kinds of mutations (disease and non-disease causing polymorphic mutations) in the onset of autoimmune diseases. - Highlights: • Protein-Protein Interaction. • Changes in Relative Solvent Accessibility (RSA). • Amino acid conversion matrix. • Polymorphic mutations. • Disease causing mutations.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Fracture and fatigue of thin crystalline SrTiO 3 membranes

The increasing availability of a variety of two-dimensional materials has generated enormous growth in the field of nanoengineering and nanomechanics. Recent developments in thin film synthesis have enabled the fabrication of freestanding functional oxide membranes that can be readily incorporated in nanomechanical devices. While many oxides are extremely brittle in bulk, recent studies have shown that, in thin membrane form, they can be much more robust to fracture as compared to their bulk counterparts. Here, we investigate the ultimate tensile strength of SrTiO 3 membranes by probing freestanding SrTiO 3 drumheads using an atomic force microscope. We demonstrate that SrTiO 3 membranes can withstand an elastic deformation with an average strain of ~6% in the sub-20 nm thickness regime, which is more than an order of magnitude beyond the bulk limit. Further, we also show that these membranes are highly resilient upon a high cycle fatigue test, surviving up to a billion cycles of force modulation at 85% of their fracture strain, demonstrating their high potential for use in nanomechanical applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Machine learning approach to pattern recognition in nuclear dynamics from the ab initio symmetry-adapted no-core shell model

A novel machine learning approach is used to provide further insight into atomic nuclei and to detect orderly patterns amid a vast data of large-scale calculations. The method utilizes a neural network that is trained on ab initio results from the symmetry-adapted no-core shell model (SA-NCSM) for light nuclei. We show that the SA-NCSM, which expands ab initio applications up to medium-mass nuclei by using dominant symmetries of nuclear dynamics, can reach heavier nuclei when coupled with the machine learning approach. In particular, we find that a neural network trained on probability amplitudes for s- and p-shell nuclear wave functions not only predicts dominant configurations for heavier nuclei but in addition, when tested for the 20 Ne ground state, accurately reproduces the probability distribution. The non-negligible configurations predicted by the network provide an important input to the SA-NCSM for reducing ultralarge model spaces to manageable sizes that can be, in turn, utilized in SA-NCSM calculations to obtain accurate observables. The neural network is capable of describing nuclear deformation and is used to track the shape evolution along the 20-42 Mg isotopic chain, suggesting a shape coexistence that is more pronounced toward the very neutron-rich isotopes. We provide first descriptions of the structure and deformation of 24 Si and 40 Mg of interest to x-ray burst nucleosynthesis, and even of the extremely heavy nuclei such as 166,168 Er and 236 U, that build on first-principles considerations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Establishing the Maximum Collectivity in Highly Deformed N = Z Nuclei

The lifetimes of the first excited 2 + states in the N = Z nuclei 80 Zr, 78 Y, and 76 Sr have been measured using the γ-ray line shape method following population via nucleon-knockout reactions from intermediate-energy rare-isotope beams. The extracted reduced electromagnetic transition strengths yield new information on where the collectivity is maximized and provide evidence for a significant, and as yet unexplained, odd-odd vs even-even staggering in the observed values. As a result, the experimental results are analyzed in the context of state-of-the-art nuclear density-functional model calculations.

59 ≤ A ≤ 89↗

DOE OSTI2018

Functional Performances of CuZnAl Shape Memory Alloy Open-Cell Foams

Shape memory alloys (SMAs) with cellular structure offer a unique mixture of thermo-physical–mechanical properties. These characteristics can be tuned by changing the pore size and make the shape memory metallic foams very attractive for developing new devices for structural and functional applications. In this work, CuZnAl SMA foams were produced through the liquid infiltration of space holder method. In comparison, a conventional CuZn brass alloy was foamed trough the same method. Functional performances were studied on both bulk and foamed SMA specimens. Calorimetric response shows similar martensitic transformation (MT) below 0 °C. Compressive response of CuZnAl revealed that mechanical behavior is strongly affected by sample morphology and that damping capacity of metallic foam is increased above the MT temperatures. The shape memory effect was detected in the CuZnAl foams. The conventional brass shows a compressive response similar to that of the martensitic CuZnAl, in which plastic deformation accumulation occurs up to the cellular structure densification after few thermal cycles.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Energy based fracture initiation criterion for strain-crystallizing rubber-like materials with pre-existing cracks

Fracture prediction is indispensable for polymers, like rubbers, which have a broad range of applications mainly due to their high extensibility. The phenomenon known as strain-induced crystallization further contributes to the fracture toughness of certain rubbers. Here in this study, a criterion based on internal bond energy, incorporating the effects of crystallization, is proposed to predict fracture initiation in rubber-like materials with pre-existing cracks. First, a multi-scale mechanical model is developed for characterizing the behavior of rubber when subjected to both uniaxial and biaxial deformation states. At the microscale, both the amorphous and crystalline chain segments are modeled as elastic in order to consider the energy contribution by the molecular bond distortions. This internal energy is considered along with the entropic and crystalline free energy for each chain. In the chain model, the effects of loading condition and the relative orientation of a chain on its crystallinity are taken into account. At the macroscopic scale, an existing crystallinity distribution function is adapted and a mixed finite element formulation with an augmented Lagrangian multiplier is utilized to impose the incompressibility constraint. A non-affine maximal advance path constraint based homogenization model is utilized for bridging the two scales. Its potential to account for anisotropy in the stretched network compels the model to be preferable due to its physical significance, for the purpose of fracture modeling. The rigidity of the crystallites is accounted for by proposing a crystallite distortion energy in addition to the critical bond dissociation energy, for fracture initiation to occur. The model is validated by comparison with existing experimental results for both crystallizing and non-crystallizing rubbers. In addition to its potential to predict the material behavior when subjected to uniaxial and biaxial loading, the capability of the model to quantitatively estimate the effect of crystallization on fracture initiation is also verified.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Local structural ordering determines the mechanical damage tolerance of amorphous grain boundary complexions

Amorphous grain boundary complexions act as toughening features within a microstructure because they can absorb dislocations more efficiently than traditional grain boundaries. This toughening effect should be a strong function of the local internal structure of the complexion, which has recently been shown to be determined by grain boundary crystallography. To test this hypothesis, molecular dynamics are used here to simulate dislocation absorption and damage nucleation for complexions with different distributions of structural short-range order. Here, the complexion with a more disordered structure away from the dislocation absorption site is actually found to better resist crack nucleation, as damage tolerance requires delocalized deformation and the operation of shear-transformation zones through the complexion thickness. The more damage tolerant complexion accommodates plastic strain efficiently within the entire complexion, providing the key mechanistic insight that local patterning and asymmetry of structural short-range order controls the toughening effect of amorphous complexions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Chiral-like doublet band structure and octupole correlations in 104 Ag

The nature of the yrast negative-parity band and its chiral-like partner band in 104 Ag is investigated experimentally and theoretically. Lifetimes of states in the negative-parity yrast band and positive-parity band based on the 4424-keV level are measured using Doppler shift attenuation technique. Lifetimes of three more states have been determined along with the upper limit for the lifetime of the highest observed yrast states. Further, lifetimes known from earlier studies are determined with better precision. The level scheme of 104 Ag has also been extended with the addition of new enhanced E1 transitions linking the positive-parity band based on the 4424-keV levels and the yrast negative-parity and its partner band. B(E1) and/or B(E1)/B(M1) values for the transitions from the positive-parity band to the yrast and its partner band have been determined for the first time; these suggest strong octupole correlation between the positive-parity and the negative-parity bands. Here, calculations based on the triaxial projected shell model (TPSM) and covariant density functional theory (CDFT) have been performed to unravel the intrinsic structures of the partner band and the excited positive-parity band. TPSM calculations predict that doublet bands have significant angular momentum contributions along the three principle axes, suggesting that bands could have chiral symmetry breaking origin. The CDFT calculations predict a π(g 9/2 ) -1 Ⓧν(h 11/2 )(g 7/2 ,d 5/2 ) 2 aligned quasiparticle configuration for the negative-parity doublet bands with deformation parameters β≈0.20 and γ≈5°. The partner band could be interpreted as a chiral vibration mode built on top of the yrast band. The excited positive-parity band is predicted to have aligned four quasiparticle configurations, namely, π(g 9/2 ) -1 Ⓧν(h 11/2 ) 2 (g 7/2 ,d 5/2 ) 1 . Further, these calculations predict significant octupole softness in 104 Ag which could be the reason for enhanced E1 transitions between the four quasiparticle positive-parity bands and the doublet negative-parity bands.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2016

Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation

Our work in this project is aimed at making fundamental advances in multiscale methods for flow and transport in highly heterogeneous porous media. The main thrust of this research is to develop a systematic multiscale analysis and efficient coarse-scale models that can capture global effects and extend existing multiscale approaches to problems with additional physics and uncertainties. A key emphasis is on problems without an apparent scale separation. Multiscale solution methods are currently under active investigation for the simulation of subsurface flow in heterogeneous formations. These procedures capture the effects of fine-scale permeability variations through the calculation of specialized coarse-scale basis functions. Most of the multiscale techniques presented to date employ localization approximations in the calculation of these basis functions. For some highly correlated (e.g., channelized) formations, however, global effects are important and these may need to be incorporated into the multiscale basis functions. Other challenging issues facing multiscale simulations are the extension of existing multiscale techniques to problems with additional physics, such as compressibility, capillary effects, etc. In our project, we explore the improvement of multiscale methods through the incorporation of additional (single-phase flow) information and the development of a general multiscale framework for flows in the presence of uncertainties, compressible flow and heterogeneous transport, and geomechanics. We have considered (1) adaptive local-global multiscale methods, (2) multiscale methods for the transport equation, (3) operator-based multiscale methods and solvers, (4) multiscale methods in the presence of uncertainties and applications, (5) multiscale finite element methods for high contrast porous media and their generalizations, and (6) multiscale methods for geomechanics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

Instability of the Human Cytochrome P450 Reductase A287P Variant Is the Major Contributor to Its Antley-Bixler Syndrome-like Phenotype

Human NADPH-cytochrome P450 oxidoreductase (POR) gene mutations are associated with severe skeletal deformities and disordered steroidogenesis. The human POR mutation A287P presents with disordered sexual development and skeletal malformations. Difficult recombinant expression and purification of this POR mutant suggested that the protein was less stable than WT. The activities of cytochrome P450 17A1, 19A1, and 21A2, critical in steroidogenesis, were similar using our purified, full-length, unmodified A287P or WT POR, as were those of several xenobiotic-metabolizing cytochromes P450, indicating that the A287P protein is functionally competent in vitro, despite its functionally deficient phenotypic behavior in vivo. Differential scanning calorimetry and limited trypsinolysis studies revealed a relatively unstable A287P compared with WT protein, leading to the hypothesis that the syndrome observed in vivo results from altered POR protein stability. The crystal structures of the soluble domains of WT and A287P reveal only subtle differences between them, but these differences are consistent with the differential scanning calorimetry results as well as the differential susceptibility of A287P and WT observed with trypsinolysis. The relative in vivo stabilities of WT and A287P proteins were also examined in an osteoblast cell line by treatment with cycloheximide, a protein synthesis inhibitor, showing that the level of A287P protein post-inhibition is lower than WT and suggesting that A287P may be degraded at a higher rate. Overall, current studies demonstrate that, unlike previously described mutations, A287P causes POR deficiency disorder due to conformational instability leading to proteolytic susceptibility in vivo, rather than through an inherent flavin-binding defect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Integrability technique for fluid flow induced deformation of a boundary hair

The deformation of a dense carpet of hair due to Stokes flow in a channel can be described by a nonlinear integrodifferential equation for the shape of a single hair, which possesses several solutions for a given choice of parameters. Although it was posed in a previous study and it bears a resemblance to the pendulum problem from mechanics, this equation has not been analytically solved until now. Despite the presence on an integral with a nonlinear functional dependence on the dependent variable, the system is integrable. We compare the analytically obtained solution to a finite-difference numerical approach, identify the physically realizable solution branch, and briefly study the solution structure through a conserved energylike quantity. Time-dependent fluid-structure interactions are a rich and complex subject to investigate, and we argue that the solution discussed herein can be used as a basis for understanding these systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Observation of the distribution of nuclear magnetization in a molecule

Precise experimental control and interrogation of molecules and calculations of their structure are enriching the investigation of nuclear and particle physics phenomena. Molecules containing heavy, octupole-deformed nuclei, such as radium, are of particular interest. Here, we report precision laser spectroscopy measurements and theoretical calculations of the structure of the radioactive radium monofluoride molecule 225 Ra 19 F. Our results reveal fine details of the short-range electron-nucleus interaction, indicating the high sensitivity of this molecule to the distribution of magnetization, within the radium nucleus. Here, these results provide a stringent test of the description of the electronic wave function inside the nuclear volume, highlighting the suitability of these molecules for investigating subatomic phenomena.

Nuclear structure↗

DOE OSTI2017

Scaling invariance of spherical projectile fragmentation upon high-velocity impact on a thin continuous shield

The problem of aluminum projectile fragmentation upon high-velocity impact on a thin aluminum shield is considered. A distinctive feature of this description is that the fragmentation has been numerically simulated using the complete system of equations of deformed solid mechanics by a method of smoothed particle hydrodynamics in three-dimensional setting. The transition from damage to fragmentation is analyzed and scaling relations are derived in terms of the impact velocity (V), ratio of shield thickness to projectile diameter (h/D), and ultimate strength (σ{sub p}) in the criterion of projectile and shield fracture. Analysis shows that the critical impact velocity V{sub c} (separating the damage and fragmentation regions) is a power function of σ{sub p} and h/D. In the supercritical region (V > V{sub c}), the weight-average fragment mass asymptotically tends to a power function of the impact velocity with exponent independent of h/D and σ{sub p}. Mean cumulative fragment mass distributions at the critical point are scale-invariant with respect to parameters h/D and σ{sub p}. Average masses of the largest fragments are also scale-invariant at V > V{sub c}, but only with respect to variable parameter σ{sub p}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗