Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative 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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

On the validity of the basis set superposition error and complete basis set limit extrapolations for the binding energy of the formic acid dimer

We report the variation of the binding energy of the Formic Acid Dimer with the size of the basis set at the Coupled Cluster with iterative Singles, Doubles and perturbatively connected Triple replacements [CCSD(T)] level of theory, estimate the Complete Basis Set (CBS) limit, and examine the validity of the Basis Set Superposition Error (BSSE)-correction for this quantity that was previously challenged by Kalescky, Kraka, and Cremer (KKC) [J. Chem. Phys. 140, 084315 (2014)]. Our results indicate that the BSSE correction, including terms that account for the substantial geometry change of the monomers due to the formation of two strong hydrogen bonds in the dimer, is indeed valid for obtaining accurate estimates for the binding energy of this system as it exhibits the expected decrease with increasing basis set size. We attribute the discrepancy between our current results and those of KKC to their use of a valence basis set in conjunction with the correlation of all electrons (i.e., including the 1s of C and O). We further show that the use of a core-valence set in conjunction with all electron correlation converges faster to the CBS limit as the BSSE correction is less than half than the valence electron/valence basis set case. The uncorrected and BSSE-corrected binding energies were found to produce the same (within 0.1 kcal/mol) CBS limits. We obtain CCSD(T)/CBS best estimates for D{sub e} = − 16.1 ± 0.1 kcal/mol and for D{sub 0} = − 14.3 ± 0.1 kcal/mol, the later in excellent agreement with the experimental value of −14.22 ± 0.12 kcal/mol.

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

Crustal structure beneath two seismic stations in the Sunda-Banda arc transition zone derived from receiver function analysis

We analyzed receiver functions to estimate the crustal thickness and velocity structure beneath two stations of Geofon (GE) network in the Sunda-Banda arc transition zone. The stations are located in two different tectonic regimes: Sumbawa Island (station PLAI) and Timor Island (station SOEI) representing the oceanic and continental characters, respectively. We analyzed teleseismic events of 80 earthquakes to calculate the receiver functions using the time-domain iterative deconvolution technique. We employed 2D grid search (H-κ) algorithm based on the Moho interaction phases to estimate crustal thickness and Vp/Vs ratio. We also derived the S-wave velocity variation with depth beneath both stations by inverting the receiver functions. We obtained that beneath station PLAI the crustal thickness is about 27.8 km with Vp/Vs ratio 2.01. As station SOEI is covered by very thick low-velocity sediment causing unstable solution for the inversion, we modified the initial velocity model by adding the sediment thickness estimated using high frequency content of receiver functions in H-κ stacking process. We obtained the crustal thickness is about 37 km with VP/Vs ratio 2.2 beneath station SOEI. We suggest that the high Vp/Vs in station PLAI may indicate the presence of fluid ascending from the subducted plate to the volcanic arc, whereas the high Vp/Vs in station SOEI could be due to the presence of sediment and rich mafic composition in the upper crust and possibly related to the serpentinization process in the lower crust. We also suggest that the difference in velocity models and crustal thicknesses between stations PLAI and SOEI are consistent with their contrasting tectonic environments.

58 GEOSCIENCES↗
DOE OSTI2015

A reanalysis of a strong-flow gyrokinetic formalism

We reanalyse an arbitrary-wavelength gyrokinetic formalism [A. M. Dimits, Phys. Plasmas 17, 055901 (2010)], which orders only the vorticity to be small and allows strong, time-varying flows on medium and long wavelengths. We obtain a simpler gyrocentre Lagrangian up to second order. In addition, the gyrokinetic Poisson equation, derived either via variation of the system Lagrangian or explicit density calculation, is consistent with that of the weak-flow gyrokinetic formalism [T. S. Hahm, Phys. Fluids 31, 2670 (1988)] at all wavelengths in the weak flow limit. The reanalysed formalism has been numerically implemented as a particle-in-cell code. An iterative scheme is described which allows for numerical solution of this system of equations, given the implicit dependence of the Euler-Lagrange equations on the time derivative of the potential.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Using an iterative eigensolver to compute vibrational energies with phase-spaced localized basis functions

