Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “configuration interaction”

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 199 records · Page 11

Ab initio Potential-Energy Surfaces and Electron-Spin-Exchange Cross Sections for H-O2 Interactions

Accurate quartet- and doublet-state potential-energy surfaces for the interaction of a hydrogen atom and an oxygen molecule in their ground states have been determined from an ab initio calculation using large-basis sets and the internally contracted multireference configuration interaction method. These potential surfaces have been used to calculate the H-O2 electron-spin-exchange cross section; the square root of the cross section (in a(sub 0)), not taking into account inelastic effects, can be obtained approximately from the expressions 2.390E(sup -1/6) and 5.266-0.708 log10(E) at low and high collision energies E (in E(sub h)), respectively. These functional forms, as well as the oscillatory structure of the cross section found at high energies, are expected from the nature of the interaction energy. The mean cross section (the cross section averaged over a Maxwellian velocity distribution) agrees reasonably well with the results of measurements.

Stallcop, James R.↗

Unusual behavior of Cooper minima of n s subshells in high- Z atoms

Here, a study of Cooper minima (CM) arising from the photoionization of 6$\textit{s}$, 5$\textit{s}$, and 4$\textit{s}$ subshells of high-$\textit{Z}$ atoms has been performed using Dirac-Fock (DF), two-channel relativistic-random-phase approximation (RRPA), and fully coupled RRPA. The results show huge splittings between $ns → εp_{3/2}$ and $ns → εp_{1/2}$ CM which increase with $\textit{Z}$ owing primarily to the relativistic interactions (spin orbit) that are attractive for the $εp_{1/2}$ final state but repulsive for the corresponding $εp_{3/2}$. In addition, it was found that correlation in the form of interchannel coupling (essentially configuration interaction in the final continuum states) plays a huge role in determining the location of the CM. For 6$\textit{s}$ photoionization, the 6$s → εp_{3/2}$ and 6$s → εp_{1/2}$ CM behave completely different as a function of $\textit{Z}$ . It was also found that for 5$\textit{s}$ and 4$\textit{s}$ photoionization, the CM move below the threshold, with increasing $\textit{Z}$ , and, at high enough $\textit{Z}$ , the 5s → εp3/2 and 4s → εp3/2 CM re-emerge into the continuum. The calculations have been carried out for the $\textit{ns}$ subshells of Hg ($\textit{Z}$ = 80), Rn ($\textit{Z}$ = 86), Ra ($\textit{Z}$ = 88), No ($\textit{Z}$ = 102), Cn ($\textit{Z}$ = 112), and Og ($\textit{Z}$ = 118).

74 ATOMIC AND MOLECULAR PHYSICS↗

System technology analysis of aeroassisted orbital transfer vehicles: Moderate lift/drag (0.75-1.5), Volume 1B, part 2, study results

A complete compilation of the results from Phase 2 of a study to identify and prioritize the technology payoffs of representative space based mid lift/drag ratio (L/D) aeroassisted orbit transfer vehicles (AOTV) and the cryofueled propulsion subsystem - configuration interactions is presented. Several combinations of basing and launch vehicle options, staging scenarios, missions (delivery, servicing, or manned round trip), and target orbits were considered. Space basing of an AOTV opens up numerous configuration opportunities. The size can exceed the launch vehicle cargo bay envelope by resorting to orbital assembly and the AOTV weight can exceed the launch vehicle capability. With the absence of fully fueled tanks during a ground based launch much lighter gossamer type structures are possible that may result in performance gains. At the Space Station, payload rearranging or manifesting may prove attractive. Several major conclusions regarding aerothermodynamics, aerodynamics, payload manifesting, the propulsion subsystem, and systems issues are discussed.

Source record↗

Diabatic Valence-Hole States in the C 2 Molecule: “Putting Humpty Dumpty Together Again”

