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 613 records · Page 34

Precision study of the massive Schwinger model near quantum criticality

We perform a numerical analysis of the massive Schwinger model in the presence of a background electric field. Using the Density Matrix Renormalization Group approach, we efficiently compute the spectrum of the Schwinger model on a staggered lattice with up to 3000 qubits. As a result, we achieve a precise computation of the critical mass of the massive Schwinger model to five digits using four different “criticality criteria,” observing perfect agreement among them Additionally, we discuss the effect of a four-fermion operator deformation of the Schwinger model and compute the critical mass for various values of the deformation parameter.

Critical phenomena↗

Beginning the 21st century with advanced Automatic Parts Identification (API)

Under the direction of the NASA George C. Marshall Space Flight Center, Huntsville, Alabama, the development and commercialization of an advanced Automated Parts Indentification (API) system is being undertaken by Rockwell International Corporation. The new API system is based on a variable sized, machine-readable, two-dimensioanl matrix symbol that can be applied directly onto most metallic and nonmetallic materials using safe, permanent marking methods. Its checkerboard-like structure is the most space efficient of all symbologies. This high data-density symbology can be applied to products of different material sizes and geometries using application-dependent, computer-driven marking devices. The high fidelity markings produced by these devices can then be captured using a specially designed camera linked to any IBM-compatible computer. Applications of compressed symbology technology will reduce costs and improve quality, productivity, and processes in a wide variety of federal and commercial applications.

Schramm, Fred↗

The least-squares finite element method for low-mach-number compressible viscous flows

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Yu, Sheng-Tao↗

Optimizing the Pyrolysis of Polypropylene and High-Density Polyethylene for Contaminated Solid Waste Treatment - 18537

Organic waste is generally characterized by high volume-to-weight ratio. Organic waste contaminated with radionuclides may further pose a hazard due to radiolysis and generation of combustible and explosive gases. Thus, minimizing waste volume as well as forming a chemically stable and conditioned matrix are of high importance in the treatment of contaminated polymers. In the present study, the pyrolytic decomposition of the commonly used polyolefins, polypropylene (PP) and high-density polyethylene (HDPE), was explored as a method for waste minimization by obtaining high gaseous and low char and liquid yields. The effects of two aluminosilicate catalysts, Zeolite Y and ZSM-5, and reaction temperature on decomposition efficiency were evaluated based on product conversions and gas composition. Under pyrolytic conditions, PP decomposition resulted in higher decomposition efficiency than HDPE. Over the catalysts, pyrolysis of the polymers over ZSM-5 resulted in higher gas conversions than over Zeolite Y, with a more significant impact on HDPE than PP. Increasing the reaction temperature from 450 deg. C to 525 deg. C also resulted in increased polymer decomposition. However, differences in gas composition, including olefinic to paraffinic ratio and carbon number distribution, support that the catalysts and the temperature impact polymer decomposition by different and opposing mechanisms. This conclusion was further supported by the finding that the decomposition efficiency was not significantly enhanced in the presence of the catalyst at 525 deg. C. Taken together, our study highlights the advantage of using ZSM-5 and optimizing the reaction temperature for treatment of contaminated organic solid waste. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preparation of the SU(3) lattice Yang-Mills vacuum with variational quantum methods

Studying QCD and other gauge theories on quantum hardware requires the preparation of physically interesting states. The variational quantum eigensolver provides a way of performing vacuum state preparation on quantum hardware. Here in this work, variational quantum eigensolver is applied to pure SU(3) lattice Yang-Mills on a single plaquette and one dimensional plaquette chains. Bayesian optimization and gradient descent were investigated for performing the classical optimization. Ansatz states for plaquette chains are constructed in a scalable manner from smaller systems using domain decomposition and a stitching procedure analogous to the density matrix renormalization group. Small examples are performed on IBM’s superconducting Manila processor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Lattice Heat Capacity of Nanostructured Materials Based on Titanium/Zirconium and Aluminum

The dynamic and thermal properties of nanostructured materials based on aluminum with the periodic inclusions of Ti or Zr clusters have been investigated by the method of molecular dynamics. The elastic moduli, spectral vibrational densities of states, and temperature dependences of heat capacity for different Ti–Al and Zr–Al systems have been obtained. The effect of specific features of the phonon spectrum on the heat capacity of the nanocomposite lattice has been examined. It is shown that the ordering type and size of Ti/Zr clusters in the aluminum matrix significantly affect the elastic properties and heat capacity. The results obtained can be used for fabricating new aluminum-, titanium-, and zirconium-based composites with the desired properties.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Plackett-Burman experimental design to facilitate syntactic foam development

