Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A transferable force field for CdS-CdSe-PbS-PbSe solid systems

A transferable force field for the PbSe-CdSe solid system using the partially charged rigid ion model has been successfully developed and was used to study the cation exchange in PbSe-CdSe heteronanocrystals [A. O. Yalcin et al., “Atomic resolution monitoring of cation exchange in CdSe-PbSe heteronanocrystals during epitaxial solid-solid-vapor growth,” Nano Lett. 14, 3661–3667 (2014)]. In this work, we extend this force field by including another two important binary semiconductors, PbS and CdS, and provide detailed information on the validation of this force field. The parameterization combines Bader charge analysis, empirical fitting, and ab initio energy surface fitting. When compared with experimental data and density functional theory calculations, it is shown that a wide range of physical properties of bulk PbS, PbSe, CdS, CdSe, and their mixed phases can be accurately reproduced using this force field. The choice of functional forms and parameterization strategy is demonstrated to be rational and effective. This transferable force field can be used in various studies on II-VI and IV-VI semiconductor materials consisting of CdS, CdSe, PbS, and PbSe. Here, we demonstrate the applicability of the force field model by molecular dynamics simulations whereby transformations are initiated by cation exchange.

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

Polychromatic sparse image reconstruction and mass attenuation spectrum estimation via B-spline basis function expansion

We develop a sparse image reconstruction method for polychromatic computed tomography (CT) measurements under the blind scenario where the material of the inspected object and the incident energy spectrum are unknown. To obtain a parsimonious measurement model parameterization, we first rewrite the measurement equation using our mass-attenuation parameterization, which has the Laplace integral form. The unknown mass-attenuation spectrum is expanded into basis functions using a B-spline basis of order one. We develop a block coordinate-descent algorithm for constrained minimization of a penalized negative log-likelihood function, where constraints and penalty terms ensure nonnegativity of the spline coefficients and sparsity of the density map image in the wavelet domain. This algorithm alternates between a Nesterov’s proximal-gradient step for estimating the density map image and an active-set step for estimating the incident spectrum parameters. Numerical simulations demonstrate the performance of the proposed scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structural stability and defect energetics of ZnO from diffusion quantum Monte Carlo

We have applied the many-body ab initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure and the energetics of the oxygen vacancy, zinc interstitial, and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory (DFT) approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O{sub 2}, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.

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

Simulation of atmospheric temperature effects on cosmic ray muon flux

The collision between a cosmic ray and an atmosphere nucleus produces a set of secondary particles, which will decay or interact with other atmosphere elements. This set of events produced a primary particle is known as an extensive air shower (EAS) and is composed by a muonic, a hadronic and an electromagnetic component. The muonic flux, produced mainly by pions and kaons decays, has a dependency with the atmosphere’s effective temperature: an increase in the effective temperature results in a lower density profile, which decreases the probability of pions and kaons to interact with the atmosphere and, consequently, resulting in a major number of meson decays. Such correlation between the muon flux and the atmosphere’s effective temperature was measured by a set of experiments, such as AMANDA, Borexino, MACRO and MINOS. This phenomena can be investigated by simulating the final muon flux produced by two different parameterizations of the isothermal atmospheric model in CORSIKA, where each parameterization is described by a depth function which can be related to the muon flux in the same way that the muon flux is related to the temperature. This research checks the agreement among different high energy hadronic interactions models and the physical expected behavior of the atmosphere temperature effect by analyzing a set of variables, such as the height of the primary interaction and the difference in the muon flux.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

CAMELOT: A machine learning approach for coarse-grained simulations of aggregation of block-copolymeric protein sequences

