Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic physics”

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,693 records · Page 94

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗

Toward Full Configuration Interaction for Transition-Metal Complexes

In this work, an efficacious approximation to full configuration interaction (FCI) is adapted to calculate singlet-triplet gaps for transition-metal complexes. This strategy, incremental FCI (iFCI), uses a many-body expansion to systematically add correlation to a simple reference wave function and therefore achieves greatly reduced computational costs compared to FCI. iFCI through the 3-body expansion is demonstrated on four model transition-metal complexes involving the metals Zn, V, and Cu. Screening techniques to increase the computational efficiency of iFCI are proposed and tested, showing reduction in the number of 3-body terms by more than 90% with controlled errors. The largest complex treated by iFCI has 142 valence electrons, all of which are correlated among the full set of 444 active orbitals. Computed spin gaps approach experimental results for the four complexes, though room for improvement remains.

74 ATOMIC AND MOLECULAR PHYSICS↗

Threshold Ionization Spectroscopy and Theoretical Calculations of LnO (Ln = La and Ce)

Mass-analyzed threshold ionization (MATI) spectroscopy was used to measure the vibronic spectra of LnO (Ln = La and Ce). Single-reference coupled cluster and relativistic multireference configuration calculations were carried out to compare with the measured vibronic energies. The spectrum of LaO displays a single vibronic band system, while that of CeO shows multiple ones. The ionization energies of LaO and CeO are measured as 5.2446 (6) and 5.3332(6) eV, respectively, which are a hundredfold improvement over the literature values. The vibrational energies of the neutral molecule and corresponding ion reveal the charge effect on the metal-oxygen bond of both species. The single band system in the spectrum of LaO arises from the transition of the ground state of the neutral molecule with the La(6s 1 )O(2p 6 ) valence configuration to the ground state of the singly charged ion with the La(6s o )O(2p 6 ) configuration. Here, the multiple band systems in the spectrum of CeO are attributed to the spin-orbit coupling for the Ce(4f 1 6s 1 )O(2p 6 ) configuration of the neutral molecule and an excited state for the Ce(4f 1 )O(2p 6 ) configuration of the ion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Neat, Simple, and Wrong: Debunking Electrostatic Fallacies Regarding Noncovalent Interactions

Multipole moments such as charge, dipole, and quadrupole are often invoked to rationalize intermolecular phenomena, but a low-order multipole expansion is rarely a valid description of electrostatics at the length scales that characterize nonbonded interactions. This is illustrated by examining several common misunderstandings rooted in erroneous electrostatic arguments. First, the notion that steric repulsion originates in Coulomb interactions is easily disproved by dissecting the interaction potential for Ar2. Second, the Hunter-Sanders model of π–π interactions, which is based on quadrupolar electrostatics, is shown to have no basis in accurate calculations. Third, curved “buckybowls” exhibit unusually large dipole moments, but these are ancillary to the forces that control their intermolecular interactions, as illustrated by two examples involving corannulene. Finally, the assumption that interactions between water and small anions are dictated by the dipole moment of H 2 O is shown to be false in the case of binary halide–water complexes. Furthermore, these examples present a compelling case that electrostatic explanations based on low-order multipole moments are very often counterfactual for nonbonded interactions at close range and should not be taken seriously in the absence of additional justification.

74 ATOMIC AND MOLECULAR PHYSICS↗

Kohn–Sham Density in a Slater Orbital Basis Set

Finite, atom-centered Slater basis sets are used to determine approximate Kohn–Sham molecular orbitals. This is achieved by minimizing the kinetic energy plus the sum-squared difference between the Kohn–Sham density and the full configuration interaction density. As a result of the finite basis, a weight factor is introduced to balance the two minimization components. Results herein show that this can be done systematically, without sensitive dependence on the choice of scaling factor. In addition, the algorithm is applied to the LiH diatomic for fractional electron counts, where stretching the bond introduces significant reorganization of the electron density. As a result, the analysis will show the correct KS orbital structure and reveal the effects of correlation and electron locality on the KS solutions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗

Accelerating Ab Initio Simulation via Nested Monte Carlo and Machine Learned Reference Potentials

