Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The unusual symmetry of hexafluoro-o-xylene—A microwave spectroscopy and computational study

The rotational constants and quartic centrifugal distortion coefficients of hexafluoro-o-xylene were precisely derived from the 8 GHz to 18 GHz gas phase microwave spectrum. In addition, the rotational constants of all singly substituted 13 C isotopologues were determined. Instead of the intuitively expected symmetry of C 2v , as in o-xylene, calculations with a variety of methods (B3LYP, CAM-B3LYP, ωB97XD, MP2, and coupled-cluster singles, doubles, and perturbative triples) predict a C 2 symmetry structure in which the two CF 3 groups rotate in opposite directions by about 16°. While the experimental results in this study are not capable of proving this unusual symmetry, they can support future microwave, gas phase electron diffraction or nuclear magnetic resonance studies. From the presented data, a preliminary r 0 structure was determined, reproducing the experimental rotational constants with deviations of no more than 15 kHz. Analysis of the interactions between the two CF 3 groups using an effective fragment potential approach identifies two major contributions to their interaction, due to exchange-repulsion and electrostatic repulsion, with electrostatic repulsion responsible for the barrier at the C 2v geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Local and global approaches to treat the torsional barriers of 4-methylacetophenone using microwave spectroscopy

The Fourier transform microwave spectrum of 4-methylacetophenone recorded from 8 GHz to 18 GHz under jet-cooled conditions has revealed large tunneling splittings arising from a low barrier to internal rotation of the ring methyl group and small splittings from a high torsional barrier of the acetyl methyl group. The large splittings are especially challenging to model, while the small splittings are difficult to analyze due to the resolution limit of 120 kHz. The combination of two methyl groups undergoing internal rotations caused each rotational transition to split into five torsional species, which were resolved and fitted using a modified version of the XIAM code and the newly devel- oped ntop code to a root-mean-square deviation close to measurement accuracy, providing an estimate of the V 3 potential barriers of about 22 cm -1 and 584–588 cm -1 for the ring and the acetyl methyl groups, respectively. The assignment was aided by separately fitting the five torsional species using odd-power order operators. Only one conformer in which all heavy atoms are located on a symmetry plane could be identified in the spectrum, in agreement with results from conformation analysis using quantum chemical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Signatures of electronic and nuclear coherences in ultrafast molecular x-ray and electron diffraction

Femtosecond x-ray and electron diffraction hold promise to image the evolving structures of single molecules. We present a unified quantum-electrodynamical formulation of diffraction signals, based on the exact many-body nuclear + electronic wavefunction that can be extracted from quantum chemistry simulations. This gives a framework for analyzing various approximate molecular dynamics simulations. We show that the complete description of ultrafast diffraction signals contains interesting contributions involving mixed elastic and inelastic scattered photons that are usually masked by other larger contributions and are neglected. These terms include overlaps of nuclear wavepackets between different electronic states that provide an electronic decoherence mechanism and are important for the time-resolved imaging of conical intersections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Laser-induced electron diffraction of the ultrafast umbrella motion in ammonia

Visualizing molecular transformations in real-time requires a structural retrieval method with Ångström spatial and femtosecond temporal atomic resolution. Imaging of hydrogen-containing molecules additionally requires an imaging method sensitive to the atomic positions of hydrogen nuclei, with most methods possessing relatively low sensitivity to hydrogen scattering. Laser-induced electron diffraction (LIED) is a table-top technique that can image ultrafast structural changes of gas-phase polyatomic molecules with sub-Ångström and femtosecond spatiotemporal resolution together with relatively high sensitivity to hydrogen scattering. Here, we image the umbrella motion of an isolated ammonia molecule (NH3) following its strong-field ionization. Upon ionization of a neutral ammonia molecule, the ammonia cation (NH3+) undergoes an ultrafast geometrical transformation from a pyramidal (ΦHNH=107°) to planar (ΦHNH=120°) structure in approximately 8 femtoseconds. Using LIED, we retrieve a near-planar (ΦHNH=117 ± 5°) field-dressed NH3+ molecular structure 7.8−9.8 femtoseconds after ionization. Our measured field-dressed NH3+ structure is in excellent agreement with our calculated equilibrium field-dressed structure using quantum chemical ab initio calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Watching a signaling protein function: What has been learned over four decades of time-resolved studies of photoactive yellow protein