We report the development and deployment of a coarse-graining method that is well suited for computer simulations of aggregation and phase separation of protein sequences with block-copolymeric architectures. Our algorithm, named CAMELOT for Coarse-grained simulations Aided by MachinE Learning Optimization and Training, leverages information from converged all atom simulations that is used to determine a suitable resolution and parameterize the coarse-grained model. To parameterize a system-specific coarse-grained model, we use a combination of Boltzmann inversion, non-linear regression, and a Gaussian process Bayesian optimization approach. The accuracy of the coarse-grained model is demonstrated through direct comparisons to results from all atom simulations. We demonstrate the utility of our coarse-graining approach using the block-copolymeric sequence from the exon 1 encoded sequence of the huntingtin protein. This sequence comprises of 17 residues from the N-terminal end of huntingtin (N17) followed by a polyglutamine (polyQ) tract. Simulations based on the CAMELOT approach are used to show that the adsorption and unfolding of the wild type N17 and its sequence variants on the surface of polyQ tracts engender a patchy colloid like architecture that promotes the formation of linear aggregates. These results provide a plausible explanation for experimental observations, which show that N17 accelerates the formation of linear aggregates in block-copolymeric N17-polyQ sequences. The CAMELOT approach is versatile and is generalizable for simulating the aggregation and phase behavior of a range of block-copolymeric protein sequences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Redesign of the DFT/MRCI Hamiltonian

The combined density functional theory and multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke [J. Chem. Phys. 111, 5645 (1999)] is a well-established semi-empirical quantum chemical method for efficiently computing excited-state properties of organic molecules. As it turns out, the method fails to treat bi-chromophores owing to the strong dependence of the parameters on the excitation class. In this work, we present an alternative form of correcting the matrix elements of a MRCI Hamiltonian which is built from a Kohn-Sham set of orbitals. It is based on the idea of constructing individual energy shifts for each of the state functions of a configuration. The new parameterization is spin-invariant and incorporates less empirism compared to the original formulation. By utilizing damping techniques together with an algorithm of selecting important configurations for treating static electron correlation, the high computational efficiency has been preserved. The robustness of the original and redesigned Hamiltonians has been tested on experimentally known vertical excitation energies of organic molecules yielding similar statistics for the two parameterizations. Besides that, our new formulation is free from artificially low-lying doubly excited states, producing qualitatively correct and consistent results for excimers. The way of modifying matrix elements of the MRCI Hamiltonian presented here shall be considered as default choice when investigating photophysical processes of bi-chromophoric systems such as singlet fission or triplet-triplet upconversion.

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

Scaling magnetized liner inertial fusion on Z and future pulsed-power accelerators

The MagLIF (Magnetized Liner Inertial Fusion) concept [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)] has demonstrated fusion–relevant plasma conditions [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] on the Z accelerator with a peak drive current of about 18 MA. We present 2D numerical simulations of the scaling of MagLIF on Z as a function of drive current, preheat energy, and applied magnetic field. The results indicate that deuterium-tritium (DT) fusion yields greater than 100 kJ could be possible on Z when all of these parameters are at the optimum values: i.e., peak current = 25 MA, deposited preheat energy = 5 kJ, and B{sub z} = 30 T. Much higher yields have been predicted [S. A. Slutz and R. A. Vesey, Phys. Rev. Lett. 108, 025003 (2012)] for MagLIF driven with larger peak currents. Two high performance pulsed-power accelerators (Z300 and Z800) based on linear-transformer-driver technology have been designed [W. A. Stygar et al., Phys. Rev. ST Accel. Beams 18, 110401 (2015)]. The Z300 design would provide 48 MA to a MagLIF load, while Z800 would provide 65 MA. Parameterized Thevenin-equivalent circuits were used to drive a series of 1D and 2D numerical MagLIF simulations with currents ranging from what Z can deliver now to what could be achieved by these conceptual future pulsed-power accelerators. 2D simulations of simple MagLIF targets containing just gaseous DT have yields of 18 MJ for Z300 and 440 MJ for Z800. The 2D simulated yield for Z800 is increased to 7 GJ by adding a layer of frozen DT ice to the inside of the liner.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Characteristics of energy exchange between inter- and intramolecular degrees of freedom in crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) with implications for coarse-grained simulations of shock waves in polyatomic molecular crystals