The use of an eight-experiment Plackett–Burman method can assess six experimental variables and eight responses in a polysiloxane-glass microsphere syntactic foam. The approach aims to decrease the time required to develop a tunable polymer composite by identifying a reduced set of variables and responses suitable for future predictive modeling. The statistical design assesses the main effects of mixing process parameters, polymer matrix composition, microsphere density and volume loading, and the blending of two grades of microspheres, using a dummy factor in statistical calculations. Responses cover rheological, physical, thermal, and mechanical properties. The cure accelerator content of the polymer matrix and the volume loading of the microspheres have the largest effects on foam properties. These factors are the most suitable for controlling the gel point of the curing foam, and the density of the cured foam. The mixing parameters introduce widespread variability and therefore should be fixed at effective levels during follow-up testing. Some responses may require greater contrast in microsphere-related factors. As a result, compared to other possible statistical approaches, the run economy of the Plackett–Burman method makes it a valuable tool for rapidly characterizing new foams.

36 MATERIALS SCIENCE↗

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermoelectric performance of nanostructured p-type Zr₀.₅Hf₀.₅Co₀.₄Rh₀.₆Sb 1–x Sn x half-Heusler alloys

Several compositions of the p-type half-Heusler alloys Zr₀.₅Hf₀.₅Co₀.₄Rh₀.₆Sb 1–x Sn x (0≤x≤0.4) were synthesized by mechanically alloying high purity elemental powders using hardened steel jars and balls on a high energy shaker mill. Powder X-ray diffraction (XRD) investigations of several aliquots taken after regularly spaced milling time suggested that single phase products with half-Heusler (HH) structure can be obtained after 10 h. However, XRD and transmission electron microscopy (TEM) studies of several specimens obtained from compacted polycrystalline powders of Z₀.₅Hf₀.₅Co₀.₄Rh₀.₆Sb 1–x Sn x alloys using a uniaxial hot press (HP) revealed the presence of CoSb inclusions with various sizes embedded inside the HH matrix. Hall effect, electrical conductivity, and thermopower data collected between 300 K and 775 K on several compositions suggested that electronic transport in the synthesized Zr₀.₅Hf₀.₅Co₀.₄Rh₀.₆Sb 1–x Sn x /CoSb composites strongly depends on the average size and/or mole fraction of the embedded CoSb inclusions rather than the fraction (x) of Sn substituting for Sb. Among the samples investigated, the nanocomposite with x=0.2, which contains nanometer-scale CoSb inclusions, showed the largest power factor (800 μW/K² m at 775 K) and the lowest lattice thermal conductivity (~2.2 W/m K at 775 K) leading to a six-fold enhancement in the figure of merit when compared to the Zr₀.₅Hf₀.₅Co₀.₄Rh₀.₆Sb₀.₉₉Sn₀.₀₁ bulk matrix. - Graphical abstract: CoSb nanoinclusions embedded into a p-type Zr₀.₅Hf₀.₅Co₀.₄Rh₀.₆Sb 1–x Sn x half-Heusler matrix simultaneously boost the thermopower and carrier mobility leading to a drastic enhancement of the power factor of the resulting bulk nanostructured materials. Highlights: • The phase composition of half-Heusler (HH) alloy is very sensitive to the synthesis method. • Mechanical alloying (MA) of elements results in bulk HH matrix with CoSb inclusions. • CoSb nanoinclusions simultaneously boost the thermopower and carrier mobility with the HH matrix. • The MA process reduces lattice thermal conductivity due to high density of grain boundaries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exciton dissociation at donor-acceptor heterojunctions: Dynamics using the collective effective mode representation of the spin-boson model

Following the recent quantum dynamics investigation of the charge transfer at an oligothiophene-fullerene heterojunction by the multi-configuration time dependent Hartree method [H. Tamura, R. Martinazzo, M. Ruckenbauer and I. Burghardt, J. Chem. Phys. 137, 22A540 (2012)], we revisit the transfer process by a perturbative non-Markovian master equation treated by the time local auxiliary density matrix approach. We compare the efficiency of the spin-boson model calibrated by quantum chemistry with the effective mode representation. A collective mode is extracted from the spin-boson spectral density. It is weakly coupled to a residual bath of vibrational modes, allowing second-order dynamics. The electron transfer is analyzed for a sampling of inter-fragment distances showing the fine interplay of the electronic coupling and energy gap on the relaxation. The electronic coherence, expected to play a role in the process, is preserved during about 200 fs.