Photoactive yellow protein (PYP) is a signaling protein whose internal p-coumaric acid chromophore undergoes reversible, light-induced trans-to-cis isomerization, which triggers a sequence of structural changes that ultimately lead to a signaling state. Since its discovery nearly 40 years ago, PYP has attracted much interest and has become one of the most extensively studied proteins found in nature. The method of time-resolved crystallography, pioneered by Keith Moffat, has successfully characterized intermediates in the PYP photocycle at near atomic resolution over 12 decades of time down to the sub-picosecond time scale, allowing one to stitch together a movie and literally watch a protein as it functions. But how close to reality is this movie? To address this question, results from numerous complementary time-resolved techniques including x-ray crystallography, x-ray scattering, and spectroscopy are discussed. Emerging from spectroscopic studies is a general consensus that three time constants are required to model the excited state relaxation, with a highly strained ground-state cis intermediate formed in less than 2.4 ps. Persistent strain drives the sequence of structural transitions that ultimately produce the signaling state. Crystal packing forces produce a restoring force that slows somewhat the rates of interconversion between the intermediates. Moreover, the solvent composition surrounding PYP can influence the number and structures of intermediates as well as the rates at which they interconvert. When chloride is present, the PYP photocycle in a crystal closely tracks that in solution, which suggests the epic movie of the PYP photocycle is indeed based in reality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Diffused morphotropic phase boundary in relaxor-PbTiO 3 crystals: High piezoelectricity with improved thermal stability

Solid solution ferroelectrics are the most widely used piezoelectric material for numerous electromechanical applications, including sensors, actuators, and transducers. A milestone in ferroelectric research is the discovery of the morphotropic phase boundary that was first reported in Pb(Zr x Ti 1–x )O 3 , which has been extensively solicited to improve the performance of various solid solution ferroelectrics, including those having the highest piezoelectricity known today. However, due to the inherent correlation between phase transition and thermodynamic imbalance, the efforts of building the phase boundary encounter the challenge that high performance materials are generally accompanied by property instability. In this paper, we report a comprehensive study on the crystalline symmetry and polar nano-regions in relaxor-PbTiO 3 crystals by synchrotron measurements and property characterizations. In contrast to the common belief that the morphotropic phase boundary is a narrow composition region, the morphotropic phase boundary in ternary Pb(In 1/2 Nb 1/2 )O 3 -Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 is demonstrated as an extensive region covering at least 9% of the PbTiO 3 content. Such a diffused morphotropic phase boundary is associated with the intensive interaction of polar nano-regions, leading to high piezoelectricity (>1500 pC/N) with greatly improved thermal stability, where the piezoelectric variation is ~90% over the temperature range of 273–373 K, which is about a factor of 3 lower compared to its binary counterpart Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 . This work sheds light on the fundamental understanding of nanoscale inhomogeneity and macroscopic symmetry in relaxor-PbTiO 3 , and successfully interlinks the structure and properties in complex solid solutions. The existence of a diffused morphotropic phase boundary is also expected to benefit other ferroic systems beyond ferroelectrics, such as ferromagnetic and multiferroic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree–Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

On the origin of ground-state vacuum-field catalysis: Equilibrium consideration

Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this “vacuum-field catalysis” remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition state configurations. In this work, we examine how VSC can affect chemical reactions rates in a cavity environment according to transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons and also assuming no charge overlap between molecules, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual micron-length cavities used to examine VSC situations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Strictly non-adiabatic quantum control of the acetylene dication using an infrared field

