Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional performance testing”

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,531 records · Page 85

Quantifying irradiation defects in zirconium alloys: a comparison between transmission electron microscopy and whole-pattern diffraction line-profile analysis - 2016-0089

Irradiation-induced dislocations significantly affect the mechanical properties of zirconium alloys, altering slip and influencing creep and growth. Thus, the quantitative characterization of irradiation defects as a function of fluence, cold work, and/or thermal treatments is important for models that attempt to predict their impact on properties. Whole-pattern diffraction line-profile analysis (DLPA) is a well-established modern tool for microstructure characterization based on first-principle physical models for dislocation density measurements in plastically deformed materials. However, applying these DLPA methods directly to irradiated materials yields higher than expected dislocation density values compared with historical transmission electron microscopy (TEM) measurements and past line-broadening analysis studies calibrated to TEM observations. In an effort to understand these differences, a new microstructural model was developed for DLPA to specifically address dislocation structures consisting of elliptical - and -component loops. To compare the refined DLPA method with TEM measurements, high-resolution neutron diffraction patterns on nonirradiated and irradiated Zr-2.5Nb samples were collected with the Neutron Powder Diffractometer instrument at the Los Alamos Neutron Science Center and were evaluated. High-resolution TEM measurements were performed at the Reactor Materials Testing Laboratory, Queen's University, for comparison with the DLPA results. The capabilities and inherent uncertainties of both the refined DLPA and TEM methods are compared and discussed in detail. We show that the differences between the density values provided by DLPA and TEM are inherent to the methods and can be reconciled with the interpretation of the data. (authors)

36 MATERIALS SCIENCE↗

An innovative heterogeneous transfer learning framework to enhance the scalability of deep reinforcement learning controllers in buildings with integrated energy systems

Deep Reinforcement Learning (DRL)-based control shows enhanced performance in the management of integrated energy systems when compared with Rule-Based Controllers (RBCs), but it still lacks scalability and generalisation due to the necessity of using tailored models for the training process. Transfer Learning (TL) is a potential solution to address this limitation. However, existing TL applications in building control have been mostly tested among buildings with similar features, not addressing the need to scale up advanced control in real-world scenarios with diverse energy systems. This paper assesses the performance of an online heterogeneous TL strategy, comparing it with RBC and offline and online DRL controllers in a simulation setup using EnergyPlus and Python. The study tests the transfer in both transductive and inductive settings of a DRL policy designed to manage a chiller coupled with a Thermal Energy Storage (TES). The control policy is pre-trained on a source building and transferred to various target buildings characterised by an integrated energy system including photovoltaic and battery energy storage systems, different building envelope features, occupancy schedule and boundary conditions (e.g., weather and price signal). The TL approach incorporates model slicing, imitation learning and fine-tuning to handle diverse state spaces and reward functions between source and target buildings. Results show that the proposed methodology leads to a reduction of 10% in electricity cost and between 10% and 40% in the mean value of the daily average temperature violation rate compared to RBC and online DRL controllers. Moreover, online TL maximises self-sufficiency and self-consumption by 9% and 11% with respect to RBC. Conversely, online TL achieves worse performance compared to offline DRL in either transductive or inductive settings. However, offline Deep Reinforcement Learning (DRL) agents should be trained at least for 15 episodes to reach the same level of performance as the online TL. Therefore, the proposed online TL methodology is effective, completely model-free and it can be directly implemented in real buildings with satisfying performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dynamical model of J/ψ photoproduction on the nucleon

