Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evaluation of Thermal Neutron Scattering Law and Cross Sections for Calcium Hydride

Presented in this paper are the calculated thermal scattering law (TSL) and thermal neutron scattering cross sections for Calcium Hydride, hereafter referred to by its chemical symbol CaH 2 . The only other such data prior to this evaluation are thermal neutron scattering libraries in the JEFF database, which suffer from nonphysical features and inaccuracies. The data in this evaluation are calculated from first principles; Density Functional Theory (DFT) is used to calculate the phonon density of states (DOS), which is the primary input required to calculate the TSL. The TSL and cross sections have been evaluated for the three non-equivalent atom cites in the CaH 2: Ca, H 1 , and H 2 . Each evaluation has been submitted to the NNDC for consideration in the next ENDF/B release.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Controlling extrapolations of nuclear properties with feature selection

Predictions of nuclear properties far from measured data are inherently inaccurate because of uncertainties in our knowledge of nuclear forces and in our treatment of quantum many-body effects in strongly-interacting systems. While the model bias can be directly calculated when experimental data is available, only an estimate can be made in the absence of such measurements. Current approaches to compute the estimated bias quickly lose predictive power when their input variables are taken far from the training region, resulting in uncontrolled uncertainties in applications such as nucleosynthesis simulations. In this letter, we present a novel technique to identify the input variables of machine learning algorithms that can provide robust estimates of model bias. Our process is based on selecting input variables, or features, based on their probability distribution functions across the entire nuclear chart. We illustrate our approach on the problem of quantifying the model bias in nuclear binding energies calculated with Density Functional Theory (DFT). We prove that building model biases with only Z and N as features leads to highly unreliable extrapolations. Conversely, we show that proper feature selection can systematically improve theoretical predictions without increasing uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

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

Electronic Structural Optimization of the Air Electrodes for Reversible Protonic Ceramic Electrochemical Cells via IIIA Cation Doping

Reversible protonic ceramic electrochemical cells (R-PCECs) have emerged as a novel technology for clean and efficient energy generation and storage. Improving the oxygen and proton conduction characteristics and stability of air electrodes at intermediate temperatures is crucial for achieving a high performance. Herein, we optimize the electronic structure of a state-of-the-art PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) air electrode via doping IIIA cations (Al 3+ , Ga 3+ , and In 3+ ). Also, it is shown that PrBa 0.8 Ca 0.2 Co 1.9 Ga 0.10 O 5+δ (PBCCGa 0.10 ) exhibits improved oxygen reaction activity and hydration capability. Density functional theory calculations confirm that Ga doping provides the most favorable electronic structure. An R-PCEC with PBCCGa 0.10 achieves a peak power density of 2.21 W cm -2 and a current density of −4.64 A cm -2 at 1.3 V at 650 °C. Additionally, the PBCCGa 0.10 electrode performs good operational stability in FC mode (for about 100 h), EC mode (for about 100 h), and reversible cyclic testing (over 200 h) at 600 °C.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Role of opacity at the 9 keV back lighter energy used in measuring the equation of state of boron at pressures up to a Gbar

Many experiments have been conducted at the National Ignition Facility to measure the equation of state Hugoniot of plastic, boron, and diamond at extreme pressures up to a Gbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. The shock front conditions are characterized using X-ray radiography, typically at energies near 9 keV. Here we examine how the opacity of boron at 9 keV changes at high pressures and temperatures. Understanding this is vital to unfolding the density in the shock front as pressures exceed 100 Mbar. We compare opacity calculations from a number of methods including the legacy XSN opacity tables, super transition array (STA), average atom (AVAT), and detailed configuration accounting (DCA) methods. We examine how the changing opacity is correlated with the K-shell occupation calculated using the opacity methods as well as other electronic structure methods such as MECCA, Purgatorio, and an all-electron Kohn-Sham density functional theory calculation that uses optimized all-electron norm-conserving Vanderbilt pseudopotentials (ONCV). We also examine the relative contribution of free-free opacity to the overall opacity, which is dominated by bound-free absorption.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Competing itinerant and local spin interactions in kagome metal FeGe

The combination of a geometrically frustrated lattice, and similar energy scales between degrees of freedom endows two-dimensional Kagome metals with a rich array of quantum phases and renders them ideal for studying strong electron correlations and band topology. The Kagome metal, FeGe is a noted example of this, exhibiting A-type collinear antiferromagnetic (AFM) order at T N ≈ 400 K, then establishes a charge density wave (CDW) phase coupled with AFM ordered moment below T CDW ≈ 110 K, and finally forms a c-axis double cone AFM structure around T Canting ≈ 60 K. Here we use neutron scattering to demonstrate the presence of gapless incommensurate spin excitations associated with the double cone AFM structure of FeGe at temperatures well above T Canting and T CDW that merge into gapped commensurate spin waves from the A-type AFM order. Commensurate spin waves follow the Bose factor and fit the Heisenberg Hamiltonian, while the incommensurate spin excitations, emerging below T N where AFM order is commensurate, start to deviate from the Bose factor around T CDW , and peaks at T Canting . This is consistent with a critical scattering of a second order magnetic phase transition with decreasing temperature. By comparing these results with density functional theory calculations, we conclude that the incommensurate magnetic structure arises from the nested Fermi surfaces of itinerant electrons and the formation of a spin density wave order.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electronic structure of reduced CeO 2 (111) surfaces interacting with hydrogen as revealed through electron energy loss spectroscopy in comparison with theoretical investigations

