SEARCH · Engineering Papers
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.
The torsion-rotation spectrum of ground states ethyl and propargyl alcohol have been successfully analyzed with an asymmetric top-asymmetric frame Hamiltonian to J=60 and J=80, respectively.
This thesis describes several new quantum algorithms. These include a polynomial time algorithm that uses a quantum fast Fourier transform to find eigenvalues and eigenvectors of a Hamiltonian operator, and that can be applied in cases for which all know classical algorithms require exponential time.
The laboratory submillimeter-wave rotational spectrum of the 13 CH 3 CN, CH 3 13 CN, and the CH 3 C 15 N isotopomers of methyl cyanide has been observed in natural abundance in the 294 to 607 GHz region. The maximum J and K values are 34 an 14, respectively. Fifteen additional CH 3 CN transitions up to K=21 were also measured. The transitions of all four species are fitted to a symmetric top Hamiltonian and the rotation and distortion constants are determined. The 14 N quadropole and spin rotation are also calculated and presented. Suggested values for many other parameters which could not be directly determined from the isotope spectra are calculated from the normal species values and isotope relationships. The determined and calculated constants should predict the spectrum of the three isotopomers to well over 1 THz accurately enough for astronomical assignments.
Molecular IR line lists computed from semi-empirically refined ab initio potential energy surface and high quality ab initio dipole moment surface may have line position accuracy typically within 0.01 - 0.05 cm−1. For high resolution spectroscopy databases, it is necessary to integrate the computed theoretical IR intensities with the line positions accurately determined from experiments and (or) reliably derived from Effective Hamiltonian models. However, it is not a trivial task, and the choices in reality are heavily contingent upon the coverage, consistency, and accuracy of the data available for a specific molecule or isotopologue. We will present several approaches applied to recent Ames IR line lists of CO 2 1, N 2 O 2 , and OCS to demonstrate how challenging such integration may become, and what we have learned From these projects. Discussions will focus on the relation between each different scenario and the corresponding choice or solution, including their advantages and limits, how the semi-empirical IR line lists can help, and what else may be needed for future improvements. We will emphasize that the “Theory+Experiment” synergy may still play significant role in the determination of the best line positions.
Abstract Quantum circuit learning is employed to simulate quantum field theories (QFTs). Typically, when simulating QFTs with quantum computers, significant challenges are encountered due to the technical limitations of quantum devices when implementing the Hamiltonian using Pauli spin matrices. To address this challenge, quantum circuit learning is leveraged, employing a compact configuration of qubits and low‐depth quantum circuits to predict real‐time dynamics in quantum field theories. The key advantage of this approach is that a single‐qubit measurement can accurately forecast various physical parameters, including fully‐connected operators. To demonstrate the effectiveness of this method, it is used to predict quench dynamics, chiral dynamics and jet production in a 1+1‐dimensional model of quantum electrodynamics. It is found that our predictions closely align with the results of rigorous classical calculations, exhibiting a high degree of accuracy. This hybrid quantum‐classical approach illustrates the feasibility of efficiently simulating large‐scale QFTs on cutting‐edge quantum devices.
Abstract The character of the electronic structure of acenes has been the subject of longstanding discussion. However, convincing experimental evidence of their open‐shell character has so far been missing. Here, we present the on‐surface synthesis of tridecacene molecules by thermal annealing of octahydrotridecacene on a Au(111) surface. We characterized the electronic structure of the tridecacene by scanning probe microscopy, which reveals the presence of an inelastic signal at 126 meV. We attribute the inelastic signal to spin excitation from the singlet diradical ground state to the triplet excited state. To rationalize the experimental findings, we carried out many‐body ab initio calculations as well as model Hamiltonians to take into account the effect of the metallic substrate. Moreover, we provide a detailed analysis of how the dynamic electron correlation and virtual charge fluctuation between the molecule and metallic surface reduces the singlet‐triplet band gap. Thus, this work provides the first experimental confirmation of the magnetic character of tridecacene.
Abstract The sustained interest in investigating magnetism in the 2D limit of insulating antiferromagnets is driven by the possibilities of discovering, or engineering, novel magnetic phases through layer stacking. However, due to the difficulty of directly measuring magnetic interactions in 2D antiferromagnets, it is not yet understood how intra layer magnetic interactions in insulating , strongly correlated, materials can be modified through layer proximity. Herein, the impact of reduced dimensionality in the model van der Waals antiferromagnet NiPS 3 is explored by measuring electronic excitations in exfoliated samples using Resonant Inelastic X‐ray Scattering (RIXS). The resulting spectra shows systematic broadening of NiS 6 multiplet excitations with decreasing layer count from bulk down to three atomic layers (3L). It is shown that these trends originate from a decrease in transition metal‐ligand and ligand–ligand hopping integrals, and by charge‐transfer energy evolving from Δ = 0.83 eV in the bulk to 0.37 eV in 3L NiPS 3 . Relevant intralayer magnetic exchange integrals computed from the electronic parameters exhibit a decrease in the average interaction strength with thickness. This study underscores the influence of inter layer electronic interactions on intra layer ones in insulating magnets, indicating that magnetic Hamiltonians in few‐layer insulating magnets can greatly deviate from their bulk counterparts.
Abstract Manipulating the magnetic ground states of 2D magnets is a focal point of recent research efforts. Various methods have demonstrated efficacy in modulating the magnetic properties inherent to van der Waals (vdW) magnetic systems. Herein, the emergence of robust anisotropic ferromagnetism within antiferromagnetic FePS 3 is unveiled via intercalation with non‐magnetic pyridinium ions. A one‐step ion exchange reaction facilitates the formation of energetically favorable B‐phase and metastable P‐phase. Notably, both B‐ and P‐phases manifest hard ferromagnetic behavior, featuring substantial unsaturated coercive fields (>7 T) and high Curie temperatures (72–87 K). First‐principles calculations elucidate the pivotal role of electron transfer from pyridinium ions to FePS 3 in engineering magnetic exchange interactions. Calculated effective spin Hamiltonian corroborates the observed hard ferromagnetism in intercalated FePS 3 . This study offers crucial insights into hard magnetism in intercalated vdW materials, thereby presenting promising avenues for 2D vdW magnet‐based magnetic devices.
Abstract Manganese (Mn) nitrides are important nitrogen (N) carriers for small‐scale intermittent ammonia (NH 3 ) synthesis. However, only 3∼8 % of lattice N are converted into NH 3 . In this study, the geometric and electronic structures of well‐defined Mn 4 N and Mn 2 N lattices were altered using transition metal heteroatoms (Cr, Fe, Co, Ni, Mo) to understand the driving force behind lattice N diffusion and extraction. Density Functional Theory (DFT) revealed that the binding of early hydrogenation product (NH) follows a linear relationship with lattice N over a wide range of close‐packed surfaces. But, the binding of NH 2 and NH 3 are more sensitive to the geometric and electronic structures. Further, the chemical bonding can be quantitatively characterized with the covalency derived from Crystal Orbital Hamiltonian Population (COHP). In the Eley Rideal‐Mars van Krevelan pathway, the overall NH 3 formation free energy ( ) and lattice N diffusion barrier ( ) are the respective thermodynamic and kinetic determining factors. Aided by a rate‐determining step (RDS) model, Mn 4 N modified with Fe, Co, Ni single‐atom dopants all show enhanced rates for NH 3 formation.
In this paper, we present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn–Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097.
Polypeptoids differ from polypeptides in that the amide bond can more frequently adopt both cis and trans conformations. The transition between the two conformations requires overcoming a large energy barrier, making it difficult for conventional molecular simulations to adequately visit the cis and trans structures. A replica-exchange method is presented that allows for easy rotations of the amide bond and also an efficient linking to a high temperature replica. The method allows for just three replicas (one at the temperature and Hamiltonian of interest, a second high temperature replica with a biased dihedral potential, and a third connecting them) to overcome the amide bond sampling problem and also enhance sampling for other coordinates. Here our results indicate that for short peptoid oligomers, the conformations can range from all cis to all trans with an average cis/trans ratio that depends on side chain and potential model.
Quantum computers are expected to outperform classical computers for specific problems in quantum chemistry. Such calculations remain expensive, but costs can be lowered through the partition of the molecular system. In the present study, partition was achieved with range-separated density functional theory (RS-DFT). The use of RS-DFT reduces both the basis set size and the active space size dependence of the ground state energy in comparison with the use of wave function theory (WFT) alone. The utilization of pair natural orbitals (PNOs) in place of canonical molecular orbitals (MOs) results in more compact qubit Hamiltonians. To test this strategy, a basis-set independent framework, known as multiresolution analysis (MRA), was employed to generate PNOs. Tests were conducted with the variational quantum eigensolver for a number of molecules. The results show that the proposed approach reduces the number of qubits needed to reach a target energy accuracy.
In this paper, the synthesis, structure, and physical properties of Y 5 Ru 2 O 12 , a new member of the Ln 5 B 2 O 12 class (Ln = lanthanide, B = Ru, Re, Mo) are reported. Single crystals of Y 5 Ru 2 O 12 are synthesized using a high-temperature, high-pressure (700 °C, 200 MPa) hydrothermal method. The structure features edge-sharing Ru—O—Ru 1D chains of distorted RuO6 octahedra. The YO n (n = 6 and 7) polyhedral groups bridge Ru—O—Ru chains across the unit cell by forming a complex 3D structure. The Ru—O—Ru chains contain alternating shorter and longer Ru—Ru distances along the chain. Single crystal neutron diffraction shows streak-like diffuse scattering along the H-direction, due to stacking faults. Temperature-dependent anisotropic magnetic measurements reveal a broad feature over a wide temperature range, with a maximum at ≈70 K indicating potential short-range ordering, with canted antiferromagnetic behavior below 30 K. The dominant magnetic interactions are antiferromagnetic with a large negative Weiss temperature θ ≈ –200 K. The compound has an effective magnetic moment that is anomalously large for this class of compounds (≈4 μ B Ru –1 ), and a phenomenological spin Hamiltonian describing the unusual anisotropic magnetization of Y 5 Ru 2 O 12 is proposed.
Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.
The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.
One of the most striking but mysterious properties of the sinh-Gordon model (ShG) is the b → 1/b self-duality of its S-matrix, of which there is no trace in its Lagrangian formulation. Here b is the coupling appearing in the model’s eponymous hyperbolic cosine present in its Lagrangian, cosh(bΦ). In this paper we develop truncated spectrum methods (TSMs) for studying the sinh-Gordon model at a finite volume as we vary the coupling constant. We obtain the expected results for b << 1 and intermediate values of b, but as the self-dual point b = 1 is approached, the basic application of the TSM to the ShG breaks down. We find that the TSM gives results with a strong cutoff Ec dependence, which disappears according only to a very slow power law in Ec. Standard renormalization group strategies — whether they be numerical or analytic — also fail to improve upon matters here. We thus explore three strategies to address the basic limitations of the TSM in the vicinity of b = 1. In the first, we focus on the small-volume spectrum. We attempt to understand how much of the physics of the ShG is encoded in the zero mode part of its Hamiltonian, in essence how ‘quantum mechanical’ vs ‘quantum field theoretic’ the problem is. In the second, we identify the divergencies present in perturbation theory and perform their resummation using a supra-Borel approximate. In the third approach, we use the exact form factors of the model to treat the ShG at one value of b as a perturbation of a ShG at a different coupling. In the light of this work, we argue that the strong coupling phase b > 1 of the Lagrangian formulation of model may be different from what is naïvely inferred from its S-matrix. In particular, we present an argument that the theory is massless for b > 1.
We obtain a conservative Hamiltonian describing the interactions of two charged bodies in Yang-Mills through $\mathcal{O}$ (α 2 ) and to all orders in velocity. Our calculation extends a recently-introduced framework based on scattering amplitudes and effective field theory (EFT) to consider color-charged objects. These results are checked against the direct integration of the observables in the Kosower-Maybee-O’Connell (KMOC) formalism. At the order we consider we find that the linear and color impulses in a scattering event can be concisely described in terms of the eikonal phase, thus extending the domain of applicability of a formula originally proposed in the context of spinning particles.
We consider a nonlocal scalar field theory inspired by the tachyon action in open string field theory. The Lorentz-covariant action is characterized by a parameter ξ 2 that quantifies the amount of nonlocality. Restricting to purely time-dependent configurations, we show that a field redefinition perturbative in ξ 2 reduces the action to a local two-derivative theory with a ξ 2 -dependent potential. This picture is supported by evidence that the redefinition maps the wildly oscillating rolling tachyon solutions of the nonlocal theory to conventional rolling in the new scalar potential. For general field configurations we exhibit an obstruction to a local Lorentz-covariant formulation, but we can still achieve a formulation local in time, as well as a light-cone formulation. These constructions provide an initial value formulation and a Hamiltonian. Their causality is consistent with a lack of superluminal behavior in the nonlocal theory.