Here, a dynamical model based on a phenomenological charm quark-nucleon (c – N) potential v cN and the Pomeronexchange mechanism is constructed to investigate the J/ψ photoproduction on the nucleon from threshold to invariant mass W = 300 GeV. The J/ψ – N potential, V J/ψN (r), is constructed by folding v cN into the wave function φ J/ψ ($c\bar{c}$) of J/ψ within a constituent quark model (CQM) of Segovia et al. [Int. J. Mod. Phys. E 22, 1330026 (2013)]. A photoproduction amplitude is also generated by v cN by a $c\bar{c}$–loop integration over the γ → $c\bar{c}$ vertex function and φ J/ψ ($c\bar{c}$). No commonly used vector meson dominance assumption is used to define this photoproduction amplitude which is needed to describe the data near the threshold. The c – N potential v cN (r) is parameterized in a form such that the predicted V J/ψN (r) at large distances has the same Yukawa potential form extracted from a lattice QCD (LQCD) calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)]. The parameters of vcN are determined by fitting the total cross-section data of Jefferson Laboratory (JLab) by performing calculations that include J/ψ – N final-state interactions (FSI). The resulting differential cross sections dσ /dt are found in good agreements with the data. It is shown that the FSI effects dominate the cross section in the very near-threshold region, allowing for sensitive testing of the predicted J/ψ – N scattering amplitudes. By imposing the constraints of J/ψ – N potential extracted from the LQCD calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)], we have obtained three J/ψ – N potentials which fit the JLab data equally well. The resulting J/ψ – N scattering lengths are in the range of a = [-0.05, -0.25] fm. With the determined v cN (r) and the wave functions generated from the same CQM, the constructed model is used to predict the cross sections of photoproduction of η c (1S) and ψ(2S) mesons for future experimental tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Strength anisotropy of shales deformed under uppermost crustal conditions

Conventional triaxial tests were performed on three sets of samples of Tournemire shale along different orientations relative to bedding (0°, 45°, and 90°). Experiments were carried out up to failure at increasing confining pressures ranging from 2.5 to 160 MPa, at strain rates ranging between 3 × 10 -7 s -1 and 3 × 10 -5 s -1 . This allowed us to determine the entire anisotropic elastic compliance matrix as a function of confining pressure. Results show that the orientation of principal stress relative to bedding plays an important role on the brittle strength, with 45° orientation being the weakest. We fit our results with a wing crack micromechanical model and an anisotropic fracture toughness. We found low values of internal friction coefficient and apparent friction coefficient in agreement with friction coefficient of clay minerals (between 0.2 and 0.3) and values of K /c comparable to that already published in the literature. We also showed that strain rate has a strong impact on peak stress and that dilatancy appears right before failure and hence highlighting the importance of plasticity mechanisms. Although brittle failure was systematically observed, stress drops and associated slips were slow and deformation always remained aseismic (no acoustic emission were detected). Lastly, this confirms that shales are good lithological candidates for shallow crust aseismic creep and slow slip events.

58 GEOSCIENCES↗

Revealing the Synergy between Oxide and Alloy Phases on the Performance of Bimetallic In–Pd Catalysts for CO 2 Hydrogenation to Methanol

In 2 O 3 has recently emerged as a promising catalyst for methanol synthesis from CO2. In this work, we present the promotional effect of Pd on this catalyst and investigate structure–performance relationships using in situ X-ray spectroscopy, ex situ characterization, and microkinetic modeling. Catalysts were synthesized with varying In:Pd ratios (1:0, 2:1, 1:1, 1:2, 0:1) and tested for methanol synthesis from CO 2 /H 2 at 40 bar and 300 °C. In:Pd(2:1)/SiO 2 shows the highest activity (5.1 μmol MeOH/g InPd s) and selectivity toward methanol (61%). While all bimetallic catalysts had enhanced catalytic performance, characterization reveals methanol synthesis was maximized when the catalyst contained both In–Pd intermetallic compounds and an indium oxide phase. Experimental results and density functional theory suggest the active phase arises from a synergy between the indium oxide phase and a bimetallic In–Pd particle with a surface enrichment of indium. We show that the promotion observed in the In–Pd system is extendable to non precious metal containing binary systems, in particular In–Ni, which displayed similar composition–activity trends to the In–Pd system. As a result, both palladium and nickel were found to form bimetallic catalysts with enhanced methanol activity and selectivity relative to that of indium oxide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimizing Binding Energies of Key Intermediates for CO 2 Hydrogenation to Methanol over Oxide-Supported Copper

