Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems
Not Available
Xiao, Jianyuan↗
SEARCH · Engineering Papers
Results for “interpolation”
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.
Thermodynamic limit of bifacial double-junction tandem solar cells
Here, a traditional single-junction solar panel cannot harness ground-scattered light (albedo reflectance, R A ), and also suffers from the fundamental sub-band-gap and the thermalization losses. In this paper, we explain how a “bifacial tandem” panel would dramatically reduce these losses, with corresponding improvement in thermodynamic performance. Our study predicts (i) the optimum combination of the band-gaps, empirically given by E opt g(t) ≈ E opt g(b) (2+R A )/3+(1–R A ) and the (ii) corresponding optimum normalized output power given by η* T(opt) ≈ R A (2η SJ(opt) ) + (1–R A )η DJ(opt) . Empirically, η* T(opt) interpolates between the thermodynamic efficiency limit of classical double-junction tandem cell (η DJ ) and twice that of a single-junction cell (η SJ ). We conclude by explaining how the fundamental loss mechanisms evolve with R A in a bifacial tandem cell.
14 SOLAR ENERGY↗
Strain effect on electronic structure and thermoelectric properties of orthorhombic SnSe: A first principles study
Strain effect on thermoelectricity of orthorhombic SnSe is studied using density function theory. The Seebeck coefficients are obtained by solving Boltzmann Transport equation (BTE) with interpolated band energies. As expected from the crystal structure, calculated Seebeck coefficients are highly anisotropic, and agree well with experiment. Changes in the Seebeck coefficients are presented, when strain is applied along b and c direction with strength from -3% to +3%, where influence by band gaps and band dispersions are significant. Moreover, for compressive strains, the sign change of Seebeck coefficients at particular direction suggests that the bipolar transport is possible for SnSe.
36 MATERIALS SCIENCE↗
A coarse-grain force field for RDX: Density dependent and energy conserving
We describe the development of a density-dependent transferable coarse-grain model of crystalline hexahydro-1,3,5-trinitro-s-triazine (RDX) that can be used with the energy conserving dissipative particle dynamics method. The model is an extension of a recently reported one-site model of RDX that was developed by using a force-matching method. The density-dependent forces in that original model are provided through an interpolation scheme that poorly conserves energy. The development of the new model presented in this work first involved a multi-objective procedure to improve the structural and thermodynamic properties of the previous model, followed by the inclusion of the density dependency via a conservative form of the force field that conserves energy. The new model accurately predicts the density, structure, pressure-volume isotherm, bulk modulus, and elastic constants of the RDX crystal at ambient pressure and exhibits transferability to a liquid phase at melt conditions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Relativistic theory of the Jahn-Teller effect: p-orbitals in tetrahedral and trigonal systems
A relativistic generalization of Jahn-Teller theory is presented which includes spin-orbit coupling effects beyond low-order Taylor expansions in vibrational coordinates. For the example of a p-electron in tetrahedral and trigonal environments, the matrix elements of the Breit-Pauli spin-orbit-coupling operator are expressed in terms of the matrix elements of the electrostatic electronic potential. Employing expansions of the latter in invariant polynomials in symmetry-adapted nuclear coordinates, the spin-orbit induced Jahn-Teller coupling terms are derived for the T{sub 2} × (t{sub 2} + e) and (E + A) × (e + a) Jahn-Teller problems up to arbitrarily high orders. The linear G{sub 3/2} × (t{sub 2} + e) Jahn-Teller Hamiltonian of Moffitt and Thorson [Phys. Rev. 108, 1251 (1957)] for tetrahedral systems is generalized to higher orders in vibrational displacements. The Jahn-Teller Hamiltonians derived in the present work are useful for the interpolation and extrapolation of Jahn-Teller distorted potential-energy surfaces of molecules and complexes with heavy elements as well as for the calculation of vibronic spectra of such systems.
74 ATOMIC AND MOLECULAR PHYSICS↗
Approximating basins of attraction for dynamical systems via stable radial bases
In applied sciences it is often required to model and supervise temporal evolution of populations via dynamical systems. In this paper, we focus on the problem of approximating the basins of attraction of such models for each stable equilibrium point. We propose to reconstruct the basins via an implicit interpolant using stable radial bases, obtaining the surfaces by partitioning the phase space into disjoint regions. An application to a competition model presenting jointly three stable equilibria is considered.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Extent of hydrogen coverage of Si(001) under chemical vapor deposition conditions from ab initio approaches
The extent of hydrogen coverage of the Si(001) c(4 × 2) surface in the presence of hydrogen gas has been studied with dispersion corrected density functional theory. Electronic energy contributions are well described using a hybrid functional. The temperature dependence of the coverage in thermodynamic equilibrium was studied computing the phonon spectrum in a supercell approach. As an approximation to these demanding computations, an interpolated phonon approach was found to give comparable accuracy. The simpler ab initio thermodynamic approach is not accurate enough for the system studied, even if corrections by the Einstein model for surface vibrations are considered. The on-set of H{sub 2} desorption from the fully hydrogenated surface is predicted to occur at temperatures around 750 K. Strong changes in hydrogen coverage are found between 1000 and 1200 K in good agreement with previous reflectance anisotropy spectroscopy experiments. These findings allow a rational choice for the surface state in the computational treatment of chemical reactions under typical metal organic vapor phase epitaxy conditions on Si(001).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The first potential energy surfaces for the C{sub 6}H{sup −}–H{sub 2} and C{sub 6}H{sup −}–He collisional systems and their corresponding inelastic cross sections
Molecular anions have recently been detected in the interstellar and circumstellar media. Accurate modeling of their abundance requires calculations of collisional data with the most abundant species that are usually He atoms and H{sub 2} molecules. In this paper, we focus on the collisional excitation of the first observed molecular anion, C{sub 6}H{sup −}, by He and H{sub 2}. Theoretical calculations of collisional cross sections rely generally on ab initio interaction potential energy surfaces (PESs). Hence, we present here the first PESs for the C{sub 6}H{sup −}–H{sub 2} and C{sub 6}H{sup −}–He van der Waals systems. The ab initio energy data for the surfaces were computed at the explicitly correlated coupled cluster with single, double, and scaled perturbative triple excitations level of theory. The method of interpolating moving least squares was used to construct 4D and 2D analytical PESs from these data. Both surfaces are characterized by deep wells and large anisotropies. Analytical models of the PESs were used in scattering calculations to obtain cross sections for low-lying rotational transitions. As could have been anticipated, important differences exist between the He and H{sub 2} cross sections. Conversely, no significant differences exist between the collisions of C{sub 6}H{sup −} with the two species of H{sub 2} (para- and ortho-H{sub 2}). We expect that these new data will help in accurately determining the abundance of the C{sub 6}H{sup −} anions in space.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Tokamak-independent software analysis suite for multi-spectral line-polarization MSE diagnostics
A tokamak-independent analysis suite has been developed to process data from Motional Stark Effect (mse) diagnostics. The software supports multi-spectral line-polarization mse diagnostics which simultaneously measure emission at the mse σ and π lines as well as at two "background" wavelengths that are displaced from the mse spectrum by a few nanometers. This analysis accurately estimates the amplitude of partially polarized background light at the σ and π wavelengths even in situations where the background light changes rapidly in time and space, a distinct improvement over traditional "time-interpolation" background estimation. The signal amplitude at many frequencies is computed using a numerical-beat algorithm which allows the retardance of the mse photo-elastic modulators (pem's) to be monitored during routine operation. It also allows the use of summed intensities at multiple frequencies in the calculation of polarization direction, which increases the effective signal strength and reduces sensitivity to pem retardance drift. The software allows the polarization angles to be corrected for calibration drift using a system that illuminates the mse diagnostic with polarized light at four known polarization angles within ten seconds of a plasma discharge. As a result, the software suite is modular, parallelized, and portable to other facilities.
97 MATHEMATICS AND COMPUTING↗
Experimental and theoretical study of rotationally inelastic diffraction of H{sub 2}(D{sub 2}) from methyl-terminated Si(111)
Fundamental details concerning the interaction between H{sub 2} and CH{sub 3}–Si(111) have been elucidated by the combination of diffractive scattering experiments and electronic structure and scattering calculations. Rotationally inelastic diffraction (RID) of H{sub 2} and D{sub 2} from this model hydrocarbon-decorated semiconductor interface has been confirmed for the first time via both time-of-flight and diffraction measurements, with modest j = 0 → 2 RID intensities for H{sub 2} compared to the strong RID features observed for D{sub 2} over a large range of kinematic scattering conditions along two high-symmetry azimuthal directions. The Debye-Waller model was applied to the thermal attenuation of diffraction peaks, allowing for precise determination of the RID probabilities by accounting for incoherent motion of the CH{sub 3}–Si(111) surface atoms. The probabilities of rotationally inelastic diffraction of H{sub 2} and D{sub 2} have been quantitatively evaluated as a function of beam energy and scattering angle, and have been compared with complementary electronic structure and scattering calculations to provide insight into the interaction potential between H{sub 2} (D{sub 2}) and hence the surface charge density distribution. Specifically, a six-dimensional potential energy surface (PES), describing the electronic structure of the H{sub 2}(D{sub 2})/CH{sub 3}−Si(111) system, has been computed based on interpolation of density functional theory energies. Quantum and classical dynamics simulations have allowed for an assessment of the accuracy of the PES, and subsequently for identification of the features of the PES that serve as classical turning points. A close scrutiny of the PES reveals the highly anisotropic character of the interaction potential at these turning points. This combination of experiment and theory provides new and important details about the interaction of H{sub 2} with a hybrid organic-semiconductor interface, which can be used to further investigate energy flow in technologically relevant systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A reactant-coordinate-based approach to state-to-state differential cross sections for tetratomic reactions
A new algorithm is proposed to compute quantum mechanically state-to-state differential cross sections for reactions involving four atoms in full dimensionality. This algorithm, which is based on the propagation of an initial state specific wave packet exclusively in reactant coordinates, extracts the S-matrix elements in the product channel by first interpolating the time-dependent wave packet using a collocation method at selected time intervals on the product coordinate grid and then projecting out the contributions of all final product states. Here, this approach is efficient and accurate, particularly for reactions that are dominated by a product well or long-range interactions. Validation of this approach is demonstrated for the H 2 + OH → H + H 2 O reaction.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
Here, an explicit high-order noncanonical symplectic algorithm for ideal two-fluid systems is developed. The fluid is discretized as particles in the Lagrangian description, while the electromagnetic fields and internal energy are treated as discrete differential form fields on a fixed mesh. With the assistance of Whitney interpolating forms, this scheme preserves the gauge symmetry of the electromagnetic field, and the pressure field is naturally derived from the discrete internal energy. The whole system is solved using the Hamiltonian splitting method discovered by He et al. [Phys. Plasmas 22, 124503 (2015)], which was been successfully adopted in constructing symplectic particle-in-cell schemes. Because of its structure preserving and explicit nature, this algorithm is especially suitable for large-scale simulations for physics problems that are multi-scale and require long-term fidelity and accuracy. The algorithm is verified via two tests: studies of the dispersion relation of waves in a two-fluid plasma system and the oscillating two-stream instability.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Efficient manifolds tracing for planar maps
Here, we present an exact calculation method and an approximation technique for tracing the invariant manifolds of unstable periodic orbits of planar maps. The exact approach relies on an adaptive refinement procedure that prevents redundant calculations occurring in other approaches, and the approximated method relies on a novel interpolation approach based on normal displacement functions. The resulting approximated manifold is precise when compared to the exact one, and its relative computational cost falls like the inverse of the manifold length. To present the tracing method, we obtain the invariant manifolds of the Chirikov-Taylor map, and as an application we illustrate the transition from homoclinic to heteroclinic chaos in the Duffing oscillator that leads from localized chaos to global chaotic motion.
74 ATOMIC AND MOLECULAR PHYSICS↗
Effects of pressure on the structure and lattice dynamics of ammonium perchlorate: A combined experimental and theoretical study
Abstract not provided
Kroonblawd, Matthew P.↗
A new method for analyzing and visualizing plasma simulations using a phase-space tessellation
We apply a novel phase-space interpolation technique referred to as the simplex-in-cell (SIC) method to analyze two- and three-dimensional particle-in-cell (PIC) simulations of electromagnetic plasmas. SIC relies on a discretization of the initial phase-space distribution function into simplices, which allows an approximation to the full, continuously defined distribution function to be constructed at any later time in the simulation. This allows densities, currents, and even full momentum distribution functions to be measured at any point in the simulation domain without averaging over control volumes. The SIC approach applies to any PIC simulation for which a tessellation of the initial particle distribution can be constructed. In this study, we use outputs from standard PIC simulations of the Weibel instability and compare physical quantities such as charge and current densities calculated in postprocessing using SIC and standard particle deposits. Using 2D simulations with 1–65 536 particles-per-cell, we find that SIC eliminates discrete particle noise and in some cases can reach a given noise level using ~1000 times fewer simulation particles than with standard particle deposition schemes. In regions of low density, such as between current filaments, SIC is able to capture small amplitude features even with fewer particles than gridpoints due to the deformable nature of the SIC volume elements. Here, by calculating momentum distributions, we show how SIC can capture low density tails in the spectrum using far fewer particles than are necessary for standard particle deposits. We calculate the charge density on spatial grids of increasing resolution to demonstrate the ability of SIC to reveal fine-scale details that are not accessible with standard particle deposits. Finally, we show how SIC can be extended to 3D and give an example of its use to calculate the charge density from 3D PIC simulations of the Weibel instability. These results motivate the future implementation of SIC directly in the simulation force calculation for a novel low-noise electromagnetic plasma simulation method.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Adaptive coupling of a deep neural network potential to a classical force field
An adaptive modeling method (AMM) that couples a deep neural network potential and a classical force field is introduced to address the accuracy-efficiency dilemma faced by the molecular simulation community. The AMM simulated system is decomposed into three types of regions. The first type captures the important phenomena in the system and requires high accuracy, for which we use the Deep Potential Molecular Dynamics (DeePMD) model in this work. The DeePMD model is trained to accurately reproduce the statistical properties of the ab initio molecular dynamics. The second type does not require high accuracy, and a classical force field is used to describe it in an efficient way. The third type is used for a smooth transition between the first and the second types of regions. By using a force interpolation scheme and imposing a thermodynamics force in the transition region, we make the DeePMD region embedded in the AMM simulated system as if it were embedded in a system that is fully described by the accurate potential. A representative example of the liquid water system is used to show the feasibility and promise of this method.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Overlapped embedded fragment stochastic density functional theory for covalently-bonded materials
The stochastic density functional theory (DFT) is a valuable linear-scaling approach to Kohn-Sham DFT that does not rely on the sparsity of the density matrix. Linear (and often sub-linear) scaling is achieved by introducing a controlled statistical error in the density, energy, and forces. The statistical error (noise) is proportional to the inverse square root of the number of stochastic orbitals and thus decreases slowly; however, by dividing the system into fragments that are embedded stochastically, the statistical error can be reduced significantly. This has been shown to provide remarkable results for non-covalently-bonded systems; however, the application to covalently bonded systems had limited success, particularly for delocalized electrons. Here, we show that the statistical error in the density correlates with both the density and the density matrix of the system and propose a new fragmentation scheme that elegantly interpolates between overlapped fragments. In conclusion, we assess the performance of the approach for bulk silicon of varying supercell sizes (up to N e = 16 384 electrons) and show that overlapped fragments reduce significantly the statistical noise even for systems with a delocalized density matrix.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Checking the Salpeter enhancement of nuclear reactions at intermediate coupling in asymmetric mixtures
In this study, we investigate the plasma enhancement of nuclear reactions in the intermediate coupling regime using orbital free molecular dynamics (OFMD) simulations. Mixtures of H-Cu and H-Ag serve as prototypes of simultaneous weak and strong couplings due to the charge asymmetry. Of particular importance is the partial ionization of Cu and Ag and the free electron polarization captured by OFMD simulations. By comparing a series of OFMD simulations at various concentrations and constant pressure to multi-component hyper-netted chain (MCHNC) calculations of effective binary ionic mixtures, we set a general procedure for computing enhancement factors. The MCHNC procedure allows extension to very low concentrations (5% or less) and to very high temperatures (few keV) unreachable by the simulations. Lastly, enhancement factors for nuclear reaction rates extracted from the MCHNC approach are compared with the Salpeter theory in the weak and strong coupling regimes, and a new interpolation is proposed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