Based on both, ab-initio cluster calculations as well as periodic density functional based loss function calculations, we have assigned the origins of the valence electron excitation regime in electron energy loss spectra of well-ordered ceria films in (111) orientation at various states of reduction as well as after exposing the reduced films to hydrogen from the gas phase (Li et al, Angew. Chem. Int. Ed., 58 (2019) 14686–14693). The explicit calculation of intensity distributions using the dipole approximation allow us to draw conclusions about the nature of oxygen vacancies, which occur upon reduction and how those interact with hydrogen to form hydride species. Herein, we also reports a brief discussion of vibrational excitations of hydrogen loaded ceria and corroborates previous interpretations, based on inelastic neutron scattering.

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

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First-principles calculations of metal surfaces. II. Properties of low-index platinum surfaces toward understanding electron emission

The stability of low-index platinum surfaces and their electronic properties is investigated with density functional theory, toward the goal of understanding the surface structure and electron emission, and identifying precursors to electrical breakdown, on nonideal platinum surfaces. Propensity for electron emission can be related to a local work function, which, in turn, is intimately dependent on the local surface structure. The (1×N) missing row reconstruction of the Pt(110) surface is systematically examined. The (1×3) missing row reconstruction is found to be the lowest in energy, with the (1×2) and (1×4) slightly less stable. In the limit of large (1×N) with wider (111) nanoterraces, the energy accurately approaches the asymptotic limit of the infinite Pt(111) surface. Furthermore, this suggests a local energetic stability of narrow (111) nanoterraces on free Pt surfaces that could be a common structural feature in the complex surface morphologies, leading to work functions consistent with those on thermally grown Pt substrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

1D Magnetic MX 3 Single-Chains (M = Cr, V and X = Cl, Br, I)

Magnetic materials in reduced dimensions are not only excellent platforms for fundamental studies of magnetism, but they play crucial roles in technological advances. The discovery of intrinsic magnetism in monolayer 2D van der Waals systems has sparked enormous interest, but the single-chain limit of 1D magnetic van der Waals materials has been largely unexplored. This paper reports on a family of 1D magnetic van der Waals materials with composition MX 3 (M = Cr, V, and X = Cl, Br, I), prepared in fully-isolated fashion within the protective cores of carbon nanotubes. Atomic-resolution scanning transmission electron microscopy identifies unique structures that differ from the well-known 2D honeycomb lattice MX 3 structure. Density functional theory calculations reveal charge-driven reversible magnetic phase transitions.

1D materials↗
DOE OSTI2024

An Adsorptive Membrane Platform for Precision Ion Separation: Membrane Design and First‐Principles Studies

One of the key challenges in separation science is the lack of precise ion separation methods and mechanistic understanding crucial for efficiently recovering critical materials from complex aqueous matrices. Herein, first-principles electronic structure calculations and in situ Raman spectroscopy are studied to elucidate the factors governing ion discrimination in an adsorptive membrane specifically designed for transition metal ion separation. Density functional theory calculations and in situ Raman data jointly reveal the thermodynamically favorable binding preferences and detailed adsorption mechanisms for competing ions. How membrane binding preferences correlate with the electronic properties of ligands is explored, such as orbital hybridization and electron localization. The findings underscore the importance of the phenolate group in oxime ligands for achieving high selectivity among competing transition metal ions. In-depth understanding on which specific atomistic site within the microenvironment of metal-ligand binding pockets governs the ion discrimination behaviors of the host will build a solid foundation to guide the rational design of next-generation materials for precision separation essential for energy technologies and environment remediation. In tandem, synthetic controllability is demonstrated to transform 3D micrometer-scale crystals to a 2D crystalline selective layer in membranes, paving the way for more precise and sustainable advances in separation science.

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

Unveiling the Role of Lithium Iodide in Stabilizing Solid Interfaces in All-Solid-State Li Metal Batteries

A critical challenge in all-solid-state lithium metal batteries (ASSLMBs) is achieving a stable interface between the lithium metal anode and the solid electrolyte. Leveraging its success in Li/I 2 batteries, lithium iodide has garnered significant attentions for its potential to enhance interfacial stability and overall cell performance in ASSLMBs. Here, we elucidate the role of lithium iodide in stabilizing the solid interface in all-solid-state Li metal batteries with a Li argyrodite electrolyte, particularly focusing on its influence on lithium deposition behavior and interfacial evolution. Through in situ optical imaging, we demonstrate more uniform lithium deposition on an iodide-contained argyrodite electrolyte compared to a chloride-based counterpart. Complementary density functional theory calculations attribute improved lithium plating behavior to the enhanced lithiophilicity and better ionic conductivity of lithium iodide at the solid interface, effectively reducing localized current density. In conclusion, these findings provide useful insights into the mechanisms through which lithium iodide enhances the interfacial stability in ASSLMBs.