Rational optimization of catalytic performance has been one of the major challenges in catalysis. In this paper, we report a bottom-up study on the ability of TiO 2 and ZrO 2 to optimize the CO 2 conversion to methanol on Cu, using combined density functional theory (DFT) calculations, kinetic Monte Carlo (KMC) simulations, in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) measurements, and steady-state flow reactor tests. Furthermore, the theoretical results from DFT and KMC agree with in situ DRIFTS measurements, showing that both TiO 2 and ZrO 2 help to promote methanol synthesis on Cu via carboxyl intermediates and the reverse water–gas-shift (RWGS) pathway; the formate intermediates, on the other hand, likely act as a spectator eventually. The origin of the superior promoting effect of ZrO 2 is associated with the fine-tuning capability of reduced Zr 3+ at the interface, being able to bind the key reaction intermediates, e.g. *CO 2 , *CO, *HCO, and *H 2 CO, moderately to facilitate methanol formation. Our study demonstrates the importance of synergy between theory and experiments to elucidate the complex reaction mechanisms of CO 2 hydrogenation for the realization of a better catalyst by design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A multi-functional testing instrument for heat assisted magnetic recording media

With recent developments in heat assisted magnetic recording (HAMR), characterization of HAMR media is becoming very important. We present a multi-functional instrument for testing HAMR media, which integrates HAMR writing, reading, and a micro-magneto-optic Kerr effect (μ-MOKE) testing function. A potential application of the present instrument is to make temperature dependent magnetic property measurement using a pump-probe configuration. In the measurement, the media is heated up by a heating (intense) beam while a testing (weak) beam is overlapped with the heating beam for MOKE measurement. By heating the media with different heating beam power, magnetic measurements by MOKE at different temperatures can be performed. Compared to traditional existing tools such as the vibrating sample magnetometer, the present instrument provides localized and efficient heating at the measurement spot. The integration of HAMR writing and μ-MOKE system can also facilitate a localized full investigation of the magnetic media by potential correlation of HAMR head independent write/read performance to localized magnetic properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vertical and adiabatic excitations in anthracene from quantum Monte Carlo: Constrained energy minimization for structural and electronic excited-state properties in the JAGP ansatz

We study the ionization energy, electron affinity, and the π → π{sup ∗} ({sup 1}L{sub a}) excitation energy of the anthracene molecule, by means of variational quantum Monte Carlo (QMC) methods based on a Jastrow correlated antisymmetrized geminal power (JAGP) wave function, developed on molecular orbitals (MOs). The MO-based JAGP ansatz allows one to rigorously treat electron transitions, such as the HOMO → LUMO one, which underlies the {sup 1}L{sub a} excited state. We present a QMC optimization scheme able to preserve the rank of the antisymmetrized geminal power matrix, thanks to a constrained minimization with projectors built upon symmetry selected MOs. We show that this approach leads to stable energy minimization and geometry relaxation of both ground and excited states, performed consistently within the correlated QMC framework. Geometry optimization of excited states is needed to make a reliable and direct comparison with experimental adiabatic excitation energies. This is particularly important in π-conjugated and polycyclic aromatic hydrocarbons, where there is a strong interplay between low-lying energy excitations and structural modifications, playing a functional role in many photochemical processes. Anthracene is an ideal benchmark to test these effects. Its geometry relaxation energies upon electron excitation are of up to 0.3 eV in the neutral {sup 1}L{sub a} excited state, while they are of the order of 0.1 eV in electron addition and removal processes. Significant modifications of the ground state bond length alternation are revealed in the QMC excited state geometry optimizations. Our QMC study yields benchmark results for both geometries and energies, with values below chemical accuracy if compared to experiments, once zero point energy effects are taken into account.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Method to Extract System-Independent Material Properties From Dual-Energy X-Ray CT

