Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics

We introduce a novel Newton–Krylov (NK)-based fully implicit algorithm for solving fluid flows in a wide range of flow conditions—from variable density nearly incompressible to supersonic shock dynamics. The key enabling feature of our all-speed solver is the ability to efficiently solve conservation laws by choosing a set of independent variables that produce a well-conditioned Jacobian matrix for the linear iterations of the global nonlinear iterative solver. In particular, instead of choosing to discretize the conservative variables (density, momentum, total energy), which is traditionally used in Eulerian high-speed compressible fluid dynamics, we demonstrate superior performance by discretizing the primitive variables—pressure–velocity–temperature in the very low-Mach flow limits or density–velocity–temperature/entropy in the shock dynamics range. Moreover, our method allows us to avoid direct inversion of the mass matrix in discrete time derivatives, which is usually an additional source for stiffness, especially pronounced when going to very high-order schemes with non-orthogonal basis functions. Here, we show robust solutions obtained for discontinuous finite element discretization up to seventh-order accuracy. Another important aspect of the solution algorithm is the Advection Upstream Splitting Method (AUSM), adopted to compute numerical fluxes within our reconstructed discontinuous Galerkin (rDG) spatial discretization scheme. The use of the low-Mach modification of the hyperbolic flux operator is found to be necessary for enabling robust simulations of very stiff liquids and metals for Mach numbers below M=10 -5 , which is well known to be very computationally challenging for compressible solvers. We demonstrate that our fully implicit rDG-NK solver with the AUSM + -up flux treatment produces efficient and high-resolution numerical solutions at all speeds, ranging from vanishing Mach numbers to transonic and supersonic, without substantial modifications of the solution procedures. (At high speed, we add limiting and use a simpler preconditioning of the Krylov solver.) Numerical examples include nearly incompressible constant-property flow past a backward-facing step with heat transfer, low-Mach variable-property channel flow of water at supercritical state, phase change and melt pool dynamics for laser spot welding and selective laser melting in additive manufacturing, and Mach 3 flow in a wind tunnel with a step.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory

The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density-matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies, relative to reference values from configuration-interaction-based calculations.

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

Computational modeling of latent-heat-storage in PCM modified interior plaster

The latent heat storage systems represent a promising way for decrease of buildings energy consumption with respect to the sustainable development principles of building industry. The presented paper is focused on the evaluation of the effect of PCM incorporation on thermal performance of cement-lime plasters. For basic characterization of the developed materials, matrix density, bulk density, and total open porosity are measured. Thermal conductivity is accessed by transient impulse method. DSC analysis is used for the identification of phase change temperature during the heating and cooling process. Using DSC data, the temperature dependent specific heat capacity is calculated. On the basis of the experiments performed, the supposed improvement of the energy efficiency of characteristic building envelope system where the designed plasters are likely to be used is evaluated by a computational analysis. Obtained experimental and computational results show a potential of PCM modified plasters for improvement of thermal stability of buildings and moderation of interior climate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ultrafast Molecular Frame Quantum Tomography

Here, we develop and experimentally demonstrate a methodology for a full molecular frame quantum tomography (MFQT) of dynamical polyatomic systems. We exemplify this approach through the complete characterization of an electronically nonadiabatic wave packet in ammonia (NH 3 ). The method exploits both energy and time-domain spectroscopic data, and yields the lab frame density matrix (LFDM) for the system, the elements of which are populations and coherences. The LFDM fully characterizes electronic and nuclear dynamics in the molecular frame, yielding the time- and orientation-angle dependent expectation values of any relevant operator. For example, the time-dependent molecular frame electronic probability density may be constructed, yielding information on electronic dynamics in the molecular frame. In NH 3 , we observe that electronic coherences are induced by nuclear dynamics which nonadiabatically drive electronic motions (charge migration) in the molecular frame. Here, the nuclear dynamics are rotational and it is nonadiabatic Coriolis coupling which drives the coherences. Interestingly, the nuclear-driven electronic coherence is preserved over longer timescales. In general, MFQT can help quantify entanglement between electronic and nuclear degrees of freedom, and provide new routes to the study of ultrafast molecular dynamics, charge migration, quantum information processing, and optimal control schemes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Development of a size reduction equation for woody biomass: The influence of branch wood properties on Rittinger's constant

