SEARCH · Engineering Papers

Results for “anharmonicity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Dimension-free path-integral molecular dynamics without preconditioning

Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics—including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods—is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann–Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of “dimension-free” numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi–Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

Korol, Roman (ORCID:0000000193076351)↗

DOE OSTI2020

ABINIT: Overview and focus on selected capabilities

ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperature-dependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. ABINIT has strong links with many other software projects that are briefly mentioned.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Ultrafast low-temperature metal–insulator interface phonon dynamics and heat transport in a Pt/Gd 3 Fe 5 O 12 heterostructure

Interfacial thermal and acoustic phenomena have an important role in quantum science and technology, including in spintronic and spincaloritronic materials and devices. Simultaneous measurements of the low-temperature thermal and acoustic properties of a metal/insulator heterostructure reveal distinct dynamics in the characteristic phonon frequency ranges of acoustic and thermal transport. The measurements probed a heterostructure consisting of a thin film of Pt on the ferrimagnetic insulator gadolinium iron garnet (Gd 3 Fe 5 O 12 , GdIG) grown epitaxially on a gadolinium gallium garnet substrate. Ultrafast structural dynamics within the Pt layer were tracked using time-resolved ultrafast x-ray diffraction and analyzed to probe interfacial acoustic and thermal properties. The rapid heating of the Pt layer by a 400 nm wavelength femtosecond-duration optical pulse produced transient structural changes that provided the stimulus for these measurements. Rapid heating produced a broadband acoustic pulse that was partially reflected by the Pt/GdIG interface. Temporal frequencies up to 740 GHz, corresponding to angular frequencies of several THz, were detected in a wavelet analysis of the acoustic oscillations of the strain in the Pt layer. The structural results were analyzed to determine (i) the acoustic damping coefficient and phonon mean free path in Pt at frequencies of hundreds of GHz and (ii) the Grüneisen anharmonicity parameter. The thermal conductance of the Pt/GdIG interface was tracked using the slower, tens-of-picosecond-scale, dynamics of the initial cooling of the heated Pt layer. Analysis using a model based on the Boltzmann transport equation shows that the phonon transmission is lower at the phonon frequencies relevant to thermal transport than for subterahertz regime acoustics.

Acoustic phenomena↗

DOE OSTI2020

High-resolution infrared spectroscopy of jet cooled CH 2 Br radicals: The symmetric CH stretch manifold and absence of nuclear spin cooling

Direct laser absorption of a slit supersonic discharge expansion provides the first high-resolution spectroscopic results on the symmetric CH stretch excitation (ν 1 ) of the bromomethyl (CH 2 Br) radical in the ground electronic state. Narrowband (<1 MHz) mid-infrared radiation is produced by difference-frequency generation of two visible laser beams, with the open shell halohydrocarbon radical generated by electron dissociative attachment of CH 2 Br 2 in a discharge and rapidly cooled to T rot = 18 ± 1 K in the subsequent slit-jet supersonic expansion. A rovibrational structure in the radical spectrum is fully resolved, as well as additional splittings due to spin-rotation effects and 79 Br/ 81 Br isotopologues in natural abundance. Spectroscopic constants and band origins are determined by fitting the transition frequencies to a non-rigid Watson Hamiltonian, yielding results consistent with a vibrationally averaged planar radical and an unpaired electron in the out-of-plane p π orbital. Additionally, extensive satellite band structure from a vibrational hot band is observed and analyzed. The hot band data is compared to CFOUR/VPT2 (CCSD(T)cc-pVQZ) ab initio anharmonic predictions of the vibration rotation alpha matrix, which permits unambiguous assignment to CH 2 symmetric-stretch excitation built on the singly excited CH 2 out-of-plane bending mode (ν 1 + ν 4 ← ν 4 ). Longitudinal cooling of the Doppler width in the slit-jet expansion geometry also reveals partially resolved hyperfine structure on transitions out of the lowest angular momentum states in excellent agreement with predictions based on microwave studies. High level ab initio MOLPRO calculations [CCSD(T)-f12b/VnZ-f12 (n = 3, 4, CBS)] are also performed with explicitly correlated f12 electron methods for the out-of-plane CH 2 bending mode over the halogen series CH 2 X (X = F, Cl, Br, I), which clearly reveals a non-planar geometry for X = F (with a ΔE ≈ 0.3 kcal/mol barrier) and yet planar equilibrium geometries for X = Cl, Br, and I. Finally, a detailed Boltzmann analysis of the transition intensities provides support for negligible collisional equilibration of the entangled H atom nuclear spin states on the few hundred microsecond time scale and high collision densities of a slit supersonic expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Experimental and theoretical studies of the doubly substituted methyl-ethyl Criegee intermediate: Infrared action spectroscopy and unimolecular decay to OH radical products

