Methane activation on Fe4 cluster : a density functional theory study

We report a comprehensive theoretical study on reaction of methane by Fe4 cluster. This Letter gains insight into the mechanism of the reaction and indicate the Fe4 cluster has strong catalytic effect on the activation reaction of methane. In detail, the results show the cleavage of the first C–H bond is both an energetically and kinetically favourable process and the breaking of the second C–H is the rate-determining step. Moreover, our Letter demonstrates that the different cluster size of iron can not only determine the catalytic activity of methane but also control the product selectivity.

A comparison of finite difference and finite volume methods for solving the space fractional advection-dispersion equation with variable coefficients

Transport processes within heterogeneous media may exhibit non-classical diffusion or dispersion; that is, not adequately described by the classical theory of Brownian motion and Fick's law. We consider a space fractional advection-dispersion equation based on a fractional Fick's law. The equation involves the Riemann-Liouville fractional derivative which arises from assuming that particles may make large jumps. Finite difference methods for solving this equation have been proposed by Meerschaert and Tadjeran. In the variable coefficient case, the product rule is first applied, and then the Riemann-Liouville fractional derivatives are discretised using standard and shifted Grunwald formulas, depending on the fractional order. In this work, we consider a finite volume method that deals directly with the equation in conservative form. Fractionally-shifted Grunwald formulas are used to discretise the fractional derivatives at control volume faces. We compare the two methods for several case studies from the literature, highlighting the convenience of the finite volume approach.

A finite element model of seat cushion indentation with a soft tissue human occupant model

Effective digital human model (DHM) simulation of automotive driver packaging ergonomics, safety and comfort depends on accurate modelling of occupant posture, which is strongly related to the mechanical interaction between human body soft tissue and flexible seat components. This paper presents a finite-element study simulating the deflection of seat cushion foam and supportive seat structures, as well as human buttock and thigh soft tissue when seated. The three-dimensional data used for modelling thigh and buttock geometry were taken on one 95th percentile male subject, representing the bivariate percentiles of the combined hip breadth (seated) and buttock-to-knee length distributions of a selected Australian and US population. A thigh-buttock surface shell based on this data was generated for the analytic model. A 6mm neoprene layer was offset from the shell to account for the compression of body tissue expected through sitting in a seat. The thigh-buttock model is therefore made of two layers, covering thin to moderate thigh and buttock proportions, but not more fleshy sizes. To replicate the effects of skin and fat, the neoprene rubber layer was modelled as a hyperelastic material with viscoelastic behaviour in a Neo-Hookean material model. Finite element (FE) analysis was performed in ANSYS V13 WB (Canonsburg, USA). It is hypothesized that the presented FE simulation delivers a valid result, compared to a standard SAE physical test and the real phenomenon of human-seat indentation. The analytical model is based on the CAD assembly of a Ford Territory seat. The optimized seat frame, suspension and foam pad CAD data were transformed and meshed into FE models and indented by the two layer, soft surface human FE model. Converging results with the least computational effort were achieved for a bonded connection between cushion and seat base as well as cushion and suspension, no separation between neoprene and indenter shell and a frictional connection between cushion pad and neoprene. The result is compared to a previous simulation of an indentation with a hard shell human finite-element model of equal geometry, and to the physical indentation result, which is approached with very high fidelity. We conclude that (a) SAE composite buttock form indentation of a suspended seat cushion can be validly simulated in a FE model of merely similar geometry, but using a two-layer hard/soft structure. (b) Human-seat indentation of a suspended seat cushion can be validly simulated with a simplified human buttock-thigh model for a selected anthropomorphism.

Comparison between analytical and 3D finite element solutions for borehole stresses in anisotropic elastic rock

We present a rigorous validation of the analytical Amadei solution for the stress concentration around an arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients b11 and b55 are not equal. It is shown from theoretical considerations and published experimental data that the b11 and b55 are not equal for realistic rocks. Second, we develop a 3D finite element elastic model within a hybrid analytical–numerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic, transverse isotropic and orthorhombic symmetries. It is concluded that the analytical Amadei solution is valid with no restriction on the borehole orientation or the symmetry of the elastic anisotropy.

How to achieve maximum charge carrier loading on heteroatom-substituted graphene nanoribbon edges: density functional theory study