Size reduction is an essential but energy-intensive process for preparing biomass for conversion processes. Three well-known scaling equations (Bond, Kick, and Rittinger) are used to estimate energy input for grinding minerals and food particles. Previous studies have shown that the Rittinger equation has the best fit to predict energy input for grinding cellulosic biomass. In the Rittinger equation, Rittinger's constant (k R ) is independent of the size of ground particles, yet we noted large variations in k R among similar particle size ranges. In this research, the dependence of k R on the physical structure and chemical composition of a number of woody materials was explored. Branches from two softwood species (Douglas fir and pine) and two hardwood species (aspen and poplar) were ground in a laboratory knife mill. The recorded data included power input, mass flow rate, and particle size before and after grinding. Nine material properties were determined: particle density, solid density (pycnometer and x-ray diffraction methods), microfibril angle, fiber coarseness, fiber length, and composition (lignin and cellulose glucan contents). The correlation matrix among the nine properties revealed high degrees of interdependence between properties. The k R value had the largest positive correlation (+0.60) with particle porosity across the species tested. As a result, particle density was strongly correlated with lignin content (0.85), microfibril angle (0.71), fiber length (0.87), and fiber coarseness (0.78). An empirical model relating k R to particle density was developed.

09 BIOMASS FUELS↗
DOE OSTI2023

Absolute electron density fluctuation reconstruction for two-dimensional hydrogen beam emission spectroscopy

Scrape-off layer (SOL) and edge plasma turbulence significantly contribute to the radial particle and heat transport, lowering the plasma confinement and increasing the heat load on the plasma facing components. SOL turbulence is predominantly intermittent, which manifests in the occurrence of isolated density filaments or blobs. Filaments propagate radially outward toward plasma facing components, limiting their lifetime by erosion and sputtering. To characterize this phenomenon in detail, few diagnostic techniques are available. Beam emission spectroscopy is a diagnostic capable of measuring plasma turbulence in both SOL and edge plasmas. Due to the finite lifetime of the excitation states during the beam–plasma interaction and the misalignment between the optics and the magnetic field, spatial smearing is introduced in the measurement. In this paper, a novel method is introduced to overcome this hindering effect by inverting the fluctuation response matrix on an optimally smoothed signal. We show that this method is fast and provides significantly more accurate absolute density fluctuation reconstruction than the direct inversion technique. Here, the presented method is usable for all types of beam emission diagnostics where the spatial resolution is higher than the combined smearing of the atomic physics and the observation.

47 OTHER INSTRUMENTATION↗
NASA NTRS1997

Low-density resin impregnated ceramic article and method for making the same

A low-density resin impregnated ceramic article advantageously employed as a structural ceramic ablator comprising a matrix of ceramic fibers. The fibers of the ceramic matrix are coated with an organic resin film. The organic resin can be a thermoplastic resin or a cured thermosetting resin. In one embodiment, the resin is uniformly distributed within the ceramic article. In a second embodiment, the resin is distributed so as to provide a density gradient along at least one direction of the ceramic article. The resin impregnated ceramic article is prepared by providing a matrix of ceramic fibers; immersing the matrix of ceramic fibers in a solution of a solvent and an organic resin infiltrant; and removing the solvent to form a resin film on the ceramic fibers.

Tran, Huy K.↗
DOE OSTI2015

Nonadiabatic dynamics with intersystem crossings: A time-dependent density functional theory implementation

