Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SECOND”

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,045 records · Page 58

Efficient linear-scaling second-order Møller-Plesset perturbation theory: The divide–expand–consolidate RI-MP2 model

The Resolution of the Identity second-order Møller-Plesset perturbation theory (RI-MP2) method is implemented within the linear-scaling Divide-Expand-Consolidate (DEC) framework. In a DEC calculation, the full molecular correlated calculation is replaced by a set of independent fragment calculations each using a subset of the total orbital space. The number of independent fragment calculations scales linearly with the system size, rendering the method linear-scaling and massively parallel. The DEC-RI-MP2 method can be viewed as an approximation to the DEC-MP2 method where the RI approximation is utilized in each fragment calculation. The individual fragment calculations scale with the fifth power of the fragment size for both methods. However, the DEC-RI-MP2 method has a reduced prefactor compared to DEC-MP2 and is well-suited for implementation on massively parallel supercomputers, as demonstrated by test calculations on a set of medium-sized molecules. The DEC error control ensures that the standard RI-MP2 energy can be obtained to the predefined precision. The errors associated with the RI and DEC approximations are compared, and it is shown that the DEC-RI-MP2 method can be applied to systems far beyond the ones that can be treated with a conventional RI-MP2 implementation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observation of turnover of spontaneous polarization in ferroelectric layer of pentacene/poly-(vinylidene-trifluoroethylene) double-layer capacitor under photo illumination by optical second-harmonic generation measurement

The details of turnover process of spontaneous polarization and associated carrier motions in indium-tin oxide/poly-(vinylidene-trifluoroethylene)/pentacene/Au capacitor were analyzed by coupling displacement current measurement (DCM) and electric-field-induced optical second-harmonic generation (EFISHG) measurement. A model was set up from DCM results to depict the relationship between electric field in semiconductor layer and applied external voltage, proving that photo illumination effect on the spontaneous polarization process lied in variation of semiconductor conductivity. The EFISHG measurement directly and selectively probed the electric field distribution in semiconductor layer, modifying the model and revealing detailed carrier behaviors involving photo illumination effect, dipole reversal, and interfacial charging in the device. A further decrease of DCM current in the low voltage region under illumination was found as the result of illumination effect, and the result was argued based on the changing of the total capacitance of the double-layer capacitors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An upper bound on the second order asymptotic expansion for the quantum communication cost of state redistribution

State redistribution is the protocol in which given an arbitrary tripartite quantum state, with two of the subsystems initially being with Alice and one being with Bob, the goal is for Alice to send one of her subsystems to Bob, possibly with the help of prior shared entanglement. We derive an upper bound on the second order asymptotic expansion for the quantum communication cost of achieving state redistribution with a given finite accuracy. In proving our result, we also obtain an upper bound on the quantum communication cost of this protocol in the one-shot setting, by using the protocol of coherent state merging as a primitive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Explicitly correlated second-order Møller-Plesset perturbation theory in a Divide-Expand-Consolidate (DEC) context

We present the DEC-RIMP2-F12 method where we have augmented the Divide Expand-Consolidate resolution-of-the-identity second-order Møller-Plesset perturbation theory method (DEC-RIMP2) [P. Baudin et al., J. Chem. Phys. 144, 054102 (2016)] with an explicitly correlated (F12) correction. Here, the new method is linear-scaling, massively parallel, and it corrects for the basis set incompleteness error in an efficient manner. In addition, we observe that the F12 contribution decreases the domain error of the DEC-RIMP2 correlation energy by roughly an order of magnitude. An important feature of the DEC scheme is the inherent error control defined by a single parameter, and this feature is also retained for the DEC-RIMP2-F12 method. In this paper we present the working equations for the DEC-RIMP2-F12 method and proof of concept numerical results for a set of test molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

G4CEP: A G4 theory modification by including pseudopotential for molecules containing first-, second- and third-row representative elements