In this work the infrared (IR) action spectrum of the doubly substituted methyl-ethyl Criegee intermediate (MECI) is observed in the CH stretch overtone region with detection of OH products. The MECI exhibits four conformers, all of which undergo unimolecular decay via a 1,4 H-atom transfer mechanism, followed by the rapid release of OH products. Conformers with different orientations of the carbonyl oxide group with respect to the methyl and ethyl substituents (i.e., anti and syn ) decay via distinct transition state barriers (16.1 kcal mol –1 and 15.4 kcal mol –1 , respectively). The observed IR action spectrum is in good agreement with the predicted anharmonic IR absorption spectrum, but exhibits significant congestion, which is attributed to couplings between spectroscopic bright states and nearby dark states. Energy-dependent OH appearance rates are measured upon IR excitation of the strongest features in the IR action spectrum and are found to be on the order of 10 6 –10 7 s –1 . The experimental rates are in good agreement with computed Rice–Ramsperger–Kassel–Marcus rates for the unimolecular decay of MECI at these energies, which incorporate quantum mechanical tunneling and sophisticated hindered rotor treatments, as well as high-level theoretical calculations of the TS barrier heights, rovibrational properties, and torsional barriers associated with the MECI conformers. Master equation modeling is used to predict thermal rates for the unimolecular decay of anti - and syn -MECI of 473 s –1 and 660 s –1 , respectively. Comparison with other previously studied Criegee intermediate systems provides insights into substituent effects on unimolecular decay under both energy-dependent and thermal conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Examining the order-of-limits problem and lattice constant performance of the Tao–Mo functional

In their recent communication a semi-local density functional derived from the density matrix expansion of the exchange hole localised by a general coordinate transformation. We show that the order-of-limits problem present in the functional, dismissed as harmless in the original publication, causes severe errors in predicted phase transition pressures. We also show that the claim that lattice volume prediction accuracy exceeds that of existing similar functionals was based on comparison to reference data that misses anharmonic zero-point expansion and consequently overestimates accuracy. Here, by highlighting these omissions, we give a more accurate assessment of the Tao-Mo functional and show a possible direction for resolving the order-of-limits problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Charge transfer rate constants for the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran: The spin-boson model vs the linearized semiclassical approximation

Charge transfer rate constants were calculated for the carotenoid-porphyrin-C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran. The calculation was based on mapping the all-atom anharmonic Hamiltonian of this system onto the spin-boson Hamiltonian. The mapping was based on discretizing the spectral density from the time correlation function of the donor–acceptor potential energy gap, as obtained from all-atom molecular dynamics simulations. Different spin-boson Hamiltonians were constructed for each of the possible transitions between the three excited electronic states in two different triad conformations. The rate constants of three possible transitions were calculated via the quantum-mechanically exact Fermi’s golden rule (FGR), as well as a progression of more approximate expressions that lead to the classical Marcus expression. The advantage of the spin-boson approach is that once the mapping is established, the quantum-mechanically exact FGR and the hierarchy of approximations are known in closed form. The classical Marcus charge transfer rate constants obtained with the spin-boson Hamiltonians were found to reproduce those obtained from all-atom simulations with the linearized semiclassical approximation, thereby confirming the equivalence of the two approaches for this system. Furthermore, within the spin-boson Hamiltonian, we also found that the quantum-mechanically exact FGR rate constants were significantly enhanced compared to the classical Marcus theory rate constants for two out of three transitions in one of the two conformations under consideration. The results confirm that mapping to the spin-boson model can yield accurate predictions for charge transfer rate constants in a system as complex as CPC 60 dissolved in tetrahydrofuran.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Understanding lead halide perovskites

It has been ten years since research on lead halide perovskites (LHPs) took off spectacularly and expanded into one of the most active fields today in chemistry, physics, materials science, and engineering. There is no doubt that the rapid growth of this field has been fueled by the exceptional performance of this class of materials in solar cells, light emitting devices, and other electronic, optoelectronic, and photonic devices. Such application oriented advance has been accompanied by major development in mechanistic understanding. Growing mechanistic studies on LHPs have revealed structure–property relationships including, among others, structural instability, unique defect chemistry, dynamic disorder and anharmonicity of phonons, local symmetry breaking, efficient screening of charge carriers in large polaron formation, and the Rashba effect due to spin–orbital coupling and the presence of local electric fields. This special issue collects the cutting edge of our understanding of why LHPs work and what lessons LHPs have taught us about structure and dynamics, in general.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Second-harmonic generation in nonlinear plasmonic lattices enhanced by quantum emitter gain medium

