Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrices”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Golden rule kinetics of transfer reactions in condensed phase: The microscopic model of electron transfer reactions in disordered solid matrices

The algorithm for a theoretical calculation of transfer reaction rates for light quantum particles (i.e., the electron and H-atom transfers) in non-polar solid matrices is formulated and justified. The mechanism postulated involves a local mode (an either intra- or inter-molecular one) serving as a mediator which accomplishes the energy exchange between the reacting high-frequency quantum mode and the phonon modes belonging to the environment. This approach uses as a background the Fermi golden rule beyond the usually applied spin-boson approximation. The dynamical treatment rests on the one-dimensional version of the standard quantum relaxation equation for the reduced density matrix, which describes the frequency fluctuation spectrum for the local mode under consideration. The temperature dependence of a reaction rate is controlled by the dimensionless parameter ξ{sub 0}=ℏω{sub 0}/k{sub B}T where ω{sub 0} is the frequency of the local mode and T is the temperature. The realization of the computational scheme is different for the high/intermediate (ξ{sub 0} < 1 − 3) and for low (ξ{sub 0}≫ 1) temperature ranges. For the first (quasi-classical) kinetic regime, the Redfield approximation to the solution of the relaxation equation proved to be sufficient and efficient in practical applications. The study of the essentially quantum-mechanical low-temperature kinetic regime in its asymptotic limit requires the implementation of the exact relaxation equation. The coherent mechanism providing a non-vanishing reaction rate has been revealed when T→ 0. An accurate computational methodology for the cross-over kinetic regime needs a further elaboration. The original model of the hopping mechanism for electronic conduction in photosensitive organic materials is considered, based on the above techniques. The electron transfer (ET) in active centers of such systems proceeds via local intra- and intermolecular modes. The active modes, as a rule, operate beyond the kinetic regimes, which are usually postulated in the existing theories of the ET. Our alternative dynamic ET model for local modes immersed in the continuum harmonic medium is formulated for both classical and quantum regimes, and accounts explicitly for the mode/medium interaction. The kinetics of the energy exchange between the local ET subsystem and the surrounding environment essentially determine the total ET rate. The efficient computer code for rate computations is elaborated on. The computations are available for a wide range of system parameters, such as the temperature, external field, local mode frequency, and characteristics of mode/medium interaction. The relation of the present approach to the Marcus ET theory and to the quantum-statistical reaction rate theory [V. G. Levich and R. R. Dogonadze, Dokl. Akad. Nauk SSSR, Ser. Fiz. Khim. 124, 213 (1959); J. Ulstrup, Charge Transfer in Condensed Media (Springer, Berlin, 1979); M. Bixon and J. Jortner, Adv. Chem. Phys. 106, 35 (1999)] underlying it is discussed and illustrated by the results of computations for practically important target systems.

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

Formation of ionic complexes in cryogenic matrices: A case study using co-deposition of Cu{sup −} with rare gas cations in solid argon

Matrix isolation spectra have been obtained for ionic species formed from a beam of mass-selected ions, with a coincident beam of externally generated counter-ions used to provide charge balance. Infrared spectra were obtained for copper carbonyl complexes formed following deposition of Cu{sup −} ions with rare-gas counter-cations into CO-doped argon matrices. Both anionic and neutral copper carbonyl complexes Cu(CO){sub n}{sup q} (n = 1–3; q = 0, −1) were observed in the spectra, with peak positions corresponding to previously reported assignments; new partially resolved bands appearing in the range 1830–1845 cm{sup −1} are assigned to larger [Cu(CO){sub 3}•(CO){sub n}]{sup −} aggregates, having additional CO ligands in the second solvation shell. The experimental geometry ensures that all Cu-centers initially arrive at the matrix as anions, so the relative abundance of anionic relative to neutral complexes is much higher than in previous studies employing alternative methods for ion deposition; this allows for monitoring of electron-transfer processes between anions and cations in the matrix. Comparison of time-dependent vs. temperature-dependent trends reveals that there are two distinct mechanisms by which the population of anionic complexes is converted into neutral complexes: short-range electron transfer between a cation-anion pair following diffusion, and long-range electron transfer involving photodetachment of an electron from the anion into the conduction band of solid argon, resulting in eventual recombination of the electron with a cation in a remote matrix site. The spectra also show a marked dependence on the deposition temperature and dopant concentration, in that 100-fold higher CO concentrations were required during deposition with the sample window at 10 K compared to that used at 20 K, in order to obtain a similar distribution of copper carbonyl complexes. Furthermore, although no carbonyl complexes are observed initially when low concentrations of CO are used at 10 K, upon warming the matrix to 15 K, the neutral di- and tricarbonyl peaks appear abruptly, which is attributed to fast diffusion of CO stimulated by the energy released upon short-range electron-transfer between Cu{sup −}:counter-cation pairs.

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