The G4CEP composite method was developed from the respective G4 all-electron version by considering the implementation of compact effective pseudopotential (CEP). The G3/05 test set was used as reference to benchmark the adaptation by treating in this work atoms and compounds from the first and second periods of the periodic table, as well as representative elements of the third period, comprising 440 thermochemical data. G4CEP has not reached a so high level of accuracy as the G4 all-electron theory. G4CEP presented a mean absolute error around 1.09 kcal mol{sup −1}, while the original method presents a deviation corresponding to 0.83 kcal mol{sup −1}. The similarity of the optimized molecular geometries between G4 and G4CEP indicates that the core-electron effects and basis set adjustments may be pointed out as a significant factor responsible for the large discrepancies between the pseudopotential results and the experimental data, or even that the all-electron calculations are more efficient either in its formulation or in the cancellation of errors. When the G4CEP mean absolute error (1.09 kcal mol{sup −1}) is compared to 1.29 kcal mol{sup −1} from G3CEP, it does not seem so efficient. However, while the G3CEP uncertainty is ±4.06 kcal mol{sup −1}, the G4CEP deviation is ±2.72 kcal mol{sup −1}. Therefore, the G4CEP theory is considerably more reliable than any previous combination of composite theory and pseudopotential, particularly for enthalpies of formation and electron affinities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Second-order asymptotics for quantum hypothesis testing in settings beyond i.i.d.—quantum lattice systems and more

Quantum Stein’s lemma is a cornerstone of quantum statistics and concerns the problem of correctly identifying a quantum state, given the knowledge that it is one of two specific states (ρ or σ). It was originally derived in the asymptotic i.i.d. setting, in which arbitrarily many (say, n) identical copies of the state (ρ{sup ⊗n} or σ{sup ⊗n}) are considered to be available. In this setting, the lemma states that, for any given upper bound on the probability α{sub n} of erroneously inferring the state to be σ, the probability β{sub n} of erroneously inferring the state to be ρ decays exponentially in n, with the rate of decay converging to the relative entropy of the two states. The second order asymptotics for quantum hypothesis testing, which establishes the speed of convergence of this rate of decay to its limiting value, was derived in the i.i.d. setting independently by Tomamichel and Hayashi, and Li. We extend this result to settings beyond i.i.d. Examples of these include Gibbs states of quantum spin systems (with finite-range, translation-invariant interactions) at high temperatures, and quasi-free states of fermionic lattice gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of first and second generation biodiesel blends on engine performance and emission

The biodiesel is a potential source of alternative fuel which can be used at different proportions with diesel fuel. This study experimentally investigated the effect of blend percentage on diesel engine performance and emission using first generation (soybean) and second generation (waste cooking) biodiesel. The characterization of the biodiesel was done according to ASTM and EN standards and compared with ultralow sulfur diesel (ULSD) fuel. A multi-cylinder test bed engine coupled with electromagnetic dynamometer and 5 gas analyzer were used for engine performance and emission test. The investigation was made using B5, B10 and B15 blends for both biodiesels. The study found that brake power (BP) and brake torque (BT) slightly decreases and brake specific fuel consumption (BSFC) slightly increases with an increase in biodiesel blends ratio. Besides, a significant reduction in exhaust emissions (except NO{sub x} emission) was found for both biodiesels compared to ULSD. Soybean biodiesel showed better engine performance and emissions reduction compared with waste cooking biodiesel. However, NO{sub x} emission for B5 waste cooking biodiesel was lower than soybean biodiesel.

09 BIOMASS FUELS↗

Analytic second derivative of the energy for density-functional tight-binding combined with the fragment molecular orbital method

The analytic second derivative of the energy is developed for the fragment molecular orbital (FMO) method combined with density-functional tight-binding (DFTB), enabling simulations of infrared and Raman spectra of large molecular systems. The accuracy of the method is established in comparison to full DFTB without fragmentation for a set of representative systems. The performance of the FMO-DFTB Hessian is discussed for molecular systems containing up to 10 041 atoms. The method is applied to the study of the binding of α-cyclodextrin to polyethylene glycol, and the calculated IR spectrum of an epoxy amine oligomer reproduces experiment reasonably well.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Second order gyrokinetic theory for particle-in-cell codes

The main idea of the gyrokinetic dynamical reduction consists in a systematical removal of the fast scale motion (the gyromotion) from the dynamics of the plasma, resulting in a considerable simplification and a significant gain of computational time. The gyrokinetic Maxwell–Vlasov equations are nowadays implemented in for modeling (both laboratory and astrophysical) strongly magnetized plasmas. Different versions of the reduced set of equations exist, depending on the construction of the gyrokinetic reduction procedure and the approximations performed in the derivation. The purpose of this article is to explicitly show the connection between the general second order gyrokinetic Maxwell–Vlasov system issued from the modern gyrokinetic theory and the model currently implemented in the global electromagnetic Particle-in-Cell code ORB5. Necessary information about the modern gyrokinetic formalism is given together with the consistent derivation of the gyrokinetic Maxwell–Vlasov equations from first principles. The variational formulation of the dynamics is used to obtain the corresponding energy conservation law, which in turn is used for the verification of energy conservation diagnostics currently implemented in ORB5. This work fits within the context of the code verification project VeriGyro currently run at IPP Max-Planck Institut in collaboration with others European institutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Note: Real-time monitoring via second-harmonic interferometry of a flow gas cell for laser wakefield acceleration