74 ATOMIC AND MOLECULAR PHYSICS↗

Transport properties in dilute UN ( X ) solid solutions ( X = Xe , Kr )

Uranium nitride (UN) is a candidate fuel for current GEN III fission reactors, for which it is investigated as an accident-tolerant fuel, as well as for future GEN IV reactors. In this study, we investigate the kinetic properties of gas fission products (Xe and Kr) in UN. Binding and migration energies are obtained using density functional theory, with an added Hubbard correlation to model f electrons, and the occupation matrix control scheme to avoid metastable states. These energies are then used as input for the self-consistent mean field method which enables to determine transport coefficients for vacancy-mediated diffusion of Xe and Kr on the U sublattice. The magnetic ordering of the UN structure is explicitly taken into account, for both energetic and transport properties. Solute diffusivities are compared with experimental measurements and the effect of various parameters on the theoretical model is carefully investigated. We find that kinetic correlations are very strong in this system, and that despite atomic migration anisotropy, macroscopic solute diffusivities show limited anisotropy. Our model indicates that the discrepancy between experimental measurements probably results from different irradiation conditions, and hence different defect concentrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Identifying Entangled Physics Relationships Through Sparse Matrix Decomposition to Inform Plasma Fusion Design

We report a sustainable burn platform through inertial confinement fusion (ICF) has been an ongoing challenge for over 50 years. Mitigating engineering limitations and improving the current design involves an understanding of the complex coupling of physical processes. While sophisticated simulation codes are used to model ICF implosions, these tools contain necessary numerical approximation but miss physical processes that limit predictive capability. Identification of relationships between controllable design inputs to ICF experiments and measurable outcomes (e.g., neutron yield, neutron velocity, areal density) from performed experiments can help guide the future design of experiments and development of simulation codes, to potentially improve the accuracy of the computational models used to simulate ICF experiments. We use sparse matrix decomposition methods to identify clusters of a few related design variables. Sparse principal component analysis (SPCA) identifies groupings that are related to the physical origin of the variables (laser, hohlraum, and capsule). A variable importance analysis finds that in addition to variables highly correlated with neutron yield, such as picket power and laser energy, variables that represent a dramatic change of the ICF design, such as number of pulse steps, are also very important. The obtained sparse components are then used to train a random forest (RF) regression surrogate for predicting total yield. The RF performance on the training and testing data compares with the performance of the RF trained using all the design variables considered. This work is intended to inform design changes in future ICF experiments by augmenting the expert intuition and simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optical characterization of free electron concentration in heteroepitaxial InN layers using Fourier transform infrared spectroscopy and a 2 Multiplication-Sign 2 transfer-matrix algebra

Fourier Transform Infrared (FTIR) reflectance spectroscopy has been implemented as a non-destructive, non-invasive, tool for the optical characterization of a set of c-plane InN single heteroepitaxial layers spanning a wide range of thicknesses (30-2000 nm). The c-plane (0001) InN epilayers were grown by plasma-assisted molecular beam epitaxy (PAMBE) on GaN(0001) buffer layers which had been grown on Al{sub 2}O{sub 3}(0001) substrates. It is shown that for arbitrary multilayers with homogeneous anisotropic layers having their principal axes coincident with the laboratory coordinates, a 2 Multiplication-Sign 2 matrix algebra based on a general transfer-matrix method (GTMM) is adequate to interpret their optical response. Analysis of optical reflectance in the far and mid infrared spectral range has been found capable to discriminate between the bulk, the surface and interface contributions of free carriers in the InN epilayers revealing the existence of electron accumulation layers with carrier concentrations in mid 10{sup 19} cm{sup -3} at both the InN surface and the InN/GaN interface. The spectra could be fitted with a three-layer model, determining the different electron concentration and mobility values of the bulk and of the surface and the interface electron accumulation layers in the InN films. The variation of these values with increasing InN thickness could be also sensitively detected by the optical measurements. The comparison between the optically determined drift mobility and the Hall mobility of the thickest sample reveals a value of r{sub H} = 1.49 for the Hall factor of InN at a carrier concentration of 1.11 Multiplication-Sign 10{sup 19} cm{sup -3} at 300 Degree-Sign {Kappa}.

36 MATERIALS SCIENCE↗

Mesh refinement in finite element analysis by minimization of the stiffness matrix trace