The practical number of charge carriers loaded is crucial to the evaluation of the capacity performance of carbon-based electrodes in service, and cannot be easily addressed experimentally. In this paper, we report a density functional theory study of charge carrier adsorption onto zigzag edge-shaped graphene nanoribbons (ZGNRs), both pristine and incorporating edge substitution with boron, nitrogen or oxygen atoms. All edge substitutions are found to be energetically favorable, especially in oxidized environments. The maximal loading of protons onto the substituted ZGNR edges obeys a rule of [8-n-1], where n is the number of valence electrons of the edge-site atom constituting the adsorption site. Hence, a maximum charge loading is achieved with boron substitution. This result correlates in a transparent manner with the electronic structure characteristics of the edge atom. The boron edge atom, characterized by the most empty p band, facilitates more than the other substitutional cases the accommodation of valence electrons transferred from the ribbon, induced by adsorption of protons. This result not only further confirms the possibility of enhancing charge storage performance of carbon-based electrochemical devices through chemical functionalization but also, more importantly, provides the physical rationale for further design strategies.

Charge carrier exchange at chemically modified graphene edges: A density functional theory study

Heteroatom doping on the edge of graphene may serve as an effective way to tune chemical activity of carbon-based electrodes with respect to charge carrier transfer in an aqueous environment. In a step towards developing mechanistic understanding of this phenomenon, we explore herein mechanisms of proton transfer from aqueous solution to pristine and doped graphene edges utilizing density functional theory. Atomic B-, N-, and O- doped edges as well as the native graphene are examined, displaying varying proton affinities and effective interaction ranges with the H3O+ charge carrier. Our study shows that the doped edges characterized by more dispersive orbitals, namely boron and nitrogen, demonstrate more energetically favourable charge carrier exchange compared with oxygen, which features more localized orbitals. Extended calculations are carried out to examine proton transfer from the hydronium ion in the presence of explicit water, with results indicating that the basic mechanistic features of the simpler model are unchanged.

Finite element modelling of load bearing cold-formed steel wall systems under fire conditions

Light Gauge Steel Framing (LSF) walls are made of cold-formed, thin-walled steel lipped channel studs with plasterboard linings on both sides. However, these thin-walled steel sections heat up quickly and lose their strength under fire conditions despite the protection provided by plasterboards. A new composite wall panel was recently proposed to improve the fire resistance rating of LSF walls, where an insulation layer was used externally between the plasterboards on both sides of the wall frame instead of using it in the cavity. A research study using both fire tests and numerical studies was undertaken to investigate the structural and thermal behaviour of load bearing LSF walls made of both conventional and the new composite panels under standard fire conditions and to determine their fire resistance rating. This paper presents the details of finite element models of LSF wall studs developed to simulate the structural performance of LSF wall panels under standard fire conditions. Finite element analyses were conducted under both steady and transient state conditions using the time-temperature profiles measured during the fire tests. The developed models were validated using the fire test results of 11 LSF wall panels with various plasterboard/insulation configurations and load ratios. They were able to predict the fire resistance rating within five minutes. The use of accurate numerical models allowed the inclusion of various complex structural and thermal effects such as local buckling, thermal bowing and neutral axis shift that occurred in thin-walled steel studs under non-uniform elevated temperature conditions. Finite element analyses also demonstrated the improvements offered by the new composite panel system over the conventional cavity insulated system.

Dosimetric effects of a high-density spinal implant

In this study, a treatment plan for a spinal lesion, with all beams transmitted though a titanium vertebral reconstruction implant, was used to investigate the potential effect of a high-density implant on a three-dimensional dose distribution for a radiotherapy treatment. The BEAMnrc/DOSXYZnrc and MCDTK Monte Carlo codes were used to simulate the treatment using both a simplified, recltilinear model and a detailed model incorporating the full complexity of the patient anatomy and treatment plan. The resulting Monte Carlo dose distributions showed that the commercial treatment planning system failed to accurately predict both the depletion of dose downstream of the implant and the increase in scattered dose adjacent to the implant. Overall, the dosimetric effect of the implant was underestimated by the commercial treatment planning system and overestimated by the simplified Monte Carlo model. The value of performing detailed Monte Carlo calculations, using the full patient and treatment geometry, was demonstrated.

Adsorption and dissociation of ammonia borane outside and inside single-walled carbon nanotubes : A density functional theory study

Amonia borane (AB) has been identified as a potential candidate highcapacity hydrogen storage material. This work probes the adsorption and dissociation of AB inside and outside single-walled carbon nanotubes (SWCNTs) within the framework of density functional theory. The dissociation barriers of AB have been calculated and compared with that of the isolated AB molecule. On the basis of the present calculations, no notable improvement results from SWCNT confinement; on the contrary, the dissociation barrier slightly increases with respect to isolated AB.

