Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 577 records · Page 32

Detailed Characterization of a Nanosecond-Lived Excited State: X-ray and Theoretical Investigation of the Quintet State in Photoexcited [Fe(terpy) 2 ] 2+

Theoretical predictions show that depending on the populations of the Fe 3dxy, 3dxz, and 3dyz orbitals two possible quintet states can exist for the high-spin state of the photoswitchable model system [Fe(terpy)2]2+. The differences in the structure and molecular properties of these 5B2 and 5E quintets are very small and pose a substantial challenge for experiments to resolve them. Yet for a better understanding of the physics of this system, which can lead to the design of novel molecules with enhanced photoswitching performance, it is vital to determine which high-spin state is reached in the transitions that follow the light excitation. The quintet state can be prepared with a short laser pulse and can be studied with cutting-edge time-resolved X-ray techniques. Here we report on the application of an extended set of X-ray spectroscopy and scattering techniques applied to investigate the quintet state of [Fe(terpy)2]2+ 80 ps after light excitation. High-quality X-ray absorption, nonresonant emission, and resonant emission spectra as well as X-ray diffuse scattering data clearly reflect the formation of the high-spin state of the [Fe(terpy)2]2+ molecule; moreover, extended X-ray absorption fine structure spectroscopy resolves the Fe-ligand bond-length variations with unprecedented bondlength accuracy in time-resolved experiments. With ab initio calculations we determine why, in contrast to most related systems, one configurational mode is insufficient for the description of the low-spin (LS)-high-spin (HS) transition. We identify the electronic structure origin of the differences between the two possible quintet modes, and finally, we unambiguously identify the formed quintet state as 5E, in agreement with our theoretical expectations.

Vankó, György↗

Testing isospin symmetry breaking in ab initio nuclear theory

In this work we present the first steps towards benchmarking isospin symmetry breaking in ab initio nuclear theory for calculations of superallowed Fermi β decay. Using the valence-space in-medium similarity renormalization group, we calculate b and c coefficients of the isobaric multiplet mass equation, starting from two different Hamiltonians constructed from chiral effective field theory. We compare results to experimental measurements for all T = 1 isobaric analog triplets of relevance to superallowed β decay for masses A = 10 to A = 74 and find an overall agreement within approximately 250 keV of experimental data for both b and c coefficients. A greater level of accuracy, however, is obtained by a phenomenological Skyrme interaction or a classical charged-sphere estimate. Lastly, we show that evolution of the valence-space operator does not meaningfully improve the quality of the coefficients with respect to experimental data, which indicates that higher-order many-body effects are likely not responsible for the observed discrepancies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

First principles electronic band structure and phonon dispersion curves for zinc blend beryllium chalcogenide

A detailed theoretical study of structural, electronic and Vibrational properties of BeX compound is presented by performing ab-initio calculations based on density-functional theory using the Espresso package. The calculated value of lattice constant and bulk modulus are compared with the available experimental and other theoretical data and agree reasonably well. BeX (X = S,Se,Te) compounds in the ZB phase are indirect wide band gap semiconductors with an ionic contribution. The phonon dispersion curves are represented which shows that these compounds are dynamically stable in ZB phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Precision Mass Measurements of Neutron-Rich Scandium Isotopes Refine the Evolution of N = 32 and N = 34 Shell Closures

In this work, we report high-precision mass measurements of 50–55 Sc isotopes performed at the LEBIT facility at NSCL and at the TITAN facility at TRIUMF. Our results provide a substantial reduction of their uncertainties and indicate significant deviations, up to 0.7 MeV, from the previously recommended mass values for 53–55 Sc. The results of this work provide an important update to the description of emerging closed-shell phenomena at neutron numbers N = 32 and N = 34 above proton-magic Z = 20. In particular, they finally enable a complete and precise characterization of the trends in ground state binding energies along the N = 32 isotone, confirming that the empirical neutron shell gap energies peak at the doubly magic 52 Ca. Moreover, our data, combined with other recent measurements, do not support the existence of a closed neutron shell in 55 Sc at N = 34. The results were compared to predictions from both ab initio and phenomenological nuclear theories, which all had success describing N = 32 neutron shell gap energies but were highly disparate in the description of the N = 34 isotone.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Atomic structures determined from digitally defined nanocrystalline regions

