Engineering Papers⌕ Search

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

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

Mode-Selective Vibrational Energy Transfer Dynamics in 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX) Thin Films

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. In this work, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N–N and axial NO2 stretching modes with the long-lived, excited phonon bath.

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

High-Accuracy Heats of Formation for Alkane Oxidation: From Small to Large via the Automated CBH-ANL Method

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce ΔH f (o K) values with 2σ uncertainties of 0.2-0.5 kcal mol -1 . This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of ΔH f (o K), in which the reference species may be up to five heavy atoms. The high-level ANLo and ANL1 reference ΔH f (o K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol -1 ) and the propagation of those uncertainties when used in the calculation of ΔH f (o K) for the 194 target species. Here we examine the sensitivity of the predicted ΔH f (o K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

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

Molecular Cage Reports on Its Contents: Spectroscopic Signatures of Cryo-Cooled K + - and Ba 2+ -Benzocryptand Complexes

UV photofragment spectroscopy and IR–UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K + or Ba 2+ imbedded inside it (labeled as K + -BzCrypt, Ba 2+ -BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH 2 -CH 2 -O-CH 2 -CH 2 -O-CH 2 -CH 2 -N chains provides a UV chromophore whose S 0 –S 1 transition we probe. K + -BzCrypt and Ba 2+ -BzCrypt have their S 0 –S 1 origin transitions at 35,925 and 36,446 cm –1 , respectively, blue-shifted by 174 and 695 cm –1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR–UV dip spectroscopy. The alkyl CH stretch region (2800–3000 cm –1 ) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K + -BzCrypt and Ba 2+ -BzCrypt, identifying two conformations (A and B) that lie below all others in energy. Here, we extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH 2 groups adjacent to these atoms. The localized CH 2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH 2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm –1 . As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH 2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.

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

Ab Initio-Based Bond Order Potential for Arsenene Polymorphs Developed via Hierarchical Reinforcement Learning

Arsenene, a less-explored two-dimensional material, holds the potential for applications in wearable electronics, memory devices, and quantum systems. This study introduces a bond-order potential model with Tersoff formalism, the ML-Tersoff, which leverages multireward hierarchical reinforcement learning (RL), trained on an ab initio data set. This data set covers a spectrum of properties for arsenene polymorphs, enhancing our understanding of its mechanical and thermal behaviors without the complexities of traditional models requiring multiple parameter sets. Our RL strategy utilizes decision trees coupled with a hierarchical reward strategy to accelerate convergence in high-dimensional continuous search spaces. Unlike the Stillinger-Weber approach, which demands separate formalisms for buckled and puckered forms, the ML-Tersoff model concurrently captures multiple properties of the two polymorphs by effectively representing the local environment, thereby avoiding the need for different atomic types. Here, we apply the ML model to understand the mechanical and thermal properties of the arsenene polymorphs and nanostructures. We observe an inverse relationship between the critical strain and temperature in arsenene. Thermal conductivity calculations in nanosheets show good agreement with ab initio data, reflecting a decrease in thermal conductivity attributable to increased anharmonic effects at higher temperatures. We also apply the model to predict the thermal behavior of arsenene nanotubes.

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

Infrared Fingerprint and Unimolecular Decay Dynamics of the Hydroperoxyalkyl Intermediate (•QOOH) in Cyclopentane Oxidation

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C 5 H 9 O 2 ), building on higher level reference calculations for ethane oxidation (C 2 H 5 O 2 ). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. Furthermore, the various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Chemical calculations↗
DOE OSTI2025

Discrepancies between Theory and Experiment in Determining the Ionization Energy of NF 3

High-accuracy ab initio thermochemical predictions for the ionization energy of NF 3 , the barrier height (to inversion) of NF 3 + , and the dissociative ionization threshold of NF 3 to NF 2 + + F are presented and incorporated into Active Thermochemical Tables. The adiabatic ionization energy of the first ionization band of NF 3 , calculated at 12.647 ± 0.010 eV, is at odds with previous experimental interpretations by nearly 0.36 eV due to unfavorable Franck-Condon factors associated with this transition. The barrier (to inversion) height is calculated to be about 0.6 eV lower in energy than the prior interpretation, which instigates a discussion of the supposed vibrational structure of the first ionization band of NF 3 . Updated assignments of the photoelectron spectrum are proposed, and the loss in vibrational spacing on the high-energy side of the experimental ionization band is discussed. Rudimentary anharmonic Franck-Condon simulations qualitatively reproduce the broad spectral features observed in experiment.

Active Thermochemical Tables↗
DOE OSTI2025

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗
DOE OSTI2025

Characterization of Two Positional Isomers of the Cs + Gly Complex Using Two-Color, IR–IR Photobleaching of the Cryogenically Cooled Ions

Metal ion binding to amino acid residues is an important interaction motif that controls the tertiary structures of oligopeptides. Analyses of the vibrational band patterns displayed by the amino acid scaffolds are commonly used to characterize the local docking motifs. Here we carry out two-color, IR-IR photobleaching measurements to obtain isomer-selective vibrational spectra of the Cs + Gly ion-molecule complex isolated in a cryogenically cooled, radiofrequency ion trap. The distinct band patterns of two non-interconverting isomers are observed and traced to different bidentate binding motifs between Cs + and the glycine scaffold. In one isomer, the ion attaches to the oxygen atoms of the carboxyl group whereas in the other it docks to the amino nitrogen and the carbonyl oxygen. Attachment to the acid head group yields a very diffuse absorption associated the OH group engaged in a strong intramolecular H-bond that closes a 5 membered ring. Furthermore, the band assignments, rearrangement pathways and electrostatic distortion of the electron density distributions in the glycine scaffold by the proximal ion are explored with electronic structure calculations and anharmonic theory.

Infrared spectroscopy↗
DOE OSTI2026

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗
DOE OSTI2022

Behavior of Long-Chain Hydrocarbons at High Pressures and Temperatures

Although long-chain aliphatic hydrocarbons are documented in meteorites, their origin is poorly understood. A key question is whether they are pristine or a byproduct of terrestrial alteration? Thus, to understand if these long-chain hydrocarbons are indigenous, it will be important to explore their thermodynamic and mechanical stability at conditions experienced by extraterrestrial objects during atmospheric entry and passage. Extreme pressures and temperatures experienced by meteorites are likely to alter the molecular organization of these long-chain hydrocarbons. These structural changes associated with extreme conditions are often documented via laboratory-based Raman spectroscopic measurements. So far, Raman spectroscopic measurements have investigated the effect of static compression on the aliphatic hydrocarbons. The effect of temperature on the structural changes remains poorly explored. To bridge this gap, in this study, we have explored the behavior of two aliphatic hydrocarbons at simultaneously high pressures and temperatures. We have used a resistively heated diamond anvil cell. On compression to moderate pressures, the appearance of new vibrational modes in the low-energy region confirms prior studies and is related to the bending of the linear chains. Upon heating to ~220 °C, we note that the new low-energy mode undergoes softening. The mode softening is likely due to the combination of unbending of the alkane chain and mode anharmonicity.

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

Correctly Assessing Defect Tolerance in Halide Perovskites

The notion of “defect tolerance” has often been invoked to explain the excellent performance of halide perovskites in optoelectronic applications. However, this concept has not been rigorously defined or assessed. A common interpretation is that all of the energetically favorable intrinsic defects are shallow. On the basis of examples in the prototypical halide perovskite CsPbI 3 , we show that this is not the case. The antisite defects Pb I and I Pb are energetically favorable, but also have levels deep in the band gap. Still, because of strong anharmonicity, they are not efficient nonradiative recombination centers. Here, our study demonstrates how to correctly evaluate the “defect tolerance” of halide perovskites.

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

Lattice Dynamics and Optoelectronic Properties of Vacancy-Ordered Double Perovskite Cs 2 TeX 6 (X = Cl – , Br – , I – ) Single Crystals

The soft, dynamic lattice of inorganic lead halide perovskite CsPbX 3 (X = Cl ⁻ , Br ⁻ , I ⁻ ) leads to the emergence of many interesting photophysical and optoelectronic phenomena. However, probing their lattice dynamics with vibrational spectroscopy remains challenging. The influence of the fundamental octahedral building block in the perovskite lattice can be better resolved in zero-dimensional (0D) vacancy-ordered double perovskites of form A 2 BX 6 . Here we study Cs 2 TeX 6 (X = Cl ⁻ , Br ⁻ , I ⁻ ) single crystals to yield detailed insight into the fundamental octahedral building block and to explore the effect that its isolation in the crystal structure has on structural and electronic properties. The isolated [TeX 6 ] 2- octahedral units serve as the vibrational, absorbing, and emitting centers within the crystal. Serving as the vibrational centers, the isolated octahedra inform the likelihood of a random distribution of 10 octahedral symmetries within the mixed-halide spaces, as well as the presence of strong exciton-phonon coupling and anharmonic lattice dynamics. Serving as the absorbing and emitting centers, the isolated octahedra exhibit compositionally tunable absorption (1.50-3.15 eV) and emission (1.31-2.11 eV) energies. Due to greater molecular orbital overlap between neighboring octahedra with increasing halide anion size, there is a transition from a more molecule-like electronic structure in Cs 2 TeCl 6 and Cs 2 TeBr 6 -as expected from the effective 0D nature of these single crystals-to a dispersive electronic structure in Cs 2 TeI 6 , typical of three-dimensional (3D) bulk single crystals.

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

Isomer-specific, Cryogenic Ion Vibrational Spectroscopy Investigation of D2-and N2- Tagged, Protonated Formic Acid Complexes Using Two Color, IR-IR Photobleaching

Here we analyze cryogenic ion vibrational spectra of tagged protonated formic acid (PFA) with electronic structure and anharmonic vibrational calculations to establish the isomers generated by electrospray ionization (ESI) followed by buffer gas cooling to ~25 K. Two isomers are identified (the trans form (E,Z) and the cis form (E,E)) and generated in comparable abundance despite the fact that the calculated E,E structure lies 6.40 kJ mol-1 above the E,Z form. A large (~60 kJ mol-1) barrier separates them such that the E,E form can be kinetically trapped upon cooling in the ion trap. The anticooperativity between the H-bonds of the OH groups is explored by measuring the shift in the D2-bound OH fundamental when a second D2 is attached. Both isomers are observed in the N2-tagged counterparts, displaying the expected red-shifted OH bands. These results indicate that ESI generates both isomers and both must be considered when analyzing cluster spectra based on the PFA core ion.

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

Signatures of Antisymmetric Vibrations in the Ultrafast Dynamics of Quadrupolar Dyes

Antisymmetric molecular vibrations are central to ultrafast, nonadiabatic photophysical and photochemical processes such as conical intersection dynamics, Herzberg–Teller couplings and, potentially, singlet fission. Direct spectroscopic identification of such vibrations is, however, challenging, since they are typically Raman inactive and affect optical transitions only weakly. Here, we report experimental signatures of vibronic couplings to a high-frequency antisymmetric vibration in the excited state dynamics of a quasi-quadrupolar molecule by ultrafast two-dimensional electronic spectroscopy (2DES). The early time, sub-50 fs 2DES maps reveal an asymmetric peak pattern with characteristic low-energy cross-peaks. We show that these peaks arise from stimulated emission transitions from an anharmonic, double-minimum excited state potential energy surface formed by vibronic coupling to a high-frequency antisymmetric mode. Simulations based on a phenomenological essential state model support the results. Our findings offer a new approach for identifying antisymmetric vibrations in ultrafast 2DES and track excited state wavepacket motion before relaxation washes out the spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Transforming Thermal Expansion from Positive to Negative: The Case of Cubic Magnetic Compounds of (Zr,Nb)Fe 2

Negative thermal expansion (NTE) is an intriguing property for not only fundamental studies but also technological applications. However, few NTE materials are available compared with the huge amount of positive thermal expansion materials. The discovery of new NTE materials remains challenging. In this paper, we report a chemical modification strategy to transform thermal expansion from positive to negative in cubic magnetic compounds of (Zr,Nb)Fe 2 by tuning the magnetic exchange interaction. Furthermore, an isotropic zero thermal expansion can be established in Zr 0.8 Nb 0.2 Fe 2 (α 1 = 1.4 x 10 -6 K -1 , 3-470 K) over a broad temperature range that is even wider than that of the prototype Invar alloy of Fe 0.64 Ni 0.36 . The NTE of (Zr,Nb)Fe 2 is originated from the weakened magnetic exchange interaction and the increased d electrons of Fe by the Nb chemical substitution, so that the magnetovolume effect overwhelms the contribution of anharmonic lattice vibration.

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

Vibrational Fermi Resonance in Atomically Thin Black Phosphorus

Fermi resonance is a phenomenon involving the hybridization of two coincidentally quasi-degenerate states that is observed in the vibrational or electronic spectra of molecules. Despite numerous examples in molecular systems, vibrational Fermi resonances in dispersive semiconducting systems remain largely unexplored due to the rarity of occurrence. Here we report a vibrational Fermi resonance in atomically thin black phosphorus. The Fermi resonance arises via anharmonic mixing of a fundamental Raman mode and a Davydov component of an infrared mode, leading to a doublet with mixed character. The extent of Fermi coupling can be modulated by the application of external biaxial strain. The consequences of Fermi hybridization are revealed by electronic resonance effects in the thickness-dependent and excitation-wavelength-dependent Raman spectrum, which is predicted by ab initio hybrid functional simulations including excitonic interactions. This work reveals new insight into electron–phonon coupling in black phosphorus and demonstrates a novel method for modulating Fermi resonances in 2D semiconductors.

74 ATOMIC AND MOLECULAR PHYSICS↗