A density functional theory study on CO2 capture and activation by graphene-like boron nitride with boron vacancy

First principle calculations for a hexagonal (graphene-like) boron nitride (g-BN) monolayer sheet in the presence of a boron-atom vacancy show promising properties for capture and activation of carbon dioxide. CO2 is found to decompose to produce an oxygen molecule via an intermediate chemisorption state on the defect g-BN sheet. The three stationary states and two transition states in the reaction pathway are confirmed by minimum energy pathway search and frequency analysis. The values computed for the two energy barriers involved in this catalytic reaction after enthalpy correction indicate that the catalytic reaction should proceed readily at room temperature.

A density functional theory study of CO2 and N2 adsorption on aluminium nitride single walled nanotubes

Strong binding of isolated carbon dioxide (CO2) on aluminium nitride (AlN) single walled nanotubes is verified using two different functionals. Two optimized configurations corresponding to physisorption and chemisorption are linked by a low energy barrier, such that the chemisorbed state is accessible and thermodynamically favored at low temperatures. In contrast, N2 is found only to form a physisorbed complex with the AlN nanotube, suggesting the potential application of aluminium nitride based materials for CO2 fixation. The effect of nanotube diameter on gas adsorption properties is also discussed. The diameter is found to have an important effect on the chemisorption of CO2, but has little effect on the physisorption of either CO2 or N2.

Adsorption of carbon dioxide and nitrogen on single-layer aluminum nitride nanostructures studied by density functional theory

The adsorption of carbon dioxide and nitrogen molecules on aluminum nitride (AlN) nanostructures has been explored using first-principle computational methods. Optimized configurations corresponding to physisorption and, subsequentially, chemisorption of CO2 are identified, in contrast to N2, for which only a physisorption structure is found. Transition-state searches imply a low energy barrier between the physisorption and chemisorption states for CO2 such that the latter is accessible and thermodynamically favored at room temperature. The effective binding energy of the optimized chemisorption structure is apparently larger than those for other CO2 adsorptive materials, suggesting the potential for application of aluminum nitride nanostructures for carbon dioxide capture and storage.

The effect of Fe doping on adsorption of CO2/N2within carbon nanotubes : a density functional theory study with dispersion corrections

An ab initio density functional theory (DFT) study with correction for dispersive interactions was performed to study the adsorption of N2 and CO2 inside an (8, 8) single-walled carbon nanotube. We find that the approach of combining DFT and van der Waals correction is very effective for describing the long-range interaction between N2/CO2 and the carbon nanotube (CNT). Surprisingly, exohedral doping of an Fe atom onto the CNT surface will only affect the adsorption energy of the quadrupolar CO2 molecule inside the CNT (20–30%), and not that of molecular N2. Our results suggest the feasibility of enhancement of CO2/N2 separation in CNT-based membranes by using exohedral doping of metal atoms.

Half metallicity in a zigzag double-walled nanotube nanodot : an Ab initio prediction

Ab initio density functional calculations were performed to study finite-length zigzag (7, 0) @ (16, 0) double-walled carbon nanotubes (DWCNTs) with H-termination at the open ends. We find that such a DWCNT nanodot displays a very large magnetic moment at the zigzag edges and the ground state displays symmetric anti-ferromagnetic coupling. When an external electric field is applied along the direction of tube axis, a gap is opened for one spin channel, whereas another spin channel remains metallic, i.e. half metallicity occurs. Our results suggest an important new avenue for the development of CNT-based spintronic materials with enhanced properties.

The catalytic role of an isolated-Ti atom in the hydrogenation of Ti-doped Al(001) surface : An ab initio density functional theory calculation

Recent work [S. Chaudhuri, J.T. Muckerman, J. Phys. Chem. B 109 (2005) 6952] reported that two Ti-substituted atoms on an Al(0 0 1) surface can form a catalytically active site for the dissociation of H2, but the diffusion barrier of atomic H away from Ti site is as high as 1.57 eV. By using ab initio density functional calculations, we found that two hydrogen molecules can dissociate on isolated-Ti atom doped Al(0 0 1) surface with small activation barriers (0.21 and 0.235 eV for first and second H2, respectively). Additionally, the diffusion barrier of atomic H away from Ti site is also moderate (0.47 eV). These results contribute further towards understanding the improved kinetics observed in recycling of hydrogen with Ti-doped NaAlH4.