In 2016, we published a method for processing dual-energy computed tomography (DECT) data called system-independent rho-e/Z-e (ρ e /Z e ) or SIRZ. Using data from multiple DECT systems and spectra, SIRZ estimated the electron density ρ e and effective atomic number Z e (based on published X-ray cross sections), for a set of known specimens. However, the decomposition process required complex spectral modeling of the DECT system, which made SIRZ difficult to implement and automate. This paper describes the subsequent work on “SIRZ-2” to simplify the spectral modeling, automate the process, and improve its range and versatility. The SIRZ-2 basis functions are more accurate for the X-ray energy range (<;200 keV) and materials that are common to DECT applications in security and nondestructive characterization. To evaluate the SIRZ-2 performance, eight well-characterized specimens (with Z e of 6-20) were scanned with multiple spectral pairs (up to 200 keV) on seven different DECT systems, including a commercial airport luggage scanner. For these varied tests, the average SIRZ-2 relative errors for ρ e estimates were four times lower (0.7% compared to 3.0% for SIRZ), while average Z e relative errors were comparably low for both methods (<;1.5%). Average precision errors were also low (<;2.3% for both). SIRZ-2 produces accurate, repeatable, and system-independent material property estimates of (ρ e , Z e ) that can be directly compared across DECT systems and over time.

42 ENGINEERING↗

Processing of Onion-like Carbon for Electrochemical Capacitors

Multi-shell fullerenes known as onion-like carbon (OLC) are especially attractive in applications relative to energy storage, such as electrochemical capacitors, due to a near-spherical shape of particles, their nanoscale diameters and high conductivity leading to fast rate performance. Because of this, onion-like carbon can be fabricated into electrodes, used as a conductive additive, and may have potential in composites and additive manufacturing. However due to agglomeration of OLC particles, creating a stable, aqueous dispersion for ink production or formulating composites proves challenging. Also, we explore how attrition milling, acid treatment, and probe sonication can be employed to decrease agglomeration and provide colloidal stability in aqueous media. We also investigate how the electrochemical performance changes with each processing step as well as the treatments in succession. When tested in electrochemical capacitors, the processing increases the capacitance by a factor of three, due to an added pseudocapacitive contribution which is not present in untreated OLC. As a result, the processing of OLC proves to be advantageous for the production of stable, aqueous solutions, which also exhibit improved electrochemical properties suitable for functional inks, conductive additives, and fabrication of composite electrodes.

25 ENERGY STORAGE↗

Disruption of Dark Matter Minihalos in the Milky Way Environment: Implications for Axion Miniclusters and Early Matter Domination

Many theories of dark matter beyond the weakly interacting massive particles paradigm feature an enhanced matter power spectrum on subparsec scales, leading to the formation of dense dark matter minihalos. Future local observations are promising to search for and constrain such substructures. The survival probability of these dense minihalos in the Milky Way environment is crucial for interpreting local observations. In this work, we investigate two environmental effects: stellar disruption and (smooth) tidal disruption. These two mechanisms are studied using semianalytic models and idealized N-body simulations. For stellar disruption, we perform a series of N-body simulations of isolated minihalo–star encounters to test and calibrate analytic models of stellar encounters before applying the model to the realistic Milky Way disk environment. For tidal disruption, we perform N-body simulations to confirm the effectiveness of the analytic treatment. Finally, we propose a framework to combine the hierarchical assembly and infall of minihalos to the Milky Way with the late-time disruption mechanisms. We make predictions for the mass functions of minihalos in the Milky Way. The mass survival fraction (at M$_{mh}$ ≥ 10$^{−12}$M$_{⊙}$) of dense dark matter minihalos, e.g., for axion miniclusters and minihalos from early matter domination, is ∼60% with the relatively low-mass, compact population surviving. The survival fraction is insensitive to the detailed model parameters. We discuss various implications of the framework and future direct detection prospects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Implementation of Augmented Teleautonomy on a Hydraulic Manipulator for Nuclear Facility D and D - 19602

Robotic and remote system have been used in nuclear applications to keep the human workers from high radiation environment and enhance task performance. However, applications in the more unstructured nuclear environments present more formidable challenges to robotics and remote systems. In particular, remote operation tasks for nuclear facility D and D (or dismantling) requires dexterous manipulation of heavy objects with simple and rugged robotic equipment. To address such challenge, this paper presents a new teleoperation method of augmented teleautonomy for remote operation of tools for D and D tasks. Rather than imposing complex functionalities on the slave robot, focused innovation is directed to providing effective operator aids in such a way as to draw more performance out of human operators. In this regard, the development entails the implementation of an enhanced telerobotic operation method by incorporating multi-modal augmented reality operator aids and autonomous motion primitives. This paper presents the relevant technology development utilizing a mixed-reality telerobotic test bed, and system integration and demonstration with a rugged, hydraulic slave robot system. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