Although phase-space localized Gaussians are themselves poor basis functions, they can be used to effectively contract a discrete variable representation basis [A. Shimshovitz and D. J. Tannor, Phys. Rev. Lett. 109, 070402 (2012)]. This works despite the fact that elements of the Hamiltonian and overlap matrices labelled by discarded Gaussians are not small. By formulating the matrix problem as a regular (i.e., not a generalized) matrix eigenvalue problem, we show that it is possible to use an iterative eigensolver to compute vibrational energy levels in the Gaussian basis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Efficient algorithms for Hirshfeld-I charges

A new viewpoint on iterative Hirshfeld charges is presented, whereby the atomic populations obtained from such a scheme are interpreted as such populations which reproduce themselves. This viewpoint yields a self-consistent requirement for the Hirshfeld-I populations rather than being understood as the result of an iterative procedure. Based on this self-consistent requirement, much faster algorithms for Hirshfeld-I charges have been developed. In addition, new atomic reference densities for the Hirshfeld-I procedure are presented. The proposed reference densities are N-representable, display proper atomic shell structure and can be computed for any charged species.

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

Granular chaos and mixing: Whirled in a grain of sand

In this paper, we overview examples of chaos in granular flows. We begin by reviewing several remarkable behaviors that have intrigued researchers over the past few decades, and we then focus on three areas in which chaos plays an intrinsic role in granular behavior. First, we discuss pattern formation in vibrated beds, which we show is a direct result of chaotic scattering combined with dynamical dissipation. Next, we consider stick-slip motion, which involves chaotic scattering on the micro-scale, and which results in complex and as yet unexplained peculiarities on the macro-scale. Finally, we examine granular mixing, which we show combines micro-scale chaotic scattering and macro-scale stick-slip motion into behaviors that are well described by dynamical systems tools, such as iterative mappings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Second-order nonlinear optical metamaterials: ABC-type nanolaminates

We demonstrate a concept for second-order nonlinear metamaterials that can be obtained from non-metallic centrosymmetric constituents with inherently low optical absorption. The concept is based on iterative atomic-layer deposition of three different materials, A = Al{sub 2}O{sub 3}, B = TiO{sub 2}, and C = HfO{sub 2}. The centrosymmetry of the resulting ABC stack is broken since the ABC and the inverted CBA sequences are not equivalent—a necessary condition for non-zero second-order nonlinearity. In our experiments, we find that the bulk second-order nonlinear susceptibility depends on the density of interfaces, leading to a nonlinear susceptibility of 0.26 pm/V at a wavelength of 800 nm. ABC-type nanolaminates can be deposited on virtually any substrate and offer a promising route towards engineering of second-order optical nonlinearities at both infrared and visible wavelengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Recent progress in plasma modelling at INFN-LNS

At Istituto Nazionale di Fisica Nucleare - Laboratori Nazionali del Sud (INFN-LNS), the development of intense ion and proton sources has been supported by a great deal of work on the modelling of microwave generated plasmas for many years. First, a stationary version of the particle-in-cell code was developed for plasma modelling starting from an iterative strategy adopted for the space charge dominated beam transport simulations. Electromagnetic properties of the plasma and full-waves simulations are now affordable for non-homogenous and non-isotropic magnetized plasma via “cold” approximation. The effects of Coulomb collisions on plasma particles dynamics was implemented with the Langevin formalism, instead of simply applying the Spitzer 90° collisions through a Monte Carlo technique. A wide database of different cross sections related to reactions occurring in a hydrogen plasma was implemented. The next step consists of merging such a variety of approaches for retrieving an “as-a-whole” picture of plasma dynamics in ion sources. The preliminary results will be summarized in the paper for a microwave discharge ion source designed for intense and high quality proton beams production, proton source for European Spallation Source project. Even if the realization of a predictive software including the complete processes involved in plasma formation is still rather far, a better comprehension of the source behavior is possible and so the simulations may support the optimization phase.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Collisional excitation of NH(X{sup 3}Σ{sup −}) by Ne: Potential energy surface, scattering calculations, and comparison with experiments