The use of a gas cell as a target for laser wakefield acceleration (LWFA) offers the possibility to obtain stable and manageable laser-plasma interaction process, a mandatory condition for practical applications of this emerging technique, especially in multi-stage accelerators. In order to obtain full control of the gas particle number density in the interaction region, thus allowing for a long term stable and manageable LWFA, real-time monitoring is necessary. In fact, the ideal gas law cannot be used to estimate the particle density inside the flow cell based on the preset backing pressure and the room temperature because the gas flow depends on several factors like tubing, regulators, and valves in the gas supply system, as well as vacuum chamber volume and vacuum pump speed/throughput. Here, second-harmonic interferometry is applied to measure the particle number density inside a flow gas cell designed for LWFA. The results demonstrate that real-time monitoring is achieved and that using low backing pressure gas (<1 bar) and different cell orifice diameters (<2 mm) it is possible to finely tune the number density up to the 10{sup 19} cm{sup −3} range well suited for LWFA.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Monte Carlo explicitly correlated second-order many-body perturbation theory

A stochastic algorithm is proposed and implemented that computes a basis-set-incompleteness (F12) correction to an ab initio second-order many-body perturbation energy as a short sum of 6- to 15-dimensional integrals of Gaussian-type orbitals, an explicit function of the electron-electron distance (geminal), and its associated excitation amplitudes held fixed at the values suggested by Ten-no. The integrals are directly evaluated (without a resolution-of-the-identity approximation or an auxiliary basis set) by the Metropolis Monte Carlo method. Applications of this method to 17 molecular correlation energies and 12 gas-phase reaction energies reveal that both the nonvariational and variational formulas for the correction give reliable correlation energies (98% or higher) and reaction energies (within 2 kJ mol -1 with a smaller statistical uncertainty) near the complete-basis-set limits by using just the aug-cc-pVDZ basis set. The nonvariational formula is found to be 2–10 times less expensive to evaluate than the variational one, though the latter yields energies that are bounded from below and is, therefore, slightly but systematically more accurate for energy differences. Being capable of using virtually any geminal form, the method confirms the best overall performance of the Slater-type geminal among 6 forms satisfying the same cusp conditions. Lastly, not having to precompute lower-dimensional integrals analytically, to store them on disk, or to transform them in a nonscalable dense-matrix-multiplication algorithm, the method scales favorably with both system size and computer size; the cost increases only as O(n 4 ) with the number of orbitals (n), and its parallel efficiency reaches 99.9% of the ideal case on going from 16 to 4096 computer processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring connections between statistical mechanics and Green’s functions for realistic systems: Temperature dependent electronic entropy and internal energy from a self-consistent second-order Green’s function

Including finite-temperature effects from the electronic degrees of freedom in electronic structure calculations of semiconductors and metals is desired; however, in practice it remains exceedingly difficult when using zero-temperature methods, since these methods require an explicit evaluation of multiple excited states in order to account for any finite-temperature effects. Using a Matsubara Green’s function formalism remains a viable alternative, since in this formalism it is easier to include thermal effects and to connect the dynamic quantities such as the self-energy with static thermodynamic quantities such as the Helmholtz energy, entropy, and internal energy. However, despite the promising properties of this formalism, little is known about the multiple solutions of the non-linear equations present in the self-consistent Matsubara formalism and only a few cases involving a full Coulomb Hamiltonian were investigated in the past. In this study, to shed some light onto the iterative nature of the Green’s function solutions, we self-consistently evaluate the thermodynamic quantities for a one-dimensional (1D) hydrogen solid at various interatomic separations and temperatures using the self-energy approximated to second-order (GF2). At many points in the phase diagram of this system, multiple phases such as a metal and an insulator exist, and we are able to determine the most stable phase from the analysis of Helmholtz energies. Additionally, we show the evolution of the spectrum of 1D boron nitride to demonstrate that GF2 is capable of qualitatively describing the temperature effects influencing the size of the band gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Random phase approximation with second-order screened exchange for current-carrying atomic states