Anions↗
DOE OSTI2022

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗
DOE OSTI2024

Corrosion Engineering of Part‐Per‐Million Single Atom Pt 1 /Ni(OH) 2 Electrocatalyst for PET Upcycling at Ampere‐Level Current Density

The plastic waste issue has posed a series of formidable challenges for the ecological environment and human health. While conventional recycling strategies often lead to plastic down-cycling, the electrochemical strategy of recovering valuable monomers enables an ideal, circular plastic economy. Here a corrosion synthesized single atom Pt 1 /Ni(OH) 2 electrocatalyst with part-per-million noble Pt loading for highly efficient and selective upcycling of polyethylene terephthalate (PET) into valuable chemicals (potassium diformate and terephthalic acid) and green hydrogen is reported. Electro-oxidation of PET hydrolysate, ethylene glycol (EG), to formate is processed with high Faraday efficiency (FE) and selectivity (>90%) at the current density close to 1000 mA cm −2 (1.444 V vs RHE). The in situ spectroscopy and density functional theory calculations provide insights into the mechanism and the understanding of the high efficiency. Remarkably, the electro-oxidation of EG at the ampere-level current density is also successfully illustrated by using a membrane-electrode assembly with high FEs to formate integrated with hydrogen production for 500 h of continuous operation. In conclusion, this process allows valuable chemical production at high space-time yield and is highly profitable (588–700 $\$$ ton −1 PET), showing an industrial perspective on single-atom catalysis of electrochemical plastic upcycling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

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

CO 2 Hydrogenation and Formic Acid Dehydrogenation Using Ir Catalysts with Amide-Based Ligands

In this work, a series of Ir catalysts bearing amide-based ligands generated by a deprotonated amide moiety was prepared with the hypotheses that the strong electron-donating ability of the coordinated anionic nitrogen atom and the proton-responsive OH group near the metal center will improve the catalytic activity for CO 2 hydrogenation and formic acid (FA) dehydrogenation. The effects of the modifications of the ligand architecture on the catalytic activity were investigated for CO 2 hydrogenation at ambient conditions (25 °C with 0.1 MPa H 2 /CO 2 (v/v = 1/1)) and under slightly harsher conditions (50 °C with 1.0 MPa H 2 /CO 2 ) in basic aqueous solutions together with deuterium kinetic isotope effects (KIEs) with selected catalysts. Cp*Ir( L12 )(H 2 O)HSO 4 ( L12 = 6-hydroxy- N -phenylpicolinamidate) that has an anionic coordinating N atom and an OH group in the second coordination sphere, exhibits a turnover frequency (TOF) of 198 h –1 based on the initial 1 h of reaction. This TOF which, to the best of our knowledge, is the highest value ever reported under ambient conditions in basic aqueous solutions. However, Cp*Ir( L10 )(H 2 O)HSO 4 ( L10 = (4-hydroxy- N -methylpicolinamidate) performs better in long-term CO 2 hydrogenation (up to a TON of 14 700 with [Ir] = 10 μM after 348 h and the final formate concentration of 0.643 M with [Ir] = 250 μM) at ambient conditions. Further, the catalytic activity for FA dehydrogenation was examined under three different conditions (pH 1.6, 2.3, and 3.5). The Cp*Ir( L12 )(H 2 O)HSO 4 complex in any of these conditions is less active compared to the picolinamidate catalysts without ortho -OH, owing to its instability. The complex without OH group, Cp*Ir( L8 )(H 2 O)HSO 4 ( L8 = N -phenyl-picolinamidate), exhibits a high TOF (up to 118 000 h -1 ) at 60 °C. Theoretical calculations were performed to examine the catalytic mechanism, and a step-by-step mechanism has been proposed for both CO 2 hydrogenation and FA dehydrogenation reactions. Density functional theory calculations of [Cp*Ir( L3 )(H 2 O)]HSO 4 ( L3 = picolinamidate) and the X-ray structure of the [Cp*Ir( L7 )(H)]·H 2 O ( L7 = N -methylpicolinamidate) complex imply a pH-dependent conformational change from N , N coordination to N , O coordination upon lowering the pH of the aqueous solution.

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

Electric field sensitivity of molecular color centers

Molecular color centers with S=1 ground states are promising candidates for quantum sensing of electric fields. These molecules have an electronic structure similar to solid state color centers, but they allow for processing modalities that permit direct interfacing with an analyte. Currently, it is unknown how sensitive these molecules are to electric fields and what molecular properties affect their sensitivity. We perform density functional theory calculations to understand the impact of electric fields on the electronic structure of five nominally tetrahedral molecular color centers exhibiting variable transition metal chemistry and ligand densities. We then extract the Stark parameters from each of these molecules and compare them to molecular properties such as the dipole moment and inner shell stiffness and find that the dipole moment of the molecule largely governs sensitivity. We predict that polar heteroleptic molecules may have electric field sensitivities comparable to solid state color centers such as nitrogen-vacancy centers in diamond.

Physics↗