The rigorous stochastic matrix multiplication scheme for the calculations of reduced equilibrium density matrices of open multilevel quantum systems

Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, i.e., how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernon influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Photochemistry of glycolaldehyde in cryogenic matrices

The photochemistry of glycolaldehyde (GA) upon irradiation at 266 nm is investigated in argon, nitrogen, neon, and para-hydrogen matrices by IR spectroscopy. Isomerization and fragmentation processes are found to compete. The hydrogen-bonded Cis-Cis form of GA is transformed mainly to the open Trans-Trans conformer and to CO and CH{sub 3}OH fragments and their mixed complexes. Different photo-induced behaviours appear depending on the matrix. In nitrogen, small amounts of Trans-Gauche and Trans-Trans conformers are detected after deposition and grow together upon irradiation. The Trans-Gauche conformer is characterized for the first time. In para-hydrogen due to a weaker cage effect additional H{sub 2}CO and HCO fragments are seen. Calculations of the potential energy surfaces of S{sub 0}, S{sub 1}, and T{sub 1} states – to analyse the torsional deformations which are involved in the isomerization process – and a kinetic analysis are presented to investigate the different relaxation pathways of GA. Fragmentation of GA under UV irradiation through the CO+CH{sub 3}OH molecular channel is a minor process, as in the gas phase.

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

Communication: Satisfying fermionic statistics in the modeling of open time-dependent quantum systems with one-electron reduced density matrices

For an open, time-dependent quantum system, Lindblad derived the most general modification of the quantum Liouville equation in the Markovian approximation that models environmental effects while preserving the non-negativity of the system’s density matrix. While Lindblad’s modification is correct for N-electron density matrices, solution of the Liouville equation with a Lindblad operator causes the one-electron reduced density matrix (1-RDM) to violate the Pauli exclusion principle. Consequently, after a short time, the 1-RDM is not representable by an ensemble N-electron density matrix (not ensemble N-representable). In this communication, we derive the necessary and sufficient constraints on the Lindbladian matrix within the Lindblad operator to ensure that the 1-RDM remains N-representable for all time. The theory is illustrated by considering the relaxation of an excitation in several molecules F{sub 2}, N{sub 2}, CO, and BeH{sub 2} subject to environmental noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

In situ formation and characterisation of singly ionised atomic europium in rare gas matrices—Luminescence spectroscopy and MP2 calculations

Irradiation of atomic europium isolated in the solid rare gases, with low intensity laser excitation of the y{sup 8}P←a{sup 8}S resonance transition at ca. 465 nm, is found to produce singly charged europium cations (Eu{sup +}) in large amounts in xenon and in smaller amounts in argon. Confirmation of the formation of matrix-isolated Eu{sup +} is obtained from characteristic absorption bands in the UV and in the visible spectral regions. The luminescence produced with excitation of the cation bands is presented in greatest detail for Eu/Xe and assigned. Excitation of the 4f{sup 7}({sup 8}S{sub 7/2})6p{sub 3/2} absorption bands of Eu{sup +} between 390 and 410 nm produces emission which is quite distinct from that resulting from excitation of the 4f{sup 7}({sup 8}S{sub 7/2})6p{sub 1/2} absorption (430 to 450 nm) features. The latter consists of narrow, resolved emission bands with Stokes shifts ten times smaller than the former. The observed spectral differences are discussed in relation to the different spatial symmetries of the p{sub 3/2} and p{sub 1/2} orbitals in these j-j coupled (7/2, 3/2){sub J} and the (7/2, 1/2){sub J} levels. Møller-Plesset calculations are conducted to obtain the molecular parameters of the neutral Eu-RG and cationic Eu{sup +}-RG diatomics (RG = Ar, Kr, Xe). From the short bond lengths and the strong binding energies obtained for the Eu{sup +}-RG species, these values suggest the isolation of the ion in small, possibly interstitial sites especially in xenon. In contrast, but consistent with previous work [O. Byrne and J. G. McCaffrey, J. Chem. Phys. 134, 124501 (2011)], the interaction potentials calculated herein for the Eu-RG diatomics suggest that the neutral Eu atom occupies tetra-vacancy (tv) and hexa-vacancy (hv) sites in the solid rare gas hosts. Possible reasons for the facile production of Eu{sup +} in the solid rare gases are discussed. The mechanism proposed is that atomic europium is also acting as an electron acceptor, providing a temporary trap for the ionised electron in the matrices.

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