Despite the long history of spectroscopic studies of the C$_2$ molecule, fundamental questions about its chemical bonding are still being hotly debated. The complex electronic structure of C$_2$ is a consequence of its dense manifold of near-degenerate, low-lying electronic states. A global multi-state diabatic model is proposed here to disentangle the numerous configuration interactions within four symmetry manifolds of C$_2$ ($^{1}\Pi_g$, $^{3}\Pi_g$, $^{1}\Sigma_u^+$, and $^{3}\Sigma_u^+$). The key concept of our model is the existence of two "valence-hole" configurations, $2\sigma_g^22\sigma_u^11\pi_{u}^33\sigma_g^2$ for $^{1,3}\Pi_g$ states and $2\sigma_g^22\sigma_u^11\pi_{u}^43\sigma_g^1$ for $^{1,3}\Sigma_u^+$ states that derive from $3\sigma_g\leftarrow2\sigma_u$ electron promotion. The lowest-energy state from each of the four C$_2$ symmetry species is dominated by this type of valence-hole configuration at its equilibrium internuclear separation. Finally, as a result of their large binding energy (nominal bond order of 3) and correlation with the 2s$^2$2p$^2$+2s2p$^3$ separated-atom configurations, the presence of these valence-hole configurations has a profound impact on the $global$ electronic structure and unimolecular dynamics of C$_2$.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗

Barrier penetration in a discrete-basis formalism

The dynamics of a many-particle system are often modeled by mapping the Hamiltonian into a Schrödinger equation. An alternative approach is to solve the Hamiltonian equations directly in a model space of many-body configurations. In a previous paper the numerical convergence of the two approaches was compared with a simplified treatment of the Hamiltonian representation. Here, we extend the comparison to the nonorthogonal model spaces that would be obtained by the generator-coordinate method. With a suitable choice of the collective-variable grid, a configuration-interaction Hamiltonian can reproduce the Schrödinger dynamics very well. Furthermore, the method as implemented here requires that the barrier height is not much larger than the zero-point energy in the collective coordinates of the configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Accelerating eigenvalue computation for nuclear structure calculations via perturbative corrections

Subspace projection methods utilizing perturbative corrections have been proposed for computing the lowest few eigenvalues and corresponding eigenvectors of large Hamiltonian matrices. In this paper, we build upon these methods and introduce the term Subspace Projection with Perturbative Corrections (SPPC) method to refer to this approach. We tailor the SPPC for nuclear many-body Hamiltonians represented in a truncated configuration interaction subspace, i.e., the no-core shell model (NCSM). We use the hierarchical structure of the NCSM Hamiltonian to partition the Hamiltonian as the sum of two matrices. The first matrix corresponds to the Hamiltonian represented in a small configuration space, whereas the second is viewed as the perturbation to the first matrix. Eigenvalues and eigenvectors of the first matrix can be computed efficiently. Because of the split, perturbative corrections to the eigenvectors of the first matrix can be obtained efficiently from the solutions of a sequence of linear systems of equations defined in the small configuration space. These correction vectors can be combined with the approximate eigenvectors of the first matrix to construct a subspace from which more accurate approximations of the desired eigenpairs can be obtained. We show by numerical examples that the SPPC method can be more efficient than conventional iterative methods for solving large-scale eigenvalue problems such as the Lanczos, block Lanczos and the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The method can also be combined with other methods to avoid convergence stagnation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of the excitation spectra with diffusion Monte Carlo on an auxiliary bosonic ground state

We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system’s ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion–boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Origin of Ferromagnetic Exchange Coupling in Donor–Acceptor Biradical Analogues of Charge-Separated Excited States