In this work, we report on a theoretical study of second-harmonic generation (SHG) in plasmonic nanostructures interacting with two-level quantum emitters (QEs) under incoherent energy pump. We generalize the driven-dissipative Tavis–Cummings model by introducing the anharmonic surface plasmon-polariton (SPP) mode coupled to QEs and examine physical properties of corresponding SPP-QE polariton states. Our calculations of the SHG efficiency for strong QE-SPP coupling demonstrate orders of magnitude enhancement facilitated by the polariton gain. We further discuss time-domain numerical simulations of SHG in a square lattice comprising Ag nanopillars coupled to QEs utilizing a fully vectorial nonperturbative nonlinear hydrodynamic model for conduction electrons coupled to Maxwell–Bloch equations for QEs. The simulations support the idea of gain enhanced SHG and show orders of magnitude increase in the SHG efficiency as the QEs are tuned in resonance with the lattice plasmon mode and brought above the population inversion threshold by incoherent pumping. By varying pump frequency and tuning QEs to a localized plasmon mode, we demonstrate further enhancement of the SHG efficiency facilitated by strong local electric fields. The incident light polarization dependence of the SHG is examined and related to the symmetries of participating plasmon modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The Li–F–H ternary system at high pressures

Evolutionary crystal structure prediction searches have been employed to explore the ternary Li-F-H system at 300 GPa. Metastable phases were uncovered within the static lattice approximation, with LiF 3 H 2 , LiF 2 H, Li 3 F 4 H, LiF 4 H 4 , Li 2 F 3 H and LiF 3 H lying within 50 meV/atom of the 0 K convex hull. All of these phases contain H n F¯ n+1 (n = 1; 2) anions, and Li + cations. Other structural motifs such as LiF slabs, $H$ $^{+}_{3}$ molecules and F δ- ions are present in some of the low enthalpy Li-F-H structures. The bonding within the H n F¯ n+1 molecules, which may be bent or linear, symmetric or asymmetric, is analyzed. The five phases closest to the hull are insulators, while LiF 3 H is metallic and predicted to have a vanishingly small superconducting critical temperature. Li 3 F 4 H is predicted to be stable at zero pressure. Furthermore, this study lays the foundation for future investigations of the role of temperature and anharmonicity on the stability and properties of compounds and alloys in the Li-F-H ternary system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Multimode two-dimensional vibronic spectroscopy. II. Simulating and extracting vibronic coupling parameters from polarization-selective spectra

We report experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang–Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles

The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

On the temperature and density dependence of dislocation drag from phonon wind

At extreme strain rates, where fast moving dislocations govern plastic deformation, anharmonic phonon scattering imparts a drag force on the dislocations. In this paper, we present calculations of the dislocation drag coefficients of aluminum and copper as functions of temperature and density. We discuss the sensitivity of the drag coefficients to changes in the third-order elastic constants with temperature and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗

DOE OSTI2022

Semi-empirical quantum optics for mid-infrared molecular nanophotonics

Nanoscale infrared (IR) resonators with sub-diffraction limited mode volumes and open geometries have emerged as new platforms for implementing cavity quantum electrodynamics at room temperature. The use of IR nanoantennas and tip nanoprobes to study strong light–matter coupling of molecular vibrations with the vacuum field can be exploited for IR quantum control with nanometer spatial and femtosecond temporal resolution. In order to advance the development of molecule-based quantum nanophotonics in the mid-IR, we propose a generally applicable semi-empirical methodology based on quantum optics to describe light–matter interaction in systems driven by mid-IR femtosecond laser pulses. The theory is shown to reproduce recent experiments on the acceleration of the vibrational relaxation rate in infrared nanostructures. It also provides physical insights on the implementation of coherent phase rotations of the near-field using broadband nanotips. We then apply the quantum framework to develop general tip-design rules for the experimental manipulation of vibrational strong coupling and Fano interference effects in open infrared resonators. We finally propose the possibility of transferring the natural anharmonicity of molecular vibrational levels to the resonator near-field in the weak coupling regime to implement intensity-dependent phase shifts of the coupled system response with strong pulses and develop a vibrational chirping model to understand the effect. The semi-empirical quantum theory is equivalent to first-principles techniques based on Maxwell’s equations, but its lower computational cost suggests its use as a rapid design tool for the development of strongly coupled infrared nanophotonic hardware for applications ranging from quantum control of materials to quantum information processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Infrared spectroscopic signature of a hydroperoxyalkyl radical (•QOOH)

Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH 2 (CH 3 ) 2 COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm –1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. Furthermore, the distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