Theoretical Insights into the Dynamics of Gas-Phase Bimolecular Reactions with Submerged Barriers
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Engineering Papers
Results for “quantum chemical dynamics”
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.
Evaluating the stability and activity of dilute Cu-based alloys for electrochemical CO 2 reduction
Cu-based catalysts currently offer the most promising route to actively and selectively produce value-added chemicals via electrochemical reduction of CO 2 (eCO 2 R); yet further improvements are required for their wide-scale deployment in carbon mitigation efforts. Here, we systematically investigate a family of dilute Cu-based alloys to explore their viability as active and selective catalysts for eCO 2 R through a combined theoretical–experimental approach. Using a quantum–classical modeling approach that accounts for dynamic solvation effects, we assess the stability and activity of model single-atom catalysts under eCO 2 R conditions. Our calculations identify that the presence of eCO2R intermediates, such as CO*, H*, and OH*, may dynamically influence the local catalyst surface composition. Additionally, we identify through binding energy descriptors of the CO*, CHO*, and OCCO* dimer intermediates that certain elements, such as group 13 elements (B, Al, and Ga), enhance the selectivity of C 2+ species relative to pure Cu by facilitating CO dimerization. The theoretical work is corroborated by preliminary testing of eCO 2 R activity and selectivity of candidate dilute Cu-based alloy catalyst films prepared by electron beam evaporation in a zero-gap gas diffusion electrode-based reactor. Of all studied alloys, dilute CuAl was found to be the most active and selective toward C 2+ products like ethylene, consistent with the theoretical predictions. Finally, we attribute the improved performance of dilute CuAl alloys to more favorable dimerization reaction energetics of bound CO species relative to that on pure Cu. In a broader context, the results presented here demonstrate the power of our simulation framework in terms of rational catalyst design.
30 DIRECT ENERGY CONVERSION↗
Trapped-ion quantum simulation of electron transfer models with tunable dissipation
Electron transfer is at the heart of many fundamental physical, chemical, and biochemical processes essential for life. The exact simulation of these reactions is often hindered by the large number of degrees of freedom and by the essential role of quantum effects. Here, we experimentally simulate a paradigmatic model of molecular electron transfer using a multispecies trapped-ion crystal, where the donor-acceptor gap, the electronic and vibronic couplings, and the bath relaxation dynamics can all be controlled independently. By manipulating both the ground-state and optical qubits, we observe the real-time dynamics of the spin excitation, measuring the transfer rate in several regimes of adiabaticity and relaxation dynamics. Our results provide a testing ground for increasingly rich models of molecular excitation transfer processes that are relevant for molecular electronics and light-harvesting systems.
Science & Technology - Other Topics↗
Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization
We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.
Jaffe, Richard L.↗
Suppression of superconductivity by anisotropic strain near a nematic quantum critical point
In most unconventional and high-temperature superconductors, superconductivity emerges as a nearby symmetry-breaking phase is suppressed by chemical doping or pressure. This has led to the belief that the fluctuations associated with the symmetry-breaking phase are beneficial, if not responsible, for the superconducting pairing. A direct test to verify this hypothesis is to observe a decrease of the superconducting critical temperature (T c ) by applying the symmetry-breaking conjugate field that suppresses the dynamic fluctuations of the competing order. However, most of the competing phases in unconventional superconductors break translational symmetry, requiring a spatially modulated conjugate field that is difficult to realize experimentally. Here, we show that anisotropic strain, the conjugate field of nematicity, reduces the T c of an iron pnictide. For optimally doped samples we show a fivefold reduction of T c with less than one per cent of strain. For underdoped samples, T c becomes zero yielding a fully metallic ground state. In addition to providing direct evidence of the role played by the nematic fluctuations in the formation of the superconducting state, these results demonstrate tunable mechanical control of a high-temperature superconductor, an important step forward for technological applications of superconductivity.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Photo‐Induced Charge State Dynamics of the Neutral and Negatively Charged Silicon Vacancy Centers in Room‐Temperature Diamond
Abstract The silicon vacancy (SiV) center in diamond is drawing much attention due to its optical and spin properties, attractive for quantum information processing and sensing. Comparatively little is known, however, about the dynamics governing SiV charge state interconversion mainly due to challenges associated with generating, stabilizing, and characterizing all possible charge states, particularly at room temperature. Here, multi‐color confocal microscopy and density functional theory are used to examine photo‐induced SiV recombination — from neutral, to single‐, to double‐negatively charged — over a broad spectral window in chemical‐vapor‐deposition (CVD) diamond under ambient conditions. For the SiV 0 to SiV ‐ transition, a linear growth of the photo‐recombination rate with laser power at all observed wavelengths is found, a hallmark of single photon dynamics. Laser excitation of SiV ‒ , on the other hand, yields only fractional recombination into SiV 2‒ , a finding that is interpreted in terms of a photo‐activated electron tunneling process from proximal nitrogen atoms.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt
Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.
42 ENGINEERING↗
MBX: A many-body energy and force calculator for data-driven many-body simulations
Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Nature of the nonequilibrium phase transition in the non-Markovian driven Dicke model
The Dicke model famously exhibits a phase transition to a superradiant phase with a macroscopic population of photons and is realized in multiple settings in open quantum systems. In this paper, we study a variant of the Dicke model where the cavity mode is lossy due to the coupling to a Markovian environment while the atomic mode is coupled to a colored bath. We analytically investigate this model by inspecting its low-frequency behavior via the Schwinger-Keldysh field theory and carefully examine the nature of the corresponding superradiant phase transition. Integrating out the fast modes, we can identify a simple effective theory allowing us to derive analytical expressions for various critical exponents including the dynamical exponent. We find excellent agreement with previous numerical results when the non-Markovian bath is at zero temperature; however, contrary to these studies, our low-frequency approach reveals that the same exponents govern the critical behavior when the colored bath is at finite temperature unless the chemical potential is zero. Furthermore, we show that the superradiant phase transition is classical in nature, while it is genuinely nonequilibrium. In this work, we derive a fractional Langevin equation and conjecture the associated fractional Fokker-Planck equation that captures the system's long-time memory as well as its nonequilibrium behavior. Finally, we consider finite-size effects at the phase transition and identify the finite-size scaling exponents, unlocking a rich behavior in both statics and dynamics of the photonic and atomic observables.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Amplifying Magnetic Field Effects on Upconversion Emission via Molecular Qubit-Driven Triplet–Triplet Annihilation
Triplet-triplet annihilation (TTA) enables photon upconversion by combining two lower-energy triplet excitons to produce a higher-energy singlet exciton. This mechanism enhances light-harvesting efficiency for solar energy conversion and enables the use of lower-energy photons in bioimaging and photoredox catalysis applications. The magnetic modulation of such high-energy excitons presents an exciting opportunity to develop molecular quantum information technologies. While the spin dynamics of triplet exciton pairs are sensitive to external magnetic fields, the magnetic field effects (MFEs) associated with these pairs are generally limited by spin statistics to at most 10% at low fields (<1 Tesla), making them challenging to apply in technological advancements. In contrast, MFEs on spin-correlated radical pairs (SCRPs) can be significantly greater, surpassing those on triplet pairs. By using SCRPs-based molecular qubits as triplet sensitizers in the sensitized TTA scheme, we can magnetically modulate TTA and consequently, the delayed fluorescence of annihilators. In our current system, we have achieved more than 70% magnetic modulation of delayed fluorescence, effectively harnessing and even amplifying magnetic modulation within SCRPs to influence high-energy excitons. In conclusion, this work opens new opportunities for advancing spin-controlled chemical reactions and molecular quantum information technologies.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Quantum algorithms for open lattice field theory
Certain aspects of some unitary quantum systems are well described by evolution via a non-Hermitian effective Hamiltonian, as in the Wigner-Weisskopf theory for spontaneous decay. Conversely, any non-Hermitian Hamiltonian evolution can be accommodated in a corresponding unitary system + environment model via a generalization of Wigner-Weisskopf theory. This demonstrates the physical relevance of novel features such as exceptional points in quantum dynamics, and opens up avenues for studying many-body systems in the complex plane of coupling constants. In the case of lattice field theory, sparsity lends these channels the promise of efficient simulation on standardized quantum hardware. We thus consider quantum operations that correspond to Suzuki-Lie-Trotter approximation of lattice field theories undergoing nonunitary time evolution, with potential applicability to studies of spin or gauge models at finite chemical potential, with topological terms, to quantum phase transitions—a range of models with sign problems. We develop non-Hermitian quantum circuits and explore their promise on a benchmark, the quantum one-dimensional Ising model with complex longitudinal magnetic field, showing that observables can probe the Lee-Yang edge singularity. The development of attractors past critical points in the space of complex couplings indicates a potential for study on near-term noisy hardware.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Ab initio Molecular Dynamics Beyond Density Functional Theory
The funding from this award supported the development of new theoretical methods that both increase the accuracy of quantum simulations beyond that of density functional theory, as well as reduce the cost of such simulations, for example, through machine learning and algorithmic improvements. The improvements lead to, in some cases, orders of magnitude increases in simulation speed as well as orders of magnitude increase in accuracy, compared to before this work was started. This translates to both greater certainty in making predictions about chemical phenomena that are currently studied computationally, as well as laying the ground-work for the study of new areas of chemistry that have not traditionally been simulated, because of the new time-scales and accuracies that can be reached.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation
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.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Dynamic aspects of graphene deformation and fracture from approximate density functional theory
Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Computing Free Energies with Fluctuation Relations on Quantum Computers
As a central thermodynamic property, free energy enables the calculation of virtually any equilibrium property of a physical system, allowing for the construction of phase diagrams and predictions about transport, chemical reactions, and biological processes. Thus, methods for efficiently computing free energies, which in general is a difficult problem, are of great interest to broad areas of physics and the natural sciences. The majority of techniques for computing free energies target classical systems, leaving the computation of free energies in quantum systems less explored. Recently developed fluctuation relations enable the computation of free energy differences in quantum systems from an ensemble of dynamic simulations. While performing such simulations is exponentially hard on classical computers, quantum computers can efficiently simulate the dynamics of quantum systems. Here, we present an algorithm utilizing a fluctuation relation known as the Jarzynski equality to approximate free energy differences of quantum systems on a quantum computer. In this work we discuss under which conditions our approximation becomes exact, and under which conditions it serves as a strict upper bound. Furthermore, we successfully demonstrate a proof of concept of our algorithm using the transverse field Ising model on a real quantum processor. As quantum hardware continues to improve, we anticipate that our algorithm will enable computation of free energy differences for a wide range of quantum systems useful across the natural sciences.
97 MATHEMATICS AND COMPUTING↗
Enhanced Emission from Bright Excitons in Asymmetrically Strained Colloidal CdSe/Cd x Zn 1– x Se Quantum Dots
Colloidal CdSe quantum dots (QDs) designed with a high degree of asymmetric internal strain have recently been shown to host a number of desirable optical properties including subthermal room-temperature line widths, suppressed spectral diffusion, and high photoluminescence (PL) quantum yields. It remains an open question, however, whether they are well-suited for applications requiring emission of identical single photons. Here we measure the low-temperature PL dynamics and the polarization-resolved fluorescence line narrowing spectra from ensembles of these strained QDs. Our spectroscopy reveals the radiative recombination rates of bright and dark excitons, the relaxation rate between the two, and the energy spectra of the quantized acoustic phonons in the QDs that can contribute to relaxation processes. In comparison to conventional colloidal CdSe/ZnS core/shell QDs, we find that in asymmetrically strained CdSe QDs over six times more light is emitted directly by the bright exciton. Furthermore, these results are therefore encouraging for the prospects of chemically synthesized colloidal QDs as emitters of single indistinguishable photons.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Formation of Carbon Nanotube Based Gears: Quantum Chemistry and Molecular Mechanics Study of the Electrophilic Addition of o-Benzyne to Fullerenes, Graphene, and Nanotubes
Considerable progress has been made in recent years in chemical functionalization of fullerene molecules. In some cases, the predominant reaction products are different from those obtained (using the same reactants) from polycyclic aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It is well established that benzyne adds across one of the rings in naphthalene, anthracene and other PAHs forming the [2+4] cycloaddition product (benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene). However, Hoke et al demonstrated that the only reaction path for o-benzyne with C60 leads to the [2+2] cycloaddition product in which benzyne adds across one of the interpentagonal bonds (forming a cyclobutene ring in the process). Either reaction product results in a loss of aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the benzyne. It is not clear, however, why different products are preferred in these cases. In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of Brenner's potential energy model to describe the structure and stability of these adducts. The Brenner potential has been widely used for describing diamondoid and graphitic carbon. Recently it has also been used for molecular mechanics and molecular dynamics simulations of fullerenes and nanotubes. However, it has not been tested for the case of functionalized fullerenes (especially with highly strained geometries). We use the Brenner potential for our companion nanogear simulations and believe that it should be calibrated to insure that those simulations are physically reasonable. In the present work, Density Functional theory (DFT) calculations are used to determine the preferred geometric structures and energetics for this calibration. The DFT method is a kind of ab initio quantum chemistry method for determining the electronic structure of molecules. For a given basis set expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less accurate than more extensive electron correlation methods such as MP4 or CCSD). However, for systems with large numbers of basis functions it more efficient than any other methods that include electron correlation effects. In this presentation we show the results of DFT calculations for the reaction of benzyne with naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] cycloaddition products. The preferred products for the naphthalene and C60 reactions have been determined by experiment and, thus, these cases serve as a validation of our quantum chemical approach. We also compare the DFT and Brenner potential results. Finally we can predict the likelihood of reaction between benzyne and nanotubes.
Jaffe, Richard↗
Improved techniques for outgoing wave variational principle calculations of converged state-to-state transition probabilities for chemical reactions
Improved techniques and well-optimized basis sets are presented for application of the outgoing wave variational principle to calculate converged quantum mechanical reaction probabilities. They are illustrated with calculations for the reactions D + H2 yields HD + H with total angular momentum J = 3 and F + H2 yields HF + H with J = 0 and 3. The optimization involves the choice of distortion potential, the grid for calculating half-integrated Green's functions, the placement, width, and number of primitive distributed Gaussians, and the computationally most efficient partition between dynamically adapted and primitive basis functions. Benchmark calculations with 224-1064 channels are presented.
Mielke, Steven L.↗