Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “developable planes”

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 613 records · Page 34

Perfect single-sided radiation and absorption without mirrors

Highly directional radiation from photonic structures is important for many applications, including high-power photonic crystal surface-emitting lasers, grating couplers, and light detection and ranging devices. However, previous dielectric, few-layer designs only achieved moderate asymmetry ratios, and a fundamental understanding of bounds on asymmetric radiation from arbitrary structures is still lacking. Here in this paper, we show that breaking the 180° rotational symmetry of the structure is crucial for achieving highly asymmetric radiation. We develop a general temporal coupled-mode theory formalism to derive bounds on the asymmetric decay rates to the top and bottom of a photonic crystal slab for a resonance with arbitrary in-plane wavevector. Guided by this formalism, we show that infinite asymmetry is still achievable even without the need for back-reflection mirrors, and we provide numerical examples of designs that achieve asymmetry ratios exceeding 10 4 . The emission direction can also be rapidly switched from top to bottom by tuning the wavevector or frequency. Furthermore, we show that with the addition of weak material absorption loss, such structures can be used to achieve perfect absorption with single-sided illumination, even for single-pass material absorption rates less than 0.5% and without back-reflection mirrors. Our work provides new design principles for achieving highly directional radiation and perfect absorption in photonics.

36 MATERIALS SCIENCE↗

Mineral Deposition on the Rough Walls of a Fracture

Modeling carbonate growth in fractures and pores is important for understanding carbon sequestration in the environment or when supersaturated solutions are injected into rocks. Here, we study the simple but nontrivial problem of calcite growth on fractures with rough walls of the same mineral using kinetic Monte Carlo simulations of attachment and detachment of molecules and scaling approaches. First, we consider wedge-shaped fracture walls whose upper terraces are in the same low-energy planes and show that the valleys are slowly filled by the propagation of parallel monolayer steps in the wedge sides. The growth ceases when the walls reach these low-energy configurations so that a gap between the walls may not be filled. Second, we consider fracture walls with equally separated monolayer steps (vicinal surfaces with roughness below 1 nm) and show that growth by step propagation will eventually clog the fracture gap. In both cases, scaling approaches predict the times to attain the final configurations as a function of the initial geometry and the step-propagation velocity, which is set by the saturation index. The same reasoning applied to a random wall geometry shows that step propagation leads to lateral filling of surface valleys until the wall reaches the low-energy crystalline plane that has the smallest initial density of molecules. Thus, the final configurations of the fracture walls are much more sensitive to the crystallography than to the roughness or the local curvature. The framework developed here may be used to determine those configurations, the times to reach them, and the mass of deposited mineral. Effects of transport limitations are discussed when the fracture gap is significantly narrowed.

calcite↗

BISON Robustness and Performance Improvements

BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SU-E-T-490: Independent Three-Dimensional (3D) Dose Verification of VMAT/SBRT Using EPID and Cloud Computing

Purpose: Dosimetric verification of VMAT/SBRT is currently performed on one or two planes in a phantom with either film or array detectors. A robust and easy-to-use 3D dosimetric tool has been sought since the advent of conformal radiation therapy. Here we present such a strategy for independent 3D VMAT/SBRT plan verification system by a combined use of EPID and cloud-based Monte Carlo (MC) dose calculation. Methods: The 3D dosimetric verification proceeds in two steps. First, the plan was delivered with a high resolution portable EPID mounted on the gantry, and the EPID-captured gantry-angle-resolved VMAT/SBRT field images were converted into fluence by using the EPID pixel response function derived from MC simulations. The fluence was resampled and used as the input for an in-house developed Amazon cloud-based MC software to reconstruct the 3D dose distribution. The accuracy of the developed 3D dosimetric tool was assessed using a Delta4 phantom with various field sizes (square, circular, rectangular, and irregular MLC fields) and different patient cases. The method was applied to validate VMAT/SBRT plans using WFF and FFF photon beams (Varian TrueBeam STX). Results: It was found that the proposed method yielded results consistent with the Delta4 measurements. For points on the two detector planes, a good agreement within 1.5% were found for all the testing fields. Patient VMAT/SBRT plan studies revealed similar level of accuracy: an average γ-index passing rate of 99.2± 0.6% (3mm/3%), 97.4± 2.4% (2mm/2%), and 72.6± 8.4 % ( 1mm/1%). Conclusion: A valuable 3D dosimetric verification strategy has been developed for VMAT/SBRT plan validation. The technique provides a viable solution for a number of intractable dosimetry problems, such as small fields and plans with high dose gradient.