Nanocrystallography has transformed our ability to interrogate the atomic structures of proteins, peptides, organic molecules and materials. By probing atomic level details in ordered sub-10 nm regions of nanocrystals, scanning nanobeam electron diffraction extends the reach of nanocrystallography and in principle obviates the need for diffraction from large portions of one or more crystals. Scanning nanobeam electron diffraction is now applied to determine atomic structures from digitally defined regions of beam-sensitive peptide nanocrystals. Using a direct electron detector, thousands of sparse diffraction patterns over multiple orientations of a given crystal are recorded. Each pattern is assigned to a specific location on a single nanocrystal with axial, lateral and angular coordinates. This approach yields a collection of patterns that represent a tilt series across an angular wedge of reciprocal space: a scanning nanobeam diffraction tomogram. Using this diffraction tomogram, intensities can be digitally extracted from any desired region of a scan in real or diffraction space, exclusive of all other scanned points. Intensities from multiple regions of a crystal or from multiple crystals can be merged to increase data completeness and mitigate missing wedges. It is demonstrated that merged intensities from digitally defined regions of two crystals of a segment from the OsPYL/RCAR5 protein produce fragment-based ab initio solutions that can be refined to atomic resolution, analogous to structures determined by selected-area electron diffraction. In allowing atomic structures to now be determined from digitally outlined regions of a nanocrystal, scanning nanobeam diffraction tomography breaks new ground in nanocrystallography.

36 MATERIALS SCIENCE↗

A Conjecture on Deducing General Relativity and the Standard Model with Its Fundamental Constants from Rational Tangles of Strands

It appears possible to deduce black holes, general relativity and the standard model of elementary particles from one-dimensional strands that fluctuate at the Planck scale. This appears possible as long as only switches of skew strand crossings are observable, but not the strands themselves. Woven fluctuating strands behave like horizons and imply black hole entropy, the field equations of general relativity and cosmological observations. Tangled fluctuating strands in flat space imply Dirac’s equation. The possible families of unknotted rational tangles produce the spectrum of elementary particles. Fluctuating rational tangles also yield the gauge groups U(1), broken SU(2), and SU(3), produce all Feynman diagrams of the standard model, and exclude any unknown elementary particle, gauge group, and Feynman diagram. The conjecture agrees with all known experimental data. Predictions for experiments arise, and the fundamental constants of the standard model can be calculated. Objections are discussed. Predictions and calculations allow testing the conjecture. As an example, an ab initio estimate of the fine structure constant is outlined.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Energy Landscape, and Structural and Spectroscopic Characterization of Diazirine and Its Cyclic Isomers

Identifying new nitrogenated hydrocarbon molecules in the interstellar medium (ISM) is challenging due to the lack of comprehensive spectroscopic data from experiments. In this computational work, we focus on investigating the structures, relative energies, spectroscopic constants, and energy landscape of the cyclic isomers of diazirine (c-CH2N2) using ab initio quantum chemical methods. Density functional theory (DFT) methods and coupled cluster theory with singles and doubles including perturbative triples [CCSD(T)] and CCSD(T) with the explicitly correlated F12 b correction [CCSD(T)-F12b] were employed for this purpose along with large correlation consistent cc-pVTZ, cc-pVQZ, and cc-pV5Z basis sets. Harmonic vibrational frequencies, infrared vibrational intensities, rotational constants, and dipole moments are reported. An harmonic vibrational fundamentals along with centrifugal distortion constants, and vibration-rotation interaction constants are also reported for all the cyclic isomers. The energies computed with the CCSD(T) and CCSD(T)-F12b methods were extrapolated to the one-particle complete basis set (CBS) limit following a three-point formula. At the CCSD(T)-F12b/CBS level of theory the 3,3H-diazirine (c-CH2N2) is the lowest energy cyclic isomer followed by 1,3H-diazirine, (E)-1,2H-Diazirine, and (Z)-1,2H-diazirine which are 20.1, 47.8, and 51.3kcal mol-1above the 3,3H-Diazirine, respectively. Accurate structures and spectroscopic constants that are reported here could be useful for future identification of these cyclic nitrogenated organic molecules in the interstellar medium or circumstellar disks.