We present a new three-dimensional potential energy surface (PES) for the NH(X{sup 3}Σ{sup −})–Ne van der Waals system, which explicitly takes into account the NH vibrational motion. Ab initio calculations of the NH–Ne PES were carried out using the open-shell single- and double-excitation coupled cluster approach with non-iterative perturbational treatment of triple excitations [RCCSD(T)]. The augmented correlation-consistent quadruple zeta (aug-cc-pVQZ) basis set was employed. Mid-bond functions were also included in order to improve the accuracy in the van der Waals well. Using this new PES, we have studied the collisional excitation of NH(X{sup 3}Σ{sup −}) by Ne. Close-coupling calculations of the collisional excitation cross sections of the fine-structure levels of NH by Ne are performed for energies up to 3000 cm{sup −1}, which yield, after thermal average, rate coefficients up to 350 K. The propensity rules between fine-structure levels are reported, and it is found that F-conserving cross sections are larger than F-changing cross sections even if the propensity rules are not as strong as for the NH–He system. The calculated rate coefficients are compared with available experimental measurements at room temperature and a fairly good agreement is found between experimental and theoretical data, confirming the good quality of the scattering calculations and also the accuracy of the potential energy surface used in this work.

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

Modeling of electron energy spectra and mobilities in semi-metallic Hg{sub 1−x}Cd{sub x}Te quantum wells

Electron mobility, energy spectra, and intrinsic carrier concentrations in the n-type Hg{sub 0.32}Cd{sub 0.68}Te/Hg{sub 1−x}Cd{sub x}Te/Hg{sub 0.32}Cd{sub 0.68}Te quantum well (QW) in semi-metallic state are numerically modeled. Energy spectra and wave functions were calculated in the framework of the 8-band k-p Hamiltonian. In our model, electron scattering on longitudinal optical phonons, charged impurities, and holes has been taken into account, and the mobility has been calculated by an iterative solution of the Boltzmann transport equation. Our results show that the increase of the electron concentration in the well enhances the screening of the 2D electron gas, decreases the hole concentration, and can ultimately lead to a high electron mobility at liquid nitrogen temperatures. The increase of the electron concentration in the QW could be achieved in situ by delta-doping of barriers or by applying the top-gate potential. Our modeling has shown that for low molar composition x the concentration of holes in the well is high in a wide range of electron concentrations; in this case, the purity of samples does not significantly influence the electron mobility. These results are important in the context of establishing optimal parameters for the fabrication of high-mobility Hg{sub 1−x}Cd{sub x}Te quantum wells able to operate at liquid nitrogen temperature and thus suitable for applications in terahertz detectors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Quantum learning of classical stochastic processes: The completely positive realization problem