A new donor–acceptor biradical complex, Tp Cum,Me Zn(SQ-VD) (Tp Cum,Me Zn + = zinc(II) hydro-tris(3-cumenyl-5-methylpyrazolyl)borate complex cation; SQ = orthosemiquinone; VD = oxoverdazyl), which is a ground-state analogue of a charge-separated excited state, has been synthesized and structurally characterized. Here, the magnetic exchange interaction between the S = 1/2 SQ and the S = 1/2 VD within the SQ-VD biradical ligand is observed to be ferromagnetic, with $J_{SQ-VD}$ = +77 cm –1 (H = -2$J_{SQ-VD}\bar{S}_{SQ}·\bar{S}_{VD}$) determined from an analysis of the variable-temperature magnetic susceptibility data. The pairwise biradical exchange interaction in Tp Cum,Me Zn(SQ-VD) can be compared with that of the related donor–acceptor biradical complex Tp Cum,Me Zn(SQ-NN) (NN = nitronyl nitroxide, S = 1/2), where $J_{SQ-NN}$ ≅ +550 cm –1 . This represents a dramatic reduction in the biradical exchange by a factor of ~7, despite the isolobal nature of the VD and NN acceptor radical SOMOs. Computations assessing the magnitude of the exchange were performed using a broken-symmetry density functional theory (DFT) approach. These computations are in good agreement with those computed at the CASSCF NEVPT2 level, which also reveals an S = 1 triplet ground state as observed in the magnetic susceptibility measurements. A combination of electronic absorption spectroscopy and CASSCF computations has been used to elucidate the electronic origin of the large difference in the magnitude of the biradical exchange coupling between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN). A Valence Bond Configuration Interaction (VBCI) model was previously employed to highlight the importance of mixing an SQ SOMO → NNLUMO charge transfer configuration into the electronic ground state to facilitate the stabilization of the high-spin triplet (S = 1) ground state in Tp Cum,Me Zn(SQ-NN). Here, CASSCF computations confirm the importance of mixing the pendant radical (e.g., VD, NN) LUMO (VD LUMO and NN LUMO ) with the SOMO of the SQ radical (SQ SOMO ) for stabilizing the triplet, in addition to spin polarization and charge transfer contributions to the exchange. An important electronic structure difference between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN), which leads to their different exchange couplings, is the reduced admixture of excited states that promote ferromagnetic exchange into the Tp Cum,Me Zn(SQ-VD) ground state, and the intrinsically weaker mixing between the VD LUMO and the SQ SOMO compared to that observed for Tp Cum,Me Zn(SQ-NN), where this orbital mixing is significant. The results of this comparative study contribute to a greater understanding of biradical exchange interactions, which are important to our understanding of excited-state singlet–triplet energy gaps, electron delocalization, and the generation of electron spin polarization in both the ground and excited states of (bpy)Pt(CAT-radical) complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Study of the Electronic Spectra of a Polycyclic Aromatic Hydrocarbon, Naphthalene, and its Derivatives

In order to preselect possible candidates for the origin of diffuse interstellar bands observed, semiempirical quantum mechanical method INDO/S was applied to the optical spectra of neutral, cationic, and anionic states of naphthalene and its hydrogen abstraction and addition derivatives. Comparison with experiment shows that the spectra of naphthalene and its ions were reliably predicted. The configuration interaction calculations with single-electron excitations provided reasonable excited state wavefunctions compared to ab initio calculations that included higher excitations. The degree of similarity of the predicted spectra of the hydrogen abstraction and derivatives to those of naphthalene and ions depends largely on the similarity of the it electron configurations. For the hydrogen addition derivatives, very little resemblance of the predicted spectra to naphthalene was found because of the disruption of the aromatic conjugation system. The relevance of these calculations to astrophysical issues is discussed within the context of these polycyclic aromatic hydrocarbon models. Comparing the calculated electronic energies to the Diffuse Interstellar Bands (DIBs), a list of possible candidates of naphthalene derivatives is established which provides selected candidates for a definitive test through laboratory studies.

Du, Ping↗

Theoretical calculations on the electron absorption spectra of selected Polycyclic Aromatic Hydrocarbons (PAH) and derivatives

