Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependant density functional theory”

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 487 records · Page 27

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure and optical properties of halide double perovskites from a Wannier-localized optimally-tuned screened range-separated hybrid functional

Halide double perovskites are a chemically diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and ab initio many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here, in this study, we apply the WOT-SRSH functional to five halide double perovskites and compare the results with those obtained from other known functionals and previous GW calculations. We also use the approach as a starting point for GW calculations and we compute the band structures and optical absorption spectrum for Cs 2 AgBiBr 6 , using both time-dependent DFT and the GW-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.

36 MATERIALS SCIENCE↗

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗

Shock Waves in a Bose-Einstein Condensate

A paper presents a theoretical study of shock waves in a trapped Bose-Einstein condensate (BEC). The mathematical model of the BEC in this study is a nonlinear Schroedinger equation (NLSE) in which (1) the role of the wave function of a single particle in the traditional Schroedinger equation is played by a space- and time-dependent complex order parameter (x,t) proportional to the square root of the density of atoms and (2) the atoms engage in a repulsive interaction characterized by a potential proportional to | (x,t)|2. Equations that describe macroscopic perturbations of the BEC at zero temperature are derived from the NLSE and simplifying assumptions are made, leading to equations for the propagation of sound waves and the transformation of sound waves into shock waves. Equations for the speeds of shock waves and the relationships between jumps of velocity and density across shock fronts are derived. Similarities and differences between this theory and the classical theory of sound waves and shocks in ordinary gases are noted. The present theory is illustrated by solving the equations for the example of a shock wave propagating in a cigar-shaped BEC.

Kulikov, Igor↗

Solution-Phase Synthesis of PdH 0.706 Nanocubes with Enhanced Stability and Activity toward Formic Acid Oxidation

We report palladium is one of the few metals capable of forming hydrides, with the catalytic properties being dependent on the elemental composition and spatial distribution of H atoms in the lattice. Herein, we report a facile method for the complete transformation of Pd nanocubes into a stable phase made of PdH 0.706 by treating them with aqueous hydrazine at a concentration as low as 9.2 mM. Using formic acid oxidation (FAO) as a model reaction, we systematically investigated the structure–catalytic property relationship of the resultant nanocubes with different degrees of hydride formation. The current density at 0.4 V was enhanced by four times when the nanocubes were completely converted from Pd to PdH 0.706 . On the basis of a set of slab models with PdH(100) overlayers on Pd(100), we conducted density functional theory calculations to demonstrate that the degree of hybrid formation could influence both the activity and selectivity toward FAO by modulating the relative stability of formate (HCOO) and carboxyl (COOH) intermediates. Furthermore, this work provides a viable strategy for augmenting the performance of Pd-based catalysts toward various reactions without altering the loading of this scarce metal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Operando probing dynamic migration of copper carbonyl during electrocatalytic CO2 reduction

Single crystals and shape-controlled nanocrystals are well known to exhibit facet-dependent catalytic properties. However, few studies have investigated how those nanocrystals evolve and (de)activate during reactions, calling for the development of nanoscale time-resolved operando methods. In this context, we have designed Cu nanocubes as a model system to elucidate the underlying driving force of dynamic nanocatalyst reconstruction during the CO2 reduction reaction (CO2RR). Operando electrochemical liquid-cell scanning transmission electron microscopy (EC-STEM) and synchrotron-based X-ray spectroscopy reveal the size- and potential-dependent complete transformation from (100)-oriented Cu@Cu2O nanocubes to polycrystalline metallic Cu nanograins under CO2RR conditions. In addition, machine learning-assisted operando four-dimensional STEM reveals that large Cu nanograins derived from nanocubes form mainly crystalline domains, while their smaller counterparts are more amorphous due to faster evolution kinetics. In situ Raman spectroscopy and density functional theory calculations suggest that CO drives the ejection of single Cu atoms, resulting in few-nanometre Cu clusters and the surface migration of highly mobile copper carbonyl (Cu–CO) species. Combined, these multimodal operando methods and theoretical approaches pave the way for understanding the complex structural evolution of energy-related nanocatalysts under electrochemical conditions.

Yang, Yao↗

Single-Atom Pt Boosting Electrochemical Nonenzymatic Glucose Sensing on Ni(OH) 2 /N-Doped Graphene