Electronic states in spherical GaN nanocrystals embedded in various dielectric matrices: The k ⋅ p-calculations

Using the envelope-function approximation, the single-particle states of electrons and holes in spherical GaN nanocrystals embedded in different amorphous dielectric matrices (SiO{sub 2}, Al{sub 2}O{sub 3}, HfO{sub 2} and Si{sub 3}N{sub 4}) have been calculated. Ground state energies of electrons and holes in GaN nanocrystals are determined using the isotropic approximation of the k ⋅ p -Hamiltonian. All the ground state energies are found to increase with lowering the nanocrystal size and are proportional to the R{sup −n}, where R is the nanocrystal radius, n =1.5-1.9 for electrons and 1.7-2.0 for holes. The optical gap of GaN nanocrystals changes from 3.8 to 5 eV for the nanocrystal radius ranging from 3 to 1 nm.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Infrared spectroscopy of water clusters isolated in methane matrices: Effects of isotope substitution and annealing

Using infrared-active solvents of CH{sub 4} and CD{sub 4} for matrix isolation, we measured infrared spectra of H{sub 2}O and D{sub 2}O clusters at 7 K. The solute-concentration dependence of the spectrum of H{sub 2}O clusters in a CH{sub 4} matrix was investigated and was used for the peak assignment. Annealing procedures were found to promote the size growth of water clusters in methane matrices for all the combinations of (H{sub 2}O, CH{sub 4}), (H{sub 2}O, CD{sub 4}), (D{sub 2}O, CH{sub 4}), and (D{sub 2}O, CD{sub 4}). We also monitored the ν{sub 3} absorption due to methane to find the annealing-induced structural change only of solid CH{sub 4}. The matrix effects on the vibrations of the clusters are discussed on the basis of “T{sub c} plots”, where their frequencies are plotted as a function of the square root of the matrix critical temperature, T{sub c}. The obtained plots assure the validity of the assignment of the cluster peaks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Permutation invariant polynomial neural network approach to fitting potential energy surfaces. IV. Coupled diabatic potential energy matrices

A machine learning method is proposed for representing the elements of diabatic potential energy matrices (PEMs) with high fidelity. This is an extension of the so-called permutation invariant polynomial-neural network (PIP-NN) method for representing adiabatic potential energy surfaces. While for one-dimensional irreducible representations the diagonal elements of a diabatic PEM are invariant under exchange of identical nuclei in a molecular system, the off-diagonal elements require special symmetry consideration, particularly in the presence of a conical intersection. A multiplicative factor is introduced to take into consideration the particular symmetry properties while maintaining the PIP-NN framework. Here, we demonstrate here that the extended PIP-NN approach is accurate in representing diabatic PEMs, as evidenced by small fitting errors and by the reproduction of absorption spectra and product branching ratios in both H 2 O$(\tilde{X} / \tilde{B})$ and NH 3 $(\tilde{X} $/Ã) non-adiabatic photodissociation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Interpretation of multiple solutions in fully iterative GF2 and GW schemes using local analysis of two-particle density matrices