As a theoretical component of the joint effort with the laboratory of Dr. Lou Allamandola to search for potential candidates for interstellar organic carbon compound that are responsible for the visible diffuse interstellar absorption bands (DIB's), quantum mechanical calculations were performed on the electron absorption spectra of selected polycyclic aromatic hydrocarbons (PAH) and derivatives. In the completed project, 15 different species of naphthalene, its hydrogen abstraction and addition derivatives, and corresponding cations and anions were studied. Using semiempirical quantum mechanical method INDO/S, the ground electronic state of each species was evaluated with restricted Hartree-Fock scheme and limited configuration interaction. The lowest energy spin state for each species was used for electron absorption calculations. Results indicate that these calculations are accurate enough to reproduce the spectra of naphthalene cation and anion observed in neon matrix. The spectral pattern of the hydrogen abstraction and addition derivatives predicted based on these results indicate that the electron configuration of the pi orbitals of these species is the dominant determinant. A combined list of 19 absorptions calculated from 4500 A to 10,400 A were compiled and suggested as potential candidates that are relevant for the DIB's absorptions. Continued studies on pyrene and derivatives revealed the ground state symmetries and multiplicities of its neutral, anionic, and cationic species. Spectral calculations show that the cation (B(sub 3g)-2) and the anion (A(sub u)-2) are more likely to have low energy absorptions in the regions between 10 kK and 20 kK, similar to naphthalene. These absorptions, together with those to be determined from the hydrogen abstraction and addition derivatives of pyrene, can be used to provide additional candidates and suggest experimental work in the search for interstellar compounds that are responsible for DIB's.

Du, Ping↗

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗

Studies of aerothermal loads generated in regions of shock/shock interaction in hypersonic flow

Experimental studies were conducted to examine the aerothermal characteristics of shock/shock/boundary layer interaction regions generated by single and multiple incident shocks. The presented experimental studies were conducted over a Mach number range from 6 to 19 for a range of Reynolds numbers to obtain both laminar and turbulent interaction regions. Detailed heat transfer and pressure measurements were made for a range of interaction types and incident shock strengths over a transverse cylinder, with emphasis on the 3 and 4 type interaction regions. The measurements were compared with the simple Edney, Keyes, and Hains models for a range of interaction configurations and freestream conditions. The complex flowfields and aerothermal loads generated by multiple-shock impingement, while not generating as large peak loads, provide important test cases for code prediction. The detailed heat transfer and pressure measurements proved a good basis for evaluating the accuracy of simple prediction methods and detailed numerical solutions for laminar and transitional regions or shock/shock interactions.

Holden, Michael S.↗

I ntera C hem : Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree–Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The dissociation energy of N2

The requirements for very accurate ab initio quantum chemical prediction of dissociation energies are examined using a detailed investigation of the nitrogen molecule. Although agreement with experiment to within 1 kcal/mol is not achieved even with the most elaborate multireference CI (configuration interaction) wave functions and largest basis sets currently feasible, it is possible to obtain agreement to within about 2 kcal/mol, or 1 percent of the dissociation energy. At this level it is necessary to account for core-valence correlation effects and to include up to h-type functions in the basis. The effect of i-type functions, the use of different reference configuration spaces, and basis set superposition error were also investigated. After discussing these results, the remaining sources of error in our best calculations are examined.

Almloef, Jan↗

The dissociation energy of N2

The requirements for very accurate ab initio quantum chemical prediction of dissociation energies are examined using a detailed investigation of the nitrogen molecule. Although agreement with experiment to within 1 kcal/mol is not achieved even with the most elaborate multireference CI (configuration interaction) wave functions and largest basis sets currently feasible, it is possible to obtain agreement to within about 2 kcal/mol, or 1 percent of the dissociation energy. At this level it is necessary to account for core-valence correlation effects and to include up to h-type functions in the basis. The effect of i-type functions, the use of different reference configuration spaces, and basis set superposition error were also investigated. After discussing these results, the remaining sources of error in our best calculations are examined.

Almlof, Jan↗

The Interactions Between Nitrogen and Oxygen Molecules

Lippincott's delta-function model for atomic interactions is analyzed, both physically and mathematically, and extended, by differentiation between K- and L-shell electrons and the introduction of a variable parameter in the expression for the delta-function strength, to cover homonuclear molecules more complex than hydrogen. In addition, modifications are made which allow treatments of diatomic, heteronuclear molecules. This theory, in conjunction with a reasonably extensive study of resonance, dispersion, and configuration interaction phenomena, as well as the use of simple quantum mechanical arguments, is then applied to the N2-N2, N2-O2, and O2-O2 interactions.

Meador, Willard E., Jr.↗