60 APPLIED LIFE SCIENCES↗

Two-Dimensional Palladium Diselenide with Strong In-Plane Optical Anisotropy and High Mobility Grown by Chemical Vapor Deposition

Two-dimensional (2D) palladium diselenide (PdSe 2 ) has strong interlayer coupling and a puckered pentagonal structure, leading to remarkable layer-dependent electronic structures and highly anisotropic in-plane optical and electronic properties. However, the lack of high-quality, 2D PdSe 2 crystals grown by bottom-up approaches limits the study of their exotic properties and practical applications. In this work, chemical vapor deposition growth of highly crystalline few-layer (≥2 layers) PdSe 2 crystals on various substrates is reported. The high quality of the PdSe 2 crystals is confirmed by low-frequency Raman spectroscopy, scanning transmission electron microscopy, and electrical characterization. In addition, strong in-plane optical anisotropy is demonstrated via polarized Raman spectroscopy and second-harmonic generation maps of the PdSe 2 flakes. A theoretical model based on kinetic Wulff construction theory and density functional theory calculations is developed and described the observed evolution of “square=like” shaped PdSe 2 crystals into rhombus due to the higher nucleation barriers for stable attachment on the (1,1) and (1,-1) edges, which results in their slower growth rates. Finally, few-layer PdSe 2 field-effect transistors reveal tunable ambipolar charge carrier conduction with an electron mobility up to ≈294 cm 2 V -1 s -1 , which is comparable to that of exfoliated PdSe 2 , indicating the promise of this anisotropic 2D material for electronics.

36 MATERIALS SCIENCE↗

Effect of Methyl Groups on Formation of Ordered or Layered Graphitic Materials from Aromatic Molecules

Abstract Developing functionally complex carbon materials from small aromatic molecules requires an understanding of how the chemistry and structure of its constituent molecules evolve and crosslink, to achieve a tailorable set of functional properties. Here, molecular dynamics (MD) simulations are used to isolate the effect of methyl groups on condensation reactions during the oxidative process and evaluate the impact on elastic modulus by considering three monodisperse pyrene‐based systems with increasing methyl group fraction. A parameter to quantify the reaction progression is designed by computing the number of new covalent bonds formed. Utilizing the previously developed MD framework, it is found that increasing methylation leads to an almost doubling of bond formation, a larger fraction of the new bonds oriented in the direction of tensile stress, and a higher basal plane alignment of the precursor molecules along the direction of tensile stress, resulting in enhanced tensile modulus. Additionally, via experiments, it is demonstrated that precursors with a higher fraction of methyl groups result in a higher alignment of molecules. Moreover, increased methylation results in the lower spread of single molecule alignment which may lead to smaller variations in tensile modulus and more consistent properties in carbon materials derived from methyl‐rich precursors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observation of a three-dimensional quasi-long-range electronic supermodulation in YBa 2 Cu 3 O 7-x /La 0.7 Ca 0.3 MnO 3 heterostructures

Recent developments in high-temperature superconductivity highlight a generic tendency of the cuprates to develop competing electronic (charge) supermodulations. While coupled with the lattice and showing different characteristics in different materials, these supermodulations themselves are generally conceived to be quasi-two-dimensional, residing mainly in individual CuO 2 planes, and poorly correlated along the c axis. Here we observed with resonant elastic X-ray scattering a distinct type of electronic supermodulation in YBa 2 Cu 3 O 7–x (YBCO) thin films grown epitaxially on La 0.7 Ca 0.3 MnO 3 (LCMO). This supermodulation has a periodicity nearly commensurate with four lattice constants in-plane, eight out of plane, with long correlation lengths in three dimensions. It sets in far above the superconducting transition temperature and competes with superconductivity below this temperature for electronic states predominantly in the CuO 2 plane. Our finding sheds light on the nature of charge ordering in cuprates as well as a reported long-range proximity effect between superconductivity and ferromagnetism in YBCO/LCMO heterostructures.

36 MATERIALS SCIENCE↗

Influence of dislocations and twin walls in BaTiO 3 on the voltage-controlled switching of perpendicular magnetization