Due to the presence of non-linear equations, iterative Green’s function methods can result in multiple different solutions even for simple molecular systems. In contrast to the wave-function methods, a detailed and careful analysis of such molecular solutions was not performed before. In this work, we use two-particle density matrices to investigate local spin and charge correlators that quantify the charge resonance and covalent characters of these solutions. When applied within the unrestricted orbital set, spin correlators elucidate the broken symmetry of the solutions, containing necessary information for building effective magnetic Hamiltonians. Based on GW and GF2 calculations of simple molecules and transition metal complexes, we construct Heisenberg Hamiltonians, four-spin-four-center corrections, and biquadratic spin–spin interactions. These Hamiltonian parameterizations are compared to previous wave-function calculations.

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

An accurate and efficient fragmentation approach via the generalized many-body expansion for density matrices

With relevant chemical space growing larger and larger by the day, the ability to extend computational tractability over that larger space is of paramount importance in virtually all fields of science. The solution we aim to provide here for this issue is in the form of the generalized many-body expansion for building density matrices (GMBE-DM) based on the set-theoretical derivation with overlapping fragments, through which the energy can be obtained by a single Fock build. In combination with the purification scheme and the truncation at the one-body level, the DM-based GMBE(1)-DM-P approach shows both highly accurate absolute and relative energies for medium-to-large size water clusters with about an order of magnitude better than the corresponding energy-based GMBE(1) scheme. Simultaneously, GMBE(1)-DM-P is about an order of magnitude faster than the previously proposed MBE-DM scheme [F. Ballesteros and K. U. Lao, J. Chem. Theory Comput. 18, 179 (2022)] and is even faster than a supersystem calculation without significant parallelization to rescue the fragmentation method. For even more challenging systems including ion–water and ion–pair clusters, GMBE(1)-DM-P also performs about 3 and 30 times better than the energy-based GMBE(1) approach, respectively. In addition, this work provides the first overlapping fragmentation algorithm with a robust and effective binning scheme implemented internally in a popular quantum chemistry software package. Thus, GMBE(1)-DM-P opens a new door to accurately and efficiently describe noncovalent clusters using quantum mechanics.

Chemistry↗
DOE OSTI2014

Lasing characteristics of difluoroborates of 2,2'-dipyrromethene derivatives in solid matrices

The luminescence-spectral, lasing and photochemical characteristics of laser media based on boron fluoride complexes of dipyrromethenes, embedded into solid bulk matrices of polymethylmethacrylate and its modifications (obtained by adding polyhedral oligomeric silsesquioxane during polymerisation) and into polymer films, in which polyhedral silsesquioxane enters the composition of monomeric unit, have been investigated. (lasers)

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

Optical fields in porous polylactide matrices

Optical fields induced by laser light inside porous polymer matrices (scaffolds) under stimulating effect on cells are investigated. The spectra of light scattering coefficients and anisotropy factors of porous polylactide scaffolds in the visible and near-IR ranges are determined by numerical simulation based on experimental data obtained using a double integrating sphere. Studies are conducted for dry and wet scaffolds fabricated by supercritical fluid foaming and surface-selective laser sintering. The Monte Carlo method is used to calculate the fluence rate distribution under the assumptions of the Henyey – Greenstein phase scattering function and the model phase function of scattering on spherical cavities. It is shown that th e fluence rate distribution mainly depends on the scattering coefficient and the mean cosine of the deflection angle (g-factor). The results obtained are of great practical importance for tissue engineering based on the use of polymeric porous scaffolds with cells to replace tissue defects; they can be used to determine the dose range for laser stimulation of cells. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2010

INFRARED SPECTROSCOPY OF NAPHTHALENE AGGREGATION AND CLUSTER FORMATION IN ARGON MATRICES

Fourier-transform, mid-infrared absorption spectra of mixed argon/naphthalene matrices at 5 K are shown with ratios of argon-to-naphthalene that vary from 1000 to 0. These spectra show the changes as naphthalene clustering and aggregation occurs, with moderate spectral shifts affecting the C-H vibrational modes and relatively small or no shifts to the C-C and C-C-C vibrational modes. The possible contribution of homogeneous naphthalene clusters to the interstellar unidentified infrared bands is discussed. The contribution of polycyclic aromatic hydrocarbon (PAH) clusters to the 7.7 {mu}m emission plateau and the blue shading of the 12.7 {mu}m emission band are identified as promising candidates for future research. In addition, since PAH clusters are model components of Jupiter and Titan's atmospheres, the information presented here may also be applicable to the spectroscopy of these objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