The direct random phase approximation (RPA) and RPA with second-order screened exchange (SOSEX) have been implemented with complex orbitals as a basis for treating open-shell atoms. Both RPA and RPA+SOSEX are natural implicit current density functionals because the paramagnetic current density implicitly is included through the use of complex orbitals. We confirm that inclusion of the SOSEX correction improves the total energy accuracy substantially compared to RPA, especially for smaller-Z atoms. Computational complexity makes post self-consistent-field (post-SCF) evaluation of RPA-type expressions commonplace, so orbital basis origins and properties become important. Sizable differences are found in correlation energies, total atomic energies, and ionization energies for RPA-type functionals evaluated in the post-SCF fashion with orbital sets obtained from different schemes. Reference orbitals from Kohn-Sham calculations with semi-local functionals are more suitable for RPA+SOSEX to generate accurate total energies, but reference orbitals from exact exchange (non-local) yield essentially energetically degenerate open-shell atom ground states. RPA+SOSEX correlation combined with exact exchange calculated from a hybrid reference orbital set (half the exchange calculated from exact-exchange orbitals, the other half of the exchange from orbitals optimized for the Perdew-Burke-Ernzerhof (PBE) exchange functional) gives the best overall performance. Finally, numerical results show that the RPA-like functional with SOSEX correction can be used as a practical implicit current density functional when current effects should be included.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Second order Møller-Plesset and coupled cluster singles and doubles methods with complex basis functions for resonances in electron-molecule scattering

The method of complex basis functions is applied to molecular resonances at correlated levels of theory. Møller-Plesset perturbation theory at second order and equation-of-motion electron attachment coupled-cluster singles and doubles (EOM-EA-CCSD) methods based on a non-Hermitian self-consistent-field reference are used to compute accurate Siegert energies for shape resonances in small molecules including N 2 - , CO - , CO 2 - , and CH 2 O - . Analytic continuation of complex θ-trajectories is used to compute Siegert energies, and the θ-trajectories of energy differences are found to yield more consistent results than those of total energies. Furthermore, the ability of such methods to accurately compute complex potential energy surfaces is investigated, and the possibility of using EOM-EA-CCSD for Feshbach resonances is explored in the context of e-helium scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Normal-mode-based analysis of electron plasma waves with second-order Hermitian formalism

The classic problem of the dynamic evolution and Landau damping of linear Langmuir electron waves in a collisionless plasma with Maxwellian background is cast as a second-order, self-adjoint problem with a continuum spectrum of real and positive squared frequencies. The corresponding complete basis of singular normal modes is obtained, along with their orthogonality relation. This yields easily the general expression of the time-reversal-invariant solution for any initial-value problem. Examples are then given for specific initial conditions that illustrate different behaviors of the Landau-damped macroscopic moments of the perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-sensitivity plasma density retrieval in a common-path second-harmonic interferometer through simultaneous group and phase velocity measurement

Precise measurements of the plasma density in ionized gas cells and discharged capillaries are critical to the design and operation of plasma-based accelerators, active plasma lenses, and plasma-based radiation sources. In this manuscript, a spectral-domain common-path second-harmonic interferometer is upgraded with the simultaneous measurement of the group and phase velocity, allowing for high-sensitivity density characterization (from the phase velocity advance) without the need for phase tracking from zero-density (enabled by the group velocity delay). The technique is applied to 1.5-cm-long plasma structures, without density ambiguity in parameter scans with >2π phase jumps. Finally, the single-shot sensitivity in phase retrieval is demonstrated at 63 mrad, equivalent to a density-length product of 1.8·10 15 cm -2 . This is an improvement of ×45 compared to group velocity analysis alone.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Second harmonic microscopy of poled x-cut thin film lithium niobate: Understanding the contrast mechanism

Thin film lithium niobate has been of great interest recently, and an understanding of periodically poled thin films is crucial for both fundamental physics and device developments. Second-harmonic (SH) microscopy allows for the noninvasive visualization and analysis of ferroelectric domain structures and walls. While the technique is well understood in bulk lithium niobate, SH microscopy in thin films is largely influenced by interfacial reflections and resonant enhancements, which depend on film thicknesses and substrate materials. We present a comprehensive analysis of SH microscopy in x-cut lithium niobate thin films, based on a full three-dimensional focus calculation and accounting for interface reflections. We show that the dominant signal in backreflection originates from a copropagating phase-matched process observed through reflections, rather than direct detection of the counterpropagating signal as in bulk samples. We simulate the SH signatures of domain structures by a simple model of the domain wall as an extensionless transition from a −χ(2) to a +χ(2) region. This allows us to explain the main observation of domain structures in the thin-film geometry, and, in particular, we show that the SH signal from thin poled films allows to unambiguously distinguish areas, which are completely or only partly inverted in depth.

Rüsing, M. (ORCID:0000000346824577)↗