FTIR and Raman Spectroscopy of Reduced Graphene Oxide Carbon Nanoribbons and Graphene Oxide on Copper (Cu)

In the past couple of decades, carbon nanomaterials produced from graphene and graphite have emerged as one of the more important recently discovered substances. These carbon nanomaterials have been found to exhibit high levels of mechanical strength, thermal resistance, and electrical conductivity, thus having a multitude of applications. Therefore, it is important to identify and analyze their traits in order to maximize their effectiveness. Two spectroscopic techniques employed for such analyses of these materials are Fourier-Transform Infrared (FTIR) and Raman Spectroscopy. The advantages of FTIR and Raman include nondestructive analysis, relatively quick measurements under normal atmospheric conditions, and enhanced sensitivity with high resolution capacity. FTIR Spectroscopy was used in this research project to determine the effects of thermal exfoliation and photoreduction on the functionality groups of a carbon nanomaterial. Raman Spectroscopy was used in this research project to determine the effects of thermal exfoliation and photoreduction on the intensity of the disorder peak, which is greater in intensity if the carbon nanomaterial exhibits greater amounts of functionalities or defects in the graphitic material. There were two sample types which were analyzed with both FTIR and Raman Spectroscopy for the research project: graphene nanoribbons and graphene oxide. The m-grade graphene nanoribbons (GNRs) were produced from m-grade carbon nanotubes (CNTs) via Tour's method. The graphene oxide (GO) was produced from graphene via Tour's method, sonicated and centrifuged, dip-coated onto Copper (Cu) foil, and then photoreduced under a UV light source. FTIR measurements and analyses were conducted using a Jasco FTIR Spectrometer. Raman measurements and analyses were conducted using a Jasco Laser Raman Spectrometer. Tasks and techniques for this project included photoreduction, dip-coating, corrosion testing, FTIR and Raman data collection and analysis, and improvement of data collection for the Raman Spectrometer. GNRs exfoliated and reduced (deoxygenation) at 270 deg. C in air showed a significant decrease in O functionalities according to FTIR. However, Raman revealed that a large increase in graphene planar and edge defects are present as indicated by the disorder band at 1340 cm-1. Further optimization could be performed by thermal exfoliation at lower temperatures assuming a need for fewer structural defects. Graphene oxide coatings showed clear reduction with an increase with photoreduction time. Raman revealed that the reduced graphene oxide films had less disruption in the planar phonon lattice vibrations which indicates that the structure is returning to sp2 hybridization of pristine graphene without adding structural defects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hardware-in-the-loop Testing of Network Protectors for Low-Voltage Networks with Distributed Energy Resources

In this report, we developed and validated a network protector relay digital twin model and interfaced a commonly used network protector relay hardware with our real-time simulation system. Hardware-in-the-loop protection studies are performed to assess the impact of distributed energy resources (DER) and benchmark a rate-of-change-based mitigation strategy. Simulation results suggest that the network protector reverse trip and auto-reclose functions are negatively impacted by the high distributed energy resource penetration. To accommodate DER backfeed while remaining secure and reliable for faults on primary feeders, we recommend options for a rate-of-change-based blocking scheme and a protection setting change. Finally, future mitigation ideas and standard revisions are discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Long-term stability of the Leksell Gamma Knife{sup ®} Perfexion™ patient positioning system (PPS)