2D-FFTLog: efficient computation of real-space covariance matrices for galaxy clustering and weak lensing

ABSTRACT Accurate covariance matrices for two-point functions are critical for inferring cosmological parameters in likelihood analyses of large-scale structure surveys. Among various approaches to obtaining the covariance, analytic computation is much faster and less noisy than estimation from data or simulations. However, the transform of covariances from Fourier space to real space involves integrals with two Bessel integrals, which are numerically slow and easily affected by numerical uncertainties. Inaccurate covariances may lead to significant errors in the inference of the cosmological parameters. In this paper, we introduce a 2D-FFTLog algorithm for efficient, accurate, and numerically stable computation of non-Gaussian real-space covariances for both 3D and projected statistics. The 2D-FFTLog algorithm is easily extended to perform real-space bin-averaging. We apply the algorithm to the covariances for galaxy clustering and weak lensing for a Dark Energy Survey Year 3-like and a Rubin Observatory’s Legacy Survey of Space and Time Year 1-like survey, and demonstrate that for both surveys, our algorithm can produce numerically stable angular bin-averaged covariances with the flat sky approximation, which are sufficiently accurate for inferring cosmological parameters. The code CosmoCov for computing the real-space covariances with or without the flat-sky approximation is released along with this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Bias on tensor-to-scalar ratio inference with estimated covariance matrices

ABSTRACT We investigate simulation-based bandpower covariance matrices commonly used in cosmological parameter inferences such as the estimation of the tensor-to-scalar ratio r. We find that upper limits on r can be biased low by tens of per cent. The underestimation of the upper limit is most severe when the number of simulation realizations is similar to the number of observables. Convergence of the covariance-matrix estimation can require a number of simulations an order of magnitude larger than the number of observables, which could mean $\mathcal {O}(10\ 000)$ simulations. This is found to be caused by an additional scatter in the posterior probability of r due to Monte Carlo noise in the estimated bandpower covariance matrix, in particular, by spurious non-zero off-diagonal elements. We show that matrix conditioning can be a viable mitigation strategy in the case that legitimate covariance assumptions can be made.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Estimating covariance matrices for two- and three-point correlation function moments in Arbitrary Survey Geometries

ABSTRACT We present configuration-space estimators for the auto- and cross-covariance of two- and three-point correlation functions (2PCF and 3PCF) in general survey geometries. These are derived in the Gaussian limit (setting higher order correlation functions to zero), but for arbitrary non-linear 2PCFs (which may be estimated from the survey itself), with a shot-noise rescaling parameter included to capture non-Gaussianity. We generalize previous approaches to include Legendre moments via a geometry-correction function calibrated from measured pair and triple counts. Making use of importance sampling and random particle catalogues, we can estimate model covariances in fractions of the time required to do so with mocks, obtaining estimates with negligible sampling noise in ∼10 (∼100) CPU-hours for the 2PCF (3PCF) autocovariance. We compare results to sample covariances from a suite of BOSS DR12 mocks and find the matrices to be in good agreement, assuming a shot-noise rescaling parameter of 1.03 (1.20) for the 2PCF (3PCF). To obtain strongest constraints on cosmological parameters, we must use multiple statistics in concert; having robust methods to measure their covariances at low computational cost is thus of great relevance to upcoming surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Analytic Gaussian covariance matrices for galaxy N-point correlation functions

Here, we derive analytic covariance matrices for the N-point correlation functions (NPCFs) of galaxies in the Gaussian limit. Our results are given for arbitrary N and projected onto the isotropic basis functions given by spherical harmonics and Wigner 3j symbols. A numerical implementation of the 4PCF covariance is compared to the sample covariance obtained from a set of lognormal simulations, Quijote dark matter halo catalogues, and MultiDark-Patchy galaxy mocks, with the latter including realistic survey geometry. The analytic formalism gives reasonable predictions for the covariances estimated from mock simulations with a periodic-box geometry. Furthermore, fitting for an effective volume and number density by maximizing a likelihood based on Kullback-Leibler divergence is shown to partially compensate for the effects of a nonuniform window function. Our result is recently shown to facilitate NPCF analysis on a realistic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