In this report, we characterize the kinetics and dynamics of energy exchange between intramolecular and intermolecular degrees of freedom (DoF) in crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). All-atom molecular dynamics (MD) simulations are used to obtain predictions for relaxation from certain limiting initial distributions of energy between the intra- and intermolecular DoF. The results are used to parameterize a coarse-grained Dissipative Particle Dynamics at constant Energy (DPDE) model for TATB. Each TATB molecule in the DPDE model is represented as an all-atom, rigid-molecule mesoparticle, with explicit external (molecular translational and rotational) DoF and coarse-grained implicit internal (vibrational) DoF. In addition to conserving linear and angular momentum, the DPDE equations of motion conserve the total system energy provided that particles can exchange energy between their external and internal DoF. The internal temperature of a TATB molecule is calculated using an internal equation of state, which we develop here, and the temperatures of the external and internal DoF are coupled using a fluctuation-dissipation relation. The DPDE force expression requires specification of the input parameter σ that determines the rate at which energy is exchanged between external and internal DoF. We adjusted σ based on the predictions for relaxation processes obtained from MD simulations. The parameterized DPDE model was employed in large-scale simulations of shock compression of TATB. We show that the rate of energy exchange governed by σ can significantly influence the transient behavior of the system behind the shock.

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

Efficient magnetic fields for supporting toroidal plasmas

The magnetic field that supports tokamak and stellarator plasmas must be produced by coils well separated from the plasma. However, the larger the separation, the more difficult it is to produce a given magnetic field in the plasma region, so plasma configurations should be chosen that can be supported as efficiently as possible by distant coils. The efficiency of an externally generated magnetic field is a measure of the field's shaping component magnitude at the plasma compared to the magnitude near the coils; the efficiency of a plasma equilibrium can be measured using the efficiency of the required external shaping field. Counterintuitively, plasma shapes with low curvature and spectral width may have low efficiency, whereas plasma shapes with sharp edges may have high efficiency. Two precise measures of magnetic field efficiency, which correctly identify such differences in difficulty, will be examined. These measures, which can be expressed as matrices, relate the externally produced normal magnetic field on the plasma surface to the either the normal field or current on a distant control surface. A singular value decomposition (SVD) of either matrix yields an efficiency ordered basis for the magnetic field distributions. Calculations are carried out for both tokamak and stellarator cases. For axisymmetric surfaces with circular cross-section, the SVD is calculated analytically, and the range of poloidal and toroidal mode numbers that can be controlled to a given desired level is determined. If formulated properly, these efficiency measures are independent of the coordinates used to parameterize the surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Coarse-grained models using local-density potentials optimized with the relative entropy: Application to implicit solvation

Bottom-up multiscale techniques are frequently used to develop coarse-grained (CG) models for simulations at extended length and time scales but are often limited by a compromise between computational efficiency and accuracy. The conventional approach to CG nonbonded interactions uses pair potentials which, while computationally efficient, can neglect the inherently multibody contributions of the local environment of a site to its energy, due to degrees of freedom that were coarse-grained out. This effect often causes the CG potential to depend strongly on the overall system density, composition, or other properties, which limits its transferability to states other than the one at which it was parameterized. Here, we propose to incorporate multibody effects into CG potentials through additional nonbonded terms, beyond pair interactions, that depend in a mean-field manner on local densities of different atomic species. This approach is analogous to embedded atom and bond-order models that seek to capture multibody electronic effects in metallic systems. We show that the relative entropy coarse-graining framework offers a systematic route to parameterizing such local density potentials. We then characterize this approach in the development of implicit solvation strategies for interactions between model hydrophobes in an aqueous environment.

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

Perspective: Ab initio force field methods derived from quantum mechanics