Conventional nanomaterials in electrochemical non-enzymatic sensing face huge challenge due to their complex size-, surface- and composition- dependent catalytic properties and low active site density. In this work, we designed a single-atom Pt supported on Ni(OH) 2 nanoplates/nitrogen-doped graphene (Pt 1 /Ni(OH) 2 /NG) as the first example for constructing a single-atom catalyst based electrochemical non-enzymatic glucose sensor. The resulting Pt 1 /Ni(OH) 2 /NG exhibited a low anode peak potential of 0.48 V and high sensitivity of 220.75 μA mM -1 cm -2 towards glucose, which are 45 mV lower and 12 times higher than those of Ni(OH) 2 , respectively. The catalyst also showed excellent selectivity for several important interferences, short response time of 4.6 s, and high stability over 4 weeks. Experimental and density functional theory (DFT) calculated results reveal that the improved performance of Pt 1 /Ni(OH) 2 /NG could be attributed to stronger binding strength of glucose on single-atom Pt active centers and their surrounding Ni atoms, combined with fast electron transfer ability by the adding of the highly conductive NG. Finally, this research sheds light on the applications of SACs in the field of electrochemical non-enzymatic sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiplicative noise and the diffusion of conserved densities

Stochastic fluid dynamics governs the long time tails of hydrodynamic correlation functions, and the critical slowing down of relaxation phenomena in the vicinity of a critical point in the phase diagram. In this work we study the role of multiplicative noise in stochastic fluid dynamics. Multiplicative noise arises from the dependence of transport coefficients, such as the diffusion constants for charge and momentum, on fluctuating hydrodynamic variables. We study long time tails and relaxation in the diffusion of a conserved density (model B), and a conserved density coupled to the transverse momentum density (model H). Careful attention is paid to fluctuation-dissipation relations. We observe that multiplicative noise contributes at the same order as non-linear interactions in model B, but is a higher order correction to the relaxation of a scalar density and the tail of the stress tensor correlation function in model H.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Terahertz parametric amplification as a reporter of exciton condensate dynamics

Condensates are a hallmark of emergence in quantum materials with superconductors and charge density wave as prominent examples. An excitonic insulator (EI) is an intriguing addition to this library, exhibiting spontaneous condensation of electron-hole pairs. However, condensate observables can be obscured through parasitic coupling to the lattice. Time-resolved terahertz (THz) spectroscopy can disentangle such obscurants through measurement of the quantum dynamics. We target Ta 2 NiSe 5 , a putative room-temperature EI where electron-lattice coupling dominates the structural transition (T c = 326 K), hindering identification of excitonic correlations. A pronounced increase in the THz reflectivity manifests following photoexcitation and exhibits a BEC-like temperature dependence. This occurs well below the T c , suggesting a novel approach to monitor exciton condensate dynamics. Nonetheless, dynamic condensate-phonon coupling remains as evidenced by peaks in the enhanced reflectivity spectrum at select infrared-active phonon frequencies. This indicates that parametric reflectivity enhancement arises from phonon squeezing, validated using Fresnel-Floquet theory and density functional calculations. In conclusion, our results highlight that coherent dynamics can drive parametric stimulated emission with concomitant possibilities, including entangled THz photon generation.

36 MATERIALS SCIENCE↗

Ultrasonic wave propagation in two-phase media: Spherical inclusions

The scattering theory, recently developed via the extended method of equivalent inclusion, is used to study the propagation of time-harmonic waves in two-phase media of elastic matrix with randomly distributed elastic spherical inclusion materials. The elastic moduli and mass density of the composite medium are determined as functions of frequencies when given properties and concentration of the spheres and the matrix. Velocity and attenuation of ultrasonic waves in two-phase media are determined for cases of distributed spheres and localized damage. An averaging theorem that requires the equivalence of the strain energy and the kinetic energy between the effective medium and the original matrix with spherical inhomogeneities is employed to derive the effective moduli and mass density. The functional dependency of these quantities upon frequencies and concentration provides a method of data analysis in ultrasonic evaluation of material properties. Numerical results or moduli, velocity and/or attenuation as functions of concentration of inclusion material, or porosity, are graphically displayed.

Fu, L. S.↗

Ultrasonic wave propagation in two-phase media - Spherical inclusions