In this work, we investigate the influence of dislocations and twin walls in BaTiO 3 on its ferroelectric response and the resulting effect on the perpendicular magnetic anisotropy (PMA) of a strain-coupled [Co\Ni] n film. A dense twinned structure in conjunction with a high dislocation density significantly reduces the converse piezoelectric effect of BaTiO 3 by hindering the propagation of newly nucleated domains with an applied electric field. This, in turn, results in a modest reduction of the PMA of the ferromagnetic layer. On the other hand, the ferroelectric polarization reorients from [100] to [001] direction in a dislocation-free BaTiO 3 , inducing the maximum achievable in-plane compressive strain of 1.1%. A large fraction of this uniaxial strain is transferred to the magnetoelastically coupled ferromagnetic layers whose magnetization switches to in plane via the inverse magnetostriction effect. This work reveals the critical role of the interplay between twin walls and dislocations within a ferroelectric substrate in the performance of multiferroic heterostructures and provides insight into the development of highly energy-efficient magnetoelectric devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

YbV 3 ⁢Sb 4 and EuV 3 ⁢Sb 4 vanadium-based kagome metals with Yb 2+ and Eu 2+ zigzag chains

Here, in this work, we present Yb⁢V 3 ⁢Sb 4 and Eu⁢V 3 ⁢Sb 4 , two compounds exhibiting slightly distorted vanadium-based kagome nets interleaved with zigzag chains of divalent Yb 2+ and Eu 2+ ions. Single crystal growth methods are reported alongside magnetic, electronic, and heat capacity measurements. Yb⁢V 3⁢ Sb 4 is a nonmagnetic metal with no collective phase transitions observed between 60 mK and 300 K. Conversely, Eu⁢V 3 ⁢Sb 4 is a magnetic kagome metal exhibiting easy-plane ferromagneticlike order below T C =32 K with hints of modulated spin texture under low field. With Yb⁢V 3 ⁢Sb 4 and Eu⁢V 3 ⁢Sb 4 we demonstrate another direction for the discovery and development of vanadium-based kagome metals while incorporating the chemical and magnetic degrees of freedom offered by a rare-earth sublattice.

36 MATERIALS SCIENCE↗

Laboratory Experiments Examining the Effect of Thermal and Mechanical Processes on Hydraulic Transmissivity Evolution

Using laboratory slide-hold-slide experiments, at temperatures from 22 to 200 degrees C, to examine effects of fracture reactivation and quasi-static loading on the evolution of fluid transport properties of simulated fractures in Westerly granite. At all temperatures, the in-plane hydraulic transmissivity consistently decays during hold periods resulting in an overall reduction in transmissivity. During the first three to fifteen hours of an experiment, transmissivity decreases rapidly due to the generation of wear products, development of a sliding surface, and compaction of the resulting gouge. Once the sliding surface has developed, the long-term transmissivity decay rate at 22 and 100 degrees C is significantly lower than the transmissivity decay rate during the initial 3-15 hours of the experiment. However, at 200 degrees C, the decay of hydraulic transmissivity remains high throughout the experiment. The long-term decay of hydraulic transmissivity can be fitted with a power law model with more rapid reduction of hydraulic transmissivity at higher temperature. Periods of sliding on the fracture surface result in transient increases in the transmissivity, due to shear dilation, as is expected for Coulomb materials. These transients are superimposed on the long-term decay. When sliding ceases and a new hold period commences, there is a rapid reduction in transmissivity and return to the long-term rate of transmissivity decay. The rate of decay of the transmissivity transients is inversely correlated with temperature, in contrast to the long-term decay and the expected behavior for processes like subcritical crack growth and indentation creep. The higher decay rates that are observed during the initial 3-15 hours of the tests and following sliding, are associated with times that the porosity of the gouge is expected to be high. The difference in decay rates suggests that when the gouge is driven far from equilibrium by active shearing, densification may be dominated by a different mechanism from long-term compaction.

15 GEOTHERMAL ENERGY↗

Kink band and shear band localization in anisotropic perfectly plastic solids