As a corollary of the rapid advances in computing, ab initio simulation is playing an increasingly important role in modeling materials at the atomic scale. Two strategies are possible, ab initio Monte Carlo (AIMC) and molecular dynamics (AIMD) simulation. The former benefits from exact sampling from the correct thermodynamic distribution, while the latter is typically more efficient with its collective all-atom coordinate updates. In this study, using a relatively simple test model comprised of helium and argon, we show that AIMC can be brought up to, and even above, the performance levels of AIMD via a hybrid nested sampling/machine learning (ML) strategy. Here, ML provides an accurate classical reference potential (up to three-body explicit interactions) that can pilot long collective Monte Carlo moves that are accepted or rejected in toto à la nested Monte Carlo (NMC); this is in contrast to the single move nature of a naive implementation. Our proposed method only requires a small up front expense from evaluating the ab initio energies and forces of (100) random configurations for training. Importantly, our method does not totally rely on the trained, assuredly imperfect, interaction. We show that high performance and exact sampling at the desired level of theory can be realized even when the trained interaction has appreciable differences from the ab initio potential. Remarkably, at the highest levels of performance realized via our approach, a pair of statistically uncorrelated atomic configurations can be generated with (1) ab initio calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Uncovering the Functional Role of Coherent Phonons during the Photoinduced Phase Transition in a Molecular Crystal

The atomic motions that make up phonons and molecular vibrations in molecular crystals influence their photophysical and electronic properties, including polaron formation, carrier mobility, and phase transitions. Discriminating between spectator and driving motions is a significant challenge hindering optimization. Unlocking this information and developing fine-tuned controls over actively participating phonon modes would not only lead to a stronger understanding of photochemistry but also provide a significant new tool in controlling solid state chemistry. We present a strategy using rationally designed double pulses to unveil the unique function of specific excited state phonon modes. Using ultrafast spectroscopy, we identified 50 and 90 cm –1 phonons involved in modulating the photoinduced spin-Peierls melting of potassium tetracyanoquinodimethane crystals. Here, we show that the 50 cm –1 phonon specifically corresponds to the coherent nuclear wavepacket involved in the charge transfer component of the overall spin-Peierls phase melting process, while the 90 cm –1 phonon facilitates the phase transition component.

74 ATOMIC AND MOLECULAR PHYSICS↗

Vibronic Excitons and Conical Intersections in Semiconductor Quantum Dots

Surface defects and organic surface-capping ligands affect the photoluminescence properties of semiconductor quantum dots (QDs) by altering the rates of competing nonradiative relaxation processes. In this study, broadband two-dimensional electronic spectroscopy reveals that absorption of light by QDs prepares vibronic excitons, excited states derived from quantum coherent mixing of the core electronic and ligand vibrational states. Rapidly damped coherent wavepacket motions of the ligands are observed during hot-carrier cooling, with vibronic coherence transferred to the photoluminescent state. Furthermore, these findings suggest a many-electron, molecular theory for the electronic structure of QDs, which is supported by calculations of the structures of conical intersections between the exciton potential surfaces of a small ammonia-passivated model CdSe nanoparticle.

36 MATERIALS SCIENCE↗

Electronic Configuration Assignments for UO from Electric Dipole Moment Measurements

Diatomic UO has more than 48 bound states within 10000 cm -1 of the ground state. This electronic state congestion has been attributed to interleaved states from the electronic configurations U 2+ ($5f^37s$)O 2- and U 2+ ($5f^27s^2$)O 2- , respectively. Ligand field theory predicts that each electronic configuration will exhibit states with distinguishable, characteristic vibrational and rotational constants. However, vibronic state mixing modifies the observed vibration–rotation constants, leading to uncertainty in the configurational assignments. The permanent electric dipole moment ($μ_e$) of an electronic state should also manifest a value that is characteristic of the parent electronic configuration. $μ_e$ and other electrostatic and magnetostatic properties should be less influenced by the vibronic state mixing, providing more robust indicators for configurational assignments. In the present study, we have measured the $μ_e$ values for four electronic states of UO. Here, the results clearly demonstrate that the ground state (X(1)4) and the first electronically excited state ((2)4) are derived from the U 2+ ($5f^37s$)O 2– and U 2+ ($5f^27s^2$)O 2– configurations, respectively.