Most finite element packages provide means to generate meshes automatically. However, the user is usually confronted with the problem of not knowing whether the mesh generated is appropriate for the problem at hand. Since the accuracy of the finite element results is mesh dependent, mesh selection forms a very important step in the analysis. Indeed, in accurate analyses, meshes need to be refined or rezoned until the solution converges to a value so that the error is below a predetermined tolerance. A-posteriori methods use error indicators, developed by using the theory of interpolation and approximation theory, for mesh refinements. Some use other criterions, such as strain energy density variation and stress contours for example, to obtain near optimal meshes. Although these methods are adaptive, they are expensive. Alternatively, a priori methods, until now available, use geometrical parameters, for example, element aspect ratio. Therefore, they are not adaptive by nature. An adaptive a-priori method is developed. The criterion is that the minimization of the trace of the stiffness matrix with respect to the nodal coordinates, leads to a minimization of the potential energy, and as a consequence provide a good starting mesh. In a few examples the method is shown to provide the optimal mesh. The method is also shown to be relatively simple and amenable to development of computer algorithms. When the procedure is used in conjunction with a-posteriori methods of grid refinement, it is shown that fewer refinement iterations and fewer degrees of freedom are required for convergence as opposed to when the procedure is not used. The mesh obtained is shown to have uniform distribution of stiffness among the nodes and elements which, as a consequence, leads to uniform error distribution. Thus the mesh obtained meets the optimality criterion of uniform error distribution.

Kittur, Madan G.↗

The all configuration mean energy multiconfiguration self-consistent-field method. I. Equal configuration weights

The All Configuration Mean Energy (ACME) conditions are a special case of state averaging for Multiconfigurational Self-Consistent-Field (MCSCF) orbital optimisation. The method is formulated using the Graphical Unitary Group Approach (GUGA) in which the Configuration State Function (CSF) basis is represented as walks within a Shavitt graph. This graphical formulation leads to efficient recursive algorithms for the energy and reduced density matrices (RDM) that are independent of the CSF dimension and that scale only as O(n2) where n is the number of occupied orbitals. The Hamiltonian matrix diagonalization step is obviated and the CSF expansion coefficients are neither referenced nor required. This allows MCSCF orbital optimisation to be performed for essentially unlimited numbers of active orbitals and arbitrarily large CSF expansions. The discussion includes various types of CSF expansion spaces, the partitioning of the essential and redundant orbital optimisation parameters, the computation of the spin-density, and the formulation of state-specific analytic gradients and nonadiabatic coupling for high-level electronic structure methods that use the ACME MCSCF orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimating the entropy and quantifying the impurity of a swarm of surface-hopping trajectories: A new perspective on decoherence

In this article, we consider the intrinsic entropy of Tully's fewest switches surface hopping (FSSH) algorithm (as estimated by the impurity of the density matrix) [J. Chem. Phys. 93, 1061 (1990)]. We show that, even for a closed system, the total impurity of a FSSH calculation increases in time (rather than stays constant). This apparent failure of the FSSH algorithm can be traced back to an incorrect, approximate treatment of the electronic coherence between wavepackets moving along different potential energy surfaces. This incorrect treatment of electronic coherence also prevents the FSSH algorithm from correctly describing wavepacket recoherences (which is a well established limitation of the FSSH method). Nevertheless, despite these limitations, the FSSH algorithm often predicts accurate observables because the electronic coherence density is modulated by a phase factor which varies rapidly in phase space and which often integrates to almost zero. Adding “decoherence” events on top of a FSSH calculation completely destroys the incorrect FSSH electronic coherence and effectively sets the Poincaré recurrence time for wavepacket recoherence to infinity; this modification usually increases FSSH accuracy (assuming there are no recoherences) while also offering long-time stability for trajectories. In practice, we show that introducing “decoherence” events does not change the total FSSH impurity significantly, but does lead to more accurate evaluations of the impurity of the electronic subsystem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stimulated X-ray Raman and Absorption Spectroscopy of Iron–Sulfur Dimers

Iron–sulfur complexes play an important role in biological processes such as metabolic electron transport. A detailed understanding of the mechanism of long-range electron transfer requires knowledge of the electronic structure of the complexes, which has traditionally been challenging to obtain, either by theory or by experiment, but the situation has begun to change with advances in quantum chemical methods and intense free electron laser light sources. We compute the spectra for stimulated X-ray Raman spectroscopy (SXRS) and absorption spectroscopy of homovalent and mixed-valence [2Fe–2S] complexes, using the ab initio density matrix renormalization group algorithm. The simulated spectra show clear signatures of the theoretically predicted dense low-lying excited states within the d–d manifold. Finally, the difference in spectral intensity between the absorption-active and Raman-active states provides a potential mechanism to selectively excite states by a proper tuning of the excitation pump, to access the electronic dynamics within this manifold.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