We demonstrate the existence of a strictly non-adiabatic control pathway in deprotonation of the acetylene dication. This pathway is identified experimentally by measuring a kinetic energy shift in an ion coincidence experiment. We use a time dependent Schrödinger equation simulation to identify which properties most strongly affect our control. We find that resonant control around conical intersections is limited by the speed of non-adiabatic dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction

We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N 3 ) operations and O(N 2 ) memory, where N is the number of basis functions. Furthermore, this limits the overall computational scaling to O(N 4 ) operations and O(N 2 ) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Multichannel dynamics in the OH+ n -butane reaction revealed by crossed-beam slice imaging and quasiclassical trajectory calculations

Multidimensional reactions present various channels that can exhibit very different dynamics and give products of varying subsequent reactivity. Here, we present a combination of experiment and theory to reveal the dynamics of hydrogen abstraction by OH radical at primary and secondary sites in n-butane at a collision energy of 8 kcal/mol. Crossed molecular beam slice imaging experiments unequivocally probe the secondary abstraction channel showing backward angular distributions with mild energy release to product translation, which are accurately captured by trajectory calculations using a specific-reaction-parameter Hamiltonian. Additionally, experiments containing both reaction channels indicate a less marked backward character in the angular distribution, whose origin is shown by trajectory calculations to appear as an evolution toward more sideways scattering from the secondary to primary channel. While the two channels have markedly different angular distributions, their energy release is largely comparable, showing ample energy release into the water product. The synergistic combination of crossed-beam imaging and trajectories opens the door to detailed reaction-dynamics studies of chemical reactions with ever-increasing complexity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations

The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H 5 O 2 + and Eigen H 9 O 4 + ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure–spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm –1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. In conclusion, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Isoenergetic two-photon excitation enhances solvent-to-solute excited-state proton transfer

Two-photon excitation (TPE) is an attractive means for controlling chemistry in both space and time. Since isoenergetic one- and twophoton excitations (OPE and TPE) in non-centrosymmetric molecules are allowed to reach the same excited state, it is usually assumed that they produce similar excited-state reactivity. We compare the solvent-to-solute excited-state proton transfer of the super photobase FR0-SB following isoenergetic OPE and TPE. We find up to 62% increased reactivity following TPE compared to OPE. From steady-state spectroscopy, we rule out the involvement of different excited states and find that OPE and TPE spectra are identical in non-polar solvents but not in polar ones. We propose that differences in the matrix elements that contribute to the two-photon absorption cross sections lead to the observed enhanced isoenergetic reactivity, consistent with the predictions of our high-level coupledcluster-based computational protocol. We find that polar solvent configurations favor greater dipole moment change between ground and excited states, which enters the probability for TPE as the absolute value squared. This, in turn, causes a difference in the Franck– Condon region reached via TPE compared to OPE. We conclude that a new method has been found for controlling chemical reactivity via the matrix elements that affect two-photon cross sections, which may be of great utility for spatial and temporal precision chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

When do short-range atomistic machine-learning models fall short?

In this work, we explore the role of long-range interactions in atomistic machine-learning models by analyzing the effects on fitting accuracy, isolated cluster properties, and bulk thermodynamic properties. Such models have become increasingly popular in molecular simulations given their ability to learn highly complex and multi-dimensional interactions within a local environment; however, many of them fundamentally lack a description of explicit long-range interactions. In order to provide a well-defined benchmark system with precisely known pairwise interactions, we chose as the reference model a flexible version of the Extended Simple Point Charge (SPC/E) water model. Our analysis shows that while local representations are sufficient for predictions of the condensed liquid phase, the short-range nature of machine-learning models falls short in representing cluster and vapor phase properties. These findings provide an improved understanding of the role of long-range interactions in machine learning models and the regimes where they are necessary.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Improved accuracy and transferability of molecular-orbital-based machine learning: Organics, transition-metal complexes, non-covalent interactions, and transition states

We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