In this work, we derive a method to perform trajectory-based nonadiabatic dynamics that is able to describe both nonadiabatic transitions and intersystem crossing events (transitions between states of different spin-multiplicity) at the same level of theory, namely, time-dependent density functional theory (TDDFT). To this end, we combined our previously developed TDDFT-based trajectory surface hopping scheme with an accurate and efficient algorithm for the calculation of the spin-orbit coupling (SOC) matrix elements. More specifically, we designed two algorithms for the calculation of intersystem crossing transitions, one based on an extended Tully’s surface hopping scheme including SOC and the second based on a Landau-Zener approximation applied to the spin sector of the electronic Hilbert space. This development allows for the design of an efficient on-the-fly nonadiabatic approach that can handle, on an equal footing, nonadiabatic and intersystem crossing transitions. The method is applied to the study of the photophysics of sulfur dioxide (SO{sub 2}) in gas and liquid phases.

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

UV-tunable laser induced phototransformations of matrix isolated anethole

A matrix isolation study of the infrared spectra and structure of anethole (1-methoxy-4-(1-propenyl)benzene) has been carried out, showing the presence of two E conformers (AE1, AE2) of the molecule in the as-deposited matrices. Irradiation using ultraviolet-tunable laser light at 308–307 nm induced conformationally selective phototransformations of these forms into two less stable Z conformers (AZ1, AZ2). The back reactions were also detected upon irradiation at 301 nm. On the whole, the obtained results allow for full assignment of the infrared spectra of all the four experimentally observed anethole isomers and showed that the narrowband UV-induced E-Z photoisomerization is an efficient and selective way to interconvert the two isomers of anethole into each other, with conformational discrimination. Photolysis of anethole was observed as well, with initial methoxyl O–C bond cleavage and formation of CH{sub 3} and p-propenylphenoxy (AR) radicals, followed by radical recombination to form 2-methyl-4-propenyl-2,4-cyclohexadienone, which subsequently undergoes ring-opening generating several conformers of long-chain conjugated ketenes. Interpretation of the experimental observations was supported by density functional theory (B3LYP and B2PLYD) calculations.

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

Electronic couplings for molecular charge transfer: Benchmarking CDFT, FODFT, and FODFTB against high-level ab initio calculations

We introduce a database (HAB11) of electronic coupling matrix elements (H{sub ab}) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n-electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute H{sub ab} values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the economic FODFT and highly cost-effective FODFTB methods underestimate couplings by 37.6% and 42.4%, respectively, due to neglect of interaction between donor and acceptor. The errors are systematic, however, and can be significantly reduced by applying a uniform scaling factor for each method. Applications to dimers outside the database, specifically rotated thiophene dimers and larger acenes up to pentacene, suggests that the same scaling procedure significantly improves the FODFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Long-range correction for tight-binding TD-DFT

We present two improvements to the tight-binding approximation of time-dependent density functional theory (TD-DFTB): First, we add an exact Hartree-Fock exchange term, which is switched on at large distances, to the ground state Hamiltonian and similarly to the coupling matrix that enters the linear response equations for the calculation of excited electronic states. We show that the excitation energies of charge transfer states are improved relative to the standard approach without the long-range correction by testing the method on a set of molecules from the database in Peach et al. [J. Chem. Phys. 128, 044118 (2008)] which are known to exhibit problematic charge transfer states. The degree of spatial overlap between occupied and virtual orbitals indicates where TD-DFTB and long-range corrected TD-DFTB (lc-TD-DFTB) can be expected to produce large errors. Second, we improve the calculation of oscillator strengths. The transition dipoles are obtained from Slater Koster files for the dipole matrix elements between valence orbitals. In particular, excitations localized on a single atom, which appear dark when using Mulliken transition charges, acquire a more realistic oscillator strength in this way. These extensions pave the way for using lc-TD-DFTB to describe the electronic structure of large chromophoric polymers, where uncorrected TD-DFTB fails to describe the high degree of conjugation and produces spurious low-lying charge transfer states.

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

Atom Probe Tomography Interlaboratory Study on Clustering Analysis in Experimental Data Using the Maximum Separation Distance Approach