The scattering theory, recently developed via the extended method of equivalent inclusion, is used to study the propagation of time-harmonic waves in two-phase media of elastic matrix with randomly distributed elastic spherical inclusion materials. The elastic moduli and mass density of the composite medium are determined as functions of frequencies when given properties and concentration of the spheres and the matrix. Velocity and attenuation of ultrasonic waves in two-phase media are determined for cases of distributed spheres and localized damage. An averaging theorem that requires the equivalence of the strain energy and the kinetic energy between the effective medium and the original matrix with spherical inhomogeneities is employed to derive the effective moduli and mass density. The functional dependency of these quantities upon frequencies and concentration provides a method of data analysis in ultrasonic evaluation of material properties. Numerical results or moduli, velocity and/or attenuation as functions of concentration of inclusion material, or porosity, are graphically displayed.

Fu, L. S.↗

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

An Impedance Humidity Sensor Based on CVD Grown WSe 2 2D Films

2D materials-based planar devices have been sought after for gas and humidity sensor applications owing to their high sensitivity due to their ultrathin nature and increased surface area. Herein, 2D WSe 2 films-based planner devices were fabricated to evaluate their performance for potential application in humidity sensors. The structure and morphology of the WSe 2 sensing films were studied by Raman spectroscopy, photoluminescence, optical microscopy, and scanning electron microscopy (SEM) techniques. The relative humidity-dependent sensing performances of the 2D WSe 2 films were evaluated using current–voltage (I–V), relative impedance, and complex total impedance (electrochemical impedance spectroscopy, EIS) techniques. This sensor showed a very stable and repeatable performance investigated over the period of 18 months. The response and recovery times of the WSe 2 -based impedance sensors were found to be 40 and 30 s, respectively. The WSe 2 sensor exhibits a 2 ± 1% minimum hysteresis at a lower humidity level (50% RH) and around 8 ± 2% maximum hysteresis at relative humidity level (70% RH). This value of hysteresis is in the acceptable range during the initial investigation of any material to explore its potential for application in humidity sensors. The density functional theory (DFT) analysis was performed with the interaction from water molecules of the WSe 2 sample for a better understanding of the sensing mechanism of devices.

2D- WSe2↗

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗

Chemical evolution in spiral and irregular galaxies

A brief review of models of chemical evolution of the interstellar medium in our galaxy and other galaxies is presented. These models predict the time variation and radial dependence of chemical composition in the gas as function of the input parameters; initial mass function, stellar birth rate, chemical composition of mass lost by stars during their evolution (yields), and the existence of large scale mass flows, like infall from the halo, outflow to the intergalactic medium or radial flows within a galaxy. At present there is a considerable wealth of observational data on the composition of HII regions in spiral and irregular galaxies to constrain the models. Comparisons are made between theory and the observed physical conditions. In particular, studies of helium, carbon, nitrogen and oxygen abundances are reviewed. In many molecular clouds the information we have on the amount of H2 is derived from the observed CO column density, and a standard CO/H2 ratio derived for the solar neighborhood. Chemical evolution models and the observed variations in O/H and N/O values, point out the need to include these results in a CO/H2 relation that should be, at least, a function of the O/H ratio. This aspect is also discussed.

Torres-Peimbert, S.↗

Identifying Strain Stacking Boundaries between Multiphase Domains in Atomically Thin Two-Dimensional Magnets

Stacking engineering of van der Waals materials is an important strategy to control the materials’ properties, such as electronic correlations, ferroelectricity, and layer-dependent two-dimensional magnetism. A timely testbed for the study of the latter is atomically thin chromium trihalides (CrX 3 , X = Cl, Br, I). Notably, by understanding the sliding mechanism between different stacking sequences, control of the stacking arrangement, and thus magnetic properties in CrX 3 , can be achieved. Such insight, however, is currently lacking. Here, in this study, advanced electron microscopy methods are used to identify multiple stacking sequences corresponding to different bulk phases in atomically thin CrX 3 (X = Cl and Br) down to bilayer thickness and with lateral domain sizes as small as tens of nanometers. Indications of nanometer scale transitions and interactions at the stacking boundaries are found, including a universally preferred sliding direction that is consistent with density functional theory calculations and the strain fields at lateral heterostructure boundaries. This study demonstrates the necessity to consider local stacking structures when interpreting averaged magnetic properties measured with macroscale probes. Additionally, the preferred sliding direction insight from this study provides a strategy to control stacking sequence in atomically thin CrX 3 samples during the exfoliation and sample fabrication process.

DFT calculations↗

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics↗