It is often desirable to accurately and efficiently model the behavior of large molecular systems in the condensed phase (thousands to tens of thousands of atoms) over long time scales (from nanoseconds to milliseconds). In these cases, ab initio methods are difficult due to the increasing computational cost with the number of electrons. A more computationally attractive alternative is to perform the simulations at the atomic level using a parameterized function to model the electronic energy. Many empirical force fields have been developed for this purpose. However, the functions that are used to model interatomic and intermolecular interactions contain many fitted parameters obtained from selected model systems, and such classical force fields cannot properly simulate important electronic effects. Furthermore, while such force fields are computationally affordable, they are not reliable when applied to systems that differ significantly from those used in their parameterization. They also cannot provide the information necessary to analyze the interactions that occur in the system, making the systematic improvement of the functional forms that are used difficult. Ab initio force field methods aim to combine the merits of both types of methods. The ideal ab initio force fields are built on first principles and require no fitted parameters. Ab initio force field methods surveyed in this perspective are based on fragmentation approaches and intermolecular perturbation theory. This perspective summarizes their theoretical foundation, key components in their formulation, and discusses key aspects of these methods such as accuracy and formal computational cost. The ab initio force fields considered here were developed for different targets, and this perspective also aims to provide a balanced presentation of their strengths and shortcomings. Finally, this perspective suggests some future directions for this actively developing area.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Rank reduced coupled cluster theory. I. Ground state energies and wavefunctions

We propose a compression of the opposite-spin coupled cluster doubles amplitudes of the form τ ab ij ≡ U V ia T VW U W jb , where U V ia are the n V -highest magnitude eigenvectors of the MP2 or MP3 doubles amplitudes. Together with a corresponding parameterization of the opposite-spin coupled cluster Lagrange multipliers of the form λ ij ab ≡ U V ia L VW U W jb , this yields a fully self-consistent parameterization of reduced-rank coupled cluster equations in terms of the Lagrangian L 0 ( T VW ,L VW ). Making this Lagrangian stationary with respect to the L VW parameters yields a perfectly determined set of equations for the T VW equations and coupled cluster energy. These equations can be solved using a Lyapunov equation for the first-order amplitude updates. We test this “rank-reduced coupled cluster” method for coupled cluster singles and doubles in medium sized molecules and find that substantial compression of the $\hat{T} $ 2 amplitudes is possible with acceptable accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

On alpha-particle transport in inertial fusion

Analytical theory and models of inertial fusion implosion use parameterized microphysics models. In this paper, we consider the DT alpha-particle transport, and report new parameterizations of the range, heating efficacy, and energy partition using modern stopping-power theory. Here, our resulting heating efficacy is lower than previously published results, which reduces the temperature and pressure generated by a dynamic implosion hot-spot evolution model, and shifts the burning-plasma regime boundary slightly farther from current experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A phaseless auxiliary-field quantum Monte Carlo perspective on the uniform electron gas at finite temperatures: Issues, observations, and benchmark study

We investigate the viability of the phaseless finite-temperature auxiliary-field quantum Monte Carlo (ph-FT-AFQMC) method for ab initio systems using the uniform electron gas as a model. Through comparisons with exact results and FT coupled cluster theory, we find that ph-FT-AFQMC is sufficiently accurate at high to intermediate electronic densities. We show, both analytically and numerically, that the phaseless constraint at FT is fundamentally different from its zero-temperature counterpart (i.e., ph-ZT-AFQMC), and generally, one should not expect ph-FT-AFQMC to agree with ph-ZT-AFQMC in the low-temperature limit. With an efficient implementation, we are able to compare exchange-correlation energies to the existing results in the thermodynamic limit and find that the existing parameterizations are highly accurate. In particular, we found that ph-FT-AFQMC exchange-correlation energies are in better agreement with a known parameterization than is restricted path-integral MC in the regime of Θ ≤ 0.5 and rs ≤ 2, which highlights the strength of ph-FT-AFQMC.

Lee, Joonho (ORCID:0000000296671081)↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2021

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