We summarize the findings from an interlaboratory study conducted between ten international research groups and investigate the use of the commonly used maximum separation distance and local concentration thresholding methods for solute clustering quantification. The study objectives are: to bring clarity to the range of applicability of the methods; identify existing and/or needed modifications; and interpretation of past published data. Participants collected experimental data from a proton-irradiated 304 stainless steel and analyzed Cu-rich and Ni–Si rich clusters. The datasets were also analyzed by one researcher to clarify variability originating from different operators. The Cu distribution fulfills the ideal requirements of the maximum separation method (MSM), namely a dilute matrix Cu concentration and concentrated Cu clusters. This enabled a relatively tight distribution of the cluster number density among the participants. By contrast, the group analysis of the Ni–Si rich clusters by the MSM was complicated by a high Ni matrix concentration and by the presence of Si-decorated dislocations, leading to larger variability among researchers. While local concentration filtering could, in principle, tighten the results, the cluster identification step inevitably maintained a high scatter. Recommendations regarding reporting, selection of analysis method, and expected variability when interpreting published data are discussed.

APT, Round Robin Studies, atom probe tomography, c↗
DOE OSTI2015

Analytic nuclear forces and molecular properties from full configuration interaction quantum Monte Carlo

Unbiased stochastic sampling of the one- and two-body reduced density matrices is achieved in full configuration interaction quantum Monte Carlo with the introduction of a second, “replica” ensemble of walkers, whose population evolves in imaginary time independently from the first and which entails only modest additional computational overheads. The matrices obtained from this approach are shown to be representative of full configuration-interaction quality and hence provide a realistic opportunity to achieve high-quality results for a range of properties whose operators do not necessarily commute with the Hamiltonian. A density-matrix formulated quasi-variational energy estimator having been already proposed and investigated, the present work extends the scope of the theory to take in studies of analytic nuclear forces, molecular dipole moments, and polarisabilities, with extensive comparison to exact results where possible. These new results confirm the suitability of the sampling technique and, where sufficiently large basis sets are available, achieve close agreement with experimental values, expanding the scope of the method to new areas of investigation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A Solubility Comparison of Cyclotetramethylenetetranitramine (HMX) and Potential Mock-HMX Candidates in Solvents Relevant to Slurry Formulating Polymer-Bonded Materials

Understanding solubility behavior of organic crystalline materials is critical in the development of slurry-formulated plastic bonded composites. In the field of high explosives (HE), a common method for obtaining high density, greater than 95 % theoretical maximum density, is to encapsulate HE in a polymer matrix resulting in the formation of plastic bonded explosive (PBX) prills. Prill formation is highly nuanced and requires understanding of particle size distributions, solvent, antisolvent, and binder system as well as interactions thereof. In this paper, in an effort to create mock-HE that mimics HE in density and mechanical strength, solubility of candidate molecules was compared to cyclotetramethylenetetranitramine (HMX) in binder solvents relevant to prill formation. The results indicate that the mock candidates chosen for this study have different solubilities compared to HMX in water and both organic solvents, 2-butanone and ethyl acetate, investigated. With these solubility differences known, prills to match PBXs based on HMX can be formulated, pressed and tested.

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

Unbiased reduced density matrices and electronic properties from full configuration interaction quantum Monte Carlo

Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the full configuration interaction quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself) and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Positron accumulation effect in particles embedded in a low-density matrix

Systematic studies of the so-called positron accumulation effect for samples with particles embedded in a matrix are reported. This effect is related to energetic positrons which penetrate inhomogeneous medium. Due to differences in the linear absorption coefficient, different amounts of positrons are accumulated and annihilate in the identical volume of both materials. Positron lifetime spectroscopy and Doppler broadening of the annihilation line using Na-22 positrons were applied to the studies of the epoxy resin samples with embedded micro-sized particles of transition metals, i.e., Ni, Sn, Mo, W, and nonmetal particles, i.e., Si and NaF. The significant difference between the determined fraction of positrons annihilating in the particles and the particle volume fraction indicates the positron accumulation effect. The simple phenomenological model and Monte Carlo simulations are able to describe the main features of the obtained dependencies. The aluminum alloy with embedded Sn nanoparticles is also considered for demonstration differences between the accumulation and another related effect, i.e., the positron affinity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

14 SOLAR ENERGY↗