Among several tasks in Machine Learning, a specially important one is the problem of inferring the latent variables of a system and their causal relations with the observed behavior. A paradigmatic instance of this is the task of inferring the hidden Markov model underlying a given stochastic process. This is known as the positive realization problem (PRP), [L. Benvenuti and L. Farina, IEEE Trans. Autom. Control 49(5), 651–664 (2004)] and constitutes a central problem in machine learning. The PRP and its solutions have far-reaching consequences in many areas of systems and control theory, and is nowadays an important piece in the broad field of positive systems theory. We consider the scenario where the latent variables are quantum (i.e., quantum states of a finite-dimensional system) and the system dynamics is constrained only by physical transformations on the quantum system. The observable dynamics is then described by a quantum instrument, and the task is to determine which quantum instrument — if any — yields the process at hand by iterative application. We take as a starting point the theory of quasi-realizations, whence a description of the dynamics of the process is given in terms of linear maps on state vectors and probabilities are given by linear functionals on the state vectors. This description, despite its remarkable resemblance with the hidden Markov model, or the iterated quantum instrument, is however devoid of any stochastic or quantum mechanical interpretation, as said maps fail to satisfy any positivity conditions. The completely positive realization problem then consists in determining whether an equivalent quantum mechanical description of the same process exists. We generalize some key results of stochastic realization theory, and show that the problem has deep connections with operator systems theory, giving possible insight to the lifting problem in quotient operator systems. Our results have potential applications in quantum machine learning, device-independent characterization and reverse-engineering of stochastic processes and quantum processors, and more generally, of dynamical processes with quantum memory [M. Guţă, Phys. Rev. A 83(6), 062324 (2011); M. Guţă and N. Yamamoto, e-print http://arxiv.org/abs/1303.3771 (2013)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Coherent diffractive imaging using randomly coded masks

We experimentally demonstrate an extension to coherent diffractive imaging that encodes additional information through the use of a series of randomly coded masks, removing the need for typical object-domain constraints while guaranteeing a unique solution to the phase retrieval problem. Phase retrieval is performed using a numerical convex relaxation routine known as “PhaseCut,” an iterative algorithm known for its stability and for its ability to find the global solution, which can be found efficiently and which is robust to noise. The experiment is performed using a laser diode at 532.2 nm, enabling rapid prototyping for future X-ray synchrotron and even free electron laser experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Polarization effects on the electric properties of urea and thiourea molecules in solid phase

We present theoretical results for the dipole moment, linear polarizability, and first hyperpolarizability of the urea and thiourea molecules in solid phase. The in-crystal electric properties were determined by applying a supermolecule approach in combination with an iterative electrostatic scheme, in which the surrounding molecules are represented by point charges. It is found for both urea and thiourea molecules that the influence of the polarization effects is mild for the linear polarizability, but it is marked for the dipole moment and first hyperpolarizability. The replacement of oxygen atoms by sulfur atoms increases, in general, the electric responses. Our second-order Møller–Plesset perturbation theory based iterative scheme predicts for the in-crystal dipole moment of urea and thiourea the values of 7.54 and 9.19 D which are, respectively, increased by 61% and 58%, in comparison with the corresponding isolated values. The result for urea is in agreement with the available experimental result of 6.56 D. In addition, we present an estimate of macroscopic quantities considering explicit unit cells of urea and thiourea crystals including environment polarization effects. These supermolecule calculations take into account partially the exchange and dispersion effects. The results illustrate the role played by the electrostatic interactions on the static second-order nonlinear susceptibility of the urea crystal.

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

The glass transition in cured epoxy thermosets: A comparative molecular dynamics study in coarse-grained and atomistic resolution

We investigate the volumetric glass transition temperature T{sub g} in epoxy thermosets by means of molecular dynamics simulations. The epoxy thermosets consist of the resin bisphenol A diglycidyl ether and the hardener diethylenetriamine. A structure based coarse-grained (CG) force field has been derived using iterative Boltzmann inversion in order to facilitate simulations of larger length scales. We observe that T{sub g} increases clearly with the degree of cross-linking for all-atomistic (AA) and CG simulations. The transition T{sub g} in CG simulations of uncured mixtures is much lower than in AA-simulations due to the soft nature of the CG potentials, but increases all the more with the formation of rigid cross-links. Additional simulations of the CG mixtures in contact with a surface show the existence of an interphase region of about 3 nm thickness in which the network properties deviate significantly from the bulk. In accordance to experimental studies, we observe that T{sub g} is reduced in this interphase region and gradually increases to its bulk value with distance from the surface. The present study shows that the glass transition is a local phenomenon that depends on the network structure in the immediate environment.

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

Iterative quantum-classical path integral with dynamically consistent state hopping

We investigate the convergence of iterative quantum-classical path integral calculations in sluggish environments strongly coupled to a quantum system. The number of classical trajectories, thus the computational cost, grows rapidly (exponentially, unless filtering techniques are employed) with the memory length included in the calculation. We argue that the choice of the (single) trajectory branch during the time preceding the memory interval can significantly affect the memory length required for convergence. At short times, the trajectory branch associated with the reactant state improves convergence by eliminating spurious memory. We also introduce an instantaneous population-based probabilistic scheme which introduces state-to-state hops in the retained pre-memory trajectory branch, and which is designed to choose primarily the trajectory branch associated with the reactant at early times, but to favor the product state more as the reaction progresses to completion. Test calculations show that the dynamically consistent state hopping scheme leads to accelerated convergence and a dramatic reduction of computational effort.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Local conditions for the Pauli potential in order to yield self-consistent electron densities exhibiting proper atomic shell structure

The local conditions for the Pauli potential that are necessary in order to yield self-consistent electron densities from orbital-free calculations are investigated for approximations that are expressed with the help of a local position variable. It is shown that those local conditions also apply when the Pauli potential is given in terms of the electron density. An explicit formula for the Ne atom is given, preserving the local conditions during the iterative procedure. The resulting orbital-free electron density exhibits proper shell structure behavior and is in close agreement with the Kohn-Sham electron density. This study demonstrates that it is possible to obtain self-consistent orbital-free electron densities with proper atomic shell structure from simple one-point approximations for the Pauli potential at local density level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Absorption spectroscopy of heavy alkaline earth metals Ba and Sr in rare gas matrices—CCSD(T) calculations and atomic site occupancies

Isolation of the heavier alkaline earth metals Ba and Sr in the solid rare gases (RGs) Ar, Kr, and Xe is analysed with absorption spectroscopy and interpreted partly with the assistance of ab initio calculations of the diatomic M ⋅ RG ground state interaction potentials. The y{sup 1}P←a{sup 1}S resonance transitions in the visible spectral region are used to compare the isolation conditions of these two metal atom systems and calcium. Complex absorption bands were recorded in all three metal atom systems even after extensive sample annealing. Coupled cluster calculations conducted on the ground states of the nine M ⋅ RG diatomics (M = Ca, Sr, and Ba; RG = Ar, Kr, and Xe) at the coupled cluster single, double, and non-iterative triple level of theory revealed long bond lengths (>5 Å) and shallow bound regions (<130 cm{sup −1}). All of the M ⋅ RG diatomics have bond lengths considerably longer than those of the rare gas dimers, with the consequence that isolation of these metal atoms in a single substitutional site of the solid rare gas is unlikely, with the possible exception of Ca/Xe. The luminescence of metal dimer bands has been recorded for Ba and Sr revealing very different behaviours. Resonance fluorescence with a lifetime of 15 ns is observed for the lowest energy transition of Sr{sub 2} while this transition is quenched in Ba{sub 2}. This behaviour is consistent with the absence of vibrational structure on the dimer absorption band in Ba{sub 2} indicating lifetime broadening arising from efficient relaxation to low-lying molecular states. More extensive 2D excitation-emission data recorded for the complex site structures present on the absorption bands of the atomic Ba and Sr systems will be presented in future publications.

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

Design of a multi beam klystron cavity from its single beam parameters

The klystron is a well-known microwave amplifier which uses kinetic energy of an electron beam for amplification of the RF signal. There are some limitations of conventional single beam klystron such as high operating voltage, low efficiency and bulky size at higher power levels, which are very effectively handled in Multi Beam Klystron (MBK) that uses multiple low purveyance electron beams for RF interaction. Each beam propagates along its individual transit path through a resonant cavity structure. Multi-Beam klystron cavity design is a critical task due to asymmetric cavity structure and can be simulated by 3D code only. The present paper shall discuss the design of multi beam RF cavities for klystrons operating at 2856 MHz (S-band) and 5 GHz (C-band) respectively. The design approach uses some scaling laws for finding the electron beam parameters of the multi beam device from their single beam counter parts. The scaled beam parameters are then used for finding the design parameters of the multi beam cavities. Design of the desired multi beam cavity can be optimized through iterative simulations in CST Microwave Studio.

42 ENGINEERING↗