Purpose: To assess the long-term mechanical stability and accuracy of the patient positioning system (PPS) of the Leksell Gamma Knife{sup ®} Perfexion™ (LGK PFX). Methods: The mechanical stability of the PPS of the LGK PFX was evaluated using measurements obtained between September 2007 and June 2011. Three methods were employed to measure the deviation of the coincidence of the radiological focus point (RFP) and the PPS calibration center point (CCP). In the first method, the onsite diode test tool with single diode detector was used together with the 4 mm collimator on a daily basis. In the second method, a service diode test tool with three diode detectors was used biannually at the time of the routine preventive maintenance. The test performed with the service diode test tool measured the deviations for all three collimators 4, 8, and 16 mm and also for three different positions of the PPS. The third method employed the conventional film pin-prick method. This test was performed annually for the 4 mm collimator at the time of the routine annual QA. To estimate the effect of the patient weight on the performance of the PPS, the focus precision tests were also conducted with varying weights on the PPS using a set of lead bricks. Results: The average deviations measured from the 641 daily focus precision tests were 0.1 ± 0.1, 0.0 ± 0.0, and 0.0 ± 0.0 mm, respectively, for the 4 mm collimator in the X (left/right of the patient), Y (anterior/posterior of the patient), and Z (superior/inferior of the patient) directions. The average of the total radial deviations as measured during ten semiannual measurements with the service diode test tool were 0.070 ± 0.029, 0.060 ± 0.022, and 0.103 ± 0.028 mm, respectively for the central, long, and short diodes for the 4 mm collimator. Similarly, the average total radial deviations measured during the semiannual measurements for the 4, 8, and 16 mm collimators and using the central diode were 0.070 ± 0.029, 0.097 ± 0.025, 0.159 ± 0.028 mm, respectively. The average values of the deviations as obtained from the five annual film pin-prick tests for the 4 mm collimator were 0.10 ± 0.06, 0.06 ± 0.09, and 0.03 ± 0.03 mm for the X, Y, Z stereotactic directions, respectively. Only a minor change was observed in the total radial deviations of the PPS as a function of the simulated patient weight up to 202 kg on the PPS. Conclusions: Excellent long-term mechanical stability and high accuracy was observed for the PPS of the LGK PFX. No PPS recalibration or any adjustment in the PPS was needed during the monitored period of time. Similarly, the weight on the PPS did not cause any significant disturbance in the performance of the PPS for up to 202 kg simulated patient weight.

61 RADIATION PROTECTION AND DOSIMETRY↗

DEVELOPMENT AND TESTING OF FAULT-DIAGNOSIS ALGORITHMS FOR REACTOR PLANT SYSTEMS

Argonne National Laboratory is further developing fault diagnosis algorithms for use by the operator of a nuclear plant to aid in improved monitoring of overall plant condition and performance. The objective is better management of plant upsets through more timely, informed decisions on control actions with the ultimate goal of improved plant safety, production, and cost management. Integration of these algorithms with visual aids for operators is taking place through a collaboration under the concept of an operator advisory system. This is a software entity whose purpose is to manage and distill the enormous amount of information an operator must process to understand the plant state, particularly in off-normal situations, and how the state trajectory will unfold in time. The fault diagnosis algorithms were exhaustively tested using computer simulations of twenty different faults introduced into the chemical and volume control system (CVCS) of a pressurized water reactor (PWR). The algorithms are unique in that each new application to a facility requires providing only the piping and instrumentation diagram (PID) and no other plant-specific information; a subject-matter expert is not needed to install and maintain each instance of an application. The testing approach followed accepted procedures for verifying and validating software. It was shown that the code satisfies its functional requirement which is to accept sensor information, identify process variable trends based on this sensor information, and then to return an accurate diagnosis based on chains of rules related to these trends. The validation and verification exercise made use of GPASS, a one-dimensional systems code, for simulating CVCS operation. Plant components were failed and the code generated the resulting plant response. Parametric studies with respect to the severity of the fault, the richness of the plant sensor set, and the accuracy of sensors were performed as part of the validation exercise. The background and overview of the software will be presented to give an overview of the approach. Following, the verification and validation effort using the GPASS code for simulation of plant transients including a sensitivity study on important parameters will be presented

Fault Diagnosis Algorithm↗

Electron dynamics in complex environments with real-time time dependent density functional theory in a QM-MM framework