Shear-driven strain localization has been observed in a wide variety of materials and may take the form of shear bands or kink bands. Based on observations of kink bands in plastically anisotropic metallic nanolaminates and single crystal metals, here we posit that, for the specific case of isochoric deformation, the kinematics of kink band formation are indistinguishable from those of plane strain shear band formation. The only distinction between shear bands and kink bands in these systems would then be that kink bands ‘lock up’ at a particular value of material rotation while shear bands may progress to arbitrarily high strains. In order to investigate whether strong material anisotropy is sufficient to arrest shear localization at a geometry that matches the classic kink band geometry, we model the development of a band of simple shear within an anisotropic perfectly plastic material. The resulting analytical model provides the stress state needed to maintain the kinematics of simple shear as a function of material anisotropy, deformation band orientation, and shear strain (or equivalently, material rotation). It is found that plastic anisotropy can promote either kink band or shear band formation depending on the loading orientation. When the deviatoric stress is positive parallel to the plane of anisotropy, shear localization may progress without bound and a shear band is produced. When the deviatoric stress is negative parallel to the plane of anisotropy, shear localization is arrested after a certain material rotation, resulting in a kink band. Examination of the requisite applied stress state during kink band formation provides an explanation for the experimentally-observed ‘lock up’ geometry. Solutions for the band boundary inclination angle are obtained and used to provide bounds on permissible band angles for both shear bands and kink bands.

36 MATERIALS SCIENCE↗

Report on Design and Failure Limits of SiC/SiC and FeCrAl ATF Cladding Concepts under RIA

Pellet-cladding mechanical interaction (PCMI) is a possible failure scenario during a postulated reactivity-initiated accident (RIA). This report summarizes fiscal year 2018 research activities that were undertaken to evaluate the PCMI-like hoop strain–driven mechanical response of accident tolerant fuel cladding candidates of silicon carbide fiber/silicon carbide matrix composite (SiC/SiC) and iron–chromium–aluminum (FeCrAl) tube specimens. To achieve various RIA-like conditions, the modified-burst test device developed during fiscal years 2016 and 2017 was used. A heating capability was also added to the system for high-temperature testing of FeCrAl samples.The first section of this report assesses the accuracy of the two-dimensional digital image correlation calculated strain values for tube samples. The speckle patterns are computer generated in a two-dimensional plane, and the generated images were wrapped into a cylindrical surface for the validation using digital image correlation technique on curved geometries. Subsequent sections address the failure behavior of SiC/SiC composites and long-term heat treated and as-drawn FeCrAl alloy tubes at various strain rates relevant to the postulated RIA events. In the last part, the effect of possible design parameters of failure strain and the strain rate are discussed for SiC/SiC composites and FeCrAl alloys during the PCMI phase of RIA.Future work is planned to develop and design safety tests for accident tolerant fuel cladding candidates and the advanced fuel forms with respect to their transient responses. A broader perspective on the mechanical behavior of the unirradiated and irradiated accident tolerant fuel cladding candidates will be developed using the modified-burst test and other mechanical tests.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SU-G-BRC-10: Feasibility of a Web-Based Monte Carlo Simulation Tool for Dynamic Electron Arc Radiotherapy (DEAR)

Purpose: DEAR is a radiation therapy technique utilizing synchronized motion of gantry and couch during delivery to optimize dose distribution homogeneity and penumbra for treatment of superficial disease. Dose calculation for DEAR is not yet supported by commercial TPSs. The purpose of this study is to demonstrate the feasibility of using a web-based Monte Carlo (MC) simulation tool (VirtuaLinac) to calculate dose distributions for a DEAR delivery. Methods: MC simulations were run through VirtuaLinac, which is based on the GEANT4 platform. VirtuaLinac utilizes detailed linac head geometry and material models, validated phase space files, and a voxelized phantom. The input was expanded to include an XML file for simulation of varying mechanical axes as a function of MU. A DEAR XML plan was generated and used in the MC simulation and delivered on a TrueBeam in Developer Mode. Radiographic film wrapped on a cylindrical phantom (12.5 cm radius) measured dose at a depth of 1.5 cm and compared to the simulation results. Results: A DEAR plan was simulated using an energy of 6 MeV and a 3×10 cm{sup 2} cut-out in a 15×15 cm{sup 2} applicator for a delivery of a 90° arc. The resulting data were found to provide qualitative and quantitative evidence that the simulation platform could be used as the basis for DEAR dose calculations. The resulting unwrapped 2D dose distributions agreed well in the cross-plane direction along the arc, with field sizes of 18.4 and 18.2 cm and penumbrae of 1.9 and 2.0 cm for measurements and simulations, respectively. Conclusion: Preliminary feasibility of a DEAR delivery using a web-based MC simulation platform has been demonstrated. This tool will benefit treatment planning for DEAR as a benchmark for developing other model based algorithms, allowing efficient optimization of trajectories, and quality assurance of plans without the need for extensive measurements.

60 APPLIED LIFE SCIENCES↗

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗

SU-F-T-230: A Simple Method to Assess Accuracy of Dynamic Wave Arc Irradiation Using An Electronic Portal Imaging Device and Log Files

Purpose: The Dynamic Wave Arc (DWA) technique, where the multi-leaf collimator (MLC) and gantry/ring move simultaneously in a predefined non-coplanar trajectory, has been developed on the Vero4DRT. The aim of this study is to develop a simple method for quality assurance of DWA delivery using an electronic portal imaging device (EPID) measurements and log files analysis. Methods: The Vero4DRT has an EPID on the beam axis, the resolution of which is 0.18 mm/pixel at the isocenter plane. EPID images were acquired automatically. To verify the detection accuracy of the MLC position by EPID images, the MLC position with intentional errors was assessed. Tests were designed considering three factors: (1) accuracy of the MLC position (2) dose output consistency with variable dose rate (160–400 MU/min), gantry speed (2.4–6°/s), ring speed (0.5–2.5°/s), and (3) MLC speed (1.6–4.2 cm/s). All the patterns were delivered to the EPID and compared with those obtained with a stationary radiation beam with a 0° gantry angle. The irradiation log, including the MLC position and gantry/ring angle, were recorded simultaneously. To perform independent checks of the machine accuracy, the MLC position and gantry/ring angle position were assessed using log files. Results: 0.1 mm intentional error can be detected by the EPID, which is smaller than the EPID pixel size. The dose outputs with different conditions of the dose rate and gantry/ring speed and MLC speed showed good agreement, with a root mean square (RMS) error of 0.76%. The RMS error between the detected and recorded data were 0.1 mm for the MLC position, 0.12° for the gantry angle, and 0.07° for the ring angle. Conclusion: The MLC position and dose outputs in variable conditions during DWA irradiation can be easily detected using EPID measurements and log file analysis. The proposed method is useful for routine verification. This research is (partially) supported by the Practical Research for Innovative Cancer Control (15Ack0106151h0001) from Japan Agency for Medical Research and development, AMED. Authors Takashi Mizowaki and Masahiro Hiraoka have consultancy agreement with Mitsubishi Heavy Industries, Ltd., Japan.

60 APPLIED LIFE SCIENCES↗

Highly Anisotropic Drude-Weight-Reduction and Enhanced Linear-Dichroism in van der Waals Weyl Semimetal T d -MoTe 2 with Coherent Interlayer Electronic Transport

Weyl semimetal (WSM) states can be achieved by breaking spatial-inversion symmetry or time reversal symmetry. However, the anisotropy of the energy reduction contributing to the emergence of WSM states has seldom been investigated by experiments. A van der Waals metal MoTe 2 exhibits a type-II WSM phase below the monoclinic-to-orthorhombic-phase-transition temperature T c ~ 250 K. Here, in this work, a combined linearly polarized optical-spectroscopy and electrical-transport study of MoTe 2 is reported at different temperatures. The Drude components in the a-axis, b-axis, and c-axis optical conductivity spectra, together with the metallic out-of-plane and in-plane electrical resistivities, indicate the coherent inter-layer and in-plane charge transports. Moreover, the Drude weight in σ 1a (ω), rather than the Drude weights in σ 1b (ω) and σ 1c (ω), decreases dramatically below T c , which exhibits a highly anisotropic decrease in its Drude weight and thus suggests a strongly anisotropic reduction of the electronic kinetic energy in the WSM phase. Furthermore, below T c , due to the in-plane anisotropic spectral-weight transfer from Drude component to high-energy region, the in-plane inter-band-absorption anisotropy increases remarkably ≈770 meV, and has the largest value (≈0.68) of normalized linear dichroism among the reported type-II WSMs. The work sheds light on seeking new WSMs and developing novel photonic devices based on WSMs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Allowing for crystalline structure effects in Geant4

In recent years, the Geant4 toolkit for the Monte Carlo simulation of radiation with matter has seen large growth in its divers user community. A fundamental aspect of a successful physics experiment is the availability of a reliable and precise simulation code. Geant4 currently does not allow for the simulation of particle interactions with anything other than amorphous matter. To overcome this limitation, the GECO (GEant4 Crystal Objects) project developed a general framework for managing solid-state structures in the Geant4 kernel and validate it against experimental data. As a result, accounting for detailed geometrical structures allows, for example, simulation of diffraction from crystal planes or the channeling of charged particle.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface

The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