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 685 records · Page 38

The segmented non-uniform dielectric module design for uniformity control of plasma profile in a capacitively coupled plasma chamber

Low-temperature plasma technique is one of the critical techniques in IC manufacturing process, such as etching and thin-film deposition, and the uniformity greatly impacts the process quality, so the design for the plasma uniformity control is very important but difficult. It is hard to finely and flexibly regulate the spatial distribution of the plasma in the chamber via controlling the discharge parameters or modifying the structure in zero-dimensional space, and it just can adjust the overall level of the process factors. In the view of this problem, a segmented non-uniform dielectric module design solution is proposed for the regulation of the plasma profile in a CCP chamber. The solution achieves refined and flexible regulation of the plasma profile in the radial direction via configuring the relative permittivity and the width of each segment. In order to solve this design problem, a novel simulation-based auto-design approach is proposed, which can automatically design the positional sequence with multi independent variables to make the output target profile in the parameterized simulation model approximate the one that users preset. This approach employs an idea of quasi-closed-loop control system, and works in an iterative mode. It starts from initial values of the design variable sequences, and predicts better sequences via the feedback of the profile error between the output target profile and the expected one. It never stops until the profile error is narrowed in the preset tolerance.

36 MATERIALS SCIENCE↗

Accurate coarse-grained models for mixtures of colloids and linear polymers under good-solvent conditions

A coarse-graining strategy, previously developed for polymer solutions, is extended here to mixtures of linear polymers and hard-sphere colloids. In this approach, groups of monomers are mapped onto a single pseudoatom (a blob) and the effective blob-blob interactions are obtained by requiring the model to reproduce some large-scale structural properties in the zero-density limit. We show that an accurate parametrization of the polymer-colloid interactions is obtained by simply introducing pair potentials between blobs and colloids. For the coarse-grained (CG) model in which polymers are modelled as four-blob chains (tetramers), the pair potentials are determined by means of the iterative Boltzmann inversion scheme, taking full-monomer (FM) pair correlation functions at zero-density as targets. For a larger number n of blobs, pair potentials are determined by using a simple transferability assumption based on the polymer self-similarity. We validate the model by comparing its predictions with full-monomer results for the interfacial properties of polymer solutions in the presence of a single colloid and for thermodynamic and structural properties in the homogeneous phase at finite polymer and colloid density. The tetramer model is quite accurate for q ≲ 1 (q=R{sup ^}{sub g}/R{sub c}, where R{sup ^}{sub g} is the zero-density polymer radius of gyration and R{sub c} is the colloid radius) and reasonably good also for q = 2. For q = 2, an accurate coarse-grained description is obtained by using the n = 10 blob model. We also compare our results with those obtained by using single-blob models with state-dependent potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectral restoration in high resolution electron energy loss spectroscopy based on iterative semi-blind Lucy-Richardson algorithm applied to rutile surfaces

A new spectral restoration algorithm of reflection electron energy loss spectra is proposed. It is based on the maximum likelihood principle as implemented in the iterative Lucy-Richardson approach. Resolution is enhanced and point spread function recovered in a semi-blind way by forcing cyclically the zero loss to converge towards a Dirac peak. Synthetic phonon spectra of TiO{sub 2} are used as a test bed to discuss resolution enhancement, convergence benefit, stability towards noise, and apparatus function recovery. Attention is focused on the interplay between spectral restoration and quasi-elastic broadening due to free carriers. A resolution enhancement by a factor up to 6 on the elastic peak width can be obtained on experimental spectra of TiO{sub 2}(110) and helps revealing mixed phonon/plasmon excitations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Fixed-node diffusion Monte Carlo potential energy curve of the fluorine molecule F{sub 2} using selected configuration interaction trial wavefunctions

The potential energy curve of the F{sub 2} molecule is calculated with Fixed-Node Diffusion Monte Carlo (FN-DMC) using Configuration Interaction (CI)-type trial wavefunctions. To keep the number of determinants reasonable and thus make FN-DMC calculations feasible in practice, the CI expansion is restricted to those determinants that contribute the most to the total energy. The selection of the determinants is made using the CIPSI approach (Configuration Interaction using a Perturbative Selection made Iteratively). The trial wavefunction used in FN-DMC is directly issued from the deterministic CI program; no Jastrow factor is used and no preliminary multi-parameter stochastic optimization of the trial wavefunction is performed. The nodes of CIPSI wavefunctions are found to reduce significantly the fixed-node error and to be systematically improved upon increasing the number of selected determinants. To reduce the non-parallelism error of the potential energy curve, a scheme based on the use of a R-dependent number of determinants is introduced. Using Dunning’s cc-pVDZ basis set, the FN-DMC energy curve of F{sub 2} is found to be of a quality similar to that obtained with full configuration interaction/cc-pVQZ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rode's iterative calculation of surface optical phonon scattering limited electron mobility in N-polar GaN devices

N-polar GaN channel mobility is important for high frequency device applications. Here, we report theoretical calculations on the surface optical (SO) phonon scattering rate of two-dimensional electron gas (2DEG) in N-polar GaN quantum well channels with high-k dielectrics. Rode's iterative calculation is used to predict the scattering rate and mobility. Coupling of the GaN plasmon modes with the SO modes is taken into account and dynamic screening is employed under linear polarization response. The effect of SO phonons on 2DEG mobility was found to be small at >5 nm channel thickness. However, the SO mobility in 3 nm N-polar GaN channels with HfO{sub 2} and ZrO{sub 2} high-k dielectrics is low and limits the total mobility. The SO scattering for SiN dielectric on GaN was found to be negligible due to its high SO phonon energy. Using Al{sub 2}O{sub 3}, the SO phonon scattering does not affect mobility significantly only except the case when the channel is too thin with a low 2DEG density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Icing modelling in NSMB with chimera overset grids

In aerospace Engineering, the accurate simulation of ice accretion is a key element to increase flight safety and avoid accidents related to icing effects. The icing code developed in the NSMB solver is based on an Eulerian formulation for droplets tracking, an iterative Messinger model using a modified water runback scheme for ice thickness calculation and mesh deformation to track the ice/air interface through time. The whole process is parallelized with MPI and applied with chimera grids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Launching the chaotic realm of iso-fractals: A short remark

In this brief note, we introduce the new, emerging sub-discipline of iso-fractals by highlighting and discussing the preliminary results of recent works. First, we note the abundance of fractal, chaotic, non-linear, and self-similar structures in nature while emphasizing the importance of studying such systems because fractal geometry is the language of chaos. Second, we outline the iso-fractal generalization of the Mandelbrot set to exemplify the newly generated Mandelbrot iso-sets. Third, we present the cutting-edge notion of dynamic iso-spaces and explain how a mathematical space can be iso-topically lifted with iso-unit functions that (continuously or discretely) change; in the discrete case examples, we mention that iteratively generated sequences like Fibonacci’s numbers and (the complex moduli of) Mandelbrot’s numbers can supply a deterministic chain of iso-units to construct an ordered series of (magnified and/or de-magnified) iso-spaces that are locally iso-morphic. Fourth, we consider the initiation of iso-fractals with Inopin’s holographic ring (IHR) topology and fractional statistics for 2D and 3D iso-spaces. In total, the reviewed iso-fractal results are a significant improvement over traditional fractals because the application of Santilli’s iso-mathematics arms us an extra degree of freedom for attacking problems in chaos. Finally, we conclude by proposing some questions and ideas for future research work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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 ↗

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↗

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↗

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↗