74 ATOMIC AND MOLECULAR PHYSICS↗

New Dimension in Ab Initio Electronic Structure Theory: Temperature, Pressure, and Chemical Potential

Ab initio electronic structure theory has transformed gas-phase molecular science with its predictive ability. In the attempt to bring such predictive ability to macroscopic systems and condensed matter, the theory must integrate quantum mechanics with statistical thermodynamics, so that thermodynamic functions such as free energy, internal energy, entropy, and chemical potentials are computed as functions of temperature in a systematically converging series of approximations. Here, a general, versatile strategy of elevating ab initio electronic structure theory to nonzero temperatures is introduced and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Femtosecond Nonadiabatic Confinement of Molecular Dication Yield

Doubly charged molecular cations often carry signatures of electronic correlation and electron–nuclear entanglement present in the parent cation. Here, we produce ethylene dications using a combination of an extreme ultraviolet pump and near-infrared probe pulses, observing a peak in the dication yield at a pump–probe delay of approximately 15 fs. Ab-initio calculations, which explicitly take into account coupled electron–nuclear dynamics induced by the pump and the multiphoton nature of the probe-induced ionization step, reproduced the observed delay in the yield. It originates from resonant enhancement of the multiphoton ionization of the electronically excited ethylene cation as the carbon–carbon double bond expands. However, this effect is tempered by rapid nonadiabatic relaxation of the excited ionic states. Our results suggest a general mechanism whereby ultrafast nonadiabatic relaxation of a molecular ion can compete with its strong-field ionization rate, confining the dication yield to a narrow temporal window of a few femtoseconds.

74 ATOMIC AND MOLECULAR PHYSICS↗

Fine-Tuning Microporosity of Crystalline Vanadomolybdate Frameworks for Selective Adsorptive Separation of Kr from Xe

Selective adsorptive capture and separation of chemically inert krypton (Kr) and xenon (Xe) noble gases with very low ppmv concentrations in air and industrial off-gases constitute an important technological challenge. Here, using a synergistic combination of experiment and theory, the microporous crystalline vanadomolybdates (MoVO x ) as highly selective Kr sorbents are studied in detail. By varying the Mo/V ratios, we show for the first time that their one-dimensional (1D) pores can be fine-tuned for the size-selective adsorption of Kr over the larger Xe with selectivities reaching >100. Using extensive electronic structure calculations and grand canonical Monte Carlo simulations, the competition between Kr uptake with CO 2 and N 2 was also investigated. As most materials reported so far are selective toward the larger, more polarizable Xe than Kr, this work constitutes an important step toward robust Kr-selective sorbent materials. Furthermore this work highlights the potential use of porous crystalline transition metal oxides as energy-efficient and selective noble gas capture sorbents for industrial applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Enhanced Ion Conductivity through Hydrated, Polyelectrolyte-Grafted Cellulose Nanocrystal Films

Here,iIonically conductive, hydrated polyelectrolyte nanocomposites are prepared from iodomethane-treated poly(2-vinylpyridine) (mPV)-grafted cellulose nanocrystals (MxG-CNC-g-mPV). These polyelectrolyte-grafted nanoparticle (PEGN) films exhibit an order-of-magnitude higher iodide ion conductivity relative to mPV films and a high in-plane/through-plane anisotropy. The PEGN architecture prevents CNC aggregation, maximizing the CNC/polyelectrolyte interface. PEGN films were prepared with varying polymer graft density (0.03-0.12 chains/nm 2 ) and molecular weight (7k-30k g/mol). The greatest ion conductivity enhancement is observed with lower molecular weight, higher density grafts: ca. 89 ± 6 mS/cm (140 ±10 mS/cm accounting for the volume of the CNC) versus 3.3 ± 0.4 mS/cm for ungrafted mPV. Poly(styrene-block-2-vinylpyridine)-grafted CNCs were prepared in which the insulating polystyrene block or the conducting mPV block was directly attached to the CNCs; only the films with the mPV block closest to the CNCs exhibited an enhancement in conductivity relative to mPV. Together, these data point to beneficial CNC/polyelectrolyte interfacial effects, resulting in significant ionic conductivity enhancement along the length of the CNCs in these films.

74 ATOMIC AND MOLECULAR PHYSICS↗