This article presents a time dependent density functional theory (TDDFT) implementation to propagate the Kohn-Sham equations in real time, including the effects of a molecular environment through a Quantum-Mechanics Molecular-Mechanics (QM-MM) hamiltonian. The code delivers an all-electron description employing Gaussian basis functions, and incorporates the Amber force-field in the QM-MM treatment. The most expensive parts of the computation, comprising the commutators between the hamiltonian and the density matrix—required to propagate the electron dynamics—, and the evaluation of the exchange-correlation energy, were migrated to the CUDA platform to run on graphics processing units, which remarkably accelerates the performance of the code. The method was validated by reproducing linear-response TDDFT results for the absorption spectra of several molecular species. Two different schemes were tested to propagate the quantum dynamics: (i) a leap-frog Verlet algorithm, and (ii) the Magnus expansion to first-order. These two approaches were confronted, to find that the Magnus scheme is more efficient by a factor of six in small molecules. Interestingly, the presence of iron was found to seriously limitate the length of the integration time step, due to the high frequencies associated with the core-electrons. This highlights the importance of pseudopotentials to alleviate the cost of the propagation of the inner states when heavy nuclei are present. Finally, the methodology was applied to investigate the shifts induced by the chemical environment on the most intense UV absorption bands of two model systems of general relevance: the formamide molecule in water solution, and the carboxy-heme group in Flavohemoglobin. In both cases, shifts of several nanometers are observed, consistently with the available experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Towards 3D printed multifunctional immobilization for proton therapy: Initial materials characterization

Purpose: 3D printing technology is investigated for the purpose of patient immobilization during proton therapy. It potentially enables a merge of patient immobilization, bolus range shifting, and other functions into one single patient-specific structure. In this first step, a set of 3D printed materials is characterized in detail, in terms of structural and radiological properties, elemental composition, directional dependence, and structural changes induced by radiation damage. These data will serve as inputs for the design of 3D printed immobilization structure prototypes. Methods: Using four different 3D printing techniques, in total eight materials were subjected to testing. Samples with a nominal dimension of 20 × 20 × 80 mm{sup 3} were 3D printed. The geometrical printing accuracy of each test sample was measured with a dial gage. To assess the mechanical response of the samples, standardized compression tests were performed to determine the Young’s modulus. To investigate the effect of radiation on the mechanical response, the mechanical tests were performed both prior and after the administration of clinically relevant dose levels (70 Gy), multiplied with a safety factor of 1.4. Dual energy computed tomography (DECT) methods were used to calculate the relative electron density to water ρ{sub e}, the effective atomic number Z{sub eff}, and the proton stopping power ratio (SPR) to water SPR. In order to validate the DECT based calculation of radiological properties, beam measurements were performed on the 3D printed samples as well. Photon irradiations were performed to measure the photon linear attenuation coefficients, while proton irradiations were performed to measure the proton range shift of the samples. The directional dependence of these properties was investigated by performing the irradiations for different orientations of the samples. Results: The printed test objects showed reduced geometric printing accuracy for 2 materials (deviation > 0.25 mm). Compression tests yielded Young’s moduli ranging from 0.6 to 2940 MPa. No deterioration in the mechanical response was observed after exposure of the samples to 100 Gy in a therapeutic MV photon beam. The DECT-based characterization yielded Z{sub eff} ranging from 5.91 to 10.43. The SPR and ρ{sub e} both ranged from 0.6 to 1.22. The measured photon attenuation coefficients at clinical energies scaled linearly with ρ{sub e}. Good agreement was seen between the DECT estimated SPR and the measured range shift, except for the higher Z{sub eff}. As opposed to the photon attenuation, the proton range shifting appeared to be printing orientation dependent for certain materials. Conclusions: In this study, the first step toward 3D printed, multifunctional immobilization was performed, by going through a candidate clinical workflow for the first time: from the material printing to DECT characterization with a verification through beam measurements. Besides a proof of concept for beam modification, the mechanical response of printed materials was also investigated to assess their capabilities for positioning functionality. For the studied set of printing techniques and materials, a wide variety of mechanical and radiological properties can be selected from for the intended purpose. Moreover the elaborated hybrid DECT methods aid in performing in-house quality assurance of 3D printed components, as these methods enable the estimation of the radiological properties relevant for use in radiation therapy.

60 APPLIED LIFE SCIENCES↗