We present a computational protocol, based on density matrix perturbation theory, to obtain non-adiabatic, frequency-dependent electron–phonon self-energies for molecules and solids. Our approach enables the evaluation of electron–phonon interaction using hybrid functionals, for spin-polarized systems, and the computational overhead to include dynamical and non-adiabatic terms in the evaluation of electron–phonon self-energies is negligible. We discuss results for molecules, as well as pristine and defective solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Time domain probabilistic seismic risk analysis using ground motion prediction equations of Fourier amplitude spectra

Modeling of Fourier amplitude spectra (FAS) of seismic motions has gained much attention in engineering seismology. In the past few years, several ground motion prediction equations (GMPEs) and inter-frequency correlation structure of FAS have been established. Due to many preferable characteristics of FAS, probabilistic seismic hazard/risk analysis is rapidly changing from ergodic, spectrum acceleration Sa(T 0 )-based approach to non-ergodic, site-specific, FAS-based approach. This paper presents time domain intrusive framework for probabilistic seismic risk analysis using GMPE of FAS. Herein, methodology for time domain stochastic ground motion modeling based on GMPEs of FAS is presented in some detail. The simulated uncertain motions are modeled as a random process and represented by polynomial chaos Karhunen-Loève expansion. The random process excitations are further propagated into the uncertain structural system using Galerkin stochastic finite element method (SFEM). Probabilistic evolution of structural response is solved, and such solution is used to develop seismic risk for any damage state. The presented framework is illustrated through seismic risk analysis of a four-story building subjected to possible earthquakes from two strike slip faults. The influences of the epistemic uncertainties in source stress drop Δσ and site attenuation κ0 on seismic risk are investigated. The need for non-ergodic seismic risk analysis with source-specific and site specific characterizations is emphasized.

58 GEOSCIENCES↗

DOE OSTI2021

Combined First-Principles Calculations of Electron–Electron and Electron–Phonon Self-Energies in Condensed Systems

In this work, we present a method to efficiently combine the computation of electron–electron and electron–phonon self-energies, which enables the evaluation of electron–phonon coupling at the G 0 W 0 level of theory for systems with hundreds of atoms. In addition, our approach, which is a generalization of a method recently proposed for molecules [J. Chem. Theory Comput.2018, 14, 6269–6275], enables the inclusion of nonadiabatic and temperature effects at no additional computational cost. We present results for diamond and defects in diamond and discuss the importance of numerically accurate G 0 W 0 band structures to obtain robust predictions of zero point renormalization (ZPR) of band gaps, and of the inclusion of nonadiabatic effects to accurately compute the ZPR of defect states in the band gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Effect of Graft Molecular Weight and Density on the Mechanical Properties of Polystyrene-Grafted Cellulose Nanocrystal Films

In this work, polymer-grafted nanoparticle (PGN) films were prepared from polystyrene (PS) grafted to rodlike cellulose nanocrystals (MxG-CNC-g-PS) with a controllable grafting density (0.03–0.25 chains/nm 2 ) and molecular weight (5–60 kg/mol). These nanorod-based PGNs are solution- and melt-processible, permitting access to one-component composite films with high nanofiller loadings (with up to 55 wt %). The impact of both grafted polymer density and molecular weight on the mechanical properties of the films was investigated and related to the polymer brush conformation: concentrated polymer brush (CPB), semidilute polymer brush (SDPB), or CPB core with SDPB corona (CPB/SDPB). The rubbery regime storage modulus (above $T_g$) showed 2 orders of magnitude increase, maximizing at a low degree of polymerization (N) and low grafting density (σ). Fracture toughness was maximized in samples with the grafted polymer in the SDPB or CPD/SDPB (higher N and relatively low σ) regime and showed enhancement relative to PS of molecular weight similar to the graft. In line with prior computation predictions, optimizing for both rubbery modulus and fracture toughness in such nanorod-based PGN films requires the polymers in the SDPB regime and CNC loading levels (ca. 50–60 wt %) that are difficult to attain in more traditional two-component CNC composites.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗

DOE OSTI2021

Quantum vibronic effects on the electronic properties of solid and molecular carbon

Here, we study the effect of quantum vibronic coupling on the electronic properties of carbon allotropes, including molecules and solids, by combining path integral first principles molecular dynamics (FPMD) with a colored noise thermostat. In addition to avoiding several approximations commonly adopted in calculations of electron-phonon coupling, our approach only adds a moderate computational cost to FPMD simulations and hence it is applicable to large supercells, such as those required to describe amorphous solids. We predict the effect of electron-phonon coupling on the fundamental gap of amorphous carbon, and we show that in diamond the zero-phonon renormalization of the band gap is larger than previously reported.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Modeling and simulation of earthquake soil structure interaction excited by inclined seismic waves

Presented is an application of wave potential formulation (WPF) together with domain reduction method (DRM) to modeling earthquake soil structure interaction (ESSI) behavior in horizontally layered ground under inclined incident seismic waves. Wave potential formulation is used to develop a spatially varying, inclined seismic wave field from incident Primary (P) and Secondary (S) waves that propagate through layered ground. Developed seismic wave field is then used to develop effective forces for Domain Reduction Method that are then used for analyzing ESSI response of a soil structure system. Developed methodology, called WPF-DRM, is verified using analytic solution for a free field response of layered ground subjected to inclined incident waves. Developed WPF-DRM methodology is illustrated through analysis of an ESSI response of a deeply embedded structure, a small modular reactor (SMR) subjected to incident S wave polarized in vertical plane (SV) with variation in inclinations and frequencies. Presented example highlights the influences of incident wave inclination and frequency on ESSI response of analyzed SMR.

58 GEOSCIENCES↗

DOE OSTI2020

Stretching-Induced Thermal Conductivity Change in Shape-Memory Polymer Composites

Explore the source record for details and available documents.

Hostler, Stephen R.↗

DOE OSTI2019

Improving the efficiency of G 0 W 0 calculations with approximate spectral decompositions of dielectric matrices

Recently, it was shown that the calculation of quasiparticle energies using the G 0 W 0 approximation can be performed without computing explicitly any virtual electronic states, by expanding the Green function and screened Coulomb interaction in terms of the eigenstates of the static dielectric matrix. Avoiding the evaluation of virtual electronic states leads to improved efficiency and ease of convergence of G 0 W 0 calculations. Here, we propose a further improvement of the efficiency of these calculations, based on an approximation of density-density response functions of molecules and solids. The approximation relies on the calculation of a subset of eigenvectors of the dielectric matrix using the kinetic operator instead of the full Hamiltonian, and it does not lead to any substantial loss of accuracy for the quasiparticle energies. The computational savings introduced by this approximation depend on the system, and they become more substantial as the number of electrons increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Energy dissipation analysis for inelastic reinforced concrete and steel beam-columns

Not Available

Yang, Han↗