Engineering Papers⌕ Search

Engineering topics

Cawkwell, Marc J.

Publications and source records attributed to Cawkwell, Marc J..

13 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dependence of the Elastic Stiffness Tensors of PETN, α‐RDX, γ‐RDX, ϵ‐RDX, ϵ‐CL‐20, DAAF, FOX‐7, and β‐HMX on Hydrostatic Compression

Abstract The dependence of the components of the elastic stiffness tensors (or elastic constants) of the organic explosives PETN, RDX, CL‐20, DAAF, FOX‐7, and HMX on hydrostatic pressure up to 10 GPa have been computed using dispersion‐corrected density functional theory. We report the evolution of lattice parameters and the non‐zero stiffnesses for the tetragonal, orthorhombic, and monoclinic crystal symmetries. Linear and quadratic dependencies of the components of the elastic stiffness tensors on volumetric compression and hydrostatic pressure are tabulated for use in single crystal plasticity models.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Linear scaling pseudo Fermi-operator expansion for fractional occupation

Recursive Fermi-operator expansion methods for the calculation of the idempotent density matrix are valid only at zero electronic temperature with integer occupation numbers. We show how such methods can be modified to include fractional occupation numbers of an approximate or pseudo Fermi-Dirac distribution and how the corresponding entropy term of the free energy is calculated. The proposed methodology is demonstrated and evaluated for different electronic structure methods including density functional tight-binding theory, Kohn-Sham density functional theory using numerical orbitals, and quantum chemistry Hartree-Fock theory using Gaussian basis functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Examining the Chemical and Structural Properties that Influence the Sensitivity of Energetic Nitrate Esters

The sensitivity of explosives is controlled by factors that span from intrinsic chemical reactivity and chemical intramolecular effects to mesoscale structure and defects, and has been a topic of extensive study for over 50 years. Due to these complex competing chemical and physical elements, a unifying relationship between molecular framework, crystal structure, and sensitivity has yet to be developed. In order to move towards this goal, ideally experimental studies should be performed on systems with small, systematic structural modifications, with modeling utilized to interpret experimental results. Pentaerythritol tetranitrate (PETN) is a common nitrate ester explosive that has been widely studied due to its use in military and commercial explosives. We have synthesized PETN derivatives with modified sensitivity characteristics by substituting the CCH 2 ONO 2 moiety with other substituents, including CH, CNH 2 , CNH3X, CCH 3 , and PO. We relate the handling sensitivity properties of each PETN derivative to its structural properties, and discuss the potential roles of thermodynamic properties such as heat capacity and heat of formation, thermal stability, crystal structure, compressibility, and inter- and intramolecular hydrogen bonding on impact sensitivity. Reactive molecular dynamics (MD) simulations of the C/H/N/O-based PETN-derivatives have been performed under cook-off conditions that mimic those accessed in impact tests. As a result, these simulations infer how changes in chemistry affect the subsequent decomposition pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Progress Version 1.0

Software library for parallel rapid order N and graph-based recursive electronic structure solvers with applications to material science, chemistry, and molecular biology.

Niklasson, Anders M.↗