Energy Landscape↗

Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature

It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. Here, we find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗

Modeling inclusive electron-nucleus scattering with Bayesian artificial neural networks

We introduce a Bayesian protocol based on artificial neural networks that is suitable for modeling inclusive electron-nucleus scattering on a variety of nuclear targets with quantified uncertainties. Unlike previous applications in the field, which directly parameterize the cross sections, our approach employs artificial neural networks to represent the longitudinal and transverse response functions. In contrast to cross sections, which depend on the incoming energy, scattering angle, and energy transfer, the response functions are determined solely by the energy and momentum transfer to the system, allowing the angular component to be treated analytically. We assess the accuracy and predictive power of our framework against the extensive data in the quasielastic inclusive electron-scattering database. Additionally, we present novel extractions of the longitudinal and transverse response functions and compare them with previous experimental analysis and nuclear ab-initio calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dispersive optical model analysis of 208 Pb generating a neutron-skin prediction beyond the mean field

In this work, a nonlocal dispersive optical model analysis is carried out for neutrons and protons in 208 Pb. Elastic-scattering angular distributions, total and reaction cross sections, single-particle energies, neutron and proton numbers, the charge distribution, and the binding energy are fitted to extract the neutron and proton self-energies both above and below the Fermi energy. From the single-particle propagator derived from these self-energies, we determine the charge and matter distributions in 208 Pb. The predicted spectroscopic factors are consistent with results from the ($\textit{e,e}'\textit{p}$) reaction and inelastic-electron-scattering data to very high-spin states. Sensible results for the high-momentum content of neutrons and protons are obtained, with protons appearing more correlated, in agreement with experiment and ab initio calculations of asymmetric matter. A neutron skin of 0.25 ± 0.05 fm is deduced. An analysis of several nuclei leads to the conclusion that finite-size effects play a nonnegligible role in the formation of the neutron skin in finite nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic Anisotropy of Graphene-Coated Thin Iron Films

The ab initio calculations of the magnetic anisotropy of thin (to six monolayers) nickel films and iron films with unusual fcc structures formed upon intercalation of graphene with Fe atoms are presented. The data have been obtained for both the pure-surface films and the films coated with graphene. The density functional theory and the pseudopotential method are used to calculate the magnetic moments of atoms of all the layers and to determine the total energies of the structures with different orientations of magnetic moments of iron and nickel atoms with respect to the film surface. A strong influence of graphene on the magnetic properties contacting iron films is demonstrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Studying electroweak few-body observables in chiral effective field theory

The use of nuclei to study electroweak probes is becoming increasingly relevant experimentally. The success of dark matter and neutrino experiments strongly depends on the ability to control nuclear effects in order to extract the fundamental parameters associated with external probes. Therefore, reliable theoretical calculations of nuclear structure and reactions, with well-controlled errors, are crucial for the success of experimental efforts. Currently, chiral effective field theory ($\chi$EFT) coupled with {\it ab-initio} methods represents one of the best approaches that fulfills these requirements. To use this approach as a tool for studying fundamental physics, it is essential to validate it against experimental data for which the calculations are well under control, such as the elastic scattering of electrons on nuclei. In this proceeding, I will present recent developments in the fitting of electromagnetic currents derived using $\chi$EFT and the calculation of electromagnetic form factors of light nuclei. The results of these calculations demonstrate the strength of the theory in describing the interaction of nuclei with electromagnetic probes over a broad range of momentum transfers and highlight the robustness of $\chi$EFT for analyzing future experimental data aimed at extracting fundamental parameters.

Gnech, Alex [Old Dominion Univ., Norfolk, VA (Unit↗

Experimental energy levels and partition function of the {sup 12}C{sub 2} molecule

The carbon dimer, the {sup 12}C{sub 2} molecule, is ubiquitous in astronomical environments. Experimental-quality rovibronic energy levels are reported for {sup 12}C{sub 2}, based on rovibronic transitions measured for and among its singlet, triplet, and quintet electronic states, reported in 42 publications. The determination utilizes the Measured Active Rotational-Vibrational Energy Levels (MARVEL) technique. The 23,343 transitions measured experimentally and validated within this study determine 5699 rovibronic energy levels, 1325, 4309, and 65 levels for the singlet, triplet, and quintet states investigated, respectively. The MARVEL analysis provides rovibronic energies for six singlet, six triplet, and two quintet electronic states. For example, the lowest measurable energy level of the a{sup 3}Π{sub u} state, corresponding to the J = 2 total angular momentum quantum number and the F {sub 1} spin-multiplet component, is 603.817(5) cm{sup −1}. This well-determined energy difference should facilitate observations of singlet–triplet intercombination lines, which are thought to occur in the interstellar medium and comets. The large number of highly accurate and clearly labeled transitions that can be derived by combining MARVEL energy levels with computed temperature-dependent intensities should help a number of astrophysical observations as well as corresponding laboratory measurements. The experimental rovibronic energy levels, augmented, where needed, with ab initio variational ones based on empirically adjusted and spin–orbit coupled potential energy curves obtained using the Duo code, are used to obtain a highly accurate partition function, and related thermodynamic data, for {sup 12}C{sub 2} up to 4000 K.

79 ASTRONOMY AND ASTROPHYSICS↗

Theoretical calculation of low-lying states of NaAr and NaXe

Potential curves as well as dipole moments and linking transition moments are calculated for the ground X 2 Sigma + and low lying excited A 2 Pi, B 2 Sigma +, C 2 Sigma +, (4) 2 Sigma +, (2) 2 Pi and (1) 2 Delta states of NaAr and NaXe. Calculations are performed using a self-consistent field plus configuration-interaction procedure with the core electrons replaced by an ab initio effective core potential. The potential curves obtained are found to be considerably less repulsive than the semiempirical curves of Pascale and Vandeplanque (1974) and to agree well with existing experimental data, although the binding energies of those states having potential minima due to van der Waals interactions are underestimated. Emission bands are also calculated for the X 2 Sigma + - C 2 Sigma + excimer transitions of NaAr and NaXe using the calculated transition moments and potential curves, and shown to agree well with experiment on the short-wavelength side of the maximum.

Laskowski, B. C.↗

Theory of collision-induced translation-rotation spectra: H2-He

An adiabatic quantal theory of spectral line shapes in collision-induced absorption and emission is presented which incorporates the induced translation-rotation and translation-vibration spectra. The generalization to account for the anisotropy of the scattering potential is given. Calculations are carried out of the collision-induced absorption spectra of He in collisions with H2 with ab initio electric dipole functions and realistic potentials. The anisotropy of the interaction potential is small and is not included in the calculations. The predicted spectra are in satisfactory agreement with experimental data though some deviations occur which may be significant. The rotational line shapes have exponential wings and are not Lorentzian. The connection between the quantal and classical theories is written out explicitly for the isotropic overlap induction.

Birnbaum, G.↗

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗

Massless and massive photons within basis light-front quantization

We compute the structure function, transverse momentum dependent parton distributions (TMDs), and generalized parton distributions (GPDs) for the photon from the light-front quantum electrodynamics (QED) Hamiltonian, determined for its constituent bare photon and electron-positron Fock components within the basis light-front quantization framework. While the structure function and GPDs align well with expectations, some oscillatory behavior in the TMDs emerges as a finite-basis artifact. We also compute the GPD and the corresponding structure function for the longitudinal polarization mode of the massive photon. We further study the photon mass dependence of the calculated observables. After performing nonperturbative renormalization, we obtain a good quality description of the perturbative QED properties of the photon, as well as good agreement with the experimental data for the massless photon’s structure function.

